An image reconstruction method is described as follows. A structural image of an object is obtained. An actual optical detected result of the object is obtained. An inhomogeneous initial state is determined based on the structural image. At least one solution converged with the actual optical detected result is determined by iteratively calculating a physical model from the inhomogeneous initial state. The image of the object is reconstructed based on the solution.