(57)【要約】【課題】被計測表面積が狭い場合でも、計測ノイズに頑健に柱状構造物のたわみ量を精度よく推定する。【解決手段】柱状構造物の中心軸から所定距離以内の３次元点群を入力し（S331）、各繰り返し処理で、３次元点群をサンプリングし（S332）、サンプリングした複数の点群に基づいて、柱状構造物の形状を規定するたわみ付き第２の円筒モデルの当てはめをし（S334）、柱状構造物のたわみ量を推定し（S334）、繰り返し処理を一定回数（Ｄloop)実行した場合（S335；Y）、サンプリングした複数の点のＤloop個の組合せの各々について推定されたたわみ量に基づいて、たわみ量の期待値を算出する（S336）。よって、たわみ付円筒モデルが、計測点にオーバーフィッティングすることを抑制でき、たわみ判定を精度よく行うことが可能である。【選択図】図１９
2015224980.tif 000002
【特許請求の範囲】
【請求項１】
柱状構造物の位置に対応する範囲内に位置する３次元点群であって、かつ、物体の表面上の位置を計測した計測結果である複数の位置を表す３次元点群から選択された複数の３次元点の組合せの各々について、前記組合せの複数の３次元点に基づいて、前記柱状構造物のたわみ量を推定し、前記複数の組合せの各々について推定されたたわみ量に基づいて、たわみ量の期待値を算出するたわみ期待値推定部
を含むたわみ推定装置。
【請求項２】
前記３次元点群に基づいて、前記柱状構造物にたわみがないと仮定した場合の前記柱状構造物の形状を規定する第１の円筒モデルのモデルパラメータを推定し、前記３次元点群と、前記推定された第１の円筒モデルのモデルパラメータとに基づいて、前記第１の円筒モデルのモデルパラメータにより形状が表される前記柱状構造物の表面に前記３次元点群が存在する度合いを示す第１のスコア値を計算し、
前記３次元点群に基づいて、前記柱状構造物にたわみがあると仮定した場合の前記柱状構造物の形状を規定する第２の円筒モデルのモデルパラメータを推定し、前記３次元点群と、前記推定された第２の円筒モデルのモデルパラメータとに基づいて、前記第２の円筒モデルのモデルパラメータにより形状が表される前記柱状構造物の表面に前記３次元点群が存在する度合いを示す第２のスコア値とを比較し、比較した結果に基づいて、前記柱状構造物にたわみがあるか否かを判断する円筒モデル比較部を更に含み、
前記たわみ期待値推定部は、前記円筒モデル比較部により、前記柱状構造物にたわみがあると判断された場合に、前記期待値を算出する
ことを特徴とする請求項１に記載のたわみ推定装置。
【請求項３】
柱状構造物の位置に対応する範囲内に位置する３次元点群であって、かつ、物体の表面上の位置を計測した計測結果である複数の位置を表す３次元点群を取得する取得部と、
前記柱状構造物の位置に対応する前記範囲に、各々前記柱状構造物の中心軸方向の隣と一部重複するように設定された複数の局所領域を設定する局所領域形状解析部と、
前記複数の局所領域の各々について、前記３次元点群から、前記局所領域に含まれる所定個の３次元点を選択し、前記複数の局所領域の各々について選択された所定個の３次元点に基づいて、前記柱状構造物のたわみ量を推定することを繰り返し、
繰り返し毎に推定されたたわみ量に基づいて、たわみ量の期待値を算出するたわみ期待値推定部と、
を含むたわみ推定装置。
【請求項４】
円筒モデル比較部を更に含み、
前記局所領域形状解析部は、前記柱状構造物の位置に対応する前記範囲に、前記複数の局所領域を設定し、前記複数の局所領域の各々について、前記３次元点群のうちの前記局所領域に含まれる３次元点に基づいて、前記局所領域における円筒形状の中心位置及び半径を含む円筒パラメータを推定し、前記複数の局所領域の各々について推定された前記円筒パラメータに基づいて、前記柱状構造物はテーパ形状を有するか否かを判断し、
前記円筒モデル比較部は、
前記３次元点群と、前記局所領域形状解析部によって判断されたテーパ形状を有するか否かの判断結果とに基づいて、前記柱状構造物にたわみがないと仮定した場合の前記柱状構造物の形状を規定する第１の円筒モデルのモデルパラメータを推定し、前記３次元点群と、前記推定された第１の円筒モデルのモデルパラメータとに基づいて、前記第１の円筒モデルのモデルパラメータにより形状が表される前記柱状構造物の表面に前記３次元点群が存在する度合いを示す第１のスコア値を計算し、
前記３次元点群と、前記局所領域形状解析部によって判断されたテーパ形状を有するか否かの判断結果とに基づいて、前記柱状構造物にたわみがあると仮定した場合の前記柱状構造物の形状を規定する第２の円筒モデルのモデルパラメータを推定し、前記３次元点群と、前記推定された第２の円筒モデルのモデルパラメータとに基づいて、前記第２の円筒モデルのモデルパラメータにより形状が表される前記柱状構造物の表面に前記３次元点群が存在する度合いを示す第２のスコア値とを比較し、比較した結果に基づいて、前記柱状構造物にたわみがあるか否かを判断し、
前記たわみ期待値推定部は、前記円筒モデル比較部により、前記柱状構造物にたわみがあると判断された場合に、前記期待値を算出する
ことを特徴とする請求項３に記載のたわみ推定装置。
【請求項５】
前記たわみ期待値推定部は、前記柱状構造物の上端と前記柱状構造物のたわみ開始位置における中心軸の接線との最短距離を、前記たわみ量として推定することを特徴とする請求項１〜請求項４の何れか１項に記載のたわみ推定装置。
【請求項６】
前記たわみ期待値推定部は、前記柱状構造物の上端と下端とを通る直線と、前記柱状構造物の中心軸との最大距離を、前記たわみ量として推定することを特徴とする請求項１〜請求項４の何れか１項に記載のたわみ推定装置。
【請求項７】
前記たわみ期待値推定部によって推定されたたわみ量の各々と、前記算出されたたわみ量の期待値とに基づいて、信頼度を算出する信頼度算出部を更に含む請求項１〜請求項６の何れか１項に記載のたわみ推定装置。
【請求項８】
請求項１〜請求項７の何れか１項記載のたわみ推定装置の各部としてコンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、たわみ推定装置、及びプログラムに係り、特に、柱状構造物のたわみ量を推定するたわみ推定装置、及びプログラムに関する。
【背景技術】
【０００２】
近年、トンネルや橋などのインフラ構造物の劣化・老朽化が問題となり、自治体において維持管理やリスク評価が課題となっている。それと同時に、インフラ設備の数の多さが問題となり、構造物の安全性の評価の効率化やコスト削減を目的とした技術開発が求められている。特に、電柱や信号などのコンクリートポールのような柱状構造物の点検業務は、その数の多さから全数点検作業の稼働コストが大きい。
【０００３】
構造物のひずみ等の形状変形に関する情報（以下、たわみと呼ぶ）は破壊につながる重要な安全指標の一つであり、従来はひずみゲージとよばれるセンサを構造物表面にとりつけることで計測を行ってきた。近年では、コストの問題やセンサ事態の耐久性を考慮して、カメラやレーザー等の光学的なセンサの開発が進み非接触な計測が可能となった。これにより、ひずみゲージ等のセンサ設置のための仮設足場などが不要となるため、短時間で複数の位置の情報を計測することが可能となった（非特許文献１）。また、柱状の構造物については、非特許文献２のような専用の製品活用により、より効率的な点検が可能となった。これにより、ある基準以上のたわみ量を有するものを不安全と判定し、その柱状構造物について保守点検を優先して行うことができる。
【０００４】
ただし、上記の非特許文献１においては、測定場所に基準点を設定して計測するときに、人手の作業が依然として必要である。非特許文献２のような専用装置を用いても、計測箇所の指定や柱状物体の中心軸位置を決定するためには人手による指示作業が必要であり、膨大な数の柱状構造物の点検作業における稼働コストは大きい。
【０００５】
一方、モバイルマッピングシステム（ＭＭＳ）と呼ばれるカメラやレーザースキャナを搭載した車が、街中を走行することで道路周辺の構造物表面の形状を計測できるシステムが普及して利用されつつある。このシステムは、GPS（全地球測位システム）やIMS(慣性計測装置)を用いて物体の表面を３次元の座標情報として記録できる。
【０００６】
また、近年では、非特許文献３のような技術により柱状構造物を検出する技術があり、この技術とＭＭＳを用いることで、屋外環境下での柱状構造物の位置を自動検出することは可能である。しかしながら、ＭＭＳで計測した点群からたわみ量は推定できないため、その柱状構造物が安全か不安全であるかの判定ができない。
【先行技術文献】
【非特許文献】
【０００７】
【非特許文献１】「光学的全視野計測法を用いた鋼・コンクリート部材の変形・ひずみ計測に関する基礎的研究」 : Nagasaki University (長崎大学), 博士(工学) (2012-09-20). 記述: 長崎大学学位論文
【非特許文献２】アイサンテクノロジ柱状物体計測ツール「バーム」、[online]、[平成２６年５月１９日検索]、インターネット＜ＵＲＬ：http://www.aisantec.co.jp/products-services/atstation/baumstation.html＞
【非特許文献３】新垣、島村、新井、谷口、「３次元点群からの局所形状と大域的な形状モデルを用いた柱状物体検出」、信学技報, vol. 112, no. 441, PRMU2012-131, pp. 7-12, 2013年2月.
【発明の概要】
【発明が解決しようとする課題】
【０００８】
従来、柱状構造物のたわみ計測をするための稼働時間が、設備の数の膨大さから問題視されてきた。現地での計測作業の時間を削減するために、点群データをもとに電柱や信号等の柱状構造物のたわみ形状の計算を、３次元ＣＡＤ装置などの専用編集ツールを利用することが考えられる。しかし、こういったＣＡＤツールでは、目視による中心軸の位置の指定作業が必須となり、手動による位置指定の精度が低い場合にはたわみ量の推定精度が低くなる。そのため、ＣＡＤ等を用いた計測には熟練の操作者が必要である。また、画面上での操作は現場への移動時間を減らせるが、画面上での細かな作業のため作業者に大きな負担をしいる。
【０００９】
例えば、図１に示すように、柱状構造物について、様々な高さの断面を作成し、その断面における中心位置を手動により指定することで、柱状構造物の中心軸の形状を推定することも可能であるが、非計測領域が狭い場合には中心位置を精度よく求めることは難しい。
【００１０】
特に、ＭＭＳのような移動する車に搭載されたセンサーシステムで屋外の構造物を計測した場合、原理的に車が走行する道路側の被写体の表面形状は測定できるが、道路反対側の表面形状は計測できないため、計測領域が狭い（表面積が小さい）ことが多い。一般的に、計算機（コンピュータアルゴリズム）によりたわみ量を自動的に推定する際にも、計測表面積が小さい場合には、中心軸の推定精度が低くなりやすい。なぜなら、非計測表面積が小さいときには、計測ノイズの影響が大きくなり、手動操作と同様に各断面における中心位置の推定精度が低くなるためである。たわみの量は微小であるため、例えば、不安全と判定される電柱でも、たわみ量が先端部でも20cmくらいであるため、中心位置の数センチの推定誤差の影響も大きい。つまり、従来技術の推定技術では、中心軸の位置推定精度の低いために、たわみ推定精度も低くなるので、柱状構造物の安全性・不安全性を判定することが難しい。
【００１１】
本発明は、このような事情を考慮してなされたものであり、その目的は、被計測表面積が狭い場合でも、計測ノイズに頑健に柱状構造物のたわみ量を精度よく推定することができるたわみ推定装置、及びプログラムを提供することにある。
【課題を解決するための手段】
【００１２】
上記の目的を達成するために本発明に係るたわみ推定装置は、柱状構造物の位置に対応する範囲内に位置する３次元点群であって、かつ、物体の表面上の位置を計測した計測結果である複数の位置を表す３次元点群から選択された複数の３次元点の組合せの各々（複数回サンプリングした点群）について、前記組合せの複数の３次元点に基づいて、前記柱状構造物のたわみ量を推定し、前記複数の組合せの各々について推定されたたわみ量に基づいて、たわみ量の期待値を算出するたわみ期待値推定部を含む。
【００１３】
本発明によれば、たわみ期待値推定部は、柱状構造物の位置に対応する範囲内に位置する３次元点群であって、かつ、物体の表面上の位置を計測した計測結果である複数の位置を表す３次元点群から選択された複数の３次元点の組合せの各々について、前記組合せの複数の３次元点に基づいて、前記柱状構造物のたわみ量を推定し、前記複数の組合せの各々について推定されたたわみ量に基づいて、たわみ量の期待値を算出する。
【００１４】
このように、本発明は、柱状構造物の位置に対応する範囲内に位置する３次元点群であって、かつ、物体の表面上の位置を計測した計測結果である複数の位置を表す３次元点群から選択された複数の３次元点の組合せの各々について推定されたたわみ量に基づいて、たわみ量の標準偏差に応じて、推定精度の信頼度を判定することができる。
【００１５】
本発明に係るプログラムは、上記たわみ推定装置としてコンピュータを機能させるためのプログラムである。
【発明の効果】
【００１６】
以上説明したように、本発明のたわみ推定装置、及びプログラムによれば、計測ノイズの影響を抑制して、たわみ量を精度よく推定することができることができる、という効果が得られる。
【図面の簡単な説明】
【００１７】
【図１】（Ａ）及び（Ｂ）は、計測点群と高さＺｈでの円筒当てはめ結果の概念図であり、（Ａ）は、測定する装置側から見た柱状構造物等の表面の点群を示し、（Ｂ）は、（Ａ）におけるＡ−Ａ断面図である。
【図２】（Ａ）〜（Ｃ）は、電柱のたわみの原因を示す概念図である。
【図３】（Ａ）は、直線上の点群にノイズを加えて曲線モデルを当てはめたときの結果の例を示し、（Ｂ）は、不安定と判定された柱状構造物の表面上の点群に曲線モデルを当てはめたときの結果の例を示す図である。
【図４】（Ａ）〜（Ｃ）は、低次元・高次元のモデルの選択の判定の説明図である。
【図５】（Ａ）及び（Ｂ）は、形状決定への計測ノイズの影響を説明し、（Ａ）は、計測ノイズが小さい場合を示し、（Ｂ）は、計測ノイズが大きい場合を示す説明図である。
【図６】（Ａ）は、たわみ量を求める第１の方法を示し、（Ｂ）は、たわみ量を求める第２の方法を説明する説明図である。
【図７】柱状構造物のたわみ推定装置１００の概要図である。
【図８】たわみ推定処理部１０５の処理の概要を示すフローチャートである。
【図９】図８のステップＳ３−１における局所領域形状解析部１２０が実行する局所領域形状解析処理の概要を示すフローチャートである。
【図１０】隣の領域と一部重複するように複数の局所領域を設定する方法を示す概念図である。
【図１１】図９のステップＳ３１３のＲＡＮＳＡＣによる円筒パラメータ登録処理の概要を示すフローチャートである。
【図１２】（Ａ）は、局所領域の登録情報を示す図であり、（Ｂ）は、局所領域ＲＳsにおける登録情報の概念図である。
【図１３】テーパの有無の判定方法を示す図である。
【図１４】図８のステップＳ３―２における円筒モデル比較部１２１が実行する円筒モデル比較処理の概要を示すフローチャートである。
【図１５】テーパ値αにおける半径の変化の説明図である。
【図１６】（Ａ）〜（Ｃ）は、柱状構造物の中心軸を表現する凸関数を求めるための方法を説明する図である。
【図１７】図１４のステップＳ３２３の円筒モデル比較部１２１が実行するＲＡＮＳＡＣによるたわみ円筒当てはめ処理の概要を示すフローチャートである。
【図１８】３次元点ｐｍと凸関数の概念図である。
【図１９】図８のステップＳ３−３におけるたわみ量期待値推定部１２２が実行するたわみ量期待値推定処理の概要を示すフローチャートである。
【図２０】（Ａ）及び（Ｂ）は、点群正規化（サンプリング）処理を説明する図である。
【図２１】推定したたわみ量の分布と信頼度との関係を示す概念図である。
【図２２】第２の実施形態におけるたわみ推定処理部１０５の処理を示すフローチャートである。
【図２３】ランダムサンプリングした点から推定するパラメータを示す図である。
【発明を実施するための形態】
【００１８】
以下、図面を参照して本発明の実施形態を詳細に説明する。
［発明の概要］
（第１の概要）
本発明では、３次元点群に対してたわみ付きの円筒モデルを当てはめることで、柱状構造物のたわみ量を推定する。ここで、当てはめるとは、最適な形状モデルパラメータを探索することであり、評価関数の値が高くなるような繰り返し処理を行い、最も評価値が良いときのモデルパラメータを出力することを意味する。
【００１９】
電柱がたわむ原因は、電力線や通信線などのケーブルの張力の影響や変圧器などの付属品の荷重である。図２に示すように、ケーブル等の張力および付属品の荷重等の合力により生じた柱状構造物への曲げモーメントによりたわみが生じる。
【００２０】
このとき、柱状構造物の中心軸は直線と仮定するよりも、図２に示すように、たわみが発生した位置に設定したローカル座標系（x’―y’座標）において凸関数、例えば、ｙ＝ｘ2(xの二乗)のような関数のような極値が１つだけ存在する関数、であると仮定した方が精度よく形状を表現できる。本発明では３次元の円筒モデルの中心軸を凸な関数で表現した３次元円筒モデル（以下、たわみ付き円筒モデル）を用いることで、形状解析を行い、たわみ量の自動推定を実現する。
【００２１】
（第２の概要）
しかしながら、歪んだ柱状構造物についても、つまりに不安全と判定されるような柱状構造物であっても、そのたわみ量は微小であるため、たわみ付き円筒モデルを当てはめた結果だけでは推定の信頼性がきわめて低い。例えば、１０ｍの柱状構造物に対してたわみ量が２０ｃｍ程度ぐらいである。以下 、柱状構造物に「たわみがある」と記載した場合には、上記のたわみ量以上の形状変形を有しているという意味とする。
【００２２】
そのため、中心軸を凸関数（例えば2次関数）で表現したときに、推定した係数の大きさだけでは判定が難しい。図３に示すように、真っ直ぐな円筒表面を計測して取得した点群について、その点群に当てはめた２次関数の曲線モデルの係数と、不安全と判定される柱状構造物の表面の点群に対して、その点群から求めた２次曲線モデルの係数について差が生じない。理由は、計測時に推定誤差（ノイズ）が生じるためである。
【００２３】
そこで本発明では、たわみの無い円筒モデル（以下、円筒モデルと記載する）とたわみ付き円筒モデルのスコア値と比較をして、たわみ円筒モデルのスコア値が大きいとき、計測ノイズではなくて、柱状構造物の形状がたわんでいることによりたわみが発生していると考える。これにより、計測ノイズや付属品等の影響により、誤ってたわみが発生している（不安全な設備である）と出力することを抑制できる。（図４参照）これは、機械学習や統計学の分野で、スコア値に差分がないときは過学習を避けるためになるべく単純なモデルを利用するような考え方と同じである。
【００２４】
（第３の概要）
更に、計測ノイズの影響が小さいと、柱状構造物の中心軸を表現する最適な１つの凸関数を求めやすい（図５（Ａ））。しかし、計測ノイズの影響が大きいとき、例えばＭＭＳで計測すると、柱状構造物の表面上の計測された面積が少ないため、計測点群がモデル当てはめをした際に、形状決定するのに十分な形状情報を有しておらず、ノイズの影響のある計測点群にオーバーフィッティング（過適応）して柱状構造物の中心軸を表現する凸関数が複数当てはまりやすい（図５（Ｂ））。この場合、仮にスコア値の差が大きいため、たわみ推定が存在することは確認できるが、そのたわみの量を正確に推定することが難しい状況がある。そこで、本発明では推定した結果が、信頼度の高い結果であるかも同時に推定する。
【００２５】
まず、点群のサンプリングによる複数の異なる条件で、例えば異なる計測点群を選択して、たわみ量を複数回推定し、その推定結果の期待値を算出する。これは機械学習や統計学の分野で用いられるクロスバリデーションの考え方と同じであり、計測ノイズにオーバーフィッティングすることを抑制する効果がある。
【００２６】
このとき、推定したたわみ量の分布状態から、推定したたわみ量の信頼度も算出し、上記のたわみ付円筒モデルのスコア値に優位性があり、かつたわみ量期待値が大きく信頼度が大きいときに、不安全な柱状構造物であると判定できる。この場合、設備管理者へ点検優先度の高いものであると示唆できる。
【００２７】
（第４の概要）
また、たわみが生じる箇所は、柱状構造物の部材の影響によるが、経験的に地面から数m高い位置のことが多い。そのため、コンクリート柱に特化して考えると、たわみが発生する箇所（以下、たわみ開始点）の位置を考慮して、円筒中心軸について、たわみ開始点よりもＺ軸の小さい側を直線モデル、たわみ開始点よりもＺ軸の大きい側を凸関数曲線モデルにするようなたわみ付き円筒モデルを用いることで、更に高精度にたわみ推定を行うことが可能である。
【００２８】
［実施形態の概説］
本発明は、電柱や信号などの柱状構造物をレーザースキャナで計測し、得られた３次元点群を用いてその柱状構造物のたわみ量を推定する技術である。
【００２９】
たわみ推定の前処理として、非特許文献３のような既存技術を用いて、柱状構造物を円筒モデルと見立てて検出し、位置と柱状構造物の中心軸の傾きを求める。次に、柱状構造物の中心軸周辺の点群を用いて、その構造物のたわみ量と信頼度を推定すればよい。
【００３０】
ここで、本発明における柱状構造物のたわみとは、外力によりまっすぐな柱状構造物が弓なりに形状を変化したことを意味する。具体的には、柱状構造物のたわみとは、図６に示すように、柱状物体の上端部の中心位置からたわみ開始位置Ｇ0における接線（柱状物体の傾きの方向）への垂線の足までの長さｄｈのこととする。もしくは、柱状構造物のたわみとは、上端部と下端部を結んだ直線と、凸関数の凸部（極値）の位置までの最短距離ｄｃのこととする。電柱点検における実際の現場では、上端部の移動量ｄｈを参考にするが、下端部近くに障害物、例えば車や木など、が存在するときにたわみ開始位置Ｇ0付近の点群が計測できないので接線方向の推定精度が低下しやすい。そのため、上端部移動量dhの大きさと相関の高いｄｃも参考指標として利用できる。
【００３１】
本発明は、電柱や信号、標識などの柱状構造物のたわみ量の推定および推定したたわみ量の信頼度算出を目的としているが、それ以外の柱状物体の形状推定も可能な技術である。
【００３２】
以下の実施形態では、具体的な例としてレーザーレンジファインダにより取得した点群を用いたたわみ推定方法について説明する。
【００３３】
ここで、３次元とは、緯度、経度、海抜（高さ）情報でもよいし、ユーザーが設定した特定の位置を原点とした３次元ユークリッド座標系でも極座標系でもよい。以下の例では、ユーザーが設定した原点における３次元ユークリッド座標系（各方向をＸ，Ｙ，Ｚ座標とする）を想定する。各座標の単位はメートル(m)やセンチメートル(ｃｍ)、ミリメートル(ｍｍ)で表現するが、他の単位でもよい。３次元点とは、各点に上記の３次元座標に、その点群が撮影された時刻や、レーザーの反射強度や赤・青・緑などの色情報等が付与されている点である。３次元点に付与される情報に制限はないが、少なくとも位置情報（Ｘ，Ｙ，Ｚ座標）は付与されたものであり、３次元点群とはその３次元点が2点以上集まった集合である。
【００３４】
[第１の実施の形態]
（たわみ推定装置全体の説明）
図７は、本発明の第１の実施形態による柱状構造物のたわみ推定装置１００の構成を示すブロック図である。たわみ推定装置１００は、ＣＰＵ（Central Processing Unit）、ＲＡＭ、２次記憶装置、プログラムを記憶するＲＯＭを備えたコンピュータで構成される。たわみ推定装置１００には、被写体計測部１０１、入力部１０２、及び出力部１０７が接続されている。
【００３５】
図７において、被写体計測部１０１は、レーザーレンジファインダや、赤外線センサ、または超音波センサなど、被写体とセンサとの距離を測定可能な装置であり、物体の表面上の位置を計測した計測結果である複数の位置を表す３次元点群を出力する。例えば、レーザーレンジファインダをＧＰＳ（Global Positioning System）が搭載された車の上、もしくはＧＰＳの搭載された飛行機に搭載し、移動しながら計測することで、屋外の環境の人工物、例えばワイヤ、建物、木、道路や道路以外の地面など不特定多数の被写体の３次元位置（点の座標）を計測するシステムである。本実施形態では、被写体計測部１０１として、車上にＧＰＳとレーザーレンジファインダとが搭載されているＭＭＳ（Mobile Mapping System）を想定している。ただし、被写体計測部１０１は、ある特定の位置（交差点など）１箇所からの計測部であってもよい。
【００３６】
入力部１０２は、マウスやキーボードなどのユーザーインターフェースであり、たわみ推定処理部１０５で使用するパラメータを入力するものである。また、パラメータを記憶したＵＳＢメモリなどの外部記憶媒体でもよく、記憶部１０３にパラメータを供給する。
【００３７】
たわみ推定装置１００は、上記２次記憶装置に設けられた記憶部１０３、柱状構造物検出部１０４、たわみ推定処理部１０５、及び構造物安全性判定部１０６を備えている。
【００３８】
記憶部１０３は、３次元点群記憶部１１０と演算処理用パラメータ記憶部１１１と柱状構造物パラメータ記憶部１１２とから構成される。３次元点群記憶部１１０は、被写体計測部１０１から取得した３次元点群を記憶し、柱状構造物検出部１０４とたわみ推定処理部１０５に供給する。演算処理用パラメータ記憶部１１１は、入力部１０２から取得したパラメータの値を記憶し、柱状構造物検出部１０４とたわみ推定処理部１０５に供給する。柱状構造物パラメータ記憶部１１２は、たわみ推定処理部１０５から取得した、たわみ状態パラメータ、スコア比、たわみ量の期待値、および信頼度を含む柱状構造物パラメータを記憶し、構造物安全性判定部１０６へ供給する。これら３つの手段により構成される記憶部１０３は、例えば、ＨＤＤ（Hard Disk Drive）や、ＳＳＤ（Solid State Drive）等のハードウェアによる記憶装置である。
【００３９】
柱状構造物検出部１０４は、３次元点群記憶部１１０から３次元点群を取得して、柱状構造物を検出し、柱状構造物毎に、柱状構造物の位置、中心軸の方向、半径、及び長さ等のパラメータを推定し、柱状構造物パラメータ記憶部１１２へ供給する。本発明では、非特許文献３のような既存技術を用いることとする。もしくは、既に設備管理ＤＢなどに登録されている場合には、入力部１０２から情報を得られるため、柱状構造物検出部１０４の処理を省くことも可能である。
【００４０】
たわみ推定処理部１０５は、入力処理部１１８、取得部１１９、局所領域形状解析部１２０、円筒モデル比較部１２１、たわみ期待値推定部１２２、信頼度算出部１２３、出力処理部１２４、及び繰り返し処理部１２５を備えている。たわみ推定処理部１０５は、記憶部１０３の３次元点群記憶部１１０から３次元点群を取得し、柱状構造物パラメータ記憶部１１２から柱状構造物パラメータを取得し、推定したたわみ量と信頼度を記憶部１０３へ供給する。
【００４１】
入力処理部１１８は、後述するパラメータを入力する。取得部１１９は、３次元点群と柱状構造物パラメータに基づいて、柱状構造物Iの中心軸から半径Ｒcut（=ＲI+ΔR）以内に含まれる３次元点群を取得する。
【００４２】
局所領域形状解析部１２０は、柱状構造物の中心軸周辺の３次元点群、及び柱状構造物パラメータを取得し、柱状構造物を含む領域を複数の領域に分割した局所領域を設定し、各局所領域の形状情報を表す円筒パラメータを推定して円筒モデル比較部１２１へ出力する。
【００４３】
円筒モデル比較部１２１は、局所領域形状解析部１２０で推定した円筒モデルパラメータ、および柱状構造物の中心軸周辺の３次元点群を取得し、第１の円筒モデルとたわみ付き第２の円筒モデルの当てはめ処理を行う。これらの結果を比較して、たわみ付き第２の円筒モデルの優位性の度合いを示すスコア比を算出し、算出したスコア比と共に、たわみがあるか否かを示すたわみ状態パラメータを、たわみ期待値推定部１２２へ供給する。
【００４４】
たわみ期待値推定部１２２は、柱状構造物パラメータを取得し、たわみ付き第２の円筒モデルたわみ付き第２の円筒モデルにおける上端部の中心位置の移動量dhもしくは上端部と下端部の中心位置と凸部（凸関数の極値位置）との最短距離dcとの距離を指標としてたわみ量を推定し、推定したたわみ量に基づいて、たわみ量の期待値を算出し、信頼度算出部１２３へ提供する。
【００４５】
信頼度算出部１２３は、推定したたわみ量及び算出したたわみ量の期待値を取得し、推定したたわみ量の信頼度を算出する。
出力処理部１２４は、スコア比、たわみ状態パラメータ、たわみ量の期待値、及び信頼度を、記憶部１０３の柱状構造物パラメータ記憶部１１２へ供給する。
繰り返し処理部１２５は、処理が繰り返し実行されるように制御する。
【００４６】
構造物安全性判定部１０６は、柱状構造物パラメータ記憶部１１２から柱状構造物パラメータを取得し、柱状構造物の安全性を判定し、判定結果を出力部１０７へ供給する。このとき、たわみ量および信頼度から不安全と判定されたときには、不安全と判定された柱状構造物全ての座標および警告信号を出力部１０７へ供給する。
【００４７】
出力部１０７は、ＣＲＴ（Cathode Ray Tube）、ＬＣＤ（Liquid Crystal Display）、ＰＤＰ（Plasma Display Panel）等のディスプレイやプリンタなどであり、本装置により求めたたわみ量および警告表示の有無を画面表示もしくは印刷等により保守管理担当者へ示すことに用いる。ここで、警告信号を構造物安全性判定部１０６から取得した場合には、不安全と判定された柱状構造物の座標情報と信頼度、たわみ量をリスト化したものを表示することに用いる。
【００４８】
なお、柱状構造物検出部１０４、たわみ推定処理部１０５（局所領域形状解析部１２０、円筒モデル比較部１２１、たわみ期待値推定部１２２、及び信頼度算出部１２３）、及び構造物安全性判定部１０６の各々をコンピュータで構成するようにしてもよい。
【００４９】
次に、本実施形態の作用を説明する。
（たわみ推定処理部１０５の動作説明）
まず、本実施形態のたわみ推定装置１００によるたわみ推定処理部１０５全体の入出力について説明する。ただし、たわみ推定処理部１０５の各部（１２０〜1２３）は、各部に対応する処理（Ｓ３＿１〜Ｓ３＿４）を実行する。図８は、本実施形態のたわみ推定処理部１０５が実行する処理のフローチャートである。
【００５０】
ステップＳ１で、入力処理部１１８は、３次元点群記憶部１１０から３次元点群、演算処理用パラメータ記憶部１１１からたわみ推定処理部１０５の演算処理で使用する演算処理用パラメータ、柱状構造物パラメータ記憶部１１２から柱状構造物パラメータ（位置、中心軸の方向、半径、長さ）を入力する。
【００５１】
ここで、各柱状構造物を区別する番号をＩ(∈1,2,...,ＮＩ)、柱状構造物の総数をＮＩで表示する。柱状構造物パラメータとは、柱状構造物の位置、例えば、下端部の３次元座標PI_bottom、中心軸の方向
【００５２】
【数１】
2015224980.tif 000003
【００５３】
、構造物の長さＬI、半径ＲI等の柱状構造物に関する形状情報のこととする。
【００５４】
ステップＳ７で、繰り返し処理部１２５は、番号Ｉに１をセット（代入）する。以下では、繰り返し処理部１２５は、ステップＳ４の終了条件を満たすまでＮＩ回、ステップＳ２とステップＳ３の処理が実行されるように処理する。
【００５５】
ステップＳ２において、取得部１１９は、柱状構造物Iの中心軸から半径Ｒcut（=ＲI+ΔR）以内に含まれる３次元点群を取得する。被写体計測部１０１は、柱状構造物Iの表面の３次元位置を計測するので、３次元点群は、柱状構造物Iの表面の３次元位置と、存在すれば柱状構造物Iの表面に取り付けられた付属品の表面の３次元位置とが含まれる。上記ΔRは演算処理用パラメータ記憶部１１１に登録されているパラメータである。以降では、このようなパラメータを「実験的に決めるパラメータである」と記載する。本実施形態においては、例えばΔＲ＝０．０５(ｍ)とした。
【００５６】
ここで、取得した３次元点群の数をＮｍ、各３次元点を区別する番号をｍ(∈1,2,...,Ｎm)で表示する。注目点ｍの位置座標をベクトルで表すと、ベクトルｐｍは次式によりＸ軸、Ｙ軸、Ｚ軸の成分で構成される。ただし、右上の添え字の記号「Ｔ」は転置を意味とする。
【００５７】
【数２】
2015224980.tif 000004
【００５８】
ただし、ｘｍ、ｙｍ、ｚｍはＸ、Ｙ、Ｚ軸の成分である。
【００５９】
ステップＳ３において、たわみ推定処理部１０５は、柱状構造物ＩのパラメータおよびステップＳ１で取得した３次元点群を入力し、たわみ推定処理を行い、たわみ量を推定し、たわみ量の期待値と信頼度を算出する。また、スコア比及びたわみ状態パラメータを出力する。
【００６０】
ステップＳ４で、繰り返し処理部１２５は、番号Ｉが柱状構造物の総数ＮＩ以上か否かを判断することにより、全ての柱状構造物について以上の処理（Ｓ２〜Ｓ３）が終了したか否かを判断する。全ての柱状構造物について以上の処理（Ｓ２〜Ｓ３）が終了したと判断された場合、ステップＳ５で、出力処理部１２４は、各柱状構造物について、スコア比と、たわみ状態パラメータと、推定したたわみ量と、たわみ量の期待値と、信頼度とを記憶部１０３（柱状構造物パラメータ記憶部１１２）に出力する。全ての柱状構造物について処理が終了していないと判断された場合、ステップＳ６で、繰り返し処理部１２５は、Ｉ←Ｉ＋１とし、処理は、ステップＳ２へ戻る。ただし、数式「Ｉ←Ｉ＋１」はＩの数を１つ増やしてＩに代入するという意味である。
【００６１】
以下、ステップＳ３について詳細な説明をする。以下、座標系の軸に対して、右上の添え字として記号「'」(ダッシュ)があるものは、ローカル座標系の軸であることを明示するために記載している。
【００６２】
（ステップＳ３−１：局所領域形状解析処理）
柱状構造物Ｉのパラメータとその柱状構造物Ｉから半径Ｒcut以内の３次元点の点群Ｐmを入力とし、複数（具体的には、Ｓ３１３で決定されるＮｒａｎ個）の局所領域円筒パラメータ（中心位置、半径、中心軸の方向）とテーパの有無を出力する。以下では、柱状構造物Ｉの中心軸から半径Ｒcut以内の３次元点群Ｐｍのことを柱状構造物Ｉの周辺点群Ｐmと記載する。
【００６３】
ステップＳ３−１の処理を示すフローチャートを図９に示す。
Ｓ３１１において、局所領域形状解析部１２０は、柱状構造物Ｉの周辺点群Ｐmを入力する。具体的には、下記の式を満たす３次元点が入力される点である。
【００６４】
【数３】
2015224980.tif 000005
【００６５】
ただし、ｐallは計測した全ての３次元点群の点を意味する。
【００６６】
Ｓ３１２において、局所領域形状解析部１２０は、柱状構造物Ｉの中心軸から半径Ｒcut以内の領域において、各々柱状構造物Ｉの中心軸方向の隣の領域と一部重複する複数の局所領域を設定する。周辺点群Ｐmに含まれる３次元点をｐｍ（ｍ∈１,２,３,..., Ｎｍ)、点群の総数をＮｍ、柱状構造物Ｉの中心軸方向をＺI'軸、中心軸に垂直な２つの軸をＹ’軸、Ｘ’軸、下端部の中心位置ＰI_bottomを原点としたローカル座標系を設定すると、図１０に示すように点周辺点群Ｐmを中心軸の方向
【００６７】
【数４】
2015224980.tif 000006
【００６８】
に分割する。ローカル座標原点PI_bottomから、ＺI'軸方向へ長さΔｓずつ区切った区間を設定すると、局所領域ＲＳs(ｓ∈１,２,３,...,ＮS）に含まれる３次元点は以下の式を満たす。
【００６９】
【数５】
2015224980.tif 000007
【００７０】
ただし、Δｓは実験的に決めるパラメータであり、本実施形態ではΔｓ＝１．５（ｍ）とする。
【００７１】
ここで数学的な記号の意味を記載すると、記号
【００７２】
【数６】
2015224980.tif 000008
【００７３】
は方向、つまり長さが１のベクトルであること（２ノルムが１であること）を明示するための記号とし、記号
【００７４】
【数７】
2015224980.tif 000009
【００７５】
はベクトルの２ノルムの大きさとする。記号ベクトル間の記号「・」はベクトルの内積演算、記号「×」はベクトルの外積、スカラー量とベクトル間の記号「・」は掛け算を意味する。
【００７６】
式（３−１）の物理的な意味は、柱状構造物Iの周辺点について、ローカル座標のＺ’軸に射影した点が(ｋ・Δｓ)から±Δｓの範囲に入る３次元点群である。
【００７７】
ステップＳ３１３において、局所領域形状解析部１２０は、各局所領域内の３次元点群に対して第１の円筒モデルの当てはめを行うことで、該当する局所領域における円筒の表面形状の情報を抽出し、各局所領域の円筒パラメータ登録処理をする。具体的には、柱状構図物の付属品の影響を抑えるために、ＲＡＮＳＡＣ処理（図１１）により円筒当てはめを行い、最もスコア値の高い第１の円筒モデルのパラメータから上位Ｎｒａｎ個について、局所領域における円筒パラメータとして登録を行う。
【００７８】
ＲＡＮＳＡＣによる円筒当てはめは既存手法である。ただし、本実施形態においてたわみ推定を高い性能で実施する必要があるために、図１１のフローチャートを用いて説明する。なお、局所領域形状解析部１２０は、図１１に示した処理を、上記ステップＳ３１２（図９）で設定された複数の局所領域の各々について実行する。
【００７９】
ステップＡ＿１において、局所領域形状解析部１２０は、局所領域内における柱状構造物Ｉの中心軸から半径Ｒcut以内に位置する３次元点群（周辺点群）を入力する。
【００８０】
ステップＡ＿２において、局所領域形状解析部１２０は、ステップＡ＿３からステップＡ＿６の処理をＫloop回まで繰り返す回数ｋを１にセットする。回数Ｋloopは実験的に決めるパラメータであり、本実施形態ではＫloop＝１０００とした。
【００８１】
ステップＡ＿３において、局所領域形状解析部１２０は、局所領域内における柱状構造物Ｉの中心軸から半径Ｒcut以内に位置する３次元点群の中からランダムに２点を選択（サンプリング）し、２点間の方向を処理ｋ回目における中心軸の方向
【００８２】
【数８】
2015224980.tif 000010
【００８３】
とする。具体的には２点の座標の差分をとり、ノルムを１にするように長さを正規化すればよい。
【００８４】
【数９】
2015224980.tif 000011
【００８５】
の記号「ｋ」はステップＡ＿２からステップＡ＿７の繰り返し処理における処理のｋ(∈1,2,...,Ｋloop)番目の回数を意味する。
【００８６】
ステップＡ＿４において、局所領域形状解析部１２０は、局所領域内における柱状構造物Ｉの中心軸から半径Ｒcut以内に位置する３次元点群の中からランダムにＮｒ個の点を選択し、選択したＮｒ個の各点の中心軸の方向
【００８７】
【数１０】
2015224980.tif 000012
【００８８】
に垂直でかつ局所領域の重心を通る平面へ射影した座標位置を求める。平面の基底ベクトルは、
【００８９】
【数１１】
2015224980.tif 000013
【００９０】
に直交して、かつ重心を通る２つのベクトルであればよい。例えば、グローバル座標のＸ軸とＹ軸について、
【００９１】
【数１２】
2015224980.tif 000014
【００９２】
方向に直交するように、グラム・シュミットの正規直交化法を用いて変換してもよい。この平面上での座標位置は、平面の基底の軸方向
【００９３】
【数１３】
2015224980.tif 000015
【００９４】
と注目点ｊの位置ベクトルｐｊとの内積で求まる。ランダムに選択した点の位置座標をｐｊ、ランダムに選択したNr個の点の重心位置をqgkとすると、上記平面での座標[ｘ'ｊ,ｙ'ｊ]は次式で求まる。
【００９５】
【数１４】
2015224980.tif 000016
【００９６】
ここで、Ｎｒは実験的に決めるパラメータであり、本実施形態ではＮｒ＝５とした。
【００９７】
ステップＡ＿５において、局所領域形状解析部１２０は、選択したＮｒ個を上記平面上に射影した点群について、最小二乗法により円の当てはめを行い、各々局所領域の形状を規定する半径ｒekと平面上での中心位置の座標［ｘek,yek］を円筒パラメータとして求める。具体的には、以下の評価式Ｅstraightを最小化するように求めればよい。
【００９８】
【数１５】
2015224980.tif 000017
【００９９】
ステップＡ＿６において、局所領域形状解析部１２０は、ステップＡ＿５で求めた円筒パラメータ（半径、中心位置）の評価を行う。評価方法は幾つか存在するが、例えば円筒表面上に存在する３次元点群の数や円筒表面上に存在する３次元点の範囲（面積）等が用いられる。
【０１００】
点の数をモデルの評価基準となるスコア値として求めた例を以下に示す。スコア値は、局所領域ＲＳkに含まれる合計Ｎs個の３次元点ｐsの各々を、中心軸
【０１０１】
【数１６】
2015224980.tif 000018
【０１０２】
に垂直な平面へ射影した点（ｘ's,y's)、半径ｒek、及び中心位置の座標(xek,yek）を用いて、下記の評価関数Scorestraightにより求まる。
【０１０３】
【数１７】
2015224980.tif 000019
【０１０４】
ここで、Ｔｈｅは実験的に決めるパラメータであり、本実施形態ではＴｈｅ＝０．０２（ｍ）とした。上記式のＦはスカラーの値であることを意味する。
【０１０５】
ステップA_7において、局所領域形状解析部１２０は、繰り返し回数がＫloop以上終了したか判定を行う。Ｋloop以上の処理が終了したと判断された場合には、処理は、ステップＳＡ＿８に進む。Ｋloop以上の処理が終了していないと判断された場合には、処理は、ステップＡ＿３に戻る。
【０１０６】
ステップＡ＿８において、局所領域形状解析部１２０は、最も高いスコア値からＮｒａｎ個の円筒パラメータを出力する（図１２（Ｂ）参照）。このとき出力する円筒パラメータは、図１２（Ａ）に示すように、局所領域の形状を規定する中心位置、半径、中心軸の方向、及びスコア値である。中心位置については、平面上の座標系(2次元ベクトル)からグローバル座標系へ変換をした位置座標ｑｓ(３次元ベクトル)として出力する。ステップＡ＿８の後、処理は、図９のステップＳ３１４に進む。
【０１０７】
【数１８】
2015224980.tif 000020
【０１０８】
以上が円筒当てはめによる、円筒パラメータ登録処理である。
【０１０９】
図９のステップＳ３１４において、局所領域形状解析部１２０は、テーパ有無の判定処理を行う。図１３に示すように、各局所領域ＲＳsについて最も良いスコア値（Scorestraightが最大）のときの半径（局所領域円筒半径）の大きさを縦軸、各局所領域の中心位置を横軸にプロットしたときの2次元座標点群を求める。横軸におけるＺｍｉｎは、最下位置の局所領域の中心位置である。このとき、点群に直線当てはめを行ったときの傾きの絶対値が、閾値ＴＨlean以上のときに「テーパが有る」と判定する。図１３に示すように、傾きをｂ、半径ｒの大きさをＹ軸に、Ｙt軸との切片の値をＣ、中心軸方向をＸt軸とした２軸を例で示すと、
【０１１０】
【数１９】
2015224980.tif 000021
【０１１１】
の式で表した時に、bの符号が負であり、かつ|b|の大きさが閾値ＴＨleanより小さいときは柱状構造物にテーパが無いと判定し、|b|の大きさが閾値ＴＨlean以上ならテーパが有ると判定する。物理的な意味を補足すると、「テーパが有る」とは、柱状構造物の形状が中心軸方向において一定の割合で小さく変化していることを意味する。
【０１１２】
なお、ステップＳ３１４については処理を省略してもよい。この場合、ステップＳ３−２におけるステップＳ３２１で入力されることがなく、ステップＳ３２２における第１のスコア値Scorestraightの算出に用いないこととなる。
【０１１３】
ステップＳ３１５において、局所領域形状解析部１２０は、各局所領域の円筒パラメータ（中心位置、半径、中心軸の方向）とテーパの有無とを出力する。ステップＳ３１５の後、処理は、図８のステップＳ３−２に戻る。
【０１１４】
（ステップＳ３−２：円筒モデル比較処理）
円筒モデル比較部１２１は、ステップＳ１で入力された柱状構造物Ｉのパラメータおよび、柱状構造物Ｉの周辺点群Pm、ステップＳ３−１の出力であるテーパ有無を入力とし、円筒モデル比較結果のスコア比Ｒａｔｉｏscoreを出力する。
【０１１５】
次に、ステップＳ３＿２の円筒モデル比較処理について説明する（図１４参照）。ステップＳ３＿２では、円筒モデル比較部１２１は、柱状物体Iについて、柱状構造物Ｉにたわみがないと仮定した場合の柱状構造物Ｉを規定する第１の円筒モデルと、柱状構造物Ｉにたわみがあると仮定した場合の柱状構造物Ｉを規定するたわみ付き第２の円筒モデルとの当てはめ処理を行い、２つのモデルに対して算出される第１のスコア値及び第２のスコア値を比較することで、不安全と判定されるべきたわみが発生しているかどうか判定を行う。これにより、本来はまっすぐな柱状構造物が、計測時のノイズの影響により、誤ってたわみが発生していると判定することを抑制する。
【０１１６】
ステップＳ３２１では、円筒モデル比較部１２１は、ステップＳ１で入力された柱状構造物Ｉのパラメータ、柱状構造物Iの周辺点群、及びステップＳ３＿１の出力であるテーパ有無の判定結果を入力する。
【０１１７】
ステップＳ３２１の後、処理は、ステップＳ３２２、Ｓ３２３に移行する。なお、円筒モデル比較部１２１は、ステップＳ３２２の後に、ステップＳ３２３を実行してもよく、ステップＳ３２３の後に、ステップＳ３２２を実行してもよい。
【０１１８】
ステップＳ３２２では、円筒モデル比較部１２１は、柱状構造物Iの周辺の３次元点群を用いて、柱状構造物Ｉにたわみがないと仮定した場合の柱状構造物Ｉを規定する第１の円筒モデルの当てはめによる第１のスコア値Scorestraightを算出する。この処理はＲＡＮＳＡＣによる円筒当てはめ処理（図９のステップＳ３１３（図１１））と同じである。ただし、ステップＳ３においてテーパが有りと判定された場合のみ、図１５のようにテーパ値により一定の割合で半径の大きさが変化する。
【０１１９】
そのため、円筒当てはめ処理ステップＡにおける入力する点が、局所領域ＲＳｋに含まれる３次元点群ではなくて、柱状構造物Iの周辺の３次元点群全てであることと、ステップＡ＿６のスコア値の算出式のみが異なる。
【０１２０】
ステップＡ＿６のスコア算出時における、繰り返し処理ｋ回目の中心軸Ｚ’kに垂直な平面の基底ベクトルを
【０１２１】
【数２０】
2015224980.tif 000022
【０１２２】
とし、テーパ値をα、重心位置をｑgk、注目点の番号をｍ、注目点の位置ベクトルをｐｍとすると、第１のスコア値Scorestraightは、次の通りである。なお、本ステップ３２２では、柱状構造物Ｉにはたわみがないと仮定している。よって、各局所領域の中心位置は、直線状に位置する。下の式の中心位置（ｘek,ｙek)及び重心位置ｑgkは、複数の局所領域の中から選択した１つの局所領域の中心位置及び重心位置である。
【０１２３】
【数２１】
2015224980.tif 000023
【０１２４】
ここで、テーパ値αは柱状構造物の規格値を使用してもよいし、実験的に決定してもよい。本実施形態では、α＝１／１５０とした。閾値ＴｈeはステップＡ＿６と同じである。
【０１２５】
ここで、上記中心位置及び半径は、柱状構造物Ｉにたわみがないと仮定した場合の柱状構造物Ｉを規定する第１の円筒モデルのモデルパラメータであり、当該柱状構造物Ｉの形状を表現する。実際の柱状構造物Ｉに付属物がなくたわみもなければ、第１のスコア値Scorestraightは、３次元点の総数Ｎmと一致致し又は近い値となる。即ち、上記第１の円筒モデルのモデルパラメータにより表現される柱状構造物の形状は、実際の柱状構造物Ｉと一致する。上記中心位置及び半径は、実際の柱状構造物Ｉの形状を正しく規定していると評価できる。
しかし、実際の柱状構造物Ｉにたわみがあると、上記中心位置及び注目点間の距離と上記半径との差Ｆの絶対値は閾値Ｔｈe以上と判定され、第１のスコア値Scorestraightは、３次元点の総数Ｎmより少なくなる。即ち、上記第１の円筒モデルのモデルパラメータにより表現される柱状構造物の形状は、実際の柱状構造物Ｉと一致しない。実際の柱状構造物Ｉにあるたわみが大きくなるに従って、第１のスコア値Scorestraightは３次元点の総数Ｎmより小さくなる。
【０１２６】
ステップＳ３２３では、円筒モデル比較部１２１は、柱状構造物Iの周辺の３次元点群を用いて、ＲＡＮＳＡＣによるたわみ付き第２の円筒モデルの当てはめ処理を行う。ステップＳ３２２の第１の円筒モデルの当てはめ処理との大きな違いは、この中心軸の設定方法であり、スコア値の算出は同様に行うことができる。中心軸は、ステップＳ３＿１の局所領域形状解析処理で求めた円筒モデルパラメータを用いて、N次の凸関数により推定する。このとき、凸関数は２次関数でも３次関数でもよいが、本実施形態では２次関数を用いた例を示す。
【０１２７】
第２の円筒モデル当てはめの詳細な説明の前に、図１６（Ａ）〜図１６（Ｃ）に中心軸推定処理の概要を説明する。まず、円筒モデル比較部１２１は、局所領域形状解析処理（図８のステップ３−１（図９））で推定した円筒モデルパラメータをランダムに２つ選び、たわみ付き第２の円筒モデルの中心軸(凸関数)が存在する平面を推定する。次に、円筒モデル比較部１２１は、この平面上に、局所領域形状解析処理で推定した円筒パラメータからランダムに数点選択した中心位置の点を射影し、平面上（ローカル座標系）での座標を求める。円筒モデル比較部１２１は、この射影した点群を用いて、凸関数のパラメータを算出する。円筒モデル比較部１２１は、最後の凸関数のパラメータから、たわみ開始位置を算出する。
【０１２８】
以下、ステップＳ３２３を実現する方法として、円筒モデル比較部１２１は、ステップＢのＲＡＮＳＡＣによるたわみ付き第２の円筒モデルの当てはめ処理を説明する（図１７参照）。
【０１２９】
ステップＢ＿１において、円筒モデル比較部１２１は、ステップＳ１で入力された柱状構造物Iのパラメータ、ステップＳ２で取得した柱状構造物Iの周辺点群Pm、Ｓ３１４で出力したテーパの有無の判定結果、及び局所領域円筒パラメータ登録処理で出力した円筒パラメータを入力する。
【０１３０】
ステップＢ_２において、円筒モデル比較部１２１は、回数ｋを１にセットする。
【０１３１】
ステップＢ_３において（図１６（Ａ）参照）、円筒モデル比較部１２１は、第２の円筒モデルの中心軸が存在する平面を推定する。まず、たわみ開始位置Ｇ0での接線方向u軸を推定する。局所領域で推定した円筒パラメータ（円の中心位置及び半径、中心軸方向）を区別する番号c(c∈1,2,3,...,Nc)、円筒パラメータcの中心軸方向をuc、局所領域において推定した円筒パラメータの数をNcとすると、ランダムに選んだ２つの局所領域の中心軸の方向ｕc1、ｕc2を用いて以下の式で接線方向ｕkを求める。
【０１３２】
【数２２】
2015224980.tif 000024
【０１３３】
ただし、frandは０〜１の値をとる乱数とする。上記式は、ランダムに選んだ２つの局所領域の中心軸の方向ｕc1、ｕc2を用いて接線方向ｕkを合成（計算）することを意味する。
【０１３４】
また、後処理のために平面の基底ベクトルであるｕkの第三成分が正となるように向きを調整する。ｕkの第三成分、つまりグローバル座標系でのZ軸の値が負の値であるとき、ｕkの向きを変更する。
【０１３５】
【数２３】
2015224980.tif 000025
【０１３６】
たわみ方向を規定するｖ軸は、ｕc1からｕkの方向成分を、グラム・シュミットの直交化法により0にして、ノルムが１となるように正規化することで求める。
【０１３７】
【数２４】
2015224980.tif 000026
【０１３８】
上記平面の中心位置は、選択した２つの局所領域の円筒パラメータの中心位置（ｑc1,ｑc2)の重心とする。
【０１３９】
【数２５】
2015224980.tif 000027
【０１４０】
ここで、添え字の「ｋ」は、繰り返し処理がｋ回目のときの推定結果であることを意味する。
【０１４１】
ステップＢ＿４において、円筒モデル比較部１２１は、たわみ付き第２の円筒モデルの中心軸の存在する平面をローカル座標系とし、ランダムに選択した円筒パラメータにおけるNj個の中心位置を中心軸が存在する平面へ射影して、その平面で基底される空間（たわみ円筒ローカル座標系）での座標を求める。
【０１４２】
以下、このたわみ円筒ローカル座標系の軸はｕ軸、ｖ軸で表現し、その座標値（スカラー値）は、それぞれｘとｙで表現する。また位置ベクトルの右上の添え字として記号「'」(ダッシュ)があるものは、ローカル座標に射影した点の位置ベクトルもしくはローカル座標系での座標値であることを明示するためにつけている。
【０１４３】
局所領域で推定した円筒パラメータについて、ランダムに選択した円筒パラメータを区別する番号をｊ、その円筒パラメータの中心位置をｐjとすると、たわみ円筒ローカル座標系での座標ｐj（＝［ｘj,ｙj］Ｔ）への変換式は次式で求まる。
【０１４４】
【数２６】
2015224980.tif 000028
【０１４５】
ステップＢ＿５において（図１６（Ｂ）参照）、円筒モデル比較部１２１は、たわみ付き第２の円筒モデルの中心軸を表現する凸関数を、たわみ付き第２の円筒モデルのモデルパラメータとして求める。凸関数は２次関数でも３次関数やそれ以上の次元の数でもよいが、本実施形態では２次関数とする。経験的に、安全・不安全の評価を行う柱状構造物は、形状の変形量が微小であるため、ノイズに影響されにくい２次か３次などの低い次元の関数の方が好ましい。本実施形態では２次関数を例として示す。２次数の係数をa、b、cとすると、下記の評価関数を最小化するように係数として求めればよい。
【０１４６】
【数２７】
2015224980.tif 000029
【０１４７】
ここで、評価関数Ｅbentについて、最小二乗法を用いて係数ａ、ｂ、ｃを求めてもよいし、非線形の数値解析手法を用いて求めてもよい。
【０１４８】
求めた係数により、たわみ開始位置ｇkは次式で求まる（図１６（Ｃ）参照）。
【０１４９】
【数２８】
2015224980.tif 000030
【０１５０】
次に、たわみ開始位置での半径の大きさｒkを設定する。本来中心軸において下端からの距離に応じて半径は変化するが、たわみ量が微小であると考えて、たわみ開始位置での接線方向の距離に比例して半径が変化する、つまり近似できると考える。
【０１５１】
選択したＮｊ個の円筒パラメータの中からランダムに１つ選んだ半径ｒselectに対して、たわみ開始位置と選択した円筒パラメータの中心位置のu軸上での位置xselectとの距離に比例した値とする。
【０１５２】
【数２９】
2015224980.tif 000031
【０１５３】
ただし、テーパがないと判定された柱状構造物Iについてはｒk＝ｒselectとする。また、Ｎｊは実験的に決めるパラメータであり、本実施形態ではＮｊ＝５とした。
【０１５４】
ステップＢ＿６において、円筒モデル比較部１２１は、第２のスコア値ScoreBentを算出する。柱状構造物Iの周辺の各３次元点ｐｍについて、たわみ付き第２の円筒モデルのモデルパラメータにより表現される柱状構造物の表面上に存在するかどうか判定し、モデル表面上に存在するときに第２のスコア値として加算する。
【０１５５】
モデル表面上に存在するかどうかは、注目点から、凸関数で表現された中心軸への最短距離により判定する。具体的には、注目点ｐｍからの垂線の足ｑhを求め、線分ｐｍｑhがたわみ付円筒モデルの半径ｒhと近い値であれば、モデル表面上に存在すると判定する。
【０１５６】
３次元空間中の曲線モデルと３次元点との最短距離を算出する方法は幾つかあるが、本実施形態では、図１８に示すようにステップB_5で設定した凸関数の基底軸（u軸、v軸）へ射影した点ｐ'ｍと凸関数の中心軸との最短距離から三平方の定理により算出する。射影点ｐ'ｍ（＝［ｘ’ｍ,ｙ’ｍ］）は、たわみ円筒ローカル座標系の座標として次式により求まる。
【０１５７】
【数３０】
2015224980.tif 000032
【０１５８】
次に、射影点のu軸上での座標により場合分けして最短点を求める。座標位置ｘ’ｍの値が、たわみ開始位置ｇkのｕ軸上での位置ｕgの値より小さいときは、ステップＢ＿６＿１を行い、u軸上での位置が大きいときはステップＢ＿６＿２を行う。
【０１５９】
ステップＢ＿６＿１において、円筒モデル比較部１２１は、３次元点ｐｍとたわみ開始点ｇkを通るｕ軸と平行な直線との距離として最短距離ｄｉｓｔ(ｐｍ)を求める。
【０１６０】
【数３１】
2015224980.tif 000033
【０１６１】
ステップＢ＿６＿２において、円筒モデル比較部１２１は、まず３次元点ｐｍを射影した射影点ｐ’ｍとたわみ開始点ｇkを通るステップＢ＿５で求めた2次関数について、最短距離dist1を求める。
【０１６２】
ここで、幾何的な性質から、点からある曲線までの最短距離の位置まで引いた線分と、最短距離の位置での接線は直交するという関係性がある。よって、最短位置ｑｈ（＝ｑ’ｈ＝［ｘ’ｈ，ｙ’ｈ］）での接線方向
【０１６３】
【数３２】
2015224980.tif 000034
【０１６４】
と射影点ｐ’ｍから最短位置ｑｈまでの方向
【０１６５】
【数３３】
2015224980.tif 000035
【０１６６】
の内積は０となり、次式が成立する。
【０１６７】
【数３４】
2015224980.tif 000036
【０１６８】
未知数（変数）はｘ’ｈのみであるから、ｘ’ｈの方程式として解くことができる。この方程式はニュートン法や修正ニュートン法などの数値解析手法で解くことが可能である。
【０１６９】
ただし、この方程式は高次の方程式のため、一般的に複数の解が求まる。そのため、求まる複数の解について、射影点ｐ’ｍから最も近い位置が最短位置ｑｈの位置ｘ’ｈとする。
【０１７０】
最短距離dist1は次式で求まる。
【０１７１】
【数３５】
2015224980.tif 000037
【０１７２】
次に、円筒モデル比較部１２１は、点ｐｍと射影した射影点ｐ’ｍとの距離dist2を求める。
【０１７３】
【数３６】
2015224980.tif 000038
【０１７４】
点ｐｍと垂線の足の位置ｑｈとの距離distは、幾何的な関係から次式で求まる。
【０１７５】
【数３７】
2015224980.tif 000039
【０１７６】
ステップＢ＿６＿２において、上記では幾何的な関係から最短距離を求めたが、この方法以外にも例えば、ラグランジュの未定乗数法を用いて解いてもよい。制約条件は、最短位置の点が２次関数上に存在することとし、射影点ｐ’ｍと凸関数曲線上の点において距離が最小となる解（極値）を求めればよい。
【０１７７】
ステップＢ＿６＿１もしくはステップＢ＿６＿２により求めた点ｐｍにおける最短距離ｄｉｓｔ（ｐｍ）を用いて、第２のスコア値は次式により求まる。
【０１７８】
【数３８】
2015224980.tif 000040
【０１７９】
ここでＴｈｅは実験的に決まるパラメータであり、本実施形態ではＴｈｅ＝０．０２とした。
【０１８０】
ここで、上記中心軸は、柱状構造物Ｉにたわみがあると仮定した場合の柱状構造物Ｉを規定するたわみ付き第２の円筒モデルのモデルパラメータであり、当該柱状構造物Ｉの形状を表現する。実際の柱状構造物Ｉに付属品がなくかつ当該中心軸に対応するたわみがあると、第２のスコア値Scorebentは、３次元点の総数Ｎmと一致致し又は近い値となる。よって、第２のスコア値Scorebentが３次元点の総数Ｎmと一致し又は近い値の場合には、実際の柱状構造物Ｉに上記中心軸に対応するたわみが存在すると評価できる。一方、実際の柱状構造物Ｉにたわみがないと、ｄｉｓｔ（ｐｍ）がＴｈe以上となる注目点が多くなり、第２のスコア値Scorebentは、３次元点の総数Ｎmより少なくなる。
【０１８１】
ステップＢ＿７において、円筒モデル比較部１２１は、回数ｋがＴloop以上か否かを判断することにより、繰り返し処理の回数ｋがＴloop回以上か否かを判断し、繰り返し処理の回数ｋがＴloop回以上と判断された場合、処理は、ステップＢ−８に進む。繰り返し処理の回数ｋがＴloop回以上と判断されなかった場合、ステップＢ−９において、円筒モデル比較部１２１は、ｋ←ｋ＋１とし、処理は、ステップＢ＿３へ戻る。
【０１８２】
ステップＢ＿８において、円筒モデル比較部１２１は、スコア値Scorebentが最も高いモデルパラメータについて、柱状構造物Iの柱状構造物パラメータとして出力する。ここで、出力するモデルパラメータとは、たわみ中心軸の存在するローカル座標系および凸関数の係数および、たわみ開始位置、たわみ開始位置での半径、たわみ開始位置での凸関数の接線方向である。
【０１８３】
以降の処理のために、出力するパラメータの記号を記載する。
【０１８４】
たわみ開始位置Ｇ０←ｇｋ、たわみ開始位置での半径Ｒ０←ｒｋ、たわみ開始位置での凸関数の接線方向
【０１８５】
【数３９】
2015224980.tif 000041
【０１８６】
として出力する。ステップＢ−８の後、処理は、図１４のステップＳ３２４に進む。
【０１８７】
ステップＳ３２４において、円筒モデル比較部１２１は、円筒モデル当てはめ処理で求めた第１のスコア値Scorestraight と、第２の円筒モデルのスコア値Scorebentとを、以下の式にあるように比較して、スコア比Ｒａｔｉｏscoreを算出し、スコア比Ｒａｔｉｏscoreが、閾値ＴＨbentより大きいか否かを判定する。
【０１８８】
【数４０】
2015224980.tif 000042
【０１８９】
ここで、ＴＨbentは実験的に決めるパラメータであり、本実施形態ではＴＨbent＝１．２とした。ただし、ＴＨbentの設定は１．０以上とする。
【０１９０】
実際の柱状構造物Ｉにたわみがなければ、第１のスコア値Scorestraightは、３次元点の総数Ｎmと一致又は近い値となるのに対し、第２のスコア値Scorebentは、３次元点の総数Ｎmより少なくなる。この場合、スコア比Ｒａｔｉｏscoreは、ＴＨbent未満となる。逆に、実際の柱状構造物Ｉに、上記中心軸に対応するたわみがあると、第１のスコア値Scorestraightは、３次元点の総数Ｎmより少なくなるのに対し、第２のスコア値Scorebentは、３次元点の総数Ｎmと一致致又は近い値となる。よって、スコア比Ｒａｔｉｏscoreは、ＴＨbent以上となる。
【０１９１】
スコア比ＲａｔｉｏscoreがＴＨbent以上の場合、即ち、たわみ付き第２の円筒モデルについて算出される第２のスコア値Scorebentに優位性があるときに、たわみが存在すると判定する。たわみが有りと判定されたとき、ステップＳ３２５で、円筒モデル比較部１２１は、たわみ状態パラメータβ＝１とし、それ以外はたわみ状態パラメータβ＝０として出力する。
【０１９２】
なお、上記優位性を判断するためには、第１のスコア値Scorestraight及び第２のスコア値Scorebentに代えて、第１の円筒モデル、たわみ付き第２の円筒モデルの表面上に存在する点の範囲（面積）を用いてもよい。
【０１９３】
（ステップＳ３−３：たわみ量期待値算出処理）
たわみ量期待値算出処理は、Ratioscore＞ＴＨbentの場合に実行される。Ratioscore＞ＴＨbentでない場合には、算出されたたわみ量から柱状構造物Ｉにたわみが発生していると判断しても、これは、柱状構造物Ｉに取り付けられた付属品の影響により柱状構造物Ｉにたわみが発生していると判断できる。よって、柱状構造物Ｉには、不安定であると判断できるたわみは発生していないと判断できる。逆に、Ratioscore＞ＴＨbentの場合、算出したたわみ量には付属品の影響は少ない。しかし、算出されたたわみ量は、計測ノイズにオーバーフィッティング（過剰適合）して、算出されている可能性がある。そこで、本実施形態では、たわみ量期待値算出処理が実行される。
【０１９４】
たわみ期待値推定部１２２は、ステップＳ１で入力された柱状構造物Ｉのパラメータおよび、柱状構造物Ｉの周辺点群を入力とし、ステップＳ３−４へ出力する。
【０１９５】
たわみ期待値推定部１２２は、ステップＳ３−３では、たわみ量期待値算出処理を行う。たわみ量期待値の処理フローチャートは図１９に示す。
【０１９６】
図１９のステップＳ３３１において、たわみ期待値推定部１２２は、柱状構造物Ｉのパラメータおよび柱状構造物Ｉの周辺の３次元点群を入力する。
【０１９７】
以下では、ｄｌ＝１からｄｌ＝ＤLoopになるまで、ステップＳ３３２からステップＳ３３５の処理を繰り返し行う。まず、ステップＳ３３１Ａで、たわみ期待値推定部１２２は、回数ｄｌに1をセットする。
【０１９８】
ステップＳ３３２において、たわみ期待値推定部１２２は、柱状構造物Ｉの周辺点群Pmについて点密度の正規化処理として、サンプリング処理を行う。各局所領域ＲＳs(s∈1,2,3,...,NS)について、Ｌs個の点群をランダムに選択する（図２０参照）。当該局所領域はＳ３１２で設定したものである。局所領域は重複しているため、同一点が２回選択されることもあるが、それは問題ない。
【０１９９】
なお、たわみ期待値推定部１２２は、後述する処理（ステップＳ３３２〜Ｓ３３５）を繰り返す毎に、柱状構造物Ｉの周辺点群の中からサンプリングした３次元点群の組み合わせに基づいて、たわみ量を推定している。よって、たわみ期待値推定部１２２は、周辺点群の内の複数の点（ＬS個の点）を有する複数の組合せの各々についてたわみ量を推定している。
【０２００】
ここで、Ｌsは実験的に決めるパラメータである。本実施形態では、各局所領域ＲＳsについて円筒当てはめを行った出力すべて（Ｎｒａｎ個×Ｎｓ個）のスコア値の中央値とする。ただし、小数点以下は切り捨てた値とする。なお、ＮｒａｎはＳ３１３で登録処理をした局所領域円筒パラメータの数である。
【０２０１】
ステップＳ３３３において、たわみ期待値推定部１２２は、たわみ付き第２の円筒モデルの当てはめ処理を行う。この処理はステップＢの処理（図１４のステップＳ３２３（図１７））と同一であるが、入力する３次元点群が正規化処理を行った後の３次元点群であることが異なる。また、dl回目の処理の入力点群Ｐdlランダムに選択しているため、選択された３次元点群は繰り返し処理ごとに異なるが、３次元点の数はＬｓ×Ｎｓ個と同じである。 
【０２０２】
ステップＳ３３４において、たわみ期待値推定部１２２は、上記ステップＳ３３３の当てはめ処理の結果に基づいて、たわみ量γdlを推定する。たわみ量の指標は、本実施形態では図６（Ａ）に示すように、上端部の中心位置とたわみ開始位置Ｇ0での接線との距離ｄh、もしくは、図６（Ｂ）に示すように、上端部の中心位置と下端部の中心位置を結んだ直線と中心軸を表現する凸関数との最も離れた位置の距離ｄcとする２つの指標とする。
【０２０３】
（指標としてｄhを用いた場合の算出方法）
まず、図６（Ａ）を参照して、指標としてｄhを用いた場合の算出方法を説明する。
たわみ期待値推定部１２２は、本実施形態では、たわみ付き第２の円筒モデルの当てはめ処理により求めた中心軸を量子化して算出した点と上端部と下端部の中心位置を結んだ直線との距離として求める。たわみ開始位置Ｇ0の線方向
【０２０４】
【数４１】
2015224980.tif 000043
【０２０５】
について、下端部と上端部の間の区間をΔｗずつ区切り、それぞれの区間ｗ（∈1,2,3,...,Ｎｗ)の中心位置での中心軸のローカル座標系での座標を求める。
【０２０６】
次に、上記ローカル座標系での座標値をグローバル座標系での３次元位置をqwで表し、たわみ開始位置Ｇ0での接線をＵ’方向で表すと、上端部の中心位置ｑtopと下端部の中心位置ｑbottomは次式で求まる。
【０２０７】
【数４２】
2015224980.tif 000044
【０２０８】
上端部の中心位置ｑtopから接線
【０２０９】
【数４３】
2015224980.tif 000045
【０２１０】
の距離ｄhは次式により求まる。
【０２１１】
【数４４】
2015224980.tif 000046
【０２１２】
たわみ量γdlは次式により求まる。
【０２１３】
【数４５】
2015224980.tif 000047
【０２１４】
もしくは、単位長さあたりの変形量として、次式により求めてもよい。
【０２１５】
【数４６】
2015224980.tif 000048
【０２１６】
ただし、Δｗは実験的に決める長さの量であり、本稿ではΔｗ＝０．０１(ｍ)とした。柱状構造物Iの周辺点群について、Ｕ’軸上に射影したとき、射影した座標の最大値と最小値の差分をΔｗで割った値がＮｗとなる。
【０２１７】
（指標としてｄcを用いた場合の算出方法）
次に、図６（Ｂ）参照してｄcを用いた場合の算出方法を説明する。
たわみ期待値推定部１２２は、本実施形態では、たわみ付き第２の円筒モデルの当てはめ処理により求めた中心軸を量子化して算出した点と上端部と下端部の中心位置を結んだ直線との距離として求める。
【０２１８】
たわみ開始位置Ｇ0の接線方向
【０２１９】
【数４７】
2015224980.tif 000049
【０２２０】
について、下端部と上端部の間の区間をΔｗずつ区切り、それぞれの区間ｗの中心位置での中心軸のローカル座標系での座標を求める。次に、この座標値をグローバル座標系での３次元位置をｑｗで表すと、次式により凸部量ｄcは求まる。
【０２２１】
【数４８】
2015224980.tif 000050
【０２２２】
たわみ量γdlは次式により求まる。
【０２２３】
【数４９】
2015224980.tif 000051
【０２２４】
もしくは、単位長さあたりの変形量として、次式により求めてもよい。
【０２２５】
【数５０】
2015224980.tif 000052
【０２２６】
ステップＳ３３５において、たわみ期待値推定部１２２は、回数ｄｌがＤLoop以上であるか否かを判断することにより、繰り返し処理の回数ｄｌがＤLoop回以上であるか否かを判断する。繰り返し処理の回数ｄｌがＤLoop回以上であれば、処理は、ステップ３３６に進む。繰り返し処理の回数ｄｌがＤLoop回以上でないときは、ステップ３３８で、たわみ期待値推定部１２２は、ｄｌ＝ｄｌ＋１とし、処理は、ステップＳ３３２の処理へ戻る。
【０２２７】
ステップＳ３３６において、たわみ期待値推定部１２２は、ステップＳ３３４によって推定されたたわみ量γdlに基づいて、たわみ量の期待値を算出する。
【０２２８】
処理ｄｌ回目の第２のスコア値をScorebent(ｄｌ)、処理ｄｌ回目で推定したたわみ量をγdlで表すと、たわみ量期待値μは次式で求まる。
【０２２９】
【数５１】
2015224980.tif 000053
【０２３０】
なお、Scorebent(dl)はステップＳ３３３の繰り返し処理で算出されたものである。
【０２３１】
ステップＳ３３７において、たわみ期待値推定部１２２は、推定したたわみ量の期待値と繰り返し処理で求めたたわみ量（γ1, γ2,...,γDloop)の全てを出力する。ステップS３３７の後、処理は、図８のステップＳ３＿４に進む。
【０２３２】
（ステップＳ３−４：信頼度算出処理）
信頼度算出部１２３は、ステップＳ３−３で算出した推定したＤLoop個のたわみ量とたわみ量期待値μを入力とし、信頼度τを出力する。
【０２３３】
ステップＳ３−４において、信頼度算出部１２３は、ステップＳ３−３で算出した推定したＤLoop個のたわみ量とたわみ量期待値μに基づいて、信頼度を算出する。信頼度算出部１２３は、ステップＳ３−３で算出したＤLoop個のたわみ量γdlにおける標準偏差から信頼度を求める。図２１に示す繰り返し処理dl回目のたわみ量γdl、たわみ量期待値μ、たわみ量標準偏差σにより、信頼度τは次式で求まる。
【０２３４】
【数５２】
2015224980.tif 000054
【０２３５】
ここで、信頼度τは、０から１の値をとり、τ＝１．０のときに最も信頼度が大きいと判定し、τ＝０のときに最も信頼度が低いと判定する値である。なお、信頼度τは、上式以外に、例えば、１／（１＋σ）として計算してもよい。
【０２３６】
（構造物安全性判定部１０６の動作説明）
最後に、構造物安全性判定部１０６について説明する。構造物安全性判定部１０６は、たわみ推定処理部から出力された柱状構造物パラメータを用いて、各柱状構造物について安全か不安全な状態であるかの判定を行う。
【０２３７】
各柱状構造物において、柱状構造物パラメータであるたわみ状態パラメータβとたわみ期待値μと信頼度τとスコア比Ｒａｔｉｏscoreについて、以下の３つの方程式を満たすとき不安全な状態であると判定する。
【０２３８】
【数５３】
2015224980.tif 000055
【０２３９】
ここで、閾値ＴＨexpと閾値ＴＨτは実験的に決めるパラメータであり、本実施形態ではＴＨexp＝０．２（ｍ）、ＴＨτ＝０．９５とした。ＲａｔｉｏscoreはステップＳ３２４の出力値である。
【０２４０】
上記３つの式を満たす場合、不安全と判定された柱状構造物の全ての座標および警告信号を出力部１０７へ供給する。出力部１０７は、上記ディスプレイに、不安全と判定された柱状構造物の全てについて、各不安全と判定された柱状構造物に対応して、当該柱状構造物の位置座標、たわみ量、期待値、及び信頼度をリスト化して表示（警告表示）したり印刷したりする。これにより保守管理担当者へ不安全と判定された柱状構造物を示すことができる。
（第２の実施形態）
【０２４１】
上記の第１の実施形態において、計算効率の点から図８に示すたわみ推定処理の中で局所領域形状解析処理を行っていた。しかしながら、オフライン作業で十分な時間をかけて点群解析するときには、以下に示すように局所領域形状解析処理を省いてもよい。（図２２）
【０２４２】
このとき、ステップ３のステップ３−１の円筒モデル比較処理において、図１７のステップＢ−３とＢ−４の処理が、局所領域で推定した円筒パラメータを用いている部分が変更になる。
【０２４３】
ステップＢのステップＢ−３とＢ−４はたわみ円筒モデルの中心軸を推定する前処理として必要な処理である。図１６に示してあるように、円筒パラメータの中心軸方向を用いてたわみ中心軸の存在平面（ローカル座標系）を設定し、円筒パラメータの中心位置をプロットした点群を凸なＮ次関数で近似することで中心軸を推定している。この部分について、局所領域で推定した円筒パラメータの「中心軸方向」および「中心位置」の算出方法について代替手段を示す。
【０２４４】
ステップＢ_３において（図２３参照）、たわみ付き第２の円筒モデルの中心軸が存在する平面を推定する。まず、たわみ開始位置GOでの接線方向u軸を推定する。注目する柱状構造物Ｉの周辺点群からランダムに選んだ２点（ｑ1,ｑ2）選ぶ。その選んだ２点の法線ベクトルを、
【０２４５】
【数５４】
2015224980.tif 000056
【０２４６】
と記載すると、外積ベクトルｃ１を次のように、算出する。
【０２４７】
【数５５】
2015224980.tif 000057
【０２４８】
同様にして、柱状構造物Ｉの周辺点群から、更に選択した２点から外積ベクトルｕｋを求める。この２つの外積ベクトルから、以下の式で接線方向ｕｋを求める。
【０２４９】
【数５６】
2015224980.tif 000058
【０２５０】
ただし、ｆrandは０〜１の値をとる乱数とする。上記式は、算出した２つの外積ベクトルｕc1とｕc2の加算ベクトルを用いて接線方向ｕkを算出することを意味する。
【０２５１】
平面の中心位置は、選択した２点の法線ベクトル
【０２５２】
【数５７】
2015224980.tif 000059
【０２５３】
について、接線方向ベクトルｕｋに垂直な平面上での交点ｑｋとする。ただし、この平面は選択した２点の内の１点(ｑ1)を通るとする。
【０２５４】
【数５８】
2015224980.tif 000060
【０２５５】
たわみ方向ｖ軸は、外積ベクトルｕc1から、接線方向ｕkの方向成分を、グラム・シュミットの直交化法により０にして、ノルムが１となるように正規化することで求める。
【０２５６】
【数５９】
2015224980.tif 000061
【０２５７】
ステップＢ−４においては、円筒パラメータで求めた中心位置の代わりに、上記のように算出した交点ｑkをＮj個用いれば、以降は同様に処理でたわみ量を推定できる。
【０２５８】
つまり、ステップＢ−３において、事前にＮｒａｎ×Ｎｓ個の法線の外積ベクトルおよび法線の交点位置を求めておけば、以降は明細書に記載の手順でたわみ量が推定可能である。
【０２５９】
パラメータについては、ステップＳ３３２において、点密度の正規化処理のサンプリング処理の点群選択数Ｌsは、明細書では「各局所領域ＲＳsについて円筒当てはめを行った出力すべてのスコア値の中央値とする」と記載した。しかし、第２の実施形態においては実験的に決めるパラメータである。
【０２６０】
密度が少ない場所の点群数を参考にすればよく、例えば、柱状構造物の上端部（密度が少ない場所）において１ｍあたりの点群数が１００点、柱状構造物の長さが１０ｍとしたら、Ｌｓ＝１０００(１００×１０)点とすればよい。
【０２６１】
［実施形態の効果］
第１の実施形態及び第２の実施形態の各々により、３次元点群を用いて柱状構造物のたわみ推定を自動化でき、柱状構造物の保守点検の稼働コストを減らす効果がある。また、レーザースキャナで計測したときに、被計測範囲が狭いため計測ノイズの影響が大きい状況でも、中心軸のたわみ量の誤推定結果を出力することを抑制できる。更に、計測した点群から条件を変化させて複数回たわみ量を推定することで、推定したたわみ推定量の信頼度を算出でき、設備管理者に対して保守点検の優先度を示唆することができる。
【０２６２】
［変形例］
なお、たわみ推定装置１００の各処理を実行するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、当該記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより、たわみ推定装置１００に係る上述した種々の処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものであってもよい。また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、フラッシュメモリ等の書き込み可能な不揮発性メモリ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。
【０２６３】
さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムが送信された場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリ（例えばＤＲＡＭ（Dynamic Random Access Memory））のように、一定時間プログラムを保持しているものも含むものとする。また、上記プログラムは、このプログラムを記憶装置等に格納したコンピュータシステムから、伝送媒体を介して、あるいは、伝送媒体中の伝送波により他のコンピュータシステムに伝送されてもよい。ここで、プログラムを伝送する「伝送媒体」は、インターネット等のネットワーク（通信網）や電話回線等の通信回線（通信線）のように情報を伝送する機能を有する媒体のことをいう。また、上記プログラムは、前述した機能の一部を実現するためのものであっても良い。さらに、前述した機能をコンピュータシステムにすでに記録されているプログラムとの組合せで実現できるもの、いわゆる差分ファイル（差分プログラム）であっても良い。
【０２６４】
以上、本発明者によってなされた発明を、前記実施形態に基づき図面を参照して具体的に説明したが、本発明は、前記実施形態に限定されるものではなく、その要旨を逸脱しない範囲において種々変更可能であることは勿論である。
【符号の説明】
【０２６５】
１００    たわみ推定装置
１０１    被写体計測部
１０２    入力部  
１０３    記憶部
１０４    柱状構造物検出部  
１０５    たわみ推定処理部
１０６    構造物安全性判定部
１０７    出力部
１１０    ３次元点群記憶部
１１１    演算処理用パラメータ記憶部
１１２    柱状構造物パラメータ記憶部
１１８    入力処理部
１１９    取得部
１２０    局所領域形状解析部
１２１    円筒モデル比較部
１２２    たわみ期待値推定部
１２３    信頼度算出部
１２４    出力処理部
１２５    繰り返し処理部
【図７】
2015224980.tif 000062
【図８】
2015224980.tif 000063
【図９】
2015224980.tif 000064
【図１０】
2015224980.tif 000065
【図１１】
2015224980.tif 000066
【図１４】
2015224980.tif 000067
【図１７】
2015224980.tif 000068
【図１９】
2015224980.tif 000069
【図２２】
2015224980.tif 000070
【図１】
2015224980.tif 000071
【図２】
2015224980.tif 000072
【図３】
2015224980.tif 000073
【図４】
2015224980.tif 000074
【図５】
2015224980.tif 000075
【図６】
2015224980.tif 000076
【図１２】
2015224980.tif 000077
【図１３】
2015224980.tif 000078
【図１５】
2015224980.tif 000079
【図１６】
2015224980.tif 000080
【図１８】
2015224980.tif 000081
【図２０】
2015224980.tif 000082
【図２１】
2015224980.tif 000083
【図２３】
2015224980.tif 000084