Abstract: In order to further explore the reconstruction algorithm of threedimensional natural magnetization, the profile reconstruction of a simplified threedimensional distribution, double layer structure, was studied in this paper. Magnetization reconstruction scheme was developed based on the concept of magnetic charge and the optimization algorithm of steepest descent method to cope with the illposeness of this typical inverse problem. In practice, a system of linear equations linking the magnetic charges of double layer structure with the corresponding leakage magnetic field was deduced at first. An initial value assignment method was developed to regularize the illposed inverse problem for magnetic charge reconstruction. Several magnetic charges distributions of double layer structure were reconstructed with the proposed method. The numerical results showed that the method was effective for the reconstruction of magnetic charge in double layer structure.