An electric properties tomography method for reconstructing a spatial distribution of electric conductivity (σ) from magnetic resonance image data representative of a magnetic resonance image of at least a portion of a subject of interest, the spatial distribution covering at least a portion of the area of the magnetic resonance image, and the method comprising following steps: - a step of segmenting the magnetic resonance image into voxels corresponding to a volume of interest and voxels outside the volume of interest, wherein the voxels corresponding to the volume of interest and the voxels outside the volume of interest are separated by a segmentation boundary; - a step of extrapolating acquired phase values of at least the magnetic resonance image data corresponding to voxels of the volume of interest proximal to the segmentation boundary to obtain extrapolated phase values for voxels of a volume outside and adjacent the volume of interest; - a step of replacing acquired phase values of magnetic resonance data corresponding to voxels of the volume outside and adjacent the volume of interest by the extrapolated phase values; - a step of transforming into the frequency domain a numerical kernel representing a numerical second derivative, the acquired phase values corresponding to the voxels of the volume of interest, and the extrapolated phase values corresponding to the voxels outside and adjacent the volume of interest; - a step of multiplying the frequency domain-transformed numerical second derivative by the frequency domain-transformed acquired phase values and the frequency domain-transformed numerical second derivative by the frequency domain-transformed extrapolated phase values, respectively; and - a step of transforming the result of the step of multiplying into the spatial domain. - a step of segmenting the magnetic resonance image into voxels corresponding to a volume of interest and voxels outside the volume of interest, wherein the voxels corresponding to the volume of interest and the voxels outside the volume of interest are separated by a segmentation boundary; - a step of extrapolating acquired phase values of at least the magnetic resonance image data corresponding to voxels of the volume of interest proximal to the segmentation boundary to obtain extrapolated phase values for voxels of a volume outside and adjacent the volume of interest; - a step of replacing acquired phase values of magnetic resonance data corresponding to voxels of the volume outside and adjacent the volume of interest by the extrapolated phase values; - a step of transforming into the frequency domain a numerical kernel representing a numerical second derivative, the acquired phase values corresponding to the voxels of the volume of interest, and the extrapolated phase values corresponding to the voxels outside and adjacent the volume of interest; - a step of multiplying the frequency domain-transformed numerical second derivative by the frequency domain-transformed acquired phase values and the frequency domain-transformed numerical second derivative by the frequency domain-transformed extrapolated phase values, respectively; and - a step of transforming the result of the step of multiplying into the spatial domain.