Systems and methods for determining a distribution map of susceptibility property of an object are provided. The method may include one or more of the following operations. A phase diagram corresponding to a magnetic resonance (MR) signal of the object may be obtained. A preliminary field map may be determined based on the phase diagram. Preliminary error limiting information associated with the preliminary field map may be obtained. A preliminary distribution map of susceptibility property of the object may be determined based on the preliminary field map and the preliminary error limiting information. An iteration process including at least one iteration may be performed to determine a target distribution map of susceptibility property of the object.