使用GDAL在Python中平滑/插值栅格？

我正在用Python开发，并使用OSGEO的GDAL来处理栅格和shapefile并与之交互。

我想获取一个具有点特征的shapefile并将其插入到表面栅格中。现在，我使用的是'RasterizeLayer'方法，该方法将点要素中的值刻录到栅格中（使用所有nodata值进行设置），但将所有未触及的像素保留为'nodata'值。因此，我剩下的是棋盘式栅格。

使用RasterizeLayer之后，我得到了什么：

我想要的最终产品是：

我相信我要寻找的功能在arcgisscripting导入中称为'Spline_sa（）'。

GDAL是否具有类似的功能，或者有其他方法来获得所需的输出？

我看一下NumPy和Scipy-在SciPy Cookbook中有一个很好的示例，它使用scipy.interpolate.griddata函数进行插值。显然，这要求您将数据存储在numpy数组中；

• 使用GDAL python绑定，您可以gdal.Dataset.ReadAsArray()将栅格中的数据读入Python 。
• 使用OGR，您可以遍历要素图层并从shapefile中提取点数据（或者更好的是，使用GEOMETRY=AS_XYZ[请参见OGR CSV文件格式]将shapefile写入CSV，然后将csv读取到Python中）。

获得网格化输出后，可以使用GDAL将生成的numpy数组写入栅格。

最后，如果您对Scipy插值库没有任何运气，也可以随时尝试scipy.ndimage。

看一下GDAL网格API。我不知道如果在Python绑定暴露，但如果没有，你打电话叫gdal_grid通过公用子模块。

GDAL网格API仅使用反距离权重，移动平均和最近邻，而不实现样条线。另一种选择是使用Scipy。

该线程有些陈旧，但我编写了一个简单的模块，该模块使用了来自sklearn的称为skspatial的KNN算法。

https://github.com/rosskush/skspatial

您可以使用geopandas导入shapefile并选择一列，它将对可以导出到栅格的曲面进行插值。这是非常基本的方法，可能不是最好的方法，但是它至少使所有内容保持纯python。