(57)【特許請求の範囲】
【請求項１】
凹凸のある変動面上に散在する計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を凹凸値とみなしてＸＹ平面に配列した二次元画像を作成し，特定時刻の前記画像から抽出した測定点の周囲の複数計測点の凹凸パターンを所定時間後の前記画像上で検索して凹凸パターンが最も類似する相似点を検出し，前記相似点及び測定点の凹凸値をＺ座標値に戻した対応三次元座標値の差分を求めてなる変動面の変位計測方法。
【請求項２】
請求項１の方法において，前記計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を輝度値に置換してＸＹ平面に配列した二次元画像を作成し，前記凹凸パターンを輝度パターンとし，前記相似点及び測定点の輝度値をＺ座標値に戻した対応三次元座標値の差分を求めてなる変動面の変位計測方法。
【請求項３】
請求項１又は２の方法において，前記計測点群の三次元座標値に対し前記変動面の平面展開用の座標変換を施した変換座標値（ｘ，ｙ，ｚ）のｚ座標値をｘｙ平面に配列して二次元画像を作成し，前記変換座標系における相似点及び測定点の対応三次元座標値の差分を求めてなる変動面の変位計測方法。
【請求項４】
請求項３の方法において，前記変換座標系における相似点及び測定点の対応三次元座標値に前記座標変換の逆変換を施し，前記逆変換後の対応三次元座標値の差分を求めてなる変動面の変位計測方法。
【請求項５】
請求項１から４の何れかの方法において，前記測定点を変動面上の既知三次元座標値とし，前記相似点の対応三次元座標値と前記測定点の既知三次元座標値との差分を求めてなる変動面の変位計測方法。
【請求項６】
請求項１から５の何れかの方法において，前記特定時刻の画像上の複数の測定点に対する所定時間後の画像上の相似点をそれぞれ検出し，前記複数の相似点及び複数の測定点の対応三次元座標値をそれぞれ結んだ相似線及び測定線の変化を求めてなる変動面の変位計測方法。
【請求項７】
凹凸のある変動面上に散在する計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を凹凸値とみなしてＸＹ平面に配列した二次元画像を作成する画像作成手段，特定時刻の前記画像から抽出した測定点の周囲の複数計測点の凹凸パターンを所定時間後の前記画像上で検索して凹凸パターンが最も類似する相似点を検出する相似点検出手段，及び前記相似点及び測定点の凹凸値をＺ座標値に戻した対応三次元座標値の差分を求める変位計測手段を備えてなる変動面の変位計測システム。
【請求項８】
請求項７のシステムにおいて，前記画像作成手段により計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を輝度値に置換してＸＹ平面に配列した二次元画像を作成し，前記凹凸パターンを輝度パターンとし，前記変位計測手段により前記相似点及び測定点の輝度値をＺ座標値に戻した対応三次元座標値の差分を求めてなる変動面の変位計測システム。
【請求項９】
請求項７又は８のシステムにおいて，前記計測点群の三次元座標値に対し前記変動面の平面展開用の座標変換を施す座標変換手段を含め，前記画像作成手段により座標変換後の変換座標値（ｘ，ｙ，ｚ）のｚ座標値をｘｙ平面に配列して二次元画像を作成し，前記変位計測手段により変換座標系における相似点及び測定点の対応三次元座標値の差分を求めてなる変動面の変位計測システム。
【請求項１０】
請求項９のシステムにおいて，前記変換座標系における相似点及び測定点の対応三次元座標値に前記座標変換の逆変換を施す逆変換手段を含め，前記変位計測手段により逆変換後の対応三次元座標値の差分を求めてなる変動面の変位計測システム。
【請求項１１】
請求項７から１０の何れかのシステムにおいて，前記測定点を変動面上の既知三次元座標値として記憶する記憶手段を設け，前記変位計測手段により相似点の対応三次元座標値と測定点の既知三次元座標値との差分を求めてなる変動面の変位計測システム。
【請求項１２】
請求項７から１１の何れかのシステムにおいて，前記相似点検出手段により特定時刻の画像上の複数の測定点に対する所定時間後の画像上の相似点をそれぞれ検出し，前記変位計測手段により複数の相似点及び複数の測定点の対応三次元座標値をそれぞれ結んだ相似線及び測定線の変化を求めてなる変動面の変位計測システム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は変動面の変位計測方法及びシステムに関し，とくに凹凸のある変動面の経時的な変位を計測する方法及びシステムに関する。
【背景技術】
【０００２】
土木・建築構造物（以下，単に構造物ということがある）を構築し又は管理する際に，安定性・安全性を確認すると共に設計・施工の妥当性を評価するため，周辺の変動しうる地山・地盤その他の自然物又は人工物表面（以下，変動面ということがある）の経時的な変位を計測することがある。例えば山岳トンネルを掘削する場合に，掘削直後の切羽において必要な支保や一次覆工を建て込むと共に，周囲地山の挙動や支保の変形を把握して施工の安全性・支保の妥当性を判断するため，切羽から離れた後方においてトンネル内空断面（変動面）の変位を継続的に計測する施工管理（トンネルＡ計測）が行われる。その変位速度が所定値以下（例えば１ｍｍ／週以下）に収束することで周囲地山が安定したと判断し，変位収束後の内空断面に最終的なトンネル断面となる二次覆工を打設する。
【０００３】
従来から変動面の変位を計測する方法として，例えばトンネル内空断面の天端部位，肩部位，両脚部位等の所定位置にそれぞれアンカーを打ち込んで測定点を設け，その測定点間の距離（測線）を直接スケール等で読み取る方法（測線法）が実施されている。また，各測定点にターゲット（反射板又はプリズム）を取り付け，トンネル底盤の所定観測点に位置決めしたトータルステーション（三次元光波式測距器）で各ターゲットを順次視準して水平角・鉛直角・距離を求めることにより，各測定点の観測点に対する三次元座標及び変位を計測する方法も実施されている（特許文献１，２参照）。トータルステーションによる変位計測法は，非接触式であるため足場等を組み立てる必要がなく，測線法に比べて簡便であることから，現在では標準的な変位計測法となっている。ただし，これらの計測法は何れもトンネル内空断面当たりの測定点の数が３〜５点程度に限られるため，その断面（変動面）の詳細な形状変化を把握することは困難である。
【０００４】
これに対し，例えば図７に示すように，変動面（図示例ではトンネル内空面）上の多数の計測点の三次元座標値から変動面の詳細な形状変化を計測する方法が提案されている（特許文献３参照）。図示例では，トンネル２の切羽２ａの後方底盤２ｂ上に三次元レーザスキャナ５を設置してトンネル内空面（変動面）１を走査することにより，内空面１上の多数の計測点の三次元座標値を取得する。三次元レーザスキャナ５とは，レーザ光の発光体及び受光センサと方位角及び仰角の切り替え装置とが組み合わされたレーザヘッド５ａを例えば三脚上に搭載し，方位角（水平角）θ及び仰角（鉛直角）φを切替えながらレーザ光を計測対象面上で走査すると共にその対象面からの反射光をセンサで検知することにより，レーザ光の往復時間（対象面までの距離ｄ）と方位角θ及び仰角φとから対象面上の各計測点の三次元座標値を取得する装置である。例えば経緯台式又は回転ミラー式の切り替え装置を用いた三次元レーザスキャナにより，トンネル内空面上の１４０万点〜７億点の三次元座標値を迅速に取得することができる。
【０００５】
図７の方法は，トンネル内空面（変動面）１上の３以上の既知位置（地球座標系の位置。以下，トンネル座標値ということがある）にそれぞれターゲット９を取り付け，そのターゲット９を含む内空面１を走査して多数の計測点の三次元座標値（スキャナ装置の座標系の位置。以下，スキャン座標値ということがある）を取得したのち，先ず計測点の中からターゲット９の位置を検出する。例えばターゲット９をレーザ光の高反射シート（全反射シート）又は吸収シート（低反射シート）とすることにより，多数の計測点の中からデータ抜け領域としてターゲット９の位置を検出し，そのスキャン座標値を特定することができる。次いで，特定したターゲット９のスキャン座標値とトンネル座標値との関係に基づき他の計測点のスキャン座標値をそれぞれトンネル座標値に変換し，変換後の各計測点のトンネル座標値に基づきトンネル内空断面の形状を計測する。
【０００６】
図７の方法によれば，トンネル内空面のような変動面の形状を継続的に計測して順次比較することにより，変動面の詳細な形状変化を把握することができる。また，計測した形状を設計形状と比較することにより，トンネル内空断面の特定部位における当たり取りや余掘の要否を判断することもできる。図７のように多数の計測点の三次元座標値から変動面の形状を把握する方法は，三次元レーザスキャナを用いてトンネル内空面の形状を計測する場合だけでなく，例えば図１のように調査船に搭載した音響測深機（ソナーヘッド）を用いて海底面・水底面（変動面）の地形を計測する場合にも適用されており（特許文献４〜６参照），航空機等に搭載した航空レーザスキャナ又はステレオ式デジタルカメラを用いて地表面（変動面）の地形を計測する場合にも適用されている（特許文献７〜９参照）。
【先行技術文献】
【特許文献】
【０００７】
【特許文献１】特開２００１−１６５６５６号公報
【特許文献２】特開２００５−０２４４９２号公報
【特許文献３】特開２０１０−２１７０１７号公報
【特許文献４】特開平０６−０９４４５６号公報
【特許文献５】特開平１０−３２５８７１号公報
【特許文献６】特開平１１−０８３４７９号公報
【特許文献７】特開平０８−１５９７６２号公報
【特許文献８】特開平１０−０８９９５８号公報
【特許文献９】特開２００２−２４３４４４号公報
【非特許文献】
【０００８】
【非特許文献１】ＭＶＴｅｃ  Ｓｏｆｔｗａｒｅ  Ｇｍｂｈほか著「画像処理アルゴリズムと実践アプリケーション」株式会社リンクス，２１０〜２１７頁，２００８年６月１１日発行
【発明の概要】
【発明が解決しようとする課題】
【０００９】
しかし，上述した三次元レーザスキャナ，音響測深機，航空レーザスキャナ，ステレオ式デジタルカメラ等のような多数の計測点群の三次元座標値を取得する装置（以下，これらをまとめて走査装置という）は，基本的に変動面（トンネル内空断面，水底面，地表面等）の形状を多数の計測点群の三次元座標値で把握するものであり，走査密度を高めて変動面の形状の把握精度を高め，更にその経時的観察によって変動面の詳細な形状変化を計測することはできるものの，変動面上の特定位置（測定点）の変位を計測できない問題点がある。すなわち，走査装置を用いた計測では，各計測点が走査密度で決まる離散的・散在的な位置であり，前回走査時に計測点となった変動面上の位置が今回走査時の計測点になるとは限らず，常に同じ位置の三次元座標値が得られるわけではない。これに対しトータルステーション等で計測する変位は，変動面上の同じ位置（測定点）の経時的な変化であり，前回計測と同じ位置で今回計測を行うことが前提となっている。もし離散的・散在的な計測点群の三次元座標値から変動面上の同じ位置（測定点）の変位を計測できれば，走査装置によって変位を計測することが可能になる。また，走査装置によって同じ位置の変位を計測できれば，トータルステーションによる計測時に必要とされるターゲット等が省略できるので，トータルステーションよりも更に簡便な変位計測法として活用が期待できる。
【００１０】
そこで本発明の目的は，変動面上の離散的・散在的な計測点群の三次元座標値から変動面上の特定位置の変位を計測する方法及びシステムを提供することにある。
【課題を解決するための手段】
【００１１】
本発明者は，土木・建築分野で変位の計測が求められる変動面の多くは表面に凹凸を有していることに注目した。また，そのような変動面は多くの場合，時間が経過して全体の形状が変化しても局所的な凹凸形状の特徴を残しており，時間経過の前後で共通する局所的な凹凸形状を手がかりに変動面上の変位を計測できる可能性があることに着目した。すなわち，特定時点における変動面の測定点近傍の特徴的な凹部及び凸部は，時間経過後においても特定時点からそれほど離れていない位置に測定点と一体的に存在しており，測定点と一体的に変位すると考えられる。本発明は，この知見に基づく研究開発の結果，完成に到ったものである。
【００１２】
図１の実施例及び図４の流れ図を参照するに，本発明による変動面の変位計測方法は，凹凸のある変動面１（図示例では水底面１）上に散在する計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を凹凸値とみなしてＸＹ平面に配列した二次元画像Ｉを作成し（ステップＳ１０２及び図６参照），特定時刻Ｔｏの二次元画像Ｉｏ（図６（Ａ）参照）から抽出した測定点Ｐｏの周囲の複数計測点の凹凸パターンＭを所定時間後Ｔｊの二次元画像ｊ（図６（Ｂ）参照）上で検索して凹凸パターンＭが最も類似する相似点Ｐｊを検出し，相似点Ｐｊ及び測定点Ｐｏの凹凸値をＺ座標値に戻した対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）を求めてなるものである（ステップＳ１０５〜Ｓ１０６参照）。
【００１３】
また図１のブロック図及び図２の流れ図を参照するに，本発明による変動面の変位計測システムは，凹凸のある変動面１（図示例では水底面１）上に散在する計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を凹凸値とみなしてＸＹ平面に配列した二次元画像Ｉを作成する画像作成手段２０（図２（Ｂ）参照），特定時刻Ｔｏの二次元画像Ｉｏ（図６（Ａ）参照）から抽出した測定点Ｐｏの周囲の複数計測点の凹凸パターンＭを所定時間後Ｔｊの二次元画像Ｉｊ（図６（Ｂ）参照）上で検索して凹凸パターンＭが最も類似する相似点Ｐｊを検出する相似点検出手段２３（図２（Ｄ）参照），及び相似点Ｐｊ及び測定点Ｐｏの凹凸値をＺ座標値に戻した対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）を求める変位計測手段２５（図２（Ｅ）参照）を備えてなるものである。
【００１４】
画像作成手段２０により，例えば計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を輝度値に置換してＸＹ平面に配列した二次元画像を作成することができ，その場合は凹凸パターンを輝度パターンとする。
【００１５】
好ましい実施例では，図１のブロック図及び図３の流れ図に示すように，計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）に対し変動面１の平面展開用の座標変換Ｒを施す座標変換手段２１を含め（図３（Ｂ）参照），画像作成手段２０により座標変換後の変換座標値（ｘ，ｙ，ｚ）のｚ座標値をｘｙ平面に配列して二次元画像Ｉを作成し（図３（Ｃ）及び図４のステップＳ１０２参照），変位計測手段２５により変換座標系における相似点Ｑｊ及び測定点Ｑｏの対応三次元座標値の差分（Δｘ，Δｙ，Δｚ）を求める（図３（Ｅ）参照）。この場合において必要な場合は，変換座標系における相似点Ｑｊ及び測定点Ｑｏの対応三次元座標値に座標変換Ｒの逆変換Ｖを施す逆変換手段２４を含め（図３（Ｆ）及び図４のステップＳ１０６参照），変位計測手段２５により逆変換後の対応三次元座標値（ΔＸ，ΔＹ，ΔＺ）の差分を求めることができる（図３（Ｇ）参照）。
【００１６】
特定時刻Ｔｏの測定点Ｐｏは変動面１上の既知三次元座標とすることができ，その場合は，測定点Ｐｏを変動面１上の既知三次元座標値して記憶する記憶手段１１を設け，変位計測手段２５により相似点Ｐｊの対応三次元座標値と測定点Ｐｏの既知三次元座標値との差分（ΔＸ，ΔＹ，ΔＺ）を求めることができる。
【００１７】
更に好ましい実施例では，図５に示すように，相似点検出手段２３により特定時刻Ｔｏの二次元画像Ｉｏ上の複数の測定点Ｐｏに対する所定時間後Ｔｊの二次元画像Ｉｊ上の相似点Ｐｊをそれぞれ検出し，変位計測手段２５により複数の相似点Ｐｊ及び複数の測定点Ｐｏの対応三次元座標値をそれぞれ結んだ相似線Ｆｊ及び測定線Ｆｏの変化を求めることができる。
【発明の効果】
【００１８】
本発明による変動面の変位計測方法及びシステムは，凹凸のある変動面１上に散在する計測点群の特定時刻Ｔｏ及び所定時間後Ｔｊにおける三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を凹凸値とみなしてＸＹ平面に配列した二次元画像Ｉｏ，Ｉｊを作成し，特定時刻Ｔｏの二次元画像Ｉｏから抽出した測定点Ｐｏの周囲の複数計測点の凹凸パターンＭを所定時間後Ｔｊの二次元画像Ｉｊ上で検索して凹凸パターンＭが最も類似する相似点Ｐｊを検出し，その相似点Ｐｊ及び測定点Ｐｏの凹凸値をＺ座標値に戻した対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）を求めて変位とするので，次の効果を奏する。
【００１９】
（イ）離散的・散在的な計測点群の三次元座標値から変動面１の変位を直接計測するのではなく，その三次元座標値のＺ座標値を凹凸値とみなしてＸＹ平面上に配列した二次元画像（Ｚ座標値が画素値として埋め込まれた画像）Ｉｏ，Ｉｊを作成し，その二次元画像Ｉｏ，Ｉｊ上で周囲の凹凸パターン（画素値の分布パターン）Ｍが相互に類似する測定点Ｐｏ，相似点Ｐｊを検出し，その測定点Ｐｏ，相似点Ｐｊの対応三次元座標値の差分から変位を計測するので，計測点群の離散性・散在性に拘らず，変動面１上の同じ位置（測定点Ｐｏ）の経時的な変位を計測することができる。
（ロ）計測点群が散在的であっても，二次元画像Ｉｏ，Ｉｊ上で測定点Ｐｏ及び相似点ＰｊのＸＹ座標値を特定することにより，散在する計測点のＺ座標値から測定点Ｐｏ及び相似点ＰｊのＺ座標値を定めることができる。
（ハ）また，周囲の凹凸パターンＭを尺度として検索することにより，従来のパターン認識技術ないし画像解析技術等を利用して特定時刻Ｔｏの測定点Ｐｏと実質上同じ位置とみなせる所定時間後Ｔｊの相似点Ｐｊを精度よく検出することができる。
【００２０】
（ニ）測定点Ｐｏにターゲット等を取り付けて変位を計測することも可能であるが，変動面１上の平坦でない部分であればターゲット等を取り付けていない任意の測定点Ｐｏの変位を計測することが可能であり，ターゲット等の取り付け作業を省略することで変位計測の作業の簡単化・効率化を図ることができる。
（ホ）変動面１がトンネル内空面のように湾曲している場合は，必要に応じて計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）に対し変動面１の平面展開用の座標変換Ｒを施した変換座標値（ｘ，ｙ，ｚ）から二次元画像Ｉを作成し，その変換座標系において相似点Ｑｊ及び測定点Ｑｏの対応三次元座標値の差分（Δｘ，Δｙ，Δｚ）を求めることにより，従来のトータルステーションを用いた計測方法等では困難であった変動面の表面に沿った方向（例えばトンネル内空面に沿ったトンネル軸方向及び周方向）の変位を計測することができる。
【図面の簡単な説明】
【００２１】
以下，添付図面を参照して本発明を実施するための形態及び実施例を説明する。
【図１】は，本発明による変動面の変位計測システムを示すブロック図の一例である。
【図２】は，所定時刻の走査図面Ｉｏ及び所定時間後の走査図面Ｉｊの説明図である
【図３】は，変動面の近似曲面の平面展開用の座標変換及びその逆変換の説明図である
【図４】は，本発明による変動面の変位計測方法を示す流れ図の一例である。
【図５】は，本発明による変動面の変位計測方法の一実施例の説明図である。
【図６】は，本発明による変動面の変位計測の原理を示す説明図である。
【図７】は，従来の三次元レーザスキャナを用いたトンネル内空断面の形状計測の説明図である。
【発明を実施するための形態】
【００２２】
図１は，例えば水底地盤中に地下トンネルを掘削する場合に，水底面を変動面１として本発明の変位計測システムを適用した実施例を示す。図示例のシステムは，水底面を走査して計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）を取得する走査装置５と，その三次元座標値を入力して水底面１上の特定の測定点Ｐｏの変位（ΔＸ，ΔＹ，ΔＺ）を計測するコンピュータ１０とを有する。図示例の走査装置５は，調査船上の所定位置に搭載して所定水平角θ及び鉛直角φのシングルビーム音波又は扇形のマルチビーム音波を水底面に向けて放射するソナーヘッド５ｂと，その調査船の位置Ｏを測量するＧＰＳ測量装置６と，調査船の方位及び姿勢を計測する姿勢計測装置７とを有する音響測深機である。ソナーヘッド５ｂによる音波の往復時間（水底面までの距離ｄ）・水平角θ・鉛直角φと，ＧＰＳ測量装置６による測量位置Ｏと，姿勢計測装置７による方位・姿勢とを計測しながら調査船を移動させることにより，水底面上の計測点群の地球座標系における三次元座標値を取得することができる。必要に応じて走査装置５に水中の温度及び塩分濃度を計測する計測器を含め，その計測値に基づき水中の音速度（例えば水深毎の音速度）を補正することにより三次元座標値の精度を高めることができる。
【００２３】
以下，図示例を参照して本発明を説明するが，本発明は水底面の変位計測への適用に限定されるわけではなく，地表面，トンネル内空面，その他の凹凸のある変動面（自然物又は人工物の表面）の変位計測に広く適用可能である。例えば上述した航空機等に搭載可能なＧＰＳ付きレーザスキャナを用いて地表面上の計測点群の三次元座標値を得ることにより，本発明を地表面の変位計測に適用することができる。また，計測点群の三次元座標値を取得する装置も図示例に限定されるものではなく，ＧＰＳ測量装置に代えて自動追尾トータルステーション，レーザ測距儀等を用いて調査船の位置Ｏを測量することもできる。更に，計測点群の三次元座標値の取得装置を航空測量等で用いるステレオ式デジタルカメラとし，ステレオ式デジタルカメラで撮影した写真画像から複数の計測点群の三次元座標値をステレオ画像法により取得して本発明に適用するこもできる。
【００２４】
図示例のコンピュータ１０は，キーボード・マウス等の入力装置１３と，ディスプレイ・プリンタ等の出力装置１５と，一次記憶装置又は二次記憶装置等の記憶手段１１とを有している。また内蔵プログラムとして，走査装置５から計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）を入力する入力手段１２と，その計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）から二次元画像Ｉを作成する画像作成手段２０を有する。画像作成手段２０は，例えば各計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）のＺ座標値を凹凸値とみなしてＸＹ平面に配列することにより，Ｚ座標値が画素値Ｉ（Ｘ，Ｙ）として埋め込まれた二次元画像Ｉを作成する。或いは，画像作成手段２０により各計測点のＺ座標値を，例えば輝度値，明度値，彩度値，色相値，ＲＧＢ値等（以下，これらをまとめて輝度値という）に置換したうえでＸＹ平面に配列して二次元画像Ｉを作成してもよい（図２（Ｂ）参照）。このような二次元画像Ｉ（以下，二次元輝度画像Ｉということがある）には，例えば色彩等を利用して輝度値を表した画像，単一の色であっても濃淡により輝度値を表した画像が含まれる。例えば，各計測点（Ｘ，Ｙ，Ｚ）のＺ座標値の平均値（又は最大値或いは最小値）を基準輝度値とし，各計測点のＺ座標値を平均値（又は最大値或いは最小値）との差に応じた輝度値へ置換することにより，ＸＹ平面上に散在する各計測点の凹凸高低差が輝度値で表わされた二次元画像Ｉを作成することができる。必要に応じて，Ｚ座標値を輝度値へ置換するための置換式又は置換テーブルを記憶手段１１に記憶しておくことができる。
【００２５】
また図示例のコンピュータ１０は，内蔵プログラムとして，特定時刻Ｔｏに作成した二次元画像Ｉｏ上から測定点Ｐｏの周囲の凹凸パターン（画素値の分布パターン，例えば輝度パターン）Ｍを抽出するパターン抽出手段２２（図６（Ａ）参照）と，そのパターンＭを所定時間後Ｔｊに作成した二次元画像Ｉｊ上で検索して最も類似する相似点Ｐｊを検出する相似点検出手段２３と，検出した相似点Ｐｊと測定点Ｐｏとの対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）を求める変位計測手段２５とを有している。
更に図示例のコンピュータ１０は，内蔵プログラムとして，入力手段１２の入力した三次元座標値，作成手段２０の作成した二次元画像Ｉ，抽出手段２２の抽出したパターンＭ，検出手段２３の検出した相似点Ｐｊ，計測手段２５の求めた変位（ΔＸ，ΔＹ，ΔＺ）等を出力装置１５へ適宜出力する出力手段１４を有している。
【００２６】
図４は，図１のコンピュータ１０により変動面１（この場合は水底面１）の変位を計測する方法の流れ図を示す。以下，図４の流れ図を参照して図１のコンピュータ１０の各内蔵プログラムの機能を説明する。先ずステップＳ１０１において，走査装置５により変動面（図示例では水底面）１を走査して多数の計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）を取得し，取得した計測点群の三次元座標値を入力手段１２経由でコンピュータ１０に入力する。ステップＳ１０２において，各計測点の三次元座標値（Ｘ，Ｙ，Ｚ）を画像作成手段２０へ入力して二次元画像Ｉを作成する（図２（Ｂ）参照）。
【００２７】
図４のステップＳ１０３において，変位計測に必要な凹凸パターンＭが抽出済みであるか否かを判断し，未抽出のときはステップＳ１０４へ進み，図２（Ｃ）及び図６（Ａ）に示すように，パターン抽出手段２２により特定時刻Ｔｏの二次元画像Ｉｏから測定点Ｐｏ＝（Ｘｏ，Ｙｏ）の周囲の凹凸パターンＭを抽出する。パターンＭは，二次元画像Ｉｏ上の任意位置に測定点Ｐｏを設定して抽出することができる。測定点Ｐｏは必ずしもステップＳ１０１で入力した計測点の何れかと一致させる必要はなく，二次元画像Ｉｏ上に散在する計測点の中間位置に測定点Ｐｏを設定することも可能であり，そのような測定点ＰｏのＺ座標値（凹凸値，輝度値）を隣接する計測点のＺ座標値から按分（平均操作）により定めることもできる。例えば，特定時刻Ｔｏの二次元画像Ｉｏを出力装置１５へ出力し，その二次元画像Ｉｏ上で指定された測定点Ｐｏを入力手段１２経由で入力してパターンＭを抽出してもよい。パターン抽出手段２２で抽出したパターンＭは，図１に示すように記憶手段１１に記憶しておくことができる。
【００２８】
ステップＳ１０４において，測定点Ｐｏとする変動面１上の位置に予めターゲット等を取り付けておく必要はないが，ターゲット等を取り付けた変動面１上の位置を測定点Ｐｏとして凹凸パターンＭを抽出することも可能である。その場合は，例えば図１に点線で示すように，ターゲット等を取り付けた変動面１上の測定点Ｐｏの三次元座標値（Ｘｏ，Ｙｏ，Ｚｏ）を予め求めてコンピュータ１０の記憶手段１１に記憶しておき，その測定点Ｐｏ＝（Ｘｏ，Ｙｏ）の周囲からパターン抽出手段２２によりパターンＭを抽出する。パターンＭを抽出したのちステップＳ１０１へ戻り，上述したステップＳ１０１〜Ｓ１０２を繰り返す。例えば所定時間後Ｔｊに変動面１を走査して取得した計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）を入力し，所定時間後Ｔｊの二次元画像Ｉｊを作成する。
【００２９】
ステップＳ１０３において凹凸パターンＭが抽出済であるときはステップＳ１０５へ進み，図２（Ｄ）及び図６（Ｂ）に示すように，相似点検出手段２３によって所定時間後Ｔｊの二次元画像Ｉｊ上で凹凸パターンＭを検索し，周囲のパターンが測定点ＰｏのパターンＭと最も類似する相似点Ｐｊ＝（Ｘｊ，Ｙｊ）を二次元画像Ｉｊから検出する。上述したように変動面１上の測定点Ｐｏは多くの場合に周囲の凹部及び凸部と一体的に移動すると考えられるので，所定時間後Ｔｊの二次元画像Ｉｊ上で凹凸パターンＭの最も類似する相似点Ｐｊは，その所定時間後Ｔｊにおける測定点Ｐｏの移動先位置とみなすことができる。
【００３０】
ステップＳ１０５において，凹凸パターンＭの類似度を尺度として所定時間後Ｔｊの二次元画像Ｉｊを検索することにより，従来のパターン認識技術等を利用して測定点Ｐｏの移動先位置とみなせる相似点Ｐｊを精度よく検出できる。二次元画像Ｉｊから検出される相似点Ｐｊは，必ずしもステップＳ１０１で入力した計測点の何れかと一致せず，二次元画像Ｉｊ上に散在する計測点の中間位置となることもあるが，そのような相似点ＰｊのＺ座標値（凹凸値，輝度値）は隣接する計測点のＺ座標値からの按分（平均操作）により定めることができる。変動面１上の計測点の密度に拘らず，測定点Ｐｏの移動先位置とみなせる二次元画像Ｉｊ上の相似点Ｐｊを精確に特定することにより，以下のステップＳ１０６において同一位置（測定点Ｐｏ）の変位を求めることが可能となる。
【００３１】
ステップＳ１０５の相似点検出手段２３は，例えばテンプレートマッチング法を利用して相似点Ｐｊを検出することができる。すなわち，例えば二次元輝度画像Ｉ上における凹凸パターン（輝度パターン）Ｍを所要大きさ（画素数ｎ）のテンプレートｔ（ｕ，ｖ）とし，そのテンプレートｔ（ｕ，ｖ）を（１）式に示すように二次元輝度画像Ｉｊ上の全ての点（ｒ，ｃ）上に移動させながら各点においてテンプレートｔ（ｕ，ｖ）との類似度ｓａｄ（ｒ，ｃ）（テンプレート全体にわたる輝度差の絶対値の総和）を算出し，その類似度ｓａｄが最も小さく（近く）なる相似点Ｐｊを検出することができる。また，（２）式のようにテンプレート全体にわたる輝度差の二乗和を類似度ｓｓｄとし，又は（３）式のように正規化相互相関を類似度ｎｃｃとするテンプレートマッチング法によって相似点Ｐｊを検出することも可能である（非特許文献１参照）。なお，（３）式の正規化相互相関における符合ｍはテンプレートｔ及び二次元画像Ｉｊの平均輝度を表し，符合ｓ２はテンプレートｔ及び二次元画像Ｉｊの輝度分散を表す（（４）式及び（５）式参照）。
【００３２】
0005943510.tif 000002
【００３３】
ただし，相似点検出手段２３による相似点Ｐｊの検出方法はテンプレートマッチング法に限定されるものではなく，凹凸パターンＭと二次元画像Ｉｊとの類似度を評価する他のパターン認識技術を利用して相似点Ｐｊを検出することができる。また，パターンＭの大きさ（面積）は，計測対象である変動面１の表面凹凸形状に応じて適宜選択できる。例えば変動面１の表面凹凸形状が比較的均一であるときは比較的大きなパターンＭを用いることで相似点Ｐｊの検出精度を高めることができ，変動面１の表面凹凸形状が位置によって相違しているときは比較的小さなパターンＭでも相似点Ｐｊを精度よく検出することができる。そのような変動面１に応じたパターンＭの大きさ（面積）は，予め試験的に求めて記憶手段１１に登録しておくことができる。
【００３４】
図４のステップＳ１０６は，ステップＳ１０５で検出した所定時間後Ｔｊの相似点Ｐｊ＝（Ｘｊ，Ｙｊ）に基づき，変位計測手段２５によって測定点Ｐｏ＝（Ｘｏ，Ｙｏ）からの変位を求める処理を示す（図２（Ｅ）参照）。例えば，ステップＳ１０１で入力した各計測点の三次元座標値を記憶手段１１に記憶しておき，ステップＳ１０６において測定点Ｐｏ及び相似点Ｐｊに対応する計測点の三次元座標値（Ｘｏ，Ｙｏ，Ｚｏ），（Ｘｊ，Ｙｊ，Ｚｊ）が存在するときは，その計測点の対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）＝（Ｘｊ−Ｘｏ，Ｙｊ−Ｙｏ，Ｚｊ−Ｚｏ）を測定点Ｐｏの変位とすることができる。或いは，測定点Ｐｏ及び相似点Ｐｊと一致する計測点が存在しない場合に，測定点Ｐｏ及び相似点Ｐｊに最も隣接する計測点をそれぞれ対応する計測点とし，その計測点の対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）を測定点Ｐｏの変位としてもよい。上述したように相似点Ｐｊは測定点Ｐｏの移動先とみなせるので，その対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）から所定時間後Ｔｊにおける測定点Ｐｏの三次元変位を求めることができる。
【００３５】
好ましくは，測定点Ｐｏ及び相似点Ｐｊに対応する計測点の対応三次元座標値に代えて，ステップＳ１０６において変位計測手段２５により，測定点Ｐｏ及び相似点Ｐｊの凹凸値（例えば輝度値）を隣接する計測点の凹凸値（例えば輝度値）から按分（平均操作）により定めたうえで，その凹凸値（例えば輝度値）をＺ座標値に戻して対応三次元座標値（Ｘｏ，Ｙｏ，Ｚｏ），（Ｘｊ，Ｙｊ，Ｚｊ）を求め，その対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）＝（Ｘｊ−Ｘｏ，Ｙｊ−Ｙｏ，Ｚｊ−Ｚｏ）を算出する。測定点Ｐｏ及び相似点Ｐｊの凹凸値又は輝度値をＺ座標値に戻した対応三次元座標値の差分を求めることにより，測定点Ｐｏ及び相似点Ｐｊが計測点と一致しない場合でも，変動面１上の同じ特定位置（測定点Ｐｏ）の三次元変位（ΔＸ，ΔＹ，ΔＺ）を求めることができる。
【００３６】
ステップＳ１０６において測定点Ｐｏの変位（ΔＸ，ΔＹ，ΔＺ）を求めたのち，ステップＳ１０８において変位計測を終了するか否かを判断し，計測を継続する場合はステップＳ１０１へ戻って上述したステップＳ１０１〜Ｓ１０６を繰り返す。すなわち，更に所定時間経過Ｔ（ｊ＋１）の計測点群の三次元座標値をコンピュータ１０に入力して二次元画像Ｉ（ｊ＋１）を作成し，その二次元画像Ｉ（ｊ＋１）上において凹凸パターンＭを検索して相似点Ｐ（ｊ＋１）＝（Ｘ（ｊ＋１），Ｙ（ｊ＋１），Ｚ（ｊ＋１））を検出し，その相似点Ｐ（ｊ＋１）及び測定点Ｐｏの対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）＝（Ｘ（ｊ＋１）−Ｘｏ，Ｙ（ｊ＋１）−Ｙｏ，Ｚ（ｊ＋１）−Ｚｏ）によって測定点Ｐｏの変位を求める。ステップＳ１０１〜Ｓ１０８のサイクルを繰り返すことにより，変動面１上の同じ特定位置（測定点Ｐｏ）の経時的な変位を順次計測できる。
【００３７】
本発明は，離散的・散在的な計測点群の三次元座標値から変動面１の変位を直接計測するのではなく，計測点群の三次元座標値からＺ座標値をＸＹ平面上に配列した二次元画像（Ｚ座標値が画素値として埋め込まれた画像）Ｉｏ，Ｉｊを作成し，二次元画像Ｉｏ上の測定点Ｐｏの周囲画素値のパターンＭとの類似度を尺度として二次元画像Ｉｊ上の相似点Ｐｊを検出し，その相似点Ｐｊと測定点Ｐｏとの対応三次元座標値の差分から変位を計測するので，計測点群の離散性・散在性に拘らず，変動面１上の特定位置（測定点Ｐｏ）の変位を計測することができる。また，測定点Ｐｏ及び相似点Ｐｊが計測点と一致しない場合でも，計測点の散在する二次元画像Ｉｏ，Ｉｊ上で測定点Ｐｏ及び相似点ＰｊのＸＹ座標値を特定することにより，隣接する計測点のＺ座標値から測定点Ｐｏ及び相似点ＰｊのＺ座標値を精度よく定めて対応三次元座標値の差分を求めることができる。
【００３８】
こうして本発明の目的である「変動面上の離散的・散在的な計測点群の三次元座標値から変動面上の特定位置の変位を計測する方法及びシステム」の提供を達成することができる。
【実施例１】
【００３９】
図３は，トンネル２の内空面を変動面１として本発明の変位計測システムを適用し，図７と同様の三次元レーザスキャナ５を用いて変動面１上の計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）を取得する実施例を示す。例えば変動面１に臨むトンネル内側の所定トンネル座標位置Ｏに三次元レーザスキャナ５を設置して各計測点の地球座標系における三次元座標値を取得するが，その設置位置Ｏのトンネル座標値は，例えばトンネル２内の複数の既知位置に設けた測量基準点８に基づき定めることができる。或いは，図７を参照して上述したように，三次元レーザスキャナ５のスキャン範囲内の３以上の既知位置にそれぞれターゲット９を取り付け，そのターゲット９を含む対象域を走査して各計測点の三次元座標値（スキャン座標値）を取得したのち，計測点から抽出したターゲット９のスキャン座標値と既知座標値（トンネル座標値）との関係に基づきレーザスキャナ５の設置位置Ｏの三次元座標値を算出することも可能である。
【００４０】
図３のように変動面１が湾曲している場合は，上述した図２（Ｂ）の場合と同様に計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）をＸＹ平面に配列して二次元画像Ｉを作成することも可能であるものの，湾曲部分の凹凸高低差が反映された二次元画像Ｉを作成することが難しくなる。このため，図３（Ｂ）〜（Ｃ）に示すように，計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）に対して変動面１の平面展開用の座標変換Ｒを施したうえで，座標変換後の変換座標値（ｘ，ｙ，ｚ）のｚ座標値をｘｙ平面に配列して二次元画像Ｉを作成することが望ましい。図１のコンピュータ１０は，走査装置５から入力した計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）に対し変動面１の平面展開用の座標変換Ｒを施す座標変換手段２１を有している。また図４のステップＳ１０２は，画像作成手段２０において計測点群の三次元座標値（Ｘ，Ｙ，Ｚ）から二次元画像Ｉを作成する際に，必要に応じて座標変換手段２１により三次元座標（Ｘ，Ｙ，Ｚ）に対して座標変換Ｒを施し，座標変換後の変換座標値（ｘ，ｙ，ｚ）から二次元画像Ｉを作成することを示している。
【００４１】
図３（Ｂ）で用いる座標変換Ｒは，例えばトンネル内空面１の近似曲面１ｓである円柱面又は円錐面を伸び縮みなく平面上に展開するものである。変動面１の三次元形状が複雑であっても，柱面，錐面その他の可展面（伸縮なく平面上に展開できる曲面）又はそれらを組合せた曲面１ｓによって近似することにより，可展面を展開する座標変換Ｒによって変動面１上の各計測点の三次元座標値を平面上に展開することができる。そのような変動面１に応じた座標変換式Ｒは，予め作成して記憶手段１１に登録しておくことができる（図１参照）。図３（Ｃ）において，座標変換Ｒによる座標変換後の変換座標値（ｘ，ｙ，ｚ）を画像作成手段２０へ入力し，そのｚ座標値をｘｙ平面に配列することにより，各計測点の凹凸高低差が表わされた二次元画像Ｉを作成する。
【００４２】
図３（Ｄ）は，ｘ軸方向をトンネル軸方向と一致させ，ｙ軸方向をトンネル周方向と一致させ，ｘｙ平面に対する凹凸高低差を表わした特定時刻Ｔｏの二次元画像Ｉｏを示す。図４のステップＳ１０４において， パターン抽出手段２２により二次元画像Ｉｏから測定点Ｑ１ｏ＝（ｘ１ｏ，ｙ１ｏ），測定点Ｑ２ｏ＝（ｘ２ｏ，ｙ２ｏ）の周囲の凹凸パターンＭ１，Ｍ２をそれぞれ抽出する。同図に示すように，パターン抽出手段２２で抽出するパターンＭは１個に限らず，二次元画像Ｉｏ上に複数の測定点Ｑを設定してその周囲からそれぞれパターンＭを抽出することができる。抽出したパターンＭは，測定点Ｑ毎に区別して記憶手段１１に記憶しておく。
【００４３】
また図３（Ｅ）は，図３（Ｄ）と同様にｘ軸方向をトンネル軸方向と一致させ，ｙ軸方向をトンネル周方向と一致させた所定時間後Ｔｊの二次元画像Ｉｊを示す。図４のステップＳ１０５において，相似点検出手段２３により二次元画像Ｉｊ上で凹凸パターンＭ１，Ｍ２をそれぞれ検索し，周囲パターンがパターンＭ１，Ｍ２と最も類似する相似点Ｑ１ｊ＝（ｘ１ｊ，ｙ１ｊ），Ｑ２ｊ＝（ｘ２ｊ，ｙ２ｊ）を二次元画像Ｉｊからそれぞれ検出する。上述したように二次元画像Ｉｊ上でパターンＭ１，Ｍ２の最も類似する相似点Ｑ１ｊ，Ｑ２ｊは，それぞれ所定時間後Ｔｊにおける測定点Ｑ１ｏ，Ｑ２ｏの移動先とみなすことができる。
【００４４】
そののち図４のステップＳ１０６において，変位計測手段２５により，測定点Ｑ１ｏ及び相似点Ｑ１ｊの対応三次元座標値（ｘ１ｏ，ｙ１ｏ，ｚ１ｏ），（ｘ１ｊ，ｙ１ｊ，ｚ１ｊ）を求め，変換座標系における対応三次元座標値の差分（Δｘ１，Δｙ１，Δｚ１）＝（ｘ１ｊ−ｘ１ｏ，ｙ１ｊ−ｙ１ｏ，ｚ１ｊ−ｚ１ｏ）を求める。また同様に，測定点Ｑ２ｏ及び相似点Ｑ２ｊの対応三次元座標値の差分（Δｘ２，Δｙ２，Δｚ２）＝（ｘ２ｊ−ｘ２ｏ，ｙ２ｊ−ｙ２ｏ，ｚ２ｊ−ｚ２ｏ）を求める。図３の二次元画像Ｉは変動面１を平面上に座標変換して展開したものであり，変換座標系における測定点Ｑｏ及び相似点Ｑｊの対応三次元座標値の差分（Δｘ，Δｙ，Δｚ）は変動面１の表面に沿った方向及び垂直方向の移動量を示している。従って，ｘ軸方向をトンネル軸方向と一致させ，ｙ軸方向をトンネル周方向と一致させ，ｚ軸方向を凹凸方向と一致させた図３の二次元画像Ｉｏ，Ｉｊから対応三次元座標値の差分（Δｘ１，Δｙ１，Δｚ１），（Δｘ２，Δｙ２，Δｚ２）を求めることにより，変動面１上の特定位置（測定点Ｑ１ｏ，Ｑ２ｏ）のトンネル軸方向，トンネル周方向，および凹凸方向の変位をそれぞれ計測することができる。
【００４５】
また図３（Ｆ）〜（Ｇ）に示すように，図４のステップＳ１０６において，変換座標系における相似点Ｑｊ及び測定点Ｑｏの対応三次元座標値に座標変換Ｒの逆変換Ｖを施し，逆変換後の対応三次元座標値（ΔＸ，ΔＹ，ΔＺ）の差分を求めることにより，変動面１の座標系（例えばトンネル座標系）における変位（ΔＸ，ΔＹ，ΔＺ）を求めることも可能である。図１のコンピュータ１０は，変換座標系の相似点Ｑｊ及び測定点Ｑｏの対応三次元座標値に対して座標変換Ｒの逆変換Ｖを施す逆変換手段２４を有している。また図４のステップＳ１０６は，変位計測手段２５において変動面１上の変位を求める際に，必要に応じて逆変換手段２４により相似点Ｑｊ及び測定点Ｑｏの対応三次元座標値に対して座標変換Ｒの逆変換Ｖを施し，逆変換後の対応三次元座標値（ΔＸ，ΔＹ，ΔＺ）の差分から変動面１上の変位を計測することを示している。
【００４６】
図３（Ｆ）で用いる逆座標変換Ｖは，図３（Ｂ）の座標変換Ｒと逆向きに二次元画像Ｉを例えばトンネル内空面１の近似曲面１ｓである円柱面又は円錐面に戻すものであり，上述した座標変換式Ｒと同様に予め作成して記憶手段１１に登録しておくことができる。この逆座標変換Ｖにより，二次元画像Ｉの変換座標系が変動面１のトンネル座標系に変換されるので，逆変換後の相似点Ｐｊ及び測定点Ｐｏの対応三次元座標値の差分（ΔＸ，ΔＹ，ΔＺ）を求めることにより，従来のトータルステーションを用いた場合と同様の変位計測（トンネル座標系における変位計測）が可能となる。上述したように，本発明の変位計測では測定点Ｐｏに予めターゲット等を取り付けておく必要がなく，平坦でない部分であれば変動面１上の任意の測定点Ｐｏの変位を計測できるので，図３（Ａ）〜（Ｇ）の流れ図により，ターゲット等の取り付け作業を省略して作業の簡単化・効率化を図りつつ，従来のトータルステーションと同様の変位計測を行うことができる。
【実施例２】
【００４７】
図４のステップＳ１０７は，例えば図５に示すように，特定時刻Ｔｏの複数の測定点Ｐ１ｏ，Ｐ２ｏ，……，Ｐ１２ｏの対応三次元座標値を結んだ測定線Ｆｏと，所定時間後Ｔｊの複数の相似点Ｐ１ｊ，Ｐ２ｊ，……，Ｐ１２ｊの対応三次元座標値を結んだ相似線Ｆｊとの変化を求める実施例を示す。例えば図７に示す従来のトンネル掘削工事の施工管理（トンネルＡ計測）において，切羽２ａから所定距離Ｄだけ離れた後方のトンネル内空断面（変動面）の形状変化（変状）を計測することが求められている。従来のトータルステーションを用いた方法では，計測対象のトンネル内空断面に沿った複数の測定点Ｐｏにそれぞれターゲットを取り付けて変位を計測し，所定時間Ｔｊ毎に複数の変位点Ｐｊを結ぶ断面形状から形状変化（変状）を求めているが，断面上の測定点Ｐｏの数が４点程度に限られるため，断面形状の変状を精度よく把握できない問題点が指摘されていた。
【００４８】
図５の実施例では，例えば図３（Ｄ）に示す特定時刻Ｔｏの二次元画像Ｉｏにおいて複数の測定点Ｐ１ｏ，Ｐ２ｏ，……，Ｐ１２ｏを設定してそれぞれ周囲のパターンＭ１，Ｍ２，……，Ｍ１２を抽出し，図３（Ｅ）に示す所定時間後Ｔｊの二次元画像Ｉｊ上において各パターンＭ１，Ｍ２，……，Ｍ１２と最も類似する相似点Ｐ１ｊ，Ｐ２ｊ，……，Ｐ１２ｊをそれぞれ検出する。変位計測手段２５において，複数の相似点Ｐ１ｊ，Ｐ２ｊ，……，Ｐ１２ｊの対応三次元座標値を結んだ相似線Ｆｊと，複数の測定点Ｐ１ｏ，Ｐ２ｏ，……，Ｐ１２ｏの対応三次元座標値を結んだ測定線Ｆｏをそれぞれ求め，相似線Ｆｊの測定線Ｆｏに対する三次元的な変化によってトンネル内空断面の変状を求める。上述したように，本発明において変位を追跡できる測定点Ｐｏの数にとくに制限はなく，平坦でない部分であれば二次元画像Ｉｏ上に任意数の測定点Ｐｏを設定して変位を追跡できる。従って，図５に示す実施例によれば，ターゲット等の取り付け作業の簡単化・効率化を図りつつ，同時にトンネル内空断面の変状計測の精度向上を図ることができる。また，本発明により求められる相似線Ｆｊの測定線Ｆｏに対する変状は，所定の断面における形状変化ではなく，断面における複数の測定点の変位である。つまり，トンネル軸方向の変位も含んだ変状を計測することができる。従って，従来のような断面における形状変化ではなく，より精度の高い変位計測が可能となる。
【符号の説明】
【００４９】
１…凹凸のある変動面（トンネル内空面）  １ｓ…近似曲面
２…構造物（トンネル）                  ２ａ…切羽
２ｂ…底盤                              ５…走査装置
５ａ…レーザヘッド                      ５ｂ…ソナーヘッド
６…ＧＰＳ測量装置                      ７…姿勢（方位）計測装置
８…測量基準点                          ９…ターゲット
１０…コンピュータ                      １１…記憶手段
１２…入力手段                          １３…入力装置
１４…出力手段                          １５…出力装置
２０…画像作成手段                      ２１…座標変換手段
２２…パターン抽出手段                  ２３…相似点検出手段
２４…逆変換手段                        ２５…変位計測手段
Ｆｏ…測定線                            Ｆｊ…相似線
Ｉ…走査図面                            Ｍ…凹凸パターン（輝度パターン）
Ｏ…走査装置の位置                      Ｐｏ…測定点
Ｐｊ…相似点                            Ｑｏ…（座標変換後の）測定点
Ｑｊ…（座標変換後の）相似点            Ｒ…座標変換式
Ｔ…時刻                                Ｖ…逆変換式
【図１】
0005943510.tif 000003
【図４】
0005943510.tif 000004
【図２】
0005943510.tif 000005
【図３】
0005943510.tif 000006
【図５】
0005943510.tif 000007
【図６】
0005943510.tif 000008
【図７】
0005943510.tif 000009