In order to apply Kriging methods for geostatistics of spatial data, an estimation of spatial coherency functions is required priorly based on the spatial distance between measurement points. In the study, the typical coherency functions, such as semi-variogram, homeogram, and covariance function, were estimated using the national geoid model. The test area consisting of 2°×2° and the Unified Control Points (UCPs) within the area were chosen as sampling measurements of the geoid. Based on the distance between the control points, a total of 100 sampling points were grouped into distinct pairs and assigned into a bin. Empirical values, which were calculated with each of the spatial coherency functions, resulted out as a wave model of a semi-variogram for the best quality of fit. Both of homeogram and covariance functions were better fitted into the exponential model. In the future, the methods of various Kriging and the functions of estimated spatial coherency need to be studied to verify the prediction accuracy and to calculate the Mean Squared Prediction Error (MSPE).