最好的方式插入深度测量法与定向组件?

Question

我有一小部分失踪的水深数据集。我需要填写这个部分。

这个地区的深度测量法有一些强烈定向的特性(例如一个峡谷,穿过失踪的地区)。使用一个标准的插值函数会削弱这些特性。

是否有一个适当的方法来估计缺失值,让更多的体重在一个给定的已知值轴?

目前数据GeoTIFF,我用QGIS、但我可以根据需要进行转换。我有简单的MATLAB和R,以及所有的草和传奇工具包含在QGIS。

Answer 1

我有同样的问题与民主党的表面。最好的解决方案,我发现在Matlab使用scatteredInterpolant类,它在Matlab是内置的。它允许自然的邻居插值(这是一个类的加权距离插值建议在之前的评论)。用法是这样的:

%如果X和Y是每个网格单元的坐标(你可以使用meshgrid %) invalidValue =南;%或-99999之类的是在你的情况中gapIdxs =民主党= = invalidValue;%的索引无效细胞interpolant = scatteredInterpolant (X (~ gapIdxs), Y (~ gapIdxs),民主党(~ gapIdxs),“自然”);% interpolant现在是一个函数,我们可以评估它在任何时候你想要的,在这种情况下,我们评估的空白细胞interpolatedValues = interpolant (X (gapIdxs), Y (gapIdxs));%现在我们补丁民主党patchedDEM =民主党;patchedDEM (gapIdxs) = interpolatedValues;

产生结果的方式比立方插值以我的经验。如果这还不够,在你开始做手工采摘的特性可以探索克里格插值。最后如果你迫使方向性,可以旋转你的数据(使用matlab imrotate),这样你的主要方向是垂直的,然后插入列的列(使用matlab interp1)。