Kriging yönteminin geoit modellemesinde kullanılabilirliğinin araştırılması
Geoit, ağırlık potansiyelinin sabit olduğu ortalama deniz yüzeyine yakın bir yüzeydir. GPS ile bulunan enlem ve boylam değerleri doğrudan kullanılmakta ancak elipsoidal yükseklik (h) değerlerinin ortometrik yüksekliğe (H) dönüştürülmesi gerekmektedir. Dönüşüm için yeterli doğrulukta geoit yüksekliklerinin bilinmesi gerekmektedir. Geoit yükseklikleri belirleme teknikleri içerisinde en yaygın kullanılanı GPS/Nivelman tekniğidir. BÖHHBÜY dönüşüm için yerel GPS/Nivelman geoidi kullanılmasını öngörmektedir. Bu çalışmada, İstanbul Metropoliten alanlarında GPS/Nivelman yöntemi ile geoit belirlemek için deterministik ve Kriging enterpolasyon yöntemleri karşılaştırılarak, Kriging yönteminin geoit yüzeyi modellemesinde kullanılabilirliğinin araştırılması yapılmıştır. Uygulamada, ArcGIS 8.3 Geostatstical Analyst yazılımı kullanılmıştır. Dayanak nokta yoğunluğunun geoit hesabına etkisini araştırmak için İstanbul’ da 50 (~ 100 km2’ye bir nokta), 200 (~ 25 km2’ye bir nokta), 393(~ 13 km2’ye bir nokta) ve 434 (~ 12 km2’ye bir nokta) noktalı geoit yüzeyleri oluşturulmuştur. Test için 50 nokta seçilmiş ve bu noktaların deterministik ve geoistatistik enterpolasyonla $N_{HESAP}$ değerleri hesaplanmıştır. $N_{HESAP-ÖLÇÜ}$ fark değerlerinden, farkların maksimum, minimum ve ortalama değerleri ile karesel ortalama hatalar karşılaştırılarak en uygun yüzeyi veren yöntem seçimi yapılmıştır. Çalışmalar sonunda, Kriging yöntemi ile geoit belirleme sonuçlarının deterministik yöntemlerden daha presizyonlu olduğu, Ordinary Kriging yönteminin Simple Kriging yönteminden bulunan sonuçlardan daha presizyonlu olduğu saptanmıştır. Multiquadratik yöntemle bulunan sonuçların Kriging yöntemlerinden bulunan sonuçlara çok yakın olduğu ve multiquadratik yöntemin en iyi sonucu veren deterministik yöntem olduğu, presizyonu arttırmak için nokta yoğunluğunu arttırmaktan çok veri kalitesini arttırmak gerektiği saptanmıştır.
Searching the use of Kriging method on geoid surface modeling
The geoid is a representation of the surface of earth that assumes the sea is covered the earth, also known as surface of equal gravitational attraction and mean sea level. The main function of the geoid in geodesy is to serve as a reference surface for leveling. The elevation measured by leveling is relative to the geoid. Three dimensional coordinates that are easily obtained by GPS easily; are widely used in such applications as large-scale map production and engineering applications. The latitude and longitude values obtained by GPS are used directly, but ellipsoidal height (h) values must be transformed into orthometric heights (H). In many surveying and engineering applications, orthometric heights are required. Ellipsoidal heights have geometric meanings in practical surveying, engineering, and geophysics and in other applications, and they bear no physical meanings. For the transformation from ellipsoidal heights to orthometric heights, which are used in applications, geoid heights (N) must be known with required accuracy. Several techniques can be used for determination of geoid heights. Most commonly used method for the determination of geoid is the combination of GPS data and leveling measurements. The by laws for large scale map production is required to be changed and allow to GPS applications, as parallel to the augmentation in application of the GPS technique. The new draft bylaw contains observation and calculation methods for obtaining orthometric heights by GPS, as well. One of them is determining the height of the geoid by GPS/Nivelman method. In this study, compared with classical methods by the deterministic and Kriging interpolation methods for determining the geoids by GPS/Nivelman method in Istanbul Metropolitan Surface deterministic and Kriging interpolations were calculated by known geoid undulation (NMEASURE) values performed by using ArcGIS 8.3 Geostatistical Analyst software. Geostatistical Analyst provides deterministic and geostatistical interpolation methods. Deterministic interpolation techniques (inverse distance weighted, radial basis functions, and local polynomial interpolation) should not be used for decision making, because they do not provide information on how good their predictions are. Geostatistical interpolation techniques (e.g., kriging) can be chosen based on the result of exploratory spatial data analysis and diagnostics (cross validation and validation). Deterministic methods use predefined mathematical functions for interpolation. Geostatistical methods rely on statistical features of the data. Geostatistical Analyst provides the necessary tools for data exploration and variography analysis. Kriging is based on the assumption that the parameter being interpolated can be treated as a regionalized variable. A regionalized variable is intermediate between a truly random variable and a completely deterministic variable in that it varies in a continuous manner from one location to the next and therefore points that are near each other have a certain degree of spatial correlation. Kriging is a set of linear regression routines which minimize estimation variance from a predefined covariance model. This method uses variogram to express the spatial variation, and it minimizes the error of predicted values which are estimated by spatial distribution of the predicted values. Besides, while the Geoid Model was forming; the effect of control point frequency to computed geoid height values was investigated. To study the effect of point frequency to the counting of geoid, geoid surfaces were formed with 50(~one point to 100 km2), 200(~one point to 25 km2), 393 (~one point to 13 km2) and 434(~one point to 12 km2) points in Istanbul. For testing with interpolations, 50 test points were chosen and $N_{CALCULUS}$ values were counted by various deterministic and geostatistical interpolations. $N_{CALCULUS-MEASURE}$ distinction values were found by subtracting $N_{CALCULUS}$ values from measurement values that were found by $N_{GPS/NİVELMAN}$ method. Mean square errors and maximum, minimum and mean error values were calculated for all methods. By comparing of these values, the method that gives the most suitable surface was chosen. It was found that the results of determining geoid by Kriging method were more precise than deterministic methods. The results of Ordinary Kriging method were more precisely than Simple Kriging method’s results but the maps produced by Simple Kriging method were more esthetic than the other equal geoid height maps. The results found by multiquadratic method between deterministic methods were close to the results found by Kriging method and multiquadratic method gave the best results. For increasing precession augmentation of survey quality is better than augmentation of model points.
___
- Akçın, H., (1998). GPS ölçülerinden pratik yüksekliklerin elde edilmesi üzerine bir çalışma, Doktora Tezi, YTÜ, Fen Bilimleri Enstitüsü, İstanbul.
- Aksoy, A., ve Güneş, İ.H., (1990). Jeodezi 1, İTÜ Kütüphanesi Sayı:1422, İstanbul.
- Aksoy, A.,Deniz, R., ve Ayan, T., (1999). Global, bölgesel ve ülke jeodezik ağları hakkında, Harita Kadastro Mühendisliği Dergisi, 86, Ankara.
- Ayan, T., ve Deniz, R., (2000). Fiziksel Jeodezi Ders Notu, İTÜ İnşaat Fakültesi, Jeodezi ve Fotogrametri Mühendisliği Bölümü, Jeodezi Anabilim Dalı, İstanbul.
- Boogaart KG. ve Schaeben, H., (2002). Kriging of regionalized directions Axes and oriantations, 1. Directions and axes, Mathematical Geology, 34, 479-503.
- BÖHHUY, (2005). Büyük Ölçekli Harita ve Harita Bilgileri Üretim Yönetmeliği. Harita ve Kadastro Mühendisleri Odası, Ankara.
- Burrough, P. A. ve McDonnell, R. A., (1998). Principles of Geographical Information Systems, Oxford University Pres, New York.
- Erkanlı, Y., (1986). Koordinatlandırılmış modelde ve alanda enterpolasyon, kollakasyon yöntemlerinin uygulanması ve neticeleri, Harita ve Kadastro Mühendisliği Dergisi, Ankara, 86, 56-57.
- Fogel, D.N. ve Tınney,L.R., (1996). Image registration using multiquadratic functions the finite element method bivariate mapping polynomials and thin plate spline, Technical Report, Santa Barbara.
- Güler, A., (1978). Sayısal arazi modellerinde interpoasyon yöntemi, Harita Dergisi, 85, 53-71, Ankara.
- İnal, C. ve Yiğit, C., (2003). Lokal alanlarda jeoit ondülasyonlarının belirlenmesinde kullanılan enterpolasyon yöntemlerinin karşılaştırılması, Selçuk Üniversitesi Jeodezi ve Fotogrametri Mühendisliği Öğretiminde 30. Yıl Sempozyumu, 97-106, Konya.
- Isaaks, E. H. ve Srivastava, M. R., (1989). An ıntroduction to applied geostatistics, Oxford University Pres, New York.
- Kleijnen, J., (2006). Design and Analysis of Monte Carlo Experiments,http://www. quantlet. Com/mdstat/scripts/csa/html/node139.html,(10.01.2006).
- Kılıçoğlu, A. ve Fırat, O., (2003). Büyük Ölçekli Harita üretiminde GPS ile ortometrik yükseklik belirlemeye yönelik jeoit modelleme ve uygulamalar, TUJK 2003 Yılı Bilimsel Toplantısı, Coğrafi Bilgi Sistemleri ve Jeodezik Ağlar Çalıştayı, 27-34, Konya.
- Lang, C., (2006). Kriging Interpolation. http://www.nbb.cornell.edu/neurobio/land/OldStudentProjects/cs490-94to95/clang/kriging.html, (10.01.2006).
- Matheron, G., (1963). Principles of geostatistcs, Economic Geology, 58, 1246-1266
- Mert, B. Mert, B.A., (2005). Jeoistatistiksel analiz için bir bilgisayar programının gelistirilmesi ve Antalya-Akseki-Kiziltas boksit yatagina uygulanmasi, Yüksek Lisans Tezi, Çukurova Üniversitesi, Fen Bilimleri Enstitüsü, Adana.
- Rohuani, S., ve Wackernagel, H., (1990). Multivariate Geostatistical Approach to Space-Time Data Analysis. Water Resources Research 26, (4), 585-591.
- Tercan, A.E., Saraç, C., (1998). Maden yataklarının değerlendirilmesinde jeoistatistiksel yöntemler, TMMOB. Maden Mühendisleri Odası Yayını, Ankara.
- Vieira, S. R., Hatfield, J. R., Nielsen, D. R. ve Biggar, J. W., (1983). Geostatistical theory and application to variability of some agronomical properties, Hilgardia, 51, 3, 1-75.
- Yanalak, M. ve İnce, C., (1997). GPS ile elde edilen elipsoid yüksekliklerinin yerel yükseklik sistemine dönüştürülmesi, VI. Harita Kurultayı. Ankara, 03-07 Mart 1996.
- Yiğit, C. Ö., (2003). Elipsoidal yüksekliklerin ortometrik yüksekliğe dönüşümünde kullanılan enterpolasyon yöntemlerinin karşılaştırılması. Selçuk Üniversitesi FBE, Yüksek Lisans Tezi, Konya.