问题说明:
我目前正在使用以下方法制作光栅文件:
sc1 = str.format('gdal_grid -a average:radiu1=30:radius2=30:min_points=3:nodata=-1 -ot Float32 -of GTiff -txe {0} {1} -tye {2} {3} -outsize {4} {5} -zfield "Slews" -l MFT_data ../MFT_data.shp Slews.tif',
xms.min(), xms.max(), yms.min(), yms.max(), np.round(xdist/2.), np.round(ydist/20.))
并使用os.system(sc1)
调用它。
我正在尝试直接在Python中实现它并使用GDAL/OGR Python API。
我的尝试:
gdal.Grid('Slews.tif', 'MFT_data.shp',format = 'GTiff',outputType = 'Float32', outputBounds = [xms.max(), yms.max(), xms.min(), yms.min()], width = np.round(xdist/2.), np.round(ydist/20.]',algorithm = 'average:radius1=30.0:radius2=30.0:min_points=3:nodata=-1',zfield ='Slews', layers = 'MFT_data')
我的尝试是否直接等同于我之前使用的方法?运行代码行并不适合我。我已经尝试搜索使用了许多网格选项的示例,但找不到任何网格选项。期待您的建议/帮助/提示。