Modern giant segmented mirror telescopes (GSMTs) such as the Extremely Large Telescope, which is currently under construction, depend heavily on adaptive optics (AO) systems to correct for atmospheric distortions. However, a residual blur always remains in the astronomical images corrected by single conjugate AO (SCAO) systems due to fitting and bandwidth errors, which can mathematically be described by a convolution of the true image with a point spread function (PSF). Due to the nature of the turbulent atmosphere and its correction, the PSF is spatially varying, which is known as an anisoplanatic effect. The PSF serves, e.g., as a quality measure for science images and therefore needs to be known as accurately as possible. In this paper, we present an algorithm for PSF reconstruction from pupil-plane data in directions apart from the guide star direction in an SCAO system. Our algorithm is adapted to the needs of GSMTs focused on estimating the contribution of the anisoplanatic and generalized fitting error to the PSF. Results obtained in an end-to-end simulation tool show a qualitatively good reconstruction of the PSF compared to the PSF calculated directly from the simulated incoming wavefront as well as stable performance with respect to imprecise knowledge of atmospheric parameters.