JP3501936B2 - 変位計測方法及び変位計測装置 - Google Patents

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、トンネル内若しくはそ
の周辺の岩盤、地下に掘削若しくは開削された岩盤、土
質若しくは覆工壁面、支保、地下構造物、斜面、地表露
頭、盛土、ダム、又は地上構造物等の計測対象を撮影し
た写真画像により変位を計測する変位計測方法に関す
る。

【０００２】

【従来の技術及び発明が解決しようとする課題】例えば
各種土木構造物を施工するに当たっては、設計値や設計
方法の考察、施工技術へのフィードバック、更には安全
管理を行うために、様々な物性値の計測が行われる。な
かでも変位計測は、基礎地盤や構造物の挙動を明らかに
するために重要な計測である。

【０００３】トンネル工事においてよく用いられるＮＡ
ＴＭ工法（Ｎｅｗ Ａｕｓｔｒｉａｎ Ｔｕｎｎｅｌｌ
ｉｎｇ Ｍｅｔｈｏｄ）を例に挙げると、坑内観察調
査、壁面の相対変位及び天端沈下の計測等、或いは地中
変位の計測等は、Ａ或いはＢ計測などと称され、トンネ
ルと周辺地山の安定、施工法の改善、更には地表面沈下
や近接構造物への影響を検討するためのものとして重要
である。

【０００４】また、例えば市街地でのトンネル施工時で
は、地上の構造物への影響を考慮する必要があり、地下
だけでなく地上の構造物の変位計測も土木構造物の施工
には重要なものとなる場合が多い。

【０００５】更に、山留め工事や擁壁の工事でも山留め
壁や擁壁の変形を計測することは盛んに行われており、
道路盛土工事における盛土の変形や地盤の沈下などの変
位、或いはダム施工時の構造物の変位計測など、各種地
下及び地上構造物の工事における調査、設計、施工及び
管理における変位計測の例を数多く挙げることができ
る。

【０００６】前記変位の計測には、沈下計、傾斜計、地
中ひずみ計、伸縮計等が用いられ、更にＮＡＴＭなどの
トンネル工事や岩盤斜面の壁面の変位計測には、レベ
ル、コンバージェンスメータ、エクステンションメー
タ、光波測距機、及びトータルステーションといった機
器を用いた測量などが広く用いられている。

【０００７】これらの機器と計測手法を用いた変位計測
においては、以下の項目が問題となっている。目的と
する位置にて計測が可能か否か。すなわち、機器が実際
に設置若しくは取り付け可能であるかどうかだけでな
く、機器の耐久性や計測値に与える作業環境の影響がど
の程度生じるかを検討する必要がある。例えば漏水や落
雷などの不測の事態や、断線などの機器の故障が生じる
危険性がないのかどうかを考慮に入れる必要がある。

【０００８】また、例えば落石を予測するための当該石
の動態観測や、人が近づき得ない場所の対象物を計測す
る際には、計測機器が設置できず、測定不可能であるこ
とが多い。経済性の問題はないか。計測点が多数必要
な場合は、特に計測機器のコストや計測に要する人的費
用が高くなる場合が多い。また、自動的に計測値が得ら
れるか否かも時間的、人的な経済性を左右する。ＮＡＴ
Ｍ工事などにおいて前記計測機器を用いる場合、測定作
業自体のコストも高くなるため、多数の計測点を安価・
迅速に計測するのは困難である。計測値の精度が十分
であるか否か。特定の計測点での精度が十分であって
も、その地点のみの値であり、面的な広がりをもって計
測される訳ではない。面的な広がりのある変位分布につ
いてよい精度が要求される場合には、多数の計測点の変
位を精度良く計測することが必要となり、経済性を考慮
した場合、施工・管理に反映させる上で十分なものとは
ならない場合が多い。

【０００９】これまで用いられてきた計測機器や計測方
法は、前記問題点の何れかを有している。例えば、内空
変位測定に用いるコンバージェンスメジャは測定前のボ
ルトの設置及び測定作業そのものに人的、経済的、時間
的コストがかかる。また遠方からの天端沈下測定に用い
る光波測距機は、機器費用が高く、移動などに難が多
い。すなわち、岩盤等における任意の計測点を計測し得
て面的な広がりのある変位分布を精度良く然も安価・迅
速に計測できる変位計測方法及び変位計測装置は未だ開
発されていないのが実状である。

