我如何在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")

这里是我要重新投影的.tiff文件的链接：dropbox.com/s/opb7bxmezjqniah/berlin.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)