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

Question

我想將作為 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")

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)