【００１０】本発明は、従来技術に存した上記のような
問題点に鑑み行われたものであって、その目的とすると
ころは、岩盤等の計測対象における計測点を任意に選ぶ
ことができ、従来の機器では設置し難かった場所や人が
近づくことのできない対象物でもその変位を求めること
ができ、低コストで計測対象における面的変位分布を容
易且つ迅速に精度良く求めることができる変位計測方法
を提供することにある。

【００１１】

【課題を解決するための手段】本発明の変位計測方法
は、計測対象の変位の前後それぞれにおいて、内部標定
要素が未知の又は十分に検定されていないデジタルカメ
ラを用いて複数の位置から各基準計測点を含む計測対象
を撮影して複数の写真画像を得、変位の前後における写
真画像上の基準計測点の座標に基づき、セルフキャリブ
レーション法を使用したバンドル法を利用して最小二乗
法を繰り返し実行することにより標定要素及び基準計測
点の三次元座標を算出し、それらの標定要素及び基準計
測点の三次元座標に基づき任意の変位計測点の三次元座
標をそれぞれ算出してその変位を求める変位計測方法で
あって、前記最小二乗法実行の繰り返しにおいて、写真
画像における基準計測点座標の観測値に対し、次式

【数５】 で表される、異常値を探索するためのウエイト関数によ
るウエイトを与えて最小二乗法を実行し、その計算結果
と写真画像における基準計測点座標の観測値との残差を
計算し、しきい値を設定してそのしきい値より小さいウ
エイトに対応する観測値を異常値と認定し、認定された
異常値を除去して再度最小二乗法を実行することにより
標定要素及び基準計測点の三次元座標の精度を高める工
程を有することを特徴とする。

【００１２】上記ウエイトは、ウエイト関数、すなわち
異常値を探索するためのウエイト関数、例えば下記式
［１］で表されるウエイト関数ｐｉ又は下記式［２］で
表されるウエイト関数（Weight）によるウエイトとする
ことができる。

【数４】 ・・・・［１］=（１／［近似値に予想される外的誤差／スケール／画像精
度］）^２ ・・・・［２］ 上記セルフキャリブレーション法の使用においては、最
小二乗法の正規方程式の拘束条件として既知距離を導入
することにより精度を向上させるものとすることが望ま
しい。この既知距離の導入は、距離拘束体における既知
距離の端点を含めて上記写真画像を得ることにより行う
ものとすることが好ましい。また、これらの場合におけ
る上記標定要素及び基準計測点の三次元座標の算出は、
既知距離の端点の座標に関する拘束条件に対応する行列
を並べ替え、既知距離の端点の座標と基準計測点の座標
を分けて計算することにより行うものとすることが好ま
しい。

【００１３】計測対象は、例えば、トンネル内若しくは
その周辺の岩盤、地下に掘削若しくは開削された岩盤、
土質若しくは覆工壁面、支保、地下構造物、斜面、地表
露頭、盛土、ダム、又は地上構造物等を任意に選択し得
る。

【００１４】

【発明の実施の形態】岩盤等における変位の計測は、変
位の前後における写真画像上の基準計測点の座標に基づ
き任意の変位計測点（基準計測点を含む）の三次元座標
をそれぞれ算出することにより行い得る。基準計測点
は、複数の位置（計測対象である被写体に対する角度の
異なる複数の位置）から撮影したものであることを要す
る。この写真画像は、基準計測点及び変位計測点を含む
被写体を撮影して得るものであってもよく、予め基準計
測点及び変位計測点を含む被写体が撮影されたものを用
いるものであってもよい。

【００１５】写真画像上の座標に基づく任意の変位計測
点の三次元座標の算出は、 （１） 写真画像座標の計測［写真画像における被写体
上の所要の基準計測点の座標の計測］ （２） 写真画像の近似解析［写真撮影の標定要素及び
基準計測点の三次元座標についての精度の良い近似値の
取得］ （３） 写真画像の精密解析［標定要素及び基準計測点
の三次元座標についての近似値に基づく写真画像の精密
解析による標定要素及び基準計測点の三次元座標の算
出、並びに必要な任意の変位計測点の三次元座標の算
出］ の３段階からなる。

