In the paper a transformation between two height datums (Kronstadt’60 and Kronstadt’86, the latter being a part of the present National Spatial Reference System in Poland) with the use of geostatistical method – kriging is presented. As the height differences between the two datums reveal visible trend a natural decision is to use the kind of kriging method that takes into account nonstationarity in the average behavior of the spatial process (height differences between the two datums). Hence, two methods were applied: hybrid technique (a method combining Trend Surface Analysis with ordinary kriging on least squares residuals) and universal kriging. The background of the two methods has been presented. The two methods were compared with respect to the prediction capabilities in a process of crossvalidation and additionally they were compared to the results obtained by applying a polynomial regression transformation model. The results obtained within this study prove that the structure hidden in the residual part of the model and used in kriging methods may improve prediction capabilities of the transformation model.