model op• 数据种类:shapemodel × points → measurement
• 调用:import shapestats; shapestats.shape_mahalanobis(model, shape, align: 'bool' = True, cumulative: 'float' = 0.99, n_modes: 'int' = 0)(或 opsshapestat.get("shape_mahalanobis"))
> 该算子的说明尚无译文,以下照原文给出。
モデルから見てその形がどれだけ異常か。→ float。
`sqrt(sum(score_i^2 / variance_i))` を先頭の何本かに限って足す。
★ 全成分を足してはいけない。 K 個体の群から出る主成分は K-1 本だが、
後ろのほうは分散が数値的なゼロまで落ちる。実測(K=12、N=80、真のモードは
2 本)の分散: ``3.2e-02, 5.6e-03, 4.0e-05, 6.7e-13, 9.7e-14, 9.9e-17, …,
1.0e-31``。この裾で割ると値が意味を失う:
========== =========== ==================
使う本数 群内の個体 0.5 の膨らみを注入
========== =========== ==================
1 0.428 1.821
2 0.885 2.617
3 1.035 18.134
5 1.373 1393035.060
11(全部) 46613.636 9.3e+13
========== =========== ==================
5 本で群内の個体が 1.37 なのに外れが 100 万、11 本では群内ですら 46614。
「異常度」ではなく「数値ゼロで割った回数」を測っている。
そこで既定は**累積寄与率 `cumulative`(0.99)までの成分だけ**を使う。
上の群では 2 本(0.99895)が選ばれ、群内 0.885 / 外れ 2.617 と素直に並ぶ。
`n_modes` を正の数で与えればその本数に固定する(比較のため本数を揃えたい
ときに使う)。どちらを使ったかで値が変わるので、報告するときは本数も一緒に。
• 示例数据目录(下载 URL / 许可证) —— 2-D 用 skimage.data(BSD/公有领域)加合成图,3-D 给出真实数据源(Stanford/PDS 等)的下载 URL。
• 算子来历与参考文献 —— 该算子族所依据的研究/方法出处。
• 算法的正典(作者・年份)与用途见上面的族使用指南。
• (暂无)
measurement 作为输入)—
model)shape_pca · shape_project · shape_reconstruct · shape_explained_variance · shape_synthesize
*Provenance: shapestats.py — SHAPESTAT 算子登记表。本条目由 tools/opdocs.py md 自动生成(请勿手工编辑)。*
© 2026 Kazufumi Furuse — Fullseye operator documentation. Licensed under Apache-2.0.