I want to get an XYZ-like surface from a Geotiff with GDAL in Python. Собственно, делаю я вот что:
1) Качаю Geotiff отсюда: http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp
2) Пишу вот такой код:
from osgeo import ogr, osr, gdal
file_path = os.path.join(os.path.dirname("/home/yak/dump/srtm_62_05.tif")
gdal_data = gdal.Open(file_path, gdal.GA_ReadOnly)
raster_band = gdal_data.GetRasterBand(1)
raster = raster_band.ReadAsArray()
Всё круто, только, во-первых, ни разу не понятно, в каких попугаях измеряется возвышение (написано, что в метрах, но относительно каких единиц? Да и в тех местах, где море, стоят какие-то страшные -32000 с гаком) — а во-вторых, мне бы это сконвертировать в обычные декартовы XYZ-координаты.
Собственно, имею спросить: а как бы это проделать-то? Явно же есть какие-то стандартные средства в самом GDAL?
