ESRI 프로젝션 좌표를 LAT-LNG로 변환합니다

StackOverflow https://stackoverflow.com/questions/1150446

  •  18-09-2019
  •  | 
  •  

문제

"NAD 1983 StatePlane Michigan South Fips 2113 피트"에 X, Y 좌표의 큰 데이터 세트가 있습니다 (일명 ESRI 102690). LAT-LNG 지점으로 변환하고 싶습니다.

이론적으로 이것은 무언가입니다 proj 다루기 위해 만들어졌지만 문서는 나에게 단서를주지 못했습니다. 훨씬 더 복잡한 사례를 설명하는 것 같습니다.

나는 a를 사용해 보았다 파이썬 인터페이스, 그렇게 :

from pyproj import Proj
p = Proj(init='esri:102690')
sx = 13304147.06410000000 #sample points
sy = 288651.94040000000
x2, y2 = p(sx, sy, inverse=True)

그러나 그것은 매우 잘못된 출력을 제공합니다.

JavaScript 라이브러리가 있지만 처리해야 할 ~ 50,000 포인트가 있으므로 적절하지 않은 것 같습니다.

나를 위해 일한 것 :

나는라는 파일을 만들었다 ptest 각 쌍의 자체 라인에 X 및 Y 좌표는 공간으로 분리됩니다.

13304147.06410000000 288651.94040000000
...

그런 다음 해당 파일을 명령에 공급하고 결과를 출력 파일로 파이프했습니다.

$>cs2cs -f %.16f +proj=lcc +lat_1=42.1 +lat_2=43.66666666666666 
+lat_0=41.5 +lon_0=-84.36666666666666 +x_0=4000000 +y_0=0 +ellps=GRS80 
+datum=NAD83 +to_meter=0.3048006096012192 +no_defs +zone=20N +to 
+proj=latlon ptest > out.txt
도움이 되었습니까?

해결책

텍스트 파일에서 재판매 만하면서 일부 데이터 마이닝을 할 수있는 경우 원하는대로 사용하고 사용합니다. http://spatialreference.org/ref/esri/102690/ 참조로.

예를 들어 Proj4를 사용하여 쉘/CMD 파일에 저장하고 Proj4 (Linux/Windows 버전 사용 가능한)로 입력 파일을 호출하여 데이터 세트 크기에 문제가 없습니다. 

cs2cs +proj=latlong +datum=NAD83 +to +proj=utm +zone=10  +datum=NAD27 -r <<EOF
cs2cs -f %.16f +proj=utm +zone=20N +to +proj=latlon - | awk '{print $1 " " $2}

따라서 귀하의 경우에는 다음과 같은 것입니다.

cs2cs -f %.16f +proj=lcc +lat_1=42.1 +lat_2=43.66666666666666 +lat_0=41.5 +lon_0=-84.36666666666666 +x_0=4000000 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs +zone=20N +to +proj=latlon

http://trac.osgeo.org/proj/wiki/man_cs2cs

http://trac.osgeo.org/proj/

다른 팁

TXT, CSV 또는 XLS 파일에 좌표가있는 경우 Ctrl+C를 수행하여 다음에 삽입 할 수 있습니다. http://cs2cs.mygeodata.eu 적절한 입력 및 원하는 출력 좌표 시스템을 설정할 수 있습니다. 다양한 형식으로 수천 개의 좌표를 삽입 할 수 있습니다 ...

라이센스 : CC-BY-SA ~와 함께 속성
제휴하지 않습니다 StackOverflow
scroll top