Построение .asc с правыми географическими координатами ⇐ Python

[Anonymous]

Я пытаюсь построить некоторые данные карты, к которым я имею доступ к «Ascii Grid и GML». Ось X/Y. Тем не менее, «по умолчанию» - это некоторые значения, которые я не совсем понимаю. < /P>
Код ниже: < /p>

Код: Выделить всё

# join ASCII tiles together to get area of interest
tile_list = glob.glob("data/01/input/OS/NY/NY[0-9][0-9].asc")
vrt = gdal.BuildVRT("merged.vrt", tile_list)
gdal.Translate("mergeDEMvrt.tif", vrt)
vrt = None

grid = Grid.from_raster('mergeDEMvrt.tif', data_name='dem')
dem = grid.read_raster('mergeDEMvrt.tif')

fig, ax = plt.subplots(figsize=(12,10))
fig.patch.set_alpha(0)

plt.imshow(dem, extent=grid.extent, cmap='terrain', zorder=1)
plt.colorbar(label='Elevation (m)')
plt.grid(zorder=0)
plt.title('Digital elevation map', size=14)
plt.xlabel('Longitude')
plt.ylabel('Latitude')
plt.tight_layout()

Вышеприведенное изображение показывает график с форматом оси «неправильной». Мне понадобятся значения Long/LAT USUSAL. Я попробовал еще один файл .tif и координаты GEO, протянутые совершенно нормально. Другими доступными типами файлов являются GML, GeoPackage или Vector Tiles.
Спасибо за помощь.

Подробнее здесь: https://stackoverflow.com/questions/774 ... oordinates