【００１６】標定要素としては、写真主点位置の二次元
座標及び画面距離等の内部標定要素、並びに、被写体空
間における適宜の座標系における写真撮影時の撮影機の
位置及び傾きに関する３つの移動量及び３つの回転量等
の外部標定要素を例として挙げることができる。内部標
定要素が既知であれば、それ以外の標定要素を求めるこ
ととなる。

【００１７】変位計測（三次元座標算出）の対象は、岩
盤に限られない。例えば種々の構造物を対象とすること
もできる。

【００１８】写真画像は、例えばＣＣＤ（Ｃｈａｒｇｅ
Ｃｏｕｐｌｅｄ Ｄｅｖｉｃｅ）等を撮像素子とする
デジタルスチルカメラやデジタルビデオカメラを用いた
撮影等により得られるデジタル情報であることが最も望
ましいが、例えばアナログスチルカメラやアナログビデ
オカメラを用いた撮影等により得られるアナログ情報を
デジタル情報に変換したものを用いることもできる。デ
ジタル情報として得た写真画像は、コンピュータで直接
取り扱うことができる。なお、計測精度等の必要に応
じ、前記デジタルカメラ等の撮影機の解像度、撮影距
離、撮影角度及び１つの基準計測点を撮影する数等を任
意に選択することができる。また、紙やプラスチックフ
ィルム等に形成した写真画像のハードコピーを直接用い
ることもでき、そのようなハードコピーを画像スキャナ
等によりデジタル情報化して用いることもできる。

【００１９】なお、これらのデジタル情報及びアナログ
情報は、記録媒体に記録しておくことができることは勿
論である。 （１）写真画像座標の計測［写真画像における被写体上
の所要の基準計測点の座標（二次元座標）の計測］ 写真画像における被写体上の所要の基準計測点の座標の
計測は、写真画像に関するデジタル情報によって行うこ
とが望ましい。複数の写真画像における同一の基準計測
点の座標の計測を精度良く行うには、写真撮影の際に、
例えば、種々の図形、記号及び色彩等の何れか若しくは
それらの組合せからなる標識を、被写体上の計測点に直
接描いたり、板状体やシート状体等に描いて被写体上の
基準計測点に固定することにより、写真画像における被
写体上の基準計測点の位置を識別し易くすることが好ま
しい。このような標識は、変位の前後にわたり被写体に
固定されていることが望ましい。なお、前記複数の写真
画像というのは、例えば、変位前に撮影された同一の計
測点を含む複数の写真画像、変位後に撮影された同一の
基準計測点を含む複数の写真画像、及び変位前後にそれ
ぞれ撮影された同一の基準計測点を含む複数の写真画像
などが挙げられる。

【００２０】写真画像における被写体上の基準計測点の
座標は、各写真画像において個別的に、例えば計測点を
写真画像を表すディスプレイ上で確認して、計測するこ
ともできるが、被写体上の同一の基準計測点を写した重
複写真の画像上の特徴から、各画像上で対応する基準計
測点を自動的に探してその座標を計測することが望まし
い。

【００２１】なお、紙やプラスチックフィルム等に形成
した写真画像のハードコピーから座標測定装置（コンパ
レータ）などにより直接基準計測点の座標を計測するこ
ともできる。 （２）写真画像の近似解析［写真撮影の標定要素及び基
準計測点の三次元座標についての精度の良い近似値の取
得］ 写真画像の精密解析における基礎方程式は非線形である
ため、写真画像の標定要素や基準計測点の三次元座標
は、それらの近似値を利用した繰り返し計算によって算
出する。従って、解の安定性を上げるためには、なるべ
く精度の良い初期近似値を与える必要がある。

【００２２】写真画像の精密解析に用いる標定要素の精
度の良い近似値及び基準計測点の三次元座標の精度の良
い近似値は、ステレオ写真画像（例えば、少なくとも一
部が重複した１対又はそれ以上の写真画像）において計
測された被写体上の基準計測点の座標（二次元座標）、
標定要素の近似値及び存在する場合には既知標定要素
（例えば内部標定要素）、並びに、所要の基準計測点に
ついての三次元座標の近似値に基づき求めることができ
る。すなわち、ステレオ写真画像の対応する光線の交会
条件によりステレオモデルを構成し、モデル結合条件に
より隣接のステレオモデル同士を接続すると共に、所要
の基準計測点について三次元座標の近似値に基づき対地
標定を行い、標定要素及び基準計測点についての三次元
座標の精度の良い近似値を求める。

