主成分分析やRANSACを利用した点群からの円柱フィッティングの方法

本スライドは2026年4月8日に開催されたオンライン勉強会（AIMTG）にて利用されました。

イベントページやYoutube動画は以下の通りです。

イベントページ
https://osaka-driven-dev.connpass.com/event/387387/

動画
https://www.youtube.com/live/9LBPwhveGAc?si=PL16bOZpyiAQceEz&t=2671

参考論文
Katayama, N., Lin, C., Inadomi, S., Itakura, K., Chen, Y., & Chun, P. J. (2025). Bridge assessment after 2024 Noto Peninsula earthquake and rainfall using point cloud. Proceedings of the Institution of Civil Engineers-Forensic Engineering, 1-10.

Transcript

1. 3次元点群からの円柱フィッティング方法について ImVisionLabs株式会社代表取締役 板倉健太 博士（農学）
2. 2 背景と目的 橋脚やパイプなどの円柱形状の構造物は、幾何モデルとして円柱にフィッティングする ことで、寸法計測や変状解析が可能 以下の図は入力の3次元点群データと橋脚を自動的に検出したときの結果を示す 画像は静岡県から公開されているオープンデータを解析した ときの結果です
3. 3 背景と目的 洪水や地震の前後で橋脚や電柱などの構造物の点群に対して円柱の検出やフィッティ ングを適用することで、傾きや変形といった幾何的な変化を把握することができる Katayama, N., Lin, C., Inadomi, S., Itakura, K., Chen, Y., & Chun, P. J. (2025). Bridge assessment after 2024 Noto Peninsula earthquake and rainfall using point cloud. Proceedings of the Institution of Civil Engineers-Forensic Engineering. 以下の論文では地震や洪水後の橋脚の3次元点群データから傾きやその変化を計算
4. 4 フィッティング対象の点群の選択 点群全体から関心領域（ROI：Region of Interest）を定義し、円柱状の物体を抽出 この際、対象の円柱状の物体以外の点はできるだけ取り除かれていることが望ましい 検出したい円柱を含む点群データ 自動的に円柱状の物体を検出した時の様子 ※上図では、地表面を除去し、かつ物体の塊（クラスタ） ごとに分けたのちに円柱フィッティングを行っている
5. 円柱の式について  変数の意味 ・ は円柱の 軸上の1点（円柱の中心軸を定義する基準点） ・ は円柱の 軸方向ベクトル ・ は円柱上の点  方程式の構造 ・ : : 柱軸上の点 から の距離 ・ 𝑥𝑥0 , 𝑦𝑦0 , 𝑧𝑧0 𝑥𝑥0 , 𝑦𝑦0 , 𝑧𝑧0 (𝑎𝑎, 𝑏𝑏, 𝑐𝑐) (𝑥𝑥, 𝑦𝑦, 𝑧𝑧) (𝑥𝑥, 𝑦𝑦, 𝑧𝑧) 円柱の軸上の点と円柱上の点にて三平方の定理を行っているイメージ 𝑥𝑥0 , 𝑦𝑦0 , 𝑧𝑧0 (𝑥𝑥, 𝑦𝑦, 𝑧𝑧)・ ・ ・ r 円柱の式： : と軸方向ベクトル の内積(軸に平行な成分) (𝑎𝑎, 𝑏𝑏, 𝑐𝑐) (𝑥𝑥 − 𝑥𝑥0 , 𝑦𝑦 − 𝑦𝑦0 , 𝑧𝑧 − 𝑧𝑧0 ) 5 ※図は有限長の円柱として示しているが、方程式が表すのは 無限円柱である
6. 円柱の式について  変数の意味 ・ は円柱の 軸上の1点（円柱の中心軸を定義する基準点 𝑂𝑂） ・ は円柱の 軸方向ベクトル ・ は円柱上の点𝐵𝐵  方程式の構造 ・ : : 柱軸上の点 から の距離 ・ を軸方向に射影した長さの二乗（図の 𝑂𝑂𝐻𝐻2） 𝑥𝑥0 , 𝑦𝑦0 , 𝑧𝑧0 𝑥𝑥0 , 𝑦𝑦0 , 𝑧𝑧0 (𝑎𝑎, 𝑏𝑏, 𝑐𝑐) (𝑥𝑥, 𝑦𝑦, 𝑧𝑧) (𝑥𝑥, 𝑦𝑦, 𝑧𝑧) 円柱の軸上の点と円柱上の点にて三平方の定理を行っているイメージ O 円柱の式： ：ベクトル(𝑥𝑥 − 𝑥𝑥0 , 𝑦𝑦 − 𝑦𝑦0 , 𝑧𝑧 − 𝑧𝑧0 ) 6 H
7. 9 円柱フィッティングの手順の例 まずは主成分分析により軸のベクトルを推定し、その後に半径を決定する（次スライド） 円柱のフィッティングには軸の向き、軸の位置、半径が必要 半径の決定にはRANSACを利用する（後述） 軸の位置を 推定 半径を 決定 ※本スライドでは円柱フィッティングを行う方法の一例を紹介
8. 10 円柱フィッティングの手順 まずは主成分分析により軸のベクトルを推定し、その後に半径を決定する（次スライド） 円柱のフィッティングには軸の向き、軸の位置、半径が必要 半径の決定にはRANSACを利用する（後述） 軸の位置を 推定 半径を 決定
9. 12 初期軸方向の推定 抽出された点群に対して主成分分析（PCA：Principal Component Analysis）を 実行し、最も分散の大きい方向を軸方向ベクトルの初期推定とする。 軸を推定
10. RANSACによる半径の推定 [a] ScanX2.0  RANSACとはデータの中からノイズによる影響を最小限にしながら、モデルのパラ メータを推定するために使用される統計的手法 13  初期軸が求まった後に、次はRANSAC (RANdom SAmple Consensus）によって半 径を推定する RANSACで 半径を決定 PCAで軸が決定している 円柱の形が決定
11. 14 RANSACによる半径の推定 1.ランダムに3点を抽出し、半径を仮定 2.各点から軸への距離と半径との差から残差を計算 3.残差が小さい点をインライアとし、最も多くのインライアを持つモデルを選定 4.参照軸（例：Z軸）との方向一致条件も導入可能 半径を仮定 インライアを推定
12. 15 フィッティングの実行結果 図1:元データ 図2:フィッティング結果 正確に円柱を抽出できていることがわかる 図2(フィッティング結果)は図1(元データ)に対 してPCAとRANSACによって抽出された円柱を 重ねたものである
13. 16 フィッティングの実行結果 回転させると右の動画のようになった どの方向から観察しても正確に円柱を抽出で きていることが分かった
14. 17 樹木のフィッティング  樹木の幹を対象として円柱フィッティングを行うと直径が0.74m (図1) 、CloudCompare 上で2点を選択し、直径を測る(図2)と0.73mとなった 円柱フィッティングにより直径を推定 手動で2点を選択し、直径を測定