xy gedreht, wenn ich gdal.ReadAsArray in Python 2.7 verwende

Question

Ich habe das gdal-Modul in Python verwendet, um ein DEM-Raster zu lesen. Wenn unter gdal die DEM x und y Größen sind korrekt (42689, 35622). Als ich die ReadAsArray-Funktion verwendete, wird die Form des DEM gedreht (35622, 42689). Gibt es eine Möglichkeit, dies zu beheben oder zumindest festzustellen, wie das DEM gedreht wurde?

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 hat die "Y" Achse als Achse 0 (die erste Achse), und die "x" Achse als Achse "1" (die zweite Achse), mit einem 0-basierten Index. Wenn Sie nach ta[3,1] fragen, erhalten Sie die Koordinate x = 1, y = 3 von einem Ursprung (0,0) oben links. Es gab keine Rotation, es ist nur eine Frage der Konvention.

Sie können das DEM unter Verwendung von matplotlib plotten oder es als neues Raster unter Verwendung von gdal speichern und Sie werden sehen, dass es die ursprüngliche Ausrichtung beibehielt.