如何在Python 3中使用gdal.Grid设置Raster的分辨率?

问题描述 投票:0回答:1

如何设置gdal.Grid中生成的栅格的输出分辨率?我使用下面的代码,它生成一个具有非常大的单元格 (~70m) 的 TIF,但是 vrt (filepath_vrt_in) 中引用的数据具有亚米分辨率的数据。当使用反距离加权算法和 x,y 半径为 1m 的最近邻时,我得到相同的结果。

from osgeo import gdal
temp_wd = gdal.Grid("C:/Temp/Python-Output-WD-01.TIF", filepath_vrt_in, zfield="Z", algorithm = "nearest:radius1=5.0:radius2=5.0:angle=0.0:nodata=-9999")
temp_wd = None

我还尝试在输出 TIF 上使用 gdal.Warp,如下所示,它会生成类似分辨率的 TIF。

temp_warpwd = gdal.Warp("C:/Temp/Python-Output-WD-WARPED-01.TIF", "C:/Temp/Python-Output-WD-01.TIF", xRes=5, yRes=5)
temp_warpwd = None
python raster gdal gdal-python-bindings
1个回答
0
投票

我有同样的问题并找到了这个解决方案:

要创建例如 100 x 200 像素的 GTiff,请写入:

from osgeo import gdal
raster= gdal.Grid(destName='output.tif',srcDS='input.shp',format='GTiff',options=gdal.GridOptions(width=100,height=200))
© www.soinside.com 2019 - 2024. All rights reserved.