在Python 2.7中使用gdal.ReadAsArray时旋转xy

Question

我使用Python中的gdal模块来读取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（第一轴），“x”轴为轴“1”（第二轴），基于0的索引。当你要求ta[3,1]时，从左上角的原点（0,0）得到坐标x = 1，y = 3。没有轮换，这只是一个惯例问题。

您可以使用matplotlib绘制DEM或使用gdal将其另存为新光栅，您会看到它保留原始方向。