2016-07-12 7 views
0

파이썬에서 gdal 모듈을 사용하여 DEM 래스터를 읽었습니다. gdal 이하인 경우 DEM x 및 y 크기가 정확합니다 (42689, 35622). ReadAsArray 함수를 사용하면 DEM의 모양이 (35622, 42689)로 회전됩니다. 이 문제를 해결하거나 최소한 DEM이 회전 된 방법을 식별 할 수있는 방법이 있습니까?Python 2.7에서 gdal.ReadAsArray를 사용할 때 xy가 순환되었습니다.

In [54]: t1.RasterYSize 
Out[54]: 35622 

In [55]: t1.RasterXSize 
Out[55]: 42689 

ta = t1.ReadAsArray 

In [64]: ta.shape 
Out[64]: (35622L, 42689L) 

답변

0

NumPy와는 0 기반 인덱스 축 0로서 "Y"축 (제 1 축)을 축 "1"(두 번째 축)와 "X"축을 갖는다. ta[3,1]을 요청하면 왼쪽 상단의 원점 (0,0)에서 좌표 x = 1, y = 3을 얻습니다. 로테이션은 없었으며 단지 협상의 문제 일뿐입니다.

matplotlib을 사용하여 DEM을 플롯하거나 gdal을 사용하여 새로운 래스터로 저장하면 원래 방향으로 유지됩니다.