A system for analyzing and improving accuracies of two-dimensional images and for forming three-dimensional images from the two-dimensional images, the system comprising: a processor; and a memory that includes a plurality of two-dimensional images and instructions, the plurality of two-dimensional images each including a same target area and being acquired at same or different times and at different collection angles, the instructions configured to, when executed by the processor, cause the processor to execute operations comprising; identifying a plurality of target features in the target area of each of the two-dimensional images; determining, independently for each of the plurality of two-dimensional images and based on image pointing parameters, a three-dimensional geolocation position for each of the plurality of target features, wherein each of the plurality of two-dimensional images comprises a stereo pair of images, and the three-dimensional geolocation position for each of the plurality of target features in each of the plurality of two-dimensional images is determined for each stereo pair of images by measuring three mutually orthogonal reference coordinates of each of the plurality of target features; calculating a weighted average of the three-dimensional geolocation positions for each of the plurality of target features across the plurality of two-dimensional images, wherein a weight factor of each of the plurality of two-dimensional images contribution to each of the three mutually orthogonal coordinates is an image error, in each coordinate, squared and divided by a sum of squared errors for all of the plurality of two-dimensional images, in each coordinate; adjusting, variably across each of the plurality of two-dimensional images, the image pointing parameters by providing an adjustment of the image pointing parameters to minimize a geolocation difference between the three-dimensional geolocation position of each of the plurality of target features in each of the plurality of two-dimensional images and the weighted average of the three-dimensional geolocation positions of each of the plurality of target features across the plurality of two-dimensional images, wherein the image pointing parameters comprise polynomial coefficients; and projecting each of the plurality of two-dimensional images onto a georeferenced three-dimensional surface model of the target area based on results of the adjusting to form georeferenced three-dimensional images from the plurality of two-dimensional images. a processor; and a memory that includes a plurality of two-dimensional images and instructions, the plurality of two-dimensional images each including a same target area and being acquired at same or different times and at different collection angles, the instructions configured to, when executed by the processor, cause the processor to execute operations comprising; identifying a plurality of target features in the target area of each of the two-dimensional images; determining, independently for each of the plurality of two-dimensional images and based on image pointing parameters, a three-dimensional geolocation position for each of the plurality of target features, wherein each of the plurality of two-dimensional images comprises a stereo pair of images, and the three-dimensional geolocation position for each of the plurality of target features in each of the plurality of two-dimensional images is determined for each stereo pair of images by measuring three mutually orthogonal reference coordinates of each of the plurality of target features; calculating a weighted average of the three-dimensional geolocation positions for each of the plurality of target features across the plurality of two-dimensional images, wherein a weight factor of each of the plurality of two-dimensional images contribution to each of the three mutually orthogonal coordinates is an image error, in each coordinate, squared and divided by a sum of squared errors for all of the plurality of two-dimensional images, in each coordinate; adjusting, variably across each of the plurality of two-dimensional images, the image pointing parameters by providing an adjustment of the image pointing parameters to minimize a geolocation difference between the three-dimensional geolocation position of each of the plurality of target features in each of the plurality of two-dimensional images and the weighted average of the three-dimensional geolocation positions of each of the plurality of target features across the plurality of two-dimensional images, wherein the image pointing parameters comprise polynomial coefficients; and projecting each of the plurality of two-dimensional images onto a georeferenced three-dimensional surface model of the target area based on results of the adjusting to form georeferenced three-dimensional images from the plurality of two-dimensional images.