. Imaging of superparamagnetic iron oxide nanoparticles based on their non-linear response to alternating magnetic fields shows promise for imaging cells and vasculature in healthy and diseased tissue. Such imaging can be achieved through x-space reconstruction typically along a unidirectional Cartesian trajectory, which rapidly convolutes the particle distribution with a 'anisotropic blurring' point spread function (PSF), leading to images with anisotropic resolution.. Here we propose combining the time domine-system matrix and x-space reconstruction methods into a forward model, where the output of the forward model is the PSF-blurred x-space reconstructed image. We then treat the blur as an inverse problem solved by Kaczmarz iteration.. After we have proposed the method optimization, the normal resolution of simulation and device images has been increased from 3.5 mm and 5.25 mm to 1.5 mm and 3.25 mm, which has reached the level in the tangential resolution. Quantitative indicators of image quality such as PSNR and SSIM have also been greatly improved.. Simulation and imaging of real phantoms indicate that our approach provides better isotropic resolution and image quality than the x-space method alone or other methods for removing PSF blur. Using our proposed method to optimize the image quality of x-space reconstructed images using unidirectional Cartesian trajectories, it will promote the clinical application of MPI in the future.