Abstract: A method and system for efficiently converting an image from a first map projection to a second map projection image. The method includes precomputing coordinates in the first map projection for a subset of the total number of pixels in an empty second map projection image and subsequently, finding the first map projection coordinates of each of the remaining pixels by linear interpolation using the precomputed first map projection coordinates of the nearest surrounding precomputed pixels. A color is assigned to the second map projection image pixel with bilinear interpolation using color values of surrounding first map projection pixels. If the first map projection has fiducials, only the area within the fiducials is considered. An embodiment converts a USGS DOQ in UTM format into Geodetic coordinates.