ความแตกต่างอย่างมากระหว่างผลลัพธ์สำหรับการแปลงตัวเลข 7-param

ฉันพยายามที่จะแปลง WGS84 lat / lon point

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

ในแพ็คเกจซอฟต์แวร์ที่แตกต่างกันสามชุด (Proj4, GeoTrans และ Leica GeoOffice) แต่ฉันได้รับความแตกต่างอย่างมากระหว่างผลลัพธ์ - ประมาณสองสามร้อยเมตรใน X และ Y! โปรดทราบว่ารูปแบบเหล่านี้จะปรากฏเฉพาะกับ 7 พารามิเตอร์ที่ไม่ใช่ 3

Proj4

cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + ถึง + proj = utm + zone = 32 + ellps = clrk80 + towgs84 = 178.3,316.7,131.5, -5.278, -5.278, -6.277, -10.9792, South + units = 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

ผลลัพธ์ตามลำดับ:

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

ฉันเชื่อว่า 3 แพ็คเกจทั้งหมดใช้วิธีการทางคณิตศาสตร์ที่เหมือนกันสำหรับการแปลง 7-param (วิธี Bursa-Wolf) อะไรเป็นสาเหตุของการเปลี่ยนแปลงครั้งใหญ่นี้

อันดับแรก Proj4 ใช้สิ่งที่ EPSG เรียกว่า "Position Vector" รุ่นของวิธีการ 7 พารามิเตอร์ เป็นไปได้ว่า GeoTrans และ Leica GeoOffice ใช้รุ่นอื่นที่ EPSG เรียกว่า "Coordinate Frame" ทั้งสองวิธีมีความเท่าเทียมกัน แต่เมทริกซ์การหมุนจะแตกต่างกันและสัญญาณของพารามิเตอร์เชิงมุมจะต้องเปลี่ยน

ประการที่สองขอบคุณที่ให้ภาพหน้าจอของคำนิยามการเปลี่ยนแปลงใน Leica GeoOffice นิยาม Proj4 หรือนิยามนี้ไม่ถูกต้อง ควรกำหนดพารามิเตอร์การกำหนด Ellipsoid A และ Ellipsoid B ปัจจุบันการแปลงนี้แปลงจาก WGS84 เป็น Congo60 ใน Proj4 โปรดทราบว่าตัวเลือกคือ + towgs84 ดังนั้นจึงกำหนดเป็นจาก Congo60 ถึง WGS84 ในการเปลี่ยนทิศทางในคำจำกัดความให้เปลี่ยนสัญลักษณ์ของพารามิเตอร์ทั้งหมด ตรวจสอบความช่วยเหลือสำหรับ GeoOffice และดูว่าพารามิเตอร์ SF ต้องการชิ้นส่วนต่อล้านรุ่นหรือค่าที่แปลงเป็นมาตราส่วนแล้ว

ฉันไม่รู้เกี่ยวกับ GeoTrans - คุณหมายถึงซอฟต์แวร์ NGA หรือไม่ อย่างไรก็ตามหวังว่าคุณจะสามารถรับ GeoOffice และ Proj4 ให้เข้ากันได้

สำหรับ proj4 สัญญาณของพารามิเตอร์จะต้องกลับด้าน

ดูหน้าคำนิยามนี้:

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"]]

ฉันไม่แน่ใจว่าคุณได้รับค่าของคุณจากที่ใด

QGIS กำหนด EPSG: 28232 Point Noire UTM 32S เป็น:

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

และผลลัพธ์พารามิเตอร์ 7 รายการไม่ควรอยู่ไกลจากนั้น

แก้ไข

ด้วย GDAL 1.10 ฉันได้ผลลัพธ์ดังต่อไปนี้:

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

ห่างกันน้อยกว่า 3 เมตร โปรดทราบว่า E และ S ต้องเป็นไปตามองศาโดยไม่มีช่องว่าง

และสำหรับบันทึก: การตั้งค่า towgs84 เป็นศูนย์ผลลัพธ์ใน:

800310.94   9422829.37 -109.32

ในขณะที่เปลี่ยนจาก Point Noire degrees EPSG: 4262 เป็น Point Noire UTM ให้ผลลัพธ์เป็น:

800311.21   9422892.49 0.00