【００２３】ステレオ写真画像としては、重複した部分
をもつ２つの写真画像を用いることが多いが、重複した
部分が３以上の写真画像に存在するようにした場合でも
同様にして解析を行うことができる。重複した部分が３
以上の写真画像に存在するようにした場合、１つ目と２
つ目の写真画像から作られるモデル、２つ目と３つ目の
写真画像から作られるモデルというように、同じ写真画
像を用いて異なる２つのモデルを形成し、モデルを次々
に接続していくことができる。このことを利用すると、
３つの写真画像が重なった部分でモデルを形成し接続し
ていくことで、被写体情報を使わずに統一モデルを形成
することが可能である。統一モデルが形成できれば、そ
の被写体空間内の５個の基準点によって、被写体空間の
３次元測量が可能である。

【００２４】このような精度の良い近似値は、ＤＬＴ法
（Ｄｉｒｅｃｔ Ｌｉｎｅａｒ Ｔｒａｎｓｆｏｒｍａ
ｔｉｏｎ Ｍｅｔｈｏｄ ）等の解析法を利用して算出
することができる。なお、妥当な範囲で、標定要素の一
部（例えば写真画像の線形歪等）を無視して計算するこ
ともできる。

【００２５】より具体的には、例えば次のａ乃至ｄの４
段階により、すなわちステレオモデルを基として相互標
定と対地標定の二段階に分けて解析を行うことにより、
標定要素及び基準計測点についての三次元座標の精度の
良い近似値を求めることができる。なお、この場合、写
真撮影は測定用カメラ（内部標定要素は既知であり、他
の歪は無視できる）で行ない、外部標定要素（回転角、
移動量などの求めるべき未知数）を標定して被写体座標
を求めるものとする。 ａ）ステレオ写真画像の対応する光線の交会条件による
ステレオモデルの構成 接続標定法により、ステレオ写真のうちの右写真（右写
真に限らず何れか一方の写真）の５個の外部標定要素に
ついて解いてステレオモデルを構成する。実際には、含
まれる誤差の調整を行うために、最小二乗法の未知量を
持った条件付観測による調整法を使って解く。 ｂ）モデル結合条件による隣接モデルの接続 １個の右写真の要素について最小二乗法の間接観測によ
る調整法を利用して決定する。 ｃ）地上基準点（被写体上の基準計測点の三次元座標の
近似値）を使用した、三次元相似変換による対地標定 ７個の対地標定（絶対標定）要素を求める。 ｄ）標定要素と被写体座標の近似値の計算 各写真の外部標定要素の中で、３個の移動量要素は三次
元相似変換式を用いて、３個の回転量要素は相互標定の
回転量要素と対地標定の回転量要素より、また、被写体
座標の近似値は三次元相似変換式を用いて計算する。 （２−１） 写真画像における基準計測点の座標（二次
元座標）の異常値探索とその除去 写真解析では入力データとして写真画像座標の量はかな
り多いので、写真画像における座標には、測定ミスや解
析入力ミス等により異常値が含まれる可能性がある。こ
のような異常値は、何れかの段階、好ましくは（２）
（近似解析）の段階で、自動的に探索して除去すること
が望ましい。

【００２６】写真画像における基準計測点の座標の異常
値についての（２） （近似解析）の段階での自動探索
及び除去は、例えば次のように行うことができる。 ステレオモデルの構成時に、写真画像における基準
計測点の座標の観測値うちｙ座標（二次元座標の何れか
一方でよく、必ずしもｙ座標に限らない）に含まれる異
常値の探索を行う。例えば、ロバスト推定法の中のウエ
イト関数法を用いて、対象となる座標の測定ミスや解析
入力ミス等による異常値を以下のように自動探索するこ
とができる。

【００２７】すなわち、ウエイト関数ｐ_i として例え
ば次式

【００２８】

