xy поворачивается при использовании gdal.ReadAsArray в Python 2.7

Question

Я использовал модуль gdal в Python для чтения растра DEM. Когда ниже gdal размеры DEM x и y правильны (42689, 35622). Поскольку я использовал функцию ReadAsArray, форма DEM поворачивается так, чтобы она была (35622, 42689). Есть ли способ исправить это или хотя бы определить, как DEM был повернут?

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) 

Answer 1

Numpy имеет «Y» ось в качестве оси 0 (первая ось), а «х» ось в качестве оси (второй оси) «1», с индексом 0 на основе. Когда вы запрашиваете ta[3,1], вы получаете координату x = 1, y = 3 от начала координат (0,0) в левом верхнем углу. Вращений нет, это просто вопрос конвенции.

Вы можете построить DEM с помощью matplotlib или сохранить его как новый растровый, используя gdal, и вы увидите, что он сохранил исходную ориентацию.