(半)楕円形に似たクラスターに配置された 3D ポイントのデータセットがあります。たとえば、MATLAB 関数http://www.mathworks.com/matlabcentral/fileexchange/24693-ellipsoid-fitに実装されている標準的な楕円体フィッティングを試みると、フィッティングをすべての点に適用しても、虚数半径が得られます楕円体ボリュームまたはその凸包上にあるポイントのみ。
クラスターの形状が完全な楕円体に似ているよりも半楕円体であるという事実が、これらの不安定性につながるのでしょうか? 私の問題を解決できるより堅牢な方法はありますか?
アップデート：
これは、そのような 3D ポイントのクラスターの 1 つの例です。青い円は、クラスターの凸包にあるポイントを示しています。
凸包の場合の最小限の作業例は次のとおりです。
k = convhull(x, y, z);
ind = unique(k(:));
[center,radii,evecs,~] = ellipsoid_fit([x(ind),y(ind),z(ind)], 1);