【数１】 を使用して、ア乃至カのように異常値を探索して除去す
ることができる。 ア）写真画像座標測定量のウエイトをすべて１として最
小二乗法を実行する。 イ）各写真画像座標の残差を計算する。 ウ）写真画像座標測定量に、上記式ｐ_i より与えられ
るウエイトを与え、最小二乗法を実行する。 エ）各写真画像座標の残差を再計算する。 オ）しきい値を設定し、しきい値より小さいウエイトに
対応する観測値を異常値と認定する。 カ）認定された異常値を取り除いて最小二乗法を実行す
る。 ｙ座標（二次元座標の一方）の処理後、次の作業を行
う。写真画像における基準計測点の座標の観測値うちｘ
座標（二次元座標の他方）に含まれる異常値の探索はス
テレオモデルの構成時にはできないが、隣接モデルの接
続時に行うことができる。例えば、隣接モデルにおける
対応する同一の基準計測点に対して、統一モデルにおけ
る座標値と各モデルにおける座標値が一致しない点のｘ
座標の観測値を異常値と認定し、除去する。 （３）写真画像の精密解析［標定要素及び基準計測点の
三次元座標についての近似値に基づく写真画像の精密解
析による標定要素及び基準計測点の三次元座標の算出、
並びに必要な任意の変位計測点の三次元座標の算出］ 標定要素の精度の良い近似値及び基準計測点の三次元座
標の精度の良い近似値を用いることにより、写真画像の
精密解析を行って標定要素及び基準計測点の三次元座標
をそれぞれ厳密に算出することができると共に、必要な
任意の変位計測点の三次元座標を厳密に算出することが
できる。

【００２９】写真画像の精密解析は、バンドル法を利用
して最小二乗法により行うのが好適である。その際、必
要な繰り返し計算を行う。また、バンドル法を利用して
最小二乗法により解析するにあたっては、例えば次のよ
うに、（ａ） 計算時間の短縮、並びに（ｂ） 計測作
業の容易化及び解析精度の向上を図ることが望ましい。 （ａ） 計算時間の短縮 撮影により得た写真画像の解析に要する時間はできるだ
け短いことが望ましい。斜面やトンネル等の岩盤構造物
及びその他の構造物の安全性の判断に利用する場合は特
にその要請が強い。

【００３０】しかしながら、一般には、バンドル法を用
いる場合の未知量の数はかなり多い（例えば、写真画像
数が１０、基準計測点が１００で、未知量の数は３６
８）ので、最小二乗法で用いられる係数行列は大きくな
り、直接逆行例を計算すると非常に時間がかかる。そこ
で、被写体における基準計測点の座標に関する係数行列
の性質を利用し、縮約正規方程式を構成して計算時間の
短縮を図る。具体的には、被写体における基準計測点の
座標に関する係数行列が３＊３小行列の形のバンド性を
持つので、被写体における基準計測点の座標に関する係
数行列の逆行列を求め、この逆行列を利用して正規方程
式を縮約し、この縮約正規方程式を用いて標定要素及び
基準計測点の三次元座標を算出する。これによって、計
算時間を大幅に短縮することができる。

【００３１】なお、他の各種数値計算法を利用して計算
時間の短縮を図ることも可能である。 （ｂ） 計測作業の容易化及び解析精度の向上 バンドル法の適用においてセルフキャリブレーション
（ｓｅｌｆ ｃａｌｉｂｒａｔｉｏｎ）法を使用すれ
ば、ＤＬＴ法において必要とされる被写体上の所要点の
高精度な三次元座標の既知情報を必要とせずに内部標定
要素を算出することができるので、様々な現場における
内部標定要素が未知の又は十分に検定されていない撮影
機を用いた計測作業が容易である。

【００３２】セルフキャリブレーション法の使用におい
ては、最小二乗法の正規方程式の拘束条件として既知距
離を導入することができる。これにより、解析精度を効
果的に向上させることができる。

【００３３】既知距離の導入は、間の距離が既知である
２以上の計測点を被写体に含めて撮影した写真画像を用
い、その距離を拘束条件とすることにより行い得る。

【００３４】間の距離が既知である２以上の計測点を被
写体に含む写真画像は、例えば、２以上の計測点を有
し、それらのうち少なくとも２つの計測点の間の距離が
既知である距離拘束体（例えば棒体）を被写体に含めて
撮影することにより得ることができる。

【００３５】距離拘束体が棒体である場合、その棒体の
長さが既知であるものを使用することができる。この場
合、棒体の両端が計測点となる。

