One of the key technologies of geomagnetic matching navigation is looking for high accuracy geomagnetic map. Based on introducing the principles of four kinds of interpolation methods (inverse distance to a power method, Kriging method, modified Shepard-s method and radial basis function method), the data of geomagnetic survey was used, and a region where the survey points were uniform was selected. Geomagnetic maps were created by all the interpolation methods, and cross-validation was used to estimate the methods. The results show that the geomagnetic map created by radial basis function method has little "bull-s-eyes", the contour is much smoother, and its standard deviation is much less than other methods, the standard error and the mean deviation is the minimum. The interpolation precision with radial basis function method is high, so the geomagnetic map created by the method better reflects the regional distribution of the earth-s magnetic field.