This paper presents a subspace-based two-step iterative shrinkage/thresholding method(S-TwIST) based on the Distorted Born iterative method (DBIM) to improve the performance of the original TwIST inverse algorithm. This method retrieves the deterministic part of the induced current from inhomogeneous Green's function operator leading to more accurate total field calculation at each iteration step than that of the original TwIST. Both inverse algorithms have been evaluated with a set of synthetic geometries with fine structures. Compared with TwIST, the results show that S-TwIST has superior accuracy in multiple objects profile (εerr=0.1454%) and 1/16λ resolution at 2GHz. Also, S-TwIST is more robust to initial guess, which means it is less likely to become unstable when the inversion procedure starts without initial guess.