使用 gdal 重新投影 .tiff 文件

Question

我需要将 .tiff 文件从 EPSG:4326 - WGS 84 - Geographic 重新投影到 EPSG: 32638（UTM Zone 38N）。当我运行下面的代码时，输​​出文件是一个重新投影的 .tiff，大小是原始文件的一半。

# Reprojection

from osgeo import gdal

filename = r"/home/DATA/Test.tif"

input_raster = gdal.Open(filename)
output_raster = r"/home/DATA/Test_UTM38N.tif"

gdal.Warp(output_raster,input_raster,dstSRS="+init=epsg:32638")

Answer 1

我在使用命令行版本的 gdalwarp 时也有过类似的经历。指定 x 和 y 分辨率解决了这个问题。即，在我的情况下，将“-tr 1000 1000”添加到 gdalwarp 命令会产生一个与原始大小大致相同的重新投影输出网格。

另请注意，'-tr' 选项 (https://gdal.org/programs/gdalwarp.html#cmdoption-gdalwarp-tr) 的 gdalwarp 文档说：

如果未指定（或未从 -te 和 -ts 推断），gdalwarp 将生成 xres=yres 的输出栅格，即使在不涉及重投影的场景中使用 gdalwarp 也是如此。