7- 파라미터 데이텀 변환 결과의 큰 차이

WGS84 위도 / 경도 포인트를 변환하려고합니다.

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

세 가지 소프트웨어 패키지 (Proj4, GeoTrans 및 Leica GeoOffice)에서 X와 Y에서 수백 미터 정도의 결과를 얻을 수 있습니다. 이러한 변형은 3이 아닌 7 개의 매개 변수에만 나타납니다.

Proj4

cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + to + proj = utm + zone = 32 + ellps = clrk80 + towgs84 = 178.3,316.7,131.5, -5.278, -6.077, -10.9792, -19.166 + 남쪽 + 단위 = m + no_defs

지오 트랜스

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)

라이카 지오 오피스

결과는 각각 다음과 같습니다.

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

3 패키지 모두 7 매개 변수 변환 (Bursa-Wolf 방법)에 대해 동일한 수학적 방법을 사용한다고 생각합니다. 이 거대한 변형의 원인은 무엇입니까?

먼저, Proj4는 EPSG가 7 매개 변수 방법의 "위치 벡터"버전이라고하는 것을 사용합니다. GeoTrans와 Leica GeoOffice는 EPSG가 "좌표 프레임"이라고하는 다른 버전을 사용할 수 있습니다. 두 방법 모두 동일하지만 회전 행렬이 다르므로 각도 매개 변수의 부호를 변경해야합니다.

둘째, Leica GeoOffice에서 변환 정의 스크린 샷을 제공해 주셔서 감사합니다. Proj4 정의 또는이 정의가 잘못 정의되었습니다. Ellipsoid A와 Ellipsoid B를 정의하는 파라미터를 전환해야합니다. 현재이 변환은 WGS84에서 Congo60으로 변환됩니다. Proj4에서 옵션은 + towgs84이므로 FROM Congo60 TO 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를 0으로 설정하면 다음과 같은 결과가 발생합니다.

800310.94   9422829.37 -109.32

Point Noire도 EPSG : 4262에서 Point Noire UTM으로 변환하는 동안 결과는 다음과 같습니다.

800311.21   9422892.49 0.00