如何在 python 中重新投影 geoTiff 格式的光栅图像？答案

【问题标题】：How do I reproject a raster image in geoTiff format in python?如何在 python 中重新投影 geoTiff 格式的光栅图像？
【发布时间】：2020-01-09 15:03:23
【问题描述】：

我想将作为 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")

【问题讨论】：

标签： python gis raster gdal geotiff

【解决方案1】：

我通过反复试验找到了解决方案。此代码适用于德国：

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)

【讨论】：