The invention belongs to the technical field of space-to-ground observation, and particularly discloses a method for calculating a three-dimensional intersection angle of satellite images. The method comprises the following steps of A, acquiring image point coordinates of same-name points of a three-dimensional satellite image pair by using a manual measurement or automatic satellite image matching method; B, implementing an RFM (rational function model)-based area network compensation method by utilizing the coordinates of the same-name points of the satellite images to realize the optimization of RPC (rational polynomial coefficient) parameters; C, calculating spatial rectangular coordinates corresponding to object-space points formed by projecting the same-name points onto projection surfaces at different elevations by using a projection trajectory method to finish the spatial reconstruction of same-name light; and D, calculating the three-dimensional intersection angle of the same-name light by utilizing direction vectors of the reconstructed same-name light. The method is applied to the calculation of the three-dimensional intersection angle of any three-dimensional model formed by overlapped satellite images without any ground control point, can be used for the mapping accuracy evaluation of the three-dimensional model, and has remarkable economic and social benefits.