【００３６】既知距離を導入すると、被写体に既知距離
の端点（例えば、本来の計測対象ではない棒体等の距離
拘束体における既知距離の端点）を含めることになる。
この場合、既知距離の端点に対応する正規方程式の係数
行列の形が複雑になるので、被写体における基準計測点
及び距離拘束体の端点等の計測点の座標に対する行列は
対角行列ではなくなる（対角行列を含むあるバンド幅を
持ったゼロの多い行列形式からスパース行列に変わ
る）。そのため、既知距離の端点の座標に関する拘束条
件に対応する行列を並べ替え、既知距離の端点の座標と
基準計測点の座標を分けて計算を行うことにより標定要
素及び基準計測点の三次元座標について解を得ることが
好ましい。

【００３７】なお、上記（１） 、（２） 及び（３）
の何れか又は２以上或いは上記（１） 、（２） 及
び（３） の何れかの一部又は２以上のうちの一部は、
それを実行し得るコンピュータ利用装置により実行する
ことが望ましい。このコンピュータ利用装置は、写真画
像を得るための撮影機を備えたものとすることもでき
る。また、このようなコンピュータ利用装置（撮影機を
備えたものも含む）により上記（１） 、（２） 及び
（３） の何れか又は２以上或いは上記（１）、（２）
及び（３） の何れかの一部又は２以上のうちの一部
を実行するためのプログラムは、任意の記録媒体に記録
しておくことができることは言うまでもない。

【００３８】

【実施例】実施例１ 上記写真画像解析の中心部分である（３） 写真画像の
精密解析に関し、解析精度を検証するシミュレーション
テストを実施した。

【００３９】シミュレーション条件は次のように設定し
た。 画像枚数：４ 基準計測点数：４９ 撮影距離：２００００ｍｍ 画像縮尺：１／４００ 制御点数：５ 長さ既知の棒の数：７ 画像距離：５０ｍｍ 画像誤差：０．０１ｍｍ 先ず、被写体空間に設定された各座標値に対して共線条
件式を用いて写真画像上に座標を変換した。

【００４０】次に、この写真画像上の座標値に対してレ
ンズ歪とランダムな測定誤差を加えることにより、基準
計測点についての実測値に近い写真画像座標を生成し
た。

【００４１】そうして、被写体空間中の設定座標値にラ
ンダムな誤差を与えて被写体の基準計測点についての初
期近似値とし、適切なウエイトを与えた上で、セルフキ
ャリブレーション法を使用したバンドル法を利用して最
小二乗法により被写体の基準計測点の三次元座標を厳密
に計算した。

【００４２】このシミュレーションにおいては、画像測
定誤差として平均１０μｍを与え、被写体の基準計測点
についての初期近似値として真値から平均５ｃｍのばら
つきを与えた。これらの値の大きさは標識を用いないで
写真画像における計測点の座標を計測する場合に妥当と
認められるものである。

【００４３】シミュレーション解析結果は次の通りであ
る。 写真画像座標の標準誤差：０．００１１６ｍｍ 内的誤差：５．４ｍｍ 外的誤差：３．９ｍｍ 解析精度の指標としては、計算により得られた被写体座
標（三次元座標）の平均内的誤差（解の分散。共分散行
列より計算し得、解の安定性を示す指標となる）、及び
平均外的誤差（真値とのずれの二乗平均）を利用した。
なお、このケースでの理論誤差は４ｍｍである。

【００４４】被写体座標の平均内的誤差に注目すると、
５．４ｍｍであり、理想精度（４ｍｍ）より少し大きか
った。従って、理想精度を得るためには、撮影写真画像
数を増加させる等の工夫が必要であることが分かる。

【００４５】得られた結果と真値とのずれを表す平均外
的誤差は３．９ｍｍで、ほぼ理論値と一致していた。従
って、画像測定誤差平均を現在の写真計測で可能なサブ
ミクロン程度まで低下させれば、内的誤差を上記の１０
分の１、すなわち１ｍｍ以下に減少させることができ
る。また、上記写真画像解析の（２） 写真画像の近似
解析により得られる標定要素及び基準計測点の三次元座
標についての近似値の精度を向上させれば、画像測定誤
差平均がこのシミュレーションテストのケースと同様で
あっても１ｍｍ以下の精度を達成することが可能であ
る。実施例２ 岩盤露頭上に３次元的に存在している不連続面を被写体
（計測対象）とした場合の実施例を示す。但し、実際の
岩盤露頭を用いるのではなく、建物の壁と板を用いた仮
想露頭面を作成して本発明を実施した。

