raster - gdal_grid 线性产生奇怪的结果

Question

我有一个来自 LiDAR las 文件的大点 shapefile(xyz，大约 65,000 个点)，我正在尝试使用 gdal_grid 将其插值到网格上:

gdal_grid -ot Float64 
          -txe  422306.5970  422343.9970 
          -tye 4037022.9899 4036967.3399 
          -outsize 747 1112
          -a linear:radius=0:nodata=0
          in.shp out.tif

这会无误地运行并生成如下图所示的 map 。您会注意到三角形图案，好像大多数点都被忽略了。这个有趣的图像中的值在它们应该是的 <1 以内，所以 gdal_grid 似乎在适本地读取 z 字段，它似乎只是缺少大多数点。如果我尝试 invdist 或 average，保持所有其他不变，问题就会消失，网格看起来应该是这样(参见具有相同色阶的第二张图像)。在这个例子中，变化在 1091 和 1093 之间。我尝试缩放 Z 以使 shapefile 中的变化更大，但仍然发现了同样的问题。我也尝试了 -z_multiply 和 -z_increase 无济于事。不幸的是，我需要双线性插值，所以我陷入了僵局。有什么想法吗？

我仅通过将插值方法更改为 invdist 得到以下输出:

gdal_grid -ot Float64 
          -txe  422306.5970  422343.9970 
          -tye 4037022.9899 4036967.3399 
          -outsize 747 1112
          -a invdist:radius1=1:radius2=1
          in.shp out.tif

Answer 1

抱歉来晚了，但我遇到了同样的问题，我找不到任何解释。然后我注意到 OP 和我的数据集原始参数都非常高。

我对Delaunay算法的这个具体实现一无所知，但我的猜测是对于X,Y的微小变化(与每个点的绝对空间位置相比较小，不确定比例是多少)被忽略并在更大的三角形下表示，从而失去细节和精度。

您可以清楚地体会到使用绝对值点云和相对值点云进行插值之间的区别:

我通过从原始点移除固定偏移、插值然后添加相同的偏移来解决。

希望对您有所帮助。