7-param datum dönüşümü sonuçları arasında büyük farklılıklar

Bir WGS84 lat / lon noktasını dönüştürmeye çalışıyorum

11d42'32.10629" E 5d12'56.75371" S

üç farklı yazılım paketinde (Proj4, GeoTrans ve Leica GeoOffice) ama sonuçlar arasında çok büyük farklar alıyorum - X ve Y'de birkaç yüz metre! Bu varyasyonların 3 ile değil sadece 7 parametrede göründüğünü unutmayın.

Proj4

cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + ila + proj = utm + bölge = 32 + ellps = clrk80 + towgs84 = 178.3,316.7,131.5, -5.278, -6.077, -10.9792, -19.166 + güney + birimler = m + no_defs

Geotrans

Delta X 178.3 
Delta Y 316.7 
Delta Z 131.5 
Rot X -5.278 
Rot Y -6.077 
Rot Z -10.9792 
SF = -19.166 / 0.999980834 (0.000019166)

Leica GeoOffice

Sonuçlar sırasıyla:

1. 800392 9422525
2. 800306 9422840
3. 800941 9422891

3 paketin hepsinin 7-param dönüşüm için aynı matematiksel yöntemleri kullandığına inanıyorum (Bursa-Wolf yöntemi). Bu büyük varyasyonun nedeni ne olabilir?

İlk olarak Proj4, EPSG'nin 7 parametre yönteminin "Konum Vektörü" sürümünü adlandırdığı şeyi kullanır. GeoTrans ve Leica GeoOffice'in EPSG'nin "Koordinat Çerçevesi" adlı diğer sürümünü kullanması mümkündür. Her iki yöntem de eşdeğerdir, ancak dönme matrisleri farklıdır ve açısal parametrelerin işaretleri değiştirilmelidir.

İkinci olarak, Leica GeoOffice'deki dönüşüm tanımının ekran görüntüsünü sağladığınız için teşekkür ederiz. Ya Proj4 tanımı ya da bu tanım yanlış tanımlanmış. Elipsoid A ve Elipsoid B tanımlayıcı parametreleri değiştirilmelidir. Şu anda, bu dönüşüm WGS84'ten Kongo60'a dönüşüyor. Proj4'te, seçeneğin + towgs84 olduğuna dikkat edin, bu yüzden Kongo'dan FOSA WGS84 olarak tanımlanır. Tanımdaki yönü değiştirmek için TÜM parametrelerin işaretlerini değiştirin. Ayrıca GeoOffice için yardımı kontrol edin ve SF parametresinin milyon sürüm başına parça isteyip istemediğini veya önceden ölçeklendirilmiş bir faktör katsayısı değerini isteyip istemediğini görün.

GeoTrans'ı bilmiyorum - NGA yazılımı mı demek istiyorsun? Neyse, umarım GeoOffice ve Proj4'ü eşleştirebilirsiniz.

Proj4 için parametrelerin işaretleri tersine çevrilmelidir.

Bu tanım sayfasına bakın:

http://www.spatialreference.org/ref/epsg/62826405/prettywkt/

GEOGCS["Pointe Noire (deg)",
    DATUM["Congo 1960 Pointe Noire",
        SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936269,
            AUTHORITY["EPSG","7011"]],
        TOWGS84[-178.3,-316.7,-131.5,5.278,6.077,10.979,3.953271276531849],
        AUTHORITY["EPSG","6282"]],
    PRIMEM["Greenwich",0.0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.017453292519943295],
    AXIS["Geodetic latitude",NORTH],
    AXIS["Geodetic longitude",EAST],
    AUTHORITY["EPSG","62826405"]]

Değerlerinizi nereden aldığınızdan emin değilim.

QGIS, EPSG: 28232 Point Noire UTM 32S'yi şu şekilde tanımlar:

+proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-148,51,-291,0,0,0,0 +units=m +no_defs

ve 7 parametre sonucu bundan çok uzak olmamalıdır.

DÜZENLE

GDAL 1.10 ile aşağıdaki sonuçları elde ediyorum:

4326-proj-3 parameters:
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-148,51,-291,0,0,0,0 +units=m +no_defs
800232.21   9423131.96 -1.76
4326-proj-7 parameters from spatialrefrence.org:
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-178.3,-316.7,-131.5,5.278,6.077,10.979,3.953271276531849 +units=m +no_defs
800230.13   9423133.46 91.31
4326-proj-7 parameters from proj4 datum_shift.csv (EPSG:1802):
cs2cs +init=epsg:4326 +to +proj=utm +zone=32 +south +a=6378249.2 +b=6356515 +towgs84=-178.3,-316.7,-131.5,5.278,6.077,10.979,19.166 +units=m +no_defs
800230.13   9423133.40 -5.72

Arada 3 metreden daha az. E ve S'nin boşluk bırakmadan dereceleri izlemesi gerektiğini unutmayın.

Ve sadece kayıt için: towgs84'ü sıfıra ayarlamak şu sonuçlarla sonuçlanır:

800310.94   9422829.37 -109.32

Point Noire derecelerinden EPSG: 4262'den Point Noire UTM'ye dönüştürürken:

800311.21   9422892.49 0.00