3次元点群から外れ値を除去

こんにちは。地理座標系（緯度[deg]、経度[deg]、高さ[m]）のデータから、scatter3で3次元図にしたデータがあります。

以下写真のように、赤丸で囲まれた外れ値を除去したいのですが、なにか良い方法はありますでしょうか？

以下の記事を見たのですが、地理座標系には対応していないようでした。

https://jp.mathworks.com/help/vision/ref/pcdenoise.html

よろしくお願いいたします。

Answer 1

[Hernia Baby]

外れ値がわかるのでしたらindexを活用するのはいかがでしょうか？

figure
[X,Y,Z] = sphere(16);
x = [0.5*X(:); 0.75*X(:); X(:)];
y = [0.5*Y(:); 0.75*Y(:); Y(:)];
z = [0.5*Z(:); 0.75*Z(:); Z(:)];
scatter3(x,y,z,'k.')
hold on

ここでx,y,zが正のものだけを取り出します

そして赤色で囲みます

idx = x > 0 & y >0 & z>0;
x1 = x(idx);
y1 = y(idx);
z1 = z(idx);
scatter3(x,y,z,'r')

赤色を外れ値として扱う場合、その部分を削除します

該当部分を青で囲みます

x2 = x; y2 = y; z2 = z;
x2(idx) = [];
y2(idx) = [];
z2(idx) = [];
scatter3(x2,y2,z2,'b')
xlim([-1 1])

[Yu]

コメントくださりありがとうございます。

説明不足で申し訳ありません。

ここでいう外れ値とは、特定の値が決まっているわけではなく、

クラスタリング（形状）に基づいて判断できるようにしたいです。

たとえば、

①クラスタを自動抽出

②クラスタの垂直長さが10m以上の場合に、抽出したいクラスタと定義

③それ以外のシグナルは除去

したいと思っています。

なにかよい方法をご存知でしたらご教示ください。

宜しくお願い致します。

[Hernia Baby]

決まってない場合でもクラスで分かれている時点でindexから抜き出す事が可能と思われます

こちら、やりたい事に近いと思います

https://jp.mathworks.com/matlabcentral/answers/524662-kmeans-of-3d-data-clustering

X = rand(500, 3);
idx = kmeans(X, 3);

scatter3(X(:,1), X(:,2), X(:,3), 10, idx, 'filled');

hold on
X(idx==2|idx==3,:) = [];
scatter3(X(:,1), X(:,2), X(:,3), 20, 'r');

[Yu]

なるほど、このようなやり方があるのですね。

とても参考になりました！ありがとうございました。