Koordinat Dönüşümü — ED50, ITRF, TUREF

Türkiye'de kullanılan coğrafi koordinat sistemleri tarihsel süreçte köklü bir dönüşüm yaşamıştır. Yarım asır boyunca geçerli olan ED50 (European Datum 1950) yerel datumundan, ITRF (International Terrestrial Reference Frame) tabanlı jeosantrik sisteme ve ulusal TUREF (Turkish Reference Frame) datumuna geçiş 2005 yılında yasal zorunluluk haline gelmiştir. Bu makale; datum kavramlarını, Türkiye'nin projeksiyon sistemini (UTM/Gauss-Krüger) ve koordinat dönüşüm hesaplamalarını (Helmert 7P, NTv2/Grid, epoch öteleme) açıklar. Türkiye'nin aktif levha hareketi (~2,5 cm/yıl) nedeniyle koordinatların zaman bağımlı doğası özellikle vurgulanmaktadır.

1. Datum ve Koordinat Sistemleri

1.1 ED50 (European Datum 1950)

ED50, Avrupa kıtası için 1950 yılında tanımlanmış yerel bir datumdur. Referans elipsoidi olarak International 1924 (Hayford) elipsoidini kullanır:

• Büyük yarı eksen: $a = 6\,378\,388{,}0$ m
• Basıklık: $f = 1/297{,}0$
• Başlangıç noktası: Potsdam, Almanya (jeosantrik değil; yerel datum)

Türkiye'de 1950'lerden 2005'e kadar resmi datum olarak kullanılmış; kadastro paftaları, imar haritaları ve topografik haritaların büyük bölümü bu datumda üretilmiştir. Türkiye Meşedağ Datumu ile ED50 arasındaki fark enlem yönünde 46,3 m, boylam yönünde 20,6 m olarak belirlenmiştir.

Tablo 1: ED50 (European Datum 1950)

Dikkat: ED50 jeosantrik değildir; GPS/GNSS çıktılarıyla doğrudan uyumsuzdur. Eski ED50 paftalarını GNSS verisiyle birlikte kullanan uygulamalarda datum dönüşümü zorunludur.

Saha Notu: Sahada eski ED50 paftasından alınan köşe koordinatları ile TUSAGA-Aktif'ten elde edilen TUREF koordinatları Türkiye orta bölgesinde ~11–15 m, doğu bölgelerinde ~12–18 m, batı bölgelerinde ~10–14 m ayrışabilir. Aplikasyon öncesi mutlaka dönüşüm yapılmalıdır.

1.2 WGS84 (World Geodetic System 1984)

WGS84, ABD Savunma Bakanlığı tarafından tanımlanmış jeosantrik bir datum ve küresel koordinat sistemidir. GPS alıcıları varsayılan olarak WGS84 koordinatı üretir.

• Referans elipsoid: GRS80 ile pratikte özdeş ($a = 6\,378\,137{,}0$ m, $f = 1/298{,}257$)
• TUREF ile fark: < 1 m (operasyonel uygulamalarda özdeş kabul edilir)
• ITRF gerçekleşimleriyle fark: ±1–2 cm düzeyinde

1.3 ITRF (International Terrestrial Reference Frame)

ITRF, Uluslararası Dünyanın Döndürme ve Referans Sistemleri Servisi (IERS) tarafından yayımlanan jeosantrik, 4-boyutlu referans çerçevesidir. Dünya genelinde ~500 GNSS, VLBI, SLR ve DORIS istasyonunun gözlemlerine dayanır.

Koordinatlar zaman bağımlıdır; her istasyonun hem koordinatı hem de hız vektörü (mm/yıl) tanımlanır. Farklı epoğ'lardaki koordinatlar aşağıdaki bağıntıyla ilişkilendirilir:

$$\vec{X}(t) = \vec{X}(t_0) + \dot{\vec{X}} \cdot (t - t_0)$$

Burada $\dot{\vec{X}}$ mm/yıl cinsinden istasyon hız vektörüdür.

Başlıca ITRF sürümleri ve Türkiye ile ilgisi:

