Abstract: Kriging method is an interpolating method based on the spatial statistical correlation of the samples. The core idea implicated in Kriging method is that each sample is assigned with different weight according to the spatial correlation among the sample points, and the estimated error is minimized. So, it can be summarized to be the best linear unbiased estimator of a random function. Non-dual Kriging method is a local interpolator for the interpolation at each node and the solution of a new Kriging linear system is required by which the location of interpolating node is explicitly depended on. Therefore, non-dual Kriging method is quite time consuming for the solution of a new Kriging linear system for every interpolating point. By equivalent transform, the Universal Kriging method can be transformed to Dual Kriging method, which is a global interpolator for its Kriging linear system is independent of the interpolating point. Therefore, the Kriging linear system is solved by the Dual Kriging method only once to interpolate all points, so the computational efficiency is significantly improved and is of great value in meteorology and oceanography where large data sets are to be interpolated. Furthermore, the result of the numerical experiment shows that Dual Kriging method is not only equivalent to Universal Kriging method in accuracy and comparable to or superior to the Cressman method built in GrADS, but also far superior to Universal Kriging method in computational efficiency. In order to improve the accuracy, efficiency and flexibility of Dual Kriging method, several unique techniques are adopted in the implementation of the computational scheme:the solution of Dual Kriging linear system using partial pivoting LU decomposition and iterative improvement, the fitting of the trend by SVD linear fitting method, Levenberg-Marquatdt iterative method for the non-linear parameters fitting of the semivariance model, and the flexible parameter configuration by the FORTRAN 90 modular interface. Finally, four kinds of semivariance expressions are obtained by fitting the statistical sample semivariance derived from the statistics of summer precipitation in the Eastern and Southern China, and the sensitivity of these four kinds of semivariance models to the accuracy of precipitation interpolation is analyzed using the Dual Kriging method. It is found that exponential semivariance model is resulted in the best analysis, spherical semivarance model is better, Gaussian and linear semivariance is the worst. It is reasonable since the exponential semivariance model is the smoothest and its range is the longest, then the interpolation is made more accurate by the smooth weight contributions from more samples around the interpolating point. Therefore, the longer the range is, the more smooth the semivariance varies, the more accurately the precipitation within the isotropic range is interpolated by the Dual Kriging method. In summary, the Dual Kriging method is superior to non-dual Kriging method and Cressman method, and as an efficient and accurate best linear unbiased interpolator, more applications in meteorology and oceanography will be gained.