【００４６】本実施例における板は、実際の露頭面上に
見られる、図１に実線で示すような亀裂と露出した面の
２種類の不連続形状を表現するためのもので、図１に破
線で示すようなものを用いた。板１２は、縁部によって
仮想亀裂１０を表現しており、その仮想亀裂１０の５箇
所に設定した各基準計測点に、半径５ｍｍの標識ラベル
Ｍを貼り付けた。また、これと別個の板１６は、平面部
分によって仮想不連続面１４を表現しており、その仮想
不連続面１４の５箇所に設定した各基準計測点に、前記
と同じ標識ラベルＭを貼り付けた。このような、亀裂を
模した仮想亀裂を有する板５枚と、露出した不連続面を
模した仮想不連続面を有する板２枚を、建物の壁面に配
置し、それを仮想露頭面１８（図２）とした。また、長
さが既知の棒として、両端に反射シール（レトロ）を貼
り付けた測量用スタッフ１２本（図示せず）を仮想露頭
面１８にできるだけ万遍なく配置した。

【００４７】次いで、デジタルスチルカメラを用い、５
個以上の基準計測点及び１本以上の測量用スタッフを含
む３つの写真画像を、前記仮想露頭面１８から約１１ｍ
の距離において、２ｍ４０ｃｍの間隔おきの３個所で、
カメラ高１６０ｍｍでの平行写真撮影により得た（図２
参照）。

【００４８】撮影時に何ら被写体上の座標計測作業を行
っていないので、解析計算開始時の初期値を精度の良い
ものとはできない。解の安定性は未知量の近似値に設定
されたウエイトの値に左右される。ウエイトは次式に従
って与えた。 Ｗｅｉｇｈｔ＝ （１／［近似値に予想される外的誤差
／スケール／画像精度］）² 本実施例では以下のような近似値を与えた。

【００４９】写真画像座標を計測した写真対のうち、
中央写真の写真画像座像を近似値計算の基礎データとし
て利用した。なお、一般に、２枚対の場合は左写真、３
枚対の場合は中央写真を用いれば良いが、これに限らな
い。

【００５０】適当な位置に座標原点を想定し、ｘ，ｙ
座標に関しては画像上の２次元座標を縮尺倍し、ｚ座標
（奥行き）に関してはカメラから露頭面（計測対象）ま
でのおおよその距離を与えた。すなわち、外部標定要素
について比較的に粗い近似値を与えた。

【００５１】これによって、第１段階として、外部標定
要素に関しては比較的に低いウエイトを与え、内部標定
要素やレンズ歪係数等の歪に関しては比較的に高いウエ
イトを与え、あたかも内部標定要素が既知のカメラで撮
影を行ったようにして、不安定な解に対して繰り返し計
算により強制的に計算を収束させ、標定要素及び基準計
測点についての解を得た。

【００５２】次に、第２段階として、第１段階で比較的
に高いウエイトであった内部標定要素とレンズ歪等の歪
に関して比較的に低いウエイトを与え、全ての変数（外
部標定要素及び基準計測点）に対して適正なウエイトを
与え、再度収束計算を行なって標定要素及び基準計測点
についての精度の良い近似値を得た。

【００５３】つまり、粗い近似値の次に精度の良い近似
値を取得する２段階の近似値取得を行なった。そして、
この標定要素及び基準計測点の三次元座標の精度の良い
近似値を用いて、セルフキャリブレーション法を使用し
たバンドル法を利用して最小二乗法により被写体の基準
計測点の三次元座標を厳密に計算した。

【００５４】写真画像座標の標準誤差の推定量（ウエイ
ト＝１を与える観測値の標準誤差）と被写体座標の真値
と解析値のずれの平均誤差である外的誤差は次の通りで
あった。すなわち、精度の良い計測を行うことができ
た。

【００５５】

【数２】 平均外的誤差＝３．９５ｍｍ（平均外的誤差は、各計測
点をトランシットによって測量した値との差を用い
た。）

【００５６】

