Sự khác biệt lớn giữa các kết quả cho biến đổi mốc 7-param

Tôi đang cố gắng chuyển đổi điểm lat / lon WGS84

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

trong ba gói phần mềm khác nhau (Proj4, GeoTrans và Leica GeoOffice) nhưng tôi nhận được sự khác biệt rất lớn giữa các kết quả - khoảng vài trăm mét trong X và Y! Lưu ý rằng các biến thể này chỉ xuất hiện với 7 tham số không phải với 3.

Dự án 4

cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + to + proj = utm + zone = 32 + ellps = clrk80 + Towgie84 = 178.3.316.7.131.5, -5.278, -6.077, -10. nam + đơn vị = 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 GeoScript

Các kết quả tương ứng là:

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

Tôi tin rằng cả 3 gói đều sử dụng các phương pháp toán học giống nhau cho phép biến đổi 7-param (phương pháp Bursa-Wolf). Điều gì có thể là nguyên nhân của sự thay đổi lớn này?

Đầu tiên, Proj4 sử dụng cái mà EPSG gọi là phiên bản "Vector vị trí" của phương thức 7 tham số. Có thể GeoTrans và Leica GeoOffice sử dụng phiên bản khác mà EPSG gọi là "Khung tọa độ". Cả hai phương pháp đều tương đương, nhưng ma trận xoay là khác nhau và các dấu hiệu của các tham số góc phải được thay đổi.

Thứ hai, cảm ơn bạn đã cung cấp ảnh chụp màn hình của định nghĩa chuyển đổi trong Leica GeoOffice. Định nghĩa Proj4 hoặc định nghĩa này được định nghĩa không chính xác. Các tham số xác định Ellipsoid A và Ellipsoid B nên được chuyển đổi. Hiện tại, sự chuyển đổi này chuyển đổi từ WGS84 sang Congo60. Trong Proj4, lưu ý rằng tùy chọn là + Towks84, do đó, nó được xác định là TỪ Congo60 ĐẾN WGS84. Để thay đổi hướng trong định nghĩa, thay đổi dấu của TẤT CẢ các tham số. Ngoài ra, hãy kiểm tra trợ giúp cho GeoScript và xem liệu tham số SF muốn các phần trên một triệu phiên bản hoặc giá trị yếu tố đã được chuyển đổi thành tỷ lệ.

Tôi không biết về GeoTrans - ý bạn là phần mềm NGA? Dù sao, hy vọng bạn có thể có được GeoScript và Proj4 để khớp.

Đối với proj4, các dấu hiệu của các tham số phải được đảo ngược.

Xem trang định nghĩa này:

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

Tôi không chắc bạn lấy giá trị từ đâu.

QGIS định nghĩa EPSG: 28,02 điểm Noire UTM 32S là:

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

và kết quả 7 tham số không nên xa đó.

BIÊN TẬP

Với GDAL 1.10, tôi nhận được các kết quả sau:

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

Vì vậy, cách nhau chưa đầy 3 mét. Lưu ý rằng E và S phải tuân theo độ không có khoảng trống.

Và, chỉ dành cho bản ghi: đặt kết quả84 thành 0, kết quả là:

800310.94   9422829.37 -109.32

trong khi chuyển đổi từ Point Noire độ EPSG: 4262 sang Point Noire UTM cho kết quả:

800311.21   9422892.49 0.00