When compiling a map, we know that the more precise the control points are, the more precise the map will be. This paper applies Integrated Geodesy approach, which integrates all available geodetic observations in a rigorous collocation model, to calculate the three dimensional coordinates of control points as the reference of the map compilation. At the same time, we can also get the functions of the disturbing potential-such as the geoidal undulation, the deflection of the vertical, and the gravity anomaly-to make a research in the gravity field and to compile the gravity field map.