本发明公开了一种饱和磁化状态下磁声磁粒子浓度图像重建方法，该方法包括：基于仿真模型，获取超声换能器处的声压数据和重建区域处的梯度磁场数据；利用每个超声换能器处获取的声压数据构建声压矩阵，利用所述重建区域处的梯度磁场数据构建系统矩阵，该系统矩阵用于描述声压与磁性纳米粒子的浓度偏导数之间存在直接的对应关系；再将两者用LSQR方法获取磁性纳米粒子的浓度偏导数分布，即可获取磁性纳米粒子浓度分布。本发明用LSQR‑梯形公式法，能够快速、稳定、高质量地重建出浓度分布，图像边界清晰，具有良好的应用前景。