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

Question

我想将作为geoTiff文件提供的光栅图像重新投影到另一个坐标系中。地图是WGS84伪墨卡托地图

但是当我运行下面的代码时，我只得到一个白色图像作为我的输出光栅。我怎样才能解决这个问题

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")

Answer 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)