Tablo 2: ITRF (International Terrestrial Reference Frame)

1.4 TUREF (Turkish Reference Frame)

TUREF, Türkiye'nin resmi ulusal jeosantrik referans çerçevesidir.

• Referans elipsoid: GRS80
• Dayandığı ITRF sürümü: ITRF96
• Referans epoğu: 2005.0
• Gerçekleşim ağı: TUSAGA-Aktif (158 sabit GNSS istasyonu)

Tablo 3: TUREF (Turkish Reference Frame)

Yasal statü: BÖHHBÜY 2018 Madde 10'a göre tüm büyük ölçekli harita üretimleri TUREF'e dayalı GRS80 elipsoidinde ve Transversal Mercator (TM) 3°'lik dilim esasına göre gerçekleştirilmelidir.

2. Projeksiyon Sistemi

2.1 UTM (Universal Transverse Mercator) — 6° Dilim

Küçük ölçekli haritalar (1:25.000 ve daha küçük) ile askeri haritacılıkta kullanılır.

Tablo 4: UTM (Universal Transverse Mercator) — 6° Dilim

UTM temel parametreleri: ölçek faktörü $k_0 = 0{,}9996$ (merkez meridyeninde), false easting $E_0 = 500\,000$ m. 180 km'lik dilim kenarında ölçek faktörü ~1,0004'e yükselir; dilim ortasında ise ~0,9996 ile en düşük düzeydedir.

UTM Projeksiyon temel formülleri:

Boylamsal eğrilik yarıçapı:

$$N = \frac{a}{\sqrt{1 - e^2 \sin^2\varphi}}$$

Doğu koordinatı (Easting):

$$E = E_0 + k_0 N \left[\ell + \frac{1-t^2+\eta^2}{6}\ell^3 + \frac{5-18t^2+t^4+14\eta^2-58t^2\eta^2}{120}\ell^5\right]$$

Burada $k_0 = 0{,}9996$ (UTM), $\ell = \lambda - \lambda_0$, $t = \tan\varphi$, $\eta^2 = e'^2\cos^2\varphi$.

2.2 Gauss-Krüger TM — 3° Dilim

BÖHHBÜY 2018 Madde 10 gereğince büyük ölçekli harita üretiminde (1:5.000 ve daha büyük ölçekler) zorunlu olarak kullanılır.

Tablo 5: Gauss-Krüger TM — 3° Dilim