【発明の効果】本発明においては、計測対象をカメラや
ビデオ等の撮影装置により撮影することによりその計測
対象における計測点の変位を求めることが可能である。
そのため、計測点を任意に選ぶことができ、従来の機器
では設置し難かった場所や人が近づくことのできない対
象物でも、その変位を求めることができる。また、変位
計測のために現場で必要な機器はカメラやビデオ等の撮
影装置であるから、携帯性に優れたものとすることがで
き、断線や機器の不良による計測作業の中断、或いは作
業環境が計測値に与える悪影響等の不都合が生じること
が防がれる。

【００５７】また、計測のための装置のコストは従来に
比し安価であり、多数点を撮影して計測するのが容易で
あるため、計測点を増してもコスト増大が小さい。

【００５８】然も、計測対象における変位分布（面的広
がりを有する変位分布）を容易且つ迅速に精度良く求め
ることができる。

【００５９】更に、変位算出のための写真画像は、特に
それがデジタル情報である場合、そのデータによりデー
タベースを作成して現場の状態を容易に再現し得、それ
により解折結果の信頼性を高め得るものとすることがで
きる。

【図面の簡単な説明】

【図１】仮想亀裂線及び仮想不連続面の例を示す斜視図
である。

【図２】仮想露頭面の撮影状況を示す概略図である。

【符号の説明】

１０ 仮想亀裂 １２ 板 １４ 仮想不連続面 １６ 板 １８ 仮想露頭面 Ｍ 標識ラベル

───────────────────────────────────────────────────── フロントページの続き (73)特許権者 502374119 岡本 隆明 滋賀県大津市比叡平３−３−16 (74)上記１名の法定代理人 599067352 岡本 真雅紀 （外１名） (72)発明者 大西 有三 西宮市剣谷町12−15 (72)発明者 岡本 厚 大津市比叡平３−３−16 (56)参考文献 特開 平８−219783（ＪＰ，Ａ) （社）日本写真測量学会・解析写真測 量委員会，解析写真測量 改訂版，日 本，（社）日本写真測量学会，1989年 ６月20日，126−148 中川徹 小柳義夫，最小二乗法による 実験データ解析，日本，東京大学出版 会，1982年 ５月30日，157−176 Ｗｉｌｌｉａｍ Ｈ Ｐｒｅｓｓ他, ニューメリカルレシピ・イン・シー，日 本，技術評論社，1997年 ５月25日 カール・クラウス，写真測量，日本, 株式会社技術書院，1993年 ７月31日, 204−262 日本測量学会，写真による三次元測定 −応用写真測量−，日本，共立出版, 1983年 ６月 ５日 (58)調査した分野(Int.Cl.⁷，ＤＢ名) G01C 11/00 G01B 11/00 G01C 7/06 G01D 21/00

Claims (2)

(57)【特許請求の範囲】
1. 【請求項１】計測対象の変位の前後それぞれにおいて、
   内部標定要素が未知の又は十分に検定されていないデジ
   タルカメラを用いて複数の位置から各基準計測点を含む
   計測対象を撮影して複数の写真画像を得、 変位の前後における写真画像上の基準計測点の座標に基
   づき、セルフキャリブレーション法を使用したバンドル
   法を利用して最小二乗法を繰り返し実行することにより
   標定要素及び基準計測点の三次元座標を算出し、それら
   の標定要素及び基準計測点の三次元座標に基づき任意の
   変位計測点の三次元座標をそれぞれ算出してその変位を
   求める変位計測方法であって、 前記最小二乗法実行の繰り返しにおいて、写真画像にお
   ける基準計測点座標の観測値に対し、次式 【数３】 で表される、異常値を探索するためのウエイト関数によ
   るウエイトを与えて最小二乗法を実行し、その計算結果
   と写真画像における基準計測点座標の観測値との残差を
   計算し、しきい値を設定してそのしきい値より小さいウ
   エイトに対応する観測値を異常値と認定し、認定された
   異常値を除去して再度最小二乗法を実行することにより
   標定要素及び基準計測点の三次元座標の精度を高める工
   程を有することを特徴とする変位計測方法。
2. 【請求項２】上記セルフキャリブレーション法の使用に
   おいて最小二乗法の正規方程式の拘束条件として既知距
   離を導入することにより精度を向上させる請求項１記載
   の変位計測方法。