如何在 python 中重新投影 geoTiff 格式的光栅图像？
How do I reproject a raster image in geoTiff format in python?
我想将作为 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")
好吧，我通过反复试验找到了解决方案。此代码适用于德国：
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)
我想将作为 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")
好吧，我通过反复试验找到了解决方案。此代码适用于德国：
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)