Gauss-Krüger temel parametreleri: ölçek faktörü $k_0 = 1{,}0000$ (merkez meridyeninde; UTM'den farklı!), false easting $E_0 = 500\,000$ m. En büyük ölçek bozulması ≈ 0,01% — UTM'den 4 kat daha iyi.

Saha Notu: CAD veya GIS yazılımında veri açarken 3° dilim numarasını (zon 9-15 arası) ve TUREF datumunu kontrol edin. Eski ED50 verisini yanlış dilimle açmak 500 m'ye varan yatay kaymaya yol açabilir.

3. Datum Dönüşümü

3.1 Helmert (7-Parametre) Dönüşümü

ED50 ↔ ITRF/WGS84 dönüşümü için kullanılan standart 3-boyutlu benzerlik dönüşümüdür:

$$\begin{pmatrix} X' \\ Y' \\ Z' \end{pmatrix} = (1 + \delta s) \begin{pmatrix} 1 & R_z & -R_y \\ -R_z & 1 & R_x \\ R_y & -R_x & 1 \end{pmatrix} \begin{pmatrix} X \\ Y \\ Z \end{pmatrix} + \begin{pmatrix} T_x \\ T_y \\ T_z \end{pmatrix}$$

Burada:

• $T_x, T_y, T_z$ — öteleme parametreleri (m)
• $R_x, R_y, R_z$ — rotasyon parametreleri (radyan veya arc-second)
• $\delta s$ — ölçek faktörü (ppm — parts per million)

ED50 → WGS84 yaklaşık parametreleri (Türkiye, EPSG:1315):

Tablo 6: Helmert (7-Parametre) Dönüşümü

Dikkat: Tablo değerleri tüm Türkiye için ortalama/yaklaşık değerlerdir. Kadastral işlemler ve halihazır harita üretiminde TKGM tarafından yayımlanan bölgesel (hücresel) dönüşüm parametreleri kullanılmalıdır.

Dönüşüm doğrulukları:

• Standart Helmert 7P: ±0,30–0,50 m (tek nokta başına)
• NTv2/Grid yöntemi (bölgesel hücresel): ±2–5 cm

3.2 TUREF → ED50 Dönüşüm Yöntemleri

Yöntem 1 — Grid (NTv2): TKGM'nin .gsb formatındaki grid dosyasıyla interpolasyon tabanlı dönüşüm. Doğruluk: ±2–5 cm. PROJ 9.x, QGIS, ArcGIS Pro, NetCAD tarafından desteklenir.

Yöntem 2 — Helmert 7P: Matematiksel olarak daha basit, doğruluğu ±0,30–0,50 m. Kadastral işlemler için yetersizdir.

Yöntem 3 — TKGM Online Servisi: tkgm.gov.tr üzerinden ED50 ↔ TUREF dönüşümü, UTM hesabı ve jeoit düzeltmesi; resmi ve güvenilir sonuç üretir.

4. Epoch Dönüşümü ve Levha Hareketi

Türkiye, Anadolu Plakası üzerinde yer alır. Bu plaka her yıl hareket ettiğinden ITRF96 datumundaki koordinatlar zamana göre değişir.

Anadolu Plakası hız değerleri:

• Ortalama: ~2,5 cm/yıl (batı yönünde)
• Batı Anadolu/Ege kıyıları: ~3,5–4 cm/yıl (GB yönünde)
• Kuzey Anadolu Fayı (Marmara): 25 mm/yıl

Farklı epoğ'lardaki koordinatlar arasında dönüşüm:

$$X(2005{,}0) = X(t_0) + (2005{,}0 - t_0) \cdot V_x$$ $$Y(2005{,}0) = Y(t_0) + (2005{,}0 - t_0) \cdot V_y$$ $$Z(2005{,}0) = Z(t_0) + (2005{,}0 - t_0) \cdot V_z$$

Tablo 7: Epoch Dönüşümü ve Levha Hareketi

Dikkat: 2025 yılında yapılan bir GNSS ölçümü, ölçü epoğunda (T ≈ 2025.0) koordinat üretir. Bu koordinat referans epoğu 2005.0'a dönüştürülmeden üst ölçekli ağa bağlanamaz. Türkiye ortası için bu fark: 20 yıl × 2,5 cm/yıl = ~50 cm düzeyindedir.

5. Yükseklik Sistemi

5.1 Üç Yükseklik Kavramı

Tablo 8: Üç Yükseklik Kavramı

Türkiye'de Helmert ortometrik yükseklik sistemi esas alınmaktadır. Referans: Antalya Mareografı.

Temel dönüşüm ilişkisi:

$$H = h - N_{TG09}$$

Burada $N_{TG09}$ jeoid ondülasyonudur (TG-09 modelinden elde edilir).

5.2 TG-09 Türkiye Jeoid Modeli

• THG-09 (Hibrit jeoid): TG-09 + GPS/nivelman jeoid yükseklikleri birleşimi
  • İç doğruluk (RMS): ±0,95 cm
  • Dış doğruluk: ±8,4 cm
• Türkiye'de $N_{TG09}$ değeri: +15 m ile +46 m arasında (kıyılarda düşük, iç bölgelerde yüksek)

Saha Notu: TUSAGA-Aktif Ağ-RTK ölçümlerinde alıcı yazılımı TG-09 jeoid modelini otomatik uygular. Ancak alıcı ayarlarında "jeoid modeli" seçeneğinin TG-09 olarak ayarlandığını kontrol edin; EGM96 veya başka bir modelin seçili olması ±30 cm farka yol açabilir.

6. Pratik Dönüşüm Araçları

Tablo 9: Pratik Dönüşüm Araçları

7. Örnek Problemler

Problem 1 — Epoch Öteleme Hesabı

Veri:

• GNSS çıktısı (ölçü epoğu 1998.0): $X = 4\,181\,235{,}412$ m, $Y = 2\,546\,870{,}134$ m, $Z = 4\,089\,221{,}756$ m
• Hız vektörü: $V_x = -0{,}0238$ m/yıl, $V_y = +0{,}0172$ m/yıl, $V_z = +0{,}0097$ m/yıl

İstenen: Koordinatları TUSAGA-Aktif referans epoğu 2005.0'a öteleyiniz.

Çözüm:

$\Delta t = 2005{,}0 - 1998{,}0 = 7{,}0$ yıl

$$X(2005{,}0) = 4\,181\,235{,}412 + 7 \times (-0{,}0238) = 4\,181\,235{,}412 - 0{,}167 = 4\,181\,235{,}245 \; \text{m}$$ $$Y(2005{,}0) = 2\,546\,870{,}134 + 7 \times 0{,}0172 = 2\,546\,870{,}134 + 0{,}120 = 2\,546\,870{,}254 \; \text{m}$$ $$Z(2005{,}0) = 4\,089\,221{,}756 + 7 \times 0{,}0097 = 4\,089\,221{,}756 + 0{,}068 = 4\,089\,221{,}824 \; \text{m}$$

Kontrol: Toplam yer değiştirme: $\sqrt{0{,}167^2 + 0{,}120^2 + 0{,}068^2} \approx 0{,}217$ m ≈ 21,7 cm (7 yılda ~3 cm/yıl hız; Türkiye için mantıklı).

Problem 2 — UTM ve Gauss-Krüger Dilim Belirleme

Veri: TUREF koordinatları: $\varphi = 39°30'00{,}000''$ K, $\lambda = 33°00'00{,}000''$ D, $h = 1\,050{,}000$ m.

İstenen: UTM dilimini, Gauss-Krüger 3° dilim numarasını ve EPSG kodunu belirleyiniz.

Çözüm:

Adım 1 — UTM dilimi:

$$\text{UTM Dilim No.} = \left\lfloor\frac{33{,}0 + 180}{6}\right\rfloor + 1 = \left\lfloor\frac{213}{6}\right\rfloor + 1 = 35 + 1 = 36\text{N}$$

Merkez meridyeni: $\lambda_0 = 33°$ (tam olarak merkez meridyeninde; en düşük UTM bozulması).

Adım 2 — 3° Gauss-Krüger dilimi:

$$\text{3° Zon No.} = \left\lfloor\frac{33{,}0}{3}\right\rfloor = 11 \Rightarrow \text{Merkez meridyeni: } 33°$$

Zon 11 → TUREF/TM33 → EPSG: 5255

Tablo 10: Problem 2 — UTM ve Gauss-Krüger Dilim Belirleme

Sonuç: Büyük ölçekli harita üretiminde EPSG:5255 kullanılmalıdır. Kontrol: $\lambda = 33°$ tam zon 11 orta meridyenine denk geldiğinden Gauss-Krüger koordinatında $E = 500\,000{,}000$ m olmalıdır.

Problem 3 — Tam Koordinat Dönüşümü ve Yükseklik Hesabı

Veri:

• WGS84 coğrafi koordinatları: $\varphi_{WGS84} = 40°12'35{,}782''$ K, $\lambda_{WGS84} = 29°04'48{,}331''$ D, $h_{WGS84} = 214{,}837$ m
• TG-09 jeoid ondülasyonu: $N_{TG09} = +38{,}24$ m
• Proje koordinat sistemi: TUREF / 3° TM Zon 10 (EPSG: 5254, merkez meridyeni 30°)
• ED50 kaykılık değerleri (HGM tablosu): $\Delta\varphi = +3{,}78''$, $\Delta\lambda = +5{,}92''$

İstenen: (1) WGS84 → TUREF dönüşümü; (2) Ortometrik yükseklik; (3) TUREF/TM Zon 10 Easting; (4) ED50 koordinatı.

Çözüm:

Adım 1 — WGS84 → TUREF (WGS84 ile TUREF/ITRF96 farkı < 1 m; operasyonel uygulamalarda özdeş):

$$\varphi_{TUREF} \approx 40°12'35{,}782'' \text{ K}, \quad \lambda_{TUREF} \approx 29°04'48{,}331'' \text{ D}$$

Adım 2 — Ortometrik yükseklik:

$$H = h - N_{TG09} = 214{,}837 - 38{,}24 = 176{,}597 \; \text{m} \approx 176{,}60 \; \text{m}$$

Adım 3 — TM Zon 10 Easting (yaklaşık, ilk terim; $\lambda_0 = 30°$):

$$\ell = 29°04'48{,}331'' - 30° = -0°55'11{,}669'' = -0{,}016058 \; \text{rad}$$ $$E = 500\,000 + k_0 \cdot N_{eğr} \cdot \ell \cdot \cos\varphi \approx 500\,000 + 0{,}9996 \cdot 6\,388\,124 \cdot (-0{,}016058) \cdot 0{,}7632$$ $$E \approx 500\,000 - 78\,323 = 421\,677 \; \text{m}$$

Adım 4 — ED50 koordinatı (HGM datum dönüşüm yöntemi):

$$\varphi_{ED50} = 40°12'35{,}782'' + 3{,}78'' = 40°12'39{,}562'' \text{ K}$$ $$\lambda_{ED50} = 29°04'48{,}331'' + 5{,}92'' = 29°04'54{,}251'' \text{ D}$$

Konum farkı büyüklüğü: $\Delta y \approx 116{,}8$ m, $\Delta x \approx 140{,}3$ m → vektör ~183 m (Türkiye kuzey-batı için tipik ED50–WGS84 farkı).

Sonuç: Ortometrik yükseklik $H = 176{,}60$ m; TUREF/TM Zon 10 Easting $E \approx 421\,677$ m (EPSG: 5254).

8. Sık Yapılan Hatalar

Tablo 11: Sık Yapılan Hatalar

9. Referanslar ve Standartlar

• BÖHHBÜY 2018, Madde 10 (RG 30403, 30.04.2018) — GRS80 + TM 3° dilim, TUREF zorunluluğu
• BÖHHBÜY 2005, Madde 5 (RG 26046, 15.07.2005) — 2005 sonrası resmi datum TUREF
• TKGM TUREF Teknik Raporu, 2008 — ITRF96, epoch 2005.0 teknik tanımı
• IERS Teknik Not No. 36 — ITRF2014 teknik dokümantasyonu
• 3194 İmar Kanunu, Madde 7/a, 38 — Harita üretim mevzuatına dayanak
• 3402 Kadastro Kanunu, Madde 47 — Kadastro haritalarında koordinat sistemi
• Kahveci M., Yıldız F. (2019). GPS/GNSS — Uydularla Konum Belirleme. Nobel Yayıncılık, Ankara.
• Kılıçoğlu A. vd. (2011). "Türkiye'de ED50-TUREF Koordinat Dönüşümü." HGM Jeodezi Dairesi Teknik Raporu.
• EPSG Registry — epsg.org (Projeksiyon kodları ve parametreleri)
• PROJ 9.x dokümantasyonu — proj.org

Kaynakça

1. İlgili Türk Standartları (TS) ve Avrupa Normları (EN)
2. TBDY 2018 — Türkiye Bina Deprem Yönetmeliği
3. İlgili ders kitapları ve teknik kaynaklar

Not: Bu makale eğitim amaçlıdır. Projelerde güncel yönetmelik ve standartlara başvurunuz.

Kaynaklar

1. BÖHHBÜY 2018.
2. Harita ve Ölçme.

Önemli Mühendislik Uyarısı: Bu içerik yalnızca bilgilendirme amaçlıdır; nihai tasarım, hesap ve uygulama kararları, güncel yönetmelikler ile proje koşulları çerçevesinde yetkili bir inşaat mühendisinin denetiminde alınmalıdır. Sayısal örnekler ve formüller genel mühendislik pratiğini yansıtır; her projenin kendine özgü zemin, yük ve çevre koşulları proje müellifince ayrıca değerlendirilmelidir.