Как мне перепроектировать растровое изображение в формате geoTiff в python? - PullRequest
Новая
       116

Как мне перепроектировать растровое изображение в формате geoTiff в python?

0 голосов
/ 09 января 2020

Я хочу перепроектировать растровое изображение, которое дается в виде файла geoTiff, в другую систему координат. Карта находится в WGS84-Pseudo-Mercator.

Но когда я запускаю следующий код, я получаю только белое изображение в качестве выходного растра. Как это исправить? 

import rioxarray

rds = rioxarray.open_rasterio("path_to_raster.tif")
crs = "EPSG:4978" # this depends on the exact projection you want to use
projected = rds.rio.reproject(crs)
projected.rio.to_raster("path_to_enu_raster.tif")

Вот ссылка на файл .tiff, который я хочу перепроектировать:
dropbox.com / s / opb7bxmezjqniah / berlin.tif

1 Ответ

1 голос
/ 09 января 2020

Ну, я нашел решение методом проб и ошибок. Этот код работает для Германии: 

from osgeo import gdal

filename = "berlin.tif"
input_raster = gdal.Open(filename)
for epsg in ["EPSG:4839", "EPSG:3068", "EPSG:25833"]:
    gdal.Warp('output_raster' + epsg + '.tif', input_raster, dstSRS=epsg)
...