MATLAB 使用 Delaunay 三角剖分查找点云中的异常值

如何解决MATLAB 使用 Delaunay 三角剖分查找点云中的异常值

我有代表 LIDAR 测量的 3D 点云，我想删除或插入异常值。

MWE 创建了一个带有一些异常值的平滑表面。然而，实际测量可以被认为是分段平滑的，因为点云代表了环境的“外壳”。因此，在 x 和 y 方向上进行简单的平滑操作并没有产生好的结果。

MWE

rng(42);
[X,Y] = meshgrid(linspace(-1,1,31),linspace(-1,31));

Z = peaks(X,Y);

T = delaunay(X,Y);

% create noisy samples
noise_idx = randi(numel(Z_noisy),20,1);
Z_noisy = Z;
Z_noisy(noise_idx) = Z_noisy(noise_idx)+3*sign(rand(size(noise_idx))-0.5);

% visualization
trisurf(T,X,Y,Z_noisy,'FaceColor','interp');
view([0 90]);

我有三个独立的问题，回答其中任何一个都会非常有帮助。

1. 如何估计异常值的指数？ （MWE 中的 noise_idx）
2. 如果异常值的索引已知，我是否可以更改 trisurf 参数以便不绘制异常值？我想阻止创建另一个 delaunay 三角剖分。
3. 如果异常值的索引是已知的，我如何获取连接点的索引以插入与异常值对应的 Z 值而不是删除异常值？

可视化