(57)【要約】【課題】撮影上の障害物が存在する条件下でも適用できる筒状リングの歪み形状計測方法及びプログラムを提供する。【解決手段】所定内周長さＧの筒状リングＡの内周面に、リング中心軸Ｃの直交断面Ｆと交差する交線Ｌ上に沿って所定角度位置θkの計測視標Tkの群を取り付け、リング内側の可動撮像機10により異なる位置・姿勢から計測視標Tkが写り込む画像Igの群を撮影する。各画像Ig内の計測視標Tkの像の二次元座標と各計測視標Tkの所定角度位置座標とから各画像Igの撮影位置・姿勢を検出し、各画像Ig内の計測視標Tkの像の二次元座標と各計測視標Tkの所定角度位置座標と各画像Igの撮影位置・姿勢とに基づくバンドル調整により各計測視標Tkの三次元座標（X,Y,Z）を算出する。更に、各三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱと所定内周長さＧとの縮尺Ｓに基づき補正した補正三次元座標（Xs,Ys,Zs）により、筒状リングＡの歪み形状を計測する。【選択図】図１
2007147498.tif 000002
【特許請求の範囲】
【請求項１】
所定内周形状及び内周長さの筒状リングの歪み形状を計測する方法において、リング中心軸の直交断面と交差するリング内周面の交線上に所定角度位置の計測視標群を取り付け、リング内側の可動撮像機により異なる位置・姿勢から前記視標が写り込む画像群を撮影し、各画像内の視標像の二次元座標と各視標の所定角度位置座標とから各画像の撮影位置・姿勢を検出し、各画像内の視標像の二次元座標と各視標の所定角度位置座標と各画像の撮影位置・姿勢とに基づくバンドル調整により各視標の三次元座標を算出し、各視標の三次元座標を結ぶ多角形の周囲長さと所定内周長さとの縮尺に基づき補正した各視標の補正三次元座標により筒状リングの歪み形状を計測してなる筒状リングの歪み形状計測方法。
【請求項２】
請求項１の計測方法において、前記リング内周面上に所定角度位置の計測視標と共に任意位置の補助視標を取り付け、前記各画像内の計測視標像の二次元座標と各計測視標の所定角度位置座標とから各画像の撮影位置・姿勢を検出すると共にその撮影位置・姿勢と各画像内の補助視標像の二次元座標とから各補助視標の位置座標を検出し、各画像内の視標像の二次元座標と各計測視標の所定角度位置座標と各補助座標の検出位置座標と各画像の撮影位置・姿勢とに基づくバンドル調整により各視標の三次元座標を算出してなる筒状リングの歪み形状計測方法。
【請求項３】
請求項１又は２の計測方法において、所定内周長さの大径筒状リングの内側に小径筒状リングを同心状に組み立てる場合に、大径リング内周面の前記交線上及び小径リング内周面の前記交線上にそれぞれ計測視標を取り付けて前記画像群を撮影し、大径リングの交線上の各視標の三次元座標を結ぶ多角形の周囲長さと所定内周長さとの縮尺に基づき補正した両リング内周面の各視標の補正三次元座標により両リング間のクリアランスを計測してなる筒状リングの歪み形状計測方法。
【請求項４】
請求項３の計測方法において、大径リング内周面のリング中心軸方向に所定距離隔てた一対の前記交線上及び小径リング内周面のリング中心軸方向に所定距離隔てた一対の前記交線上にそれぞれ計測視標を取り付けて前記画像群を撮影し、各リング内周面の交線対上の各視標の補正三次元座標と交線対間の所定距離とから両リング中心軸の交差角度を計測し、その交差角度に基づき両リング間の中心軸方向のクリアランス変位を計測してなる筒状リングの歪み形状計測方法。
【請求項５】
請求項１から４の何れかの計測方法において、前記各計測視標を所定形状とし且つ所定向きに取り付け、前記撮影した各画像内の計測視標像の一部に切欠けがある場合に、各撮影画像内の無切欠け視標像の二次元座標と各視標の所定角度位置座標とから各画像の撮影位置・姿勢を検出し、その撮影位置・姿勢と各視標の所定角度位置と所定形状及び取り付け向きとから切欠け視標の仮想画像を作成し且つその仮想画像と撮影画像との重ね合わせにより切欠け視標像の二次元座標を抽出し、撮影画像内の無切欠け視標像の二次元座標と切欠け視標像の抽出二次元座標と各視標の所定角度位置座標と各画像の撮影位置・姿勢とに基づくバンドル調整により各視標の三次元座標を算出してなる筒状リングの歪み形状計測方法。
【請求項６】
所定内周形状及び内周長さの筒状リングの歪み形状を計測するためコンピュータを、リング中心軸の直交断面と交差するリング内周面の交線上に取り付けた所定角度位置の計測視標群をリング内側の異なる位置・姿勢から撮影した画像群を入力する入力手段、筒状リングの所定内周形状及び内周長さと各計測視標の所定角度位置座標とを記憶する記憶手段、各画像内の視標像の二次元座標を検出し且つ検出二次元座標と各視標の所定角度位置座標とから各画像の撮影位置・姿勢を検出する検出手段、各画像内の視標像の二次元座標と各視標の所定角度位置座標と各画像の撮影位置・姿勢とに基づくバンドル調整により各視標の三次元座標を算出する座標算出手段、各視標の三次元座標を結ぶ多角形の周囲長さと所定内周長さとの縮尺に基づき各視標の三次元座標を補正する座標補正手段、並びに各視標の補正三次元座標により筒状リングの歪み形状を出力する出力手段として機能させる筒状リングの歪み形状計測プログラム。
【請求項７】
請求項６の計測プログラムにおいて、前記入力手段により前記リング内周面上に所定角度位置の計測視標と共に取り付けた任意位置の補助視標を撮影した前記画像群を入力し、前記検出手段により各画像内の計測視標像の二次元座標と各計測視標の所定角度位置座標とから各画像の撮影位置・姿勢を検出すると共にその撮影位置・姿勢と各画像内の補助視標像の二次元座標とから各補助視標の位置座標を検出し、前記座標算出手段により各画像内の視標像の二次元座標と各計測視標の所定角度位置座標と各補助座標の検出位置座標と各画像の撮影位置・姿勢とに基づくバンドル調整により各視標の三次元座標を算出してなる筒状リングの歪み形状計測プログラム。
【請求項８】
請求項６又は７の計測プログラムにおいて、所定内周長さの大径筒状リングの内側に小径筒状リングを同心状に組み立てる場合に、前記入力手段により大径リング内周面の前記交線上及び小径リング内周面の前記交線上にそれぞれ取り付けた計測視標群を撮影した前記画像群を入力し、前記補正手段により大径リングの交線上の各視標の三次元座標を結ぶ多角形の周囲長さと所定内周長さとの縮尺に基づき両リング内周面の各視標の三次元座標を補正し、両リング内周面の各視標の補正三次元座標により大径及び小径リング間のクリアランスを計測するクリアランス計測手段を設け、前記出力手段により両リング間のクリアランスを出力してなる筒状リングの歪み形状計測プログラム。
【請求項９】
請求項８の計測プログラムにおいて、前記入力手段により大径リング内周面のリング中心軸方向に所定距離隔てた一対の前記交線上及び小径リング内周面のリング中心軸方向に所定距離隔てた一対の前記交線上にそれぞれ取り付けた計測視標群を撮影した前記画像群を入力し、前記記憶手段に両リング内周面の交線対間の所定距離を記憶し、各リング内周面の交線対上の各視標の補正三次元座標と交線対間の所定距離とから両リング中心軸の交差角度を計測する角度計測手段を設け、前記クリアランス計測手段により両リング中心軸の交差角度に基づき両リング間の中心軸方向のクリアランス変位を計測してなる筒状リングの歪み形状計測プログラム。
【請求項１０】
請求項６から９の何れかの計測プログラムにおいて、前記記憶手段に各計測視標の所定形状及び取り付け向きを記憶し、前記入力手段の入力画像内の計測視標像の一部に切欠けがある場合に、前記検出手段により各入力画像内の無切欠け視標像の二次元座標と各視標の所定角度位置座標とから各画像の撮影位置・姿勢を検出し、その撮影位置・姿勢と各視標の所定角度位置と所定形状及び取り付け向きとから切欠け視標の仮想画像を作成し且つその仮想画像と撮影画像との重ね合わせにより切欠け視標像の二次元座標を抽出する座標抽出手段を設け、前記座標算出手段により入力画像内の無切欠け視標像の二次元座標と切欠け視標像の抽出二次元座標と各視標の所定角度位置座標と各画像の撮影位置・姿勢とに基づくバンドル調整により各視標の三次元座標を算出してなる筒状リングの歪み形状計測プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は筒状リングの歪み形状計測方法及びプログラムに関し、とくに筒状リングの画像からその精確な歪み形状を計測する方法及びプログラムに関する。本発明は、シールドトンネルを掘削するシールド掘削機のスキンプレート、そのスキンプレートの内側に組み立てるシールドセグメント等の筒状リングの歪み形状計測に適用することができる。
【背景技術】
【０００２】
シールド掘進機を用いるトンネル掘削工事では、図１２に示すように、地盤を一定距離だけ掘進する毎に、シールド掘進機１の筒状外周壁であるスキンプレート２のテール部（後端部）の内側に、トンネル内壁となるセグメントリング７を組み立てる。通常のセグメントリング７は複数の円弧状又は楕円弧状のセグメントピース８により構成され、テール部に設けたエレクタ６により複数のセグメントピース８を例えば真円又は楕円の筒状リングに組み立てた後、ボルト等で本締めして構築する。エレクタ６の真円保持装置等によりセグメントリング７をできるかぎり精確に組み立てることが大切であるが、組み立てたセグメントリング７を本締めする前に、その詳細な形状を計測してチェックすることが品質管理上重要とされている。
【０００３】
従来から、組み立て直後のセグメントリング７の外周面とスキンプレート２の内周面と間の距離（テールクリアランス）をコンベックスやスケール等で測定し、その測定値とスキンプレート２の内周面の断面形状（以下、内周形状という）とからセグメントリング７の形状を計測する方法（テールクリアランス計測法）が行われている。例えばスキンプレート２の内周形状を真円とすれば、その内周面の半径とテールクリアランス測定値との差としてセグメントリング７の外周面の半径を計測できる。コンベックス等による手動測定に代えて、特許文献１及び２が開示するように、自動化等を目的としたテールクリアランス測定装置等も提案されている。
【０００４】
また図１２に示すように、エレクタ６の先端部にレーザ距離計36等を取り付け、エレクタ６を回転させながらセグメントリング７の内径を測定し、その測定値からセグメントリング７の内周形状を計測する方法（内径計測法）も行われている。エレクタ６を用いる方法に代えて、セグメントリング７の内周面上に予め測定用罫線を設けておき、組み立て直後の相対する罫線間の距離（直径）をレーザ距離計等で測定してセグメントリング７のリング内径を測定する場合もある。
【０００５】
しかし、従来のテールクリアランス計測法及び内径計測法は何れも、十分な精度でセグメントリング７の形状を求めることができず、しかも計測作業に手間がかかる問題点がある。シールド掘進機１のスキンプレート２の内周形状は掘削作業によって変形し、当初は真円であってもセグメントリング７の組み立て時の正確な形状が不明であるため、テールクリアランス計測法ではスキンプレート２の変形に伴う誤差が生じる。また内径計測法では、エレクタ６の機械的な軸ずれやたわみ等による誤差が生じ、手作業で測定する場合は人為的な読み取り誤差等も生じる。更に、セグメントリング７の詳細な形状を計測するには測定箇所を増やす必要があるが、従来の計測法は測定箇所の増加に応じて作業量も増大するので計測作業に時間がかかる。とくに内径計測法は、リング全体にわたり見通しを確保する必要があり、他の作業と並行に計測作業を進めることも難しいため、組み立て直後のセグメントリング７の形状を迅速に把握することは困難であった。
【０００６】
これに対し特許文献３は、２台のカメラ（撮像機）を用いた画像計測法によりセグメントリング７の形状を計測する方法を提案している。画像計測法は、図１０に示すように、対象点Ｐとカメラ撮像面の中心（カメラ中心）Ｏとその撮像面上の像点ｐとの３点が同一直線上に存在するという幾何学的原理に基づくものである。地上座標系における対象点Ｐの三次元座標を（X,Y,Z）、カメラ中心Ｏの三次元座標を（X0,Y0,Z0）、撮像機の焦点距離をc、カメラ座標系における撮像機のx軸、y軸及びz軸の回りの回転角度を（ω,φ,κ）、カメラ座標系における像点ｐの三次元座標を（x,y,-c）とした場合、画像計測の幾何学的原理は式(1)〜(3)の共線条件式として表すことができる。また、カメラの焦点距離、主点位置、レンズ歪等を考慮する場合は、式(1)、(2)の共線条件式を変形して式(4)及び(5)とすることができる。式(4)及び(5)におけるΔx、Δyは、焦点距離、主点位置、レンズ歪係数により定まる補正項である。
【０００７】
【数１】
2007147498.tif 000003
【０００８】
共線条件式(1)及び(2)はカメラ位置（X0,Y0,Z0）及びカメラ角度（ω,φ,κ）の６つの未知数（外部標定要素）を含むが、画像計測では既知三次元座標（X,Y,Z）の複数の対象点Ｐ（以下、基準点ということがある）を写し込み、撮像面上の基準像点ｐの二次元座標（x,y）を検出することにより共線条件式の未知数を標定することができる。特許文献３のように２台のカメラを用いた場合は、１つの基準点Ｐに対し２つの像点ｐ1、ｐ2の座標が得られるので、原理的には１枚の画像中に少なくとも３つの基準点Ｐを写し込めば共線条件式(1)及び(2)の６つの未知数を標定できる。共線条件式(4)及び(5)を用いる場合は、カメラ位置（X0,Y0,Z0）及びカメラ角度（ω,φ,κ）に加えて焦点距離、主点位置、レンズ歪係数等が未知数（内部標定要素）であるため、評定のために更に多くの基準点を必要とする。共線条件式における未知数が標定できれば、撮像面上における任意像点ｐの二次元座標を共線条件式へ代入することにより、それと対応する地上座標系の対象点Ｐ（以下、計測点Ｐということがある）の三次元座標を算出することができる（前方交会法）。
【０００９】
特許文献３の画像計測法では、図１２に示すように、セグメントリング７の内周面上の複数の対象点に反射板又はカラーマーキング等の視標Ｔを取り付け、シールド掘削機１の異なる所定位置に所定姿勢（向き）で固定した一対の撮像機10により各視標Ｔを含むセグメントリング７の画像Igを撮影し、撮像機10の所定位置及び姿勢（以下、位置・姿勢と表すことがある）と画像Ig内の視標像の二次元位置とから各視標Ｔのシールド掘削機１に対する三次元位置を算出し、その三次元位置によりセグメントリング７の三次元形状を計測する。この画像計測法によれば、セグメントリング７の画像を撮影するだけで三次元形状を計測できるので計測作業時間を短縮できると共に、スキンプレート２の変形やエレクタ６の機械的精度に関係しない誤差の小さなセグメントリング７の三次元形状の計測が期待できる。
【００１０】
【特許文献１】特開２００４−０６８３６６号公報
【特許文献２】特開２００３−０３０８８８号公報
【特許文献３】特許第２６８５４０７号公報
【特許文献４】特許第３５３０９７８号公報
【非特許文献１】秋本圭一・服部進「画像計測の基礎」岡山職業能力開発短期大学校紀要、第11号、1997年3月
【発明の開示】
【発明が解決しようとする課題】
【００１１】
しかし、特許文献３の計測方法は、セグメントリング７の三次元形状を計測するために撮像機10の位置・姿勢を固定する必要があり、その位置・姿勢の測定に手間がかかる問題点がある。図１２に示すように実際のシールド掘進機１の切羽付近にはコンベアや配管、風管その他の様々な器具や装置が配置されており、撮像機10を一定の位置・姿勢に固定したままでは、器具等が障害物となってセグメントリング７の視標Ｔを写し込むことができない場合がある。そのため撮像機10の位置・姿勢を適宜移動させる必要があるが、特許文献３の方法では撮像機10を移動させる都度その固定位置・姿勢を測定する必要があるため、結果的に計測作業に要する時間が長くなってしまう。組み立て直後のセグメントリング７の形状を短時間で計測するためには、撮像機10の位置・姿勢の固定及び測定を必要としない画像計測法が必要である。
【００１２】
撮像機10の位置・姿勢の測定を必要としない画像計測法として、特許文献４は、異なる任意の位置・姿勢から撮影した複数の画像を用いたバンドル調整法による画像計測法を提案している。バンドル調整法とは、上述したように基準点Ｐの写り込んだ画像毎に共線条件式の未知数を標定する方法（単写真標定）ではなく、図１１に示すように、基準点Ｐの写り込んだ複数の画像から共線条件式の未知数の標定と計測点Ｐの三次元座標の算出とを最小二乗法により同時に行う方法である（非特許文献１参照）。すなわち、基準点及び計測点Ｐの三次元座標、カメラ中心Ｏの三次元座標、及びカメラ角度の真値を各々の近似値（X',Y',Z'）（X'0,Y'0,Z'0）及び（ω',φ',κ'）に補正量を加えたもの（X'＋ΔX,Y'＋ΔY,Z'＋ΔZ）、（X'0＋ΔX0,Y'0＋ΔY0,Z'0＋ΔZ0）及び（ω'＋Δω,φ'＋Δφ,κ'＋Δκ）とし、基準点及び計測点Ｐの像点ｐの二次元座標の真値をその二次元座標の検出値（x',y'）に誤差を加えたもの（x'＋Δx,y'＋Δy）とする。補正量を加えた近似値と誤差を含む検出値とを複数の共線条件式（式(1)及び(2)）に代入し、それらの共線条件式をテーラー展開により線形化した上で、各補正量（ΔX,ΔY,ΔZ）、（ΔX0,ΔY0,ΔZ0）、（Δω,Δφ,Δκ）及び誤差（Δx,Δy）を最小にする解を逐次繰り返し法（逐次近似解法）の収束解として求める。共線条件式(4)及び(5)を用いたバンドル調整法によれば、レンズ歪係数等の補正量を最小にする解も求めることができる。バンドル調整法には、基準点Ｐの正確な三次元座標を求める必要がないので計測作業を簡単化できると共に、補正量の最小化により通常要求される1/20,000〜1/50,000（長さ10ｍに対し200〜500μｍの誤差）以上の高精度で計測点Ｐの三次元座標を計測できる利点もある。
【００１３】
特許文献４の画像計測法を、図８の流れ図及び図９（Ａ）のブロック図を参照して、本発明の理解に必要な限度において説明する。先ずステップ401〜402において計測対象30上の計測部位にそれぞれ相互識別可能な識別視標Ta（図９（Ｂ）及び（Ｃ）参照）を取り付け、計測対象30上又は近傍の基準位置に所定相互間隔Ｌの基準視標群Tb（同図（Ｄ）参照）を固定する。ステップ403において、可動撮像機10により異なる任意の位置・姿勢から、基準視標群Tbと３以上の識別視標Taとが共通に写り込み且つ非共通の識別視標Taが含まれる第１画像Ig1及び第２画像Ig2を撮影する。ステップ404〜405において、第１画像Ig1及び第２画像Ig2を視標像座標検出手段31（図９（Ａ）参照）に入力して各画像Ig1、Ig2における各基準視標Tb及び３以上の識別視標Taの視標像の二次元座標を検出し、撮影位置・姿勢検出手段19（図９（Ａ）参照）により視標Ta、Tbの像の二次元座標と基準視標群Tbの所定相互間隔Ｌとから基準位置５を原点とする座標系における第１画像Ig1及び第２画像Ig2の撮影の位置・姿勢と３以上の識別視標Taの三次元座標とを算出する。
【００１４】
図９（Ｄ）に示す基準視標群Tbは、直交する２直線の交点上に配置した基準視標Tb1と、一方の直線上に配置した２つの基準視標Tb2、Tb3と、他方の直線上に配置した３つの基準視標Tb4、Tb5、Tb6とを含み、例えば交点の基準視標Tb1を基準位置に位置合わせして固定する。２直線上に配置される基準視標Tb1〜Tb6の数が相違するので、画像Ig1、Ig2に写し込まれた基準視標群Tb1〜Tb6から２直線を識別して座標軸とし、基準視標Tb1〜Tb6の相互間隔Ｌを各座標軸の単位長さとすることができる。撮影位置・姿勢検出手段19は、座標軸と相互間隔Ｌとで定まる基準視標Tb1〜Tb6の三次元座標とその視標像の二次元座標とを共線条件式(1)及び(2)へ代入することにより共線条件式の未知数、すなわち第１及び第２画像Ig1、Ig2の撮影位置・姿勢を算出する（ステップ405）。その後、第１及び第２画像Ig1、Ig2内の各識別視標Taの像の二次元座標を共線条件式(1)及び(2)へ代入することにより、基準位置を原点とする座標系における各識別視標Taの三次元座標を算出する（ステップ405）。
【００１５】
次いでステップ406〜407において可動撮像機10を移動させ、第ｎ画像Ign（ｎは３以上の自然数）を、その中に第（ｎ−１）画像Ig(n-1)内の座標算出済の３以上の識別視標Taと座標未算出の識別視標Taとが共通に写り込み且つ第（ｎ−１）画像Ig(n-1)と非共通の識別視標Taが含まれる位置・姿勢で撮影する。またステップ408〜409において、撮影した第（ｎ−１）画像Ig(n-1)及び第ｎ画像Ignを視標像座標検出手段31に入力して各画像Ig(n-1)、Ignにおける視標Taの像の二次元座標を検出し、撮影位置・姿勢検出手段19によりその視標Taの像の二次元座標と座標算出済視標Taの三次元座標とから、第ｎ画像Ignの撮影位置・姿勢と座標未算出の識別視標Taの三次元座標とを算出する。第ｎ画像Ignと第（ｎ−１）画像Ig(n-1)とに座標算出済の３以上の識別視標Taを共通に写し込むことにより、その識別座標の三次元座標を用いて第ｎ画像Ignの撮影位置・姿勢を算出し、第（ｎ−１）画像Ig(n-1)において座標未算出の識別視標Taの三次元座標を第ｎ画像Ignの利用により算出する。更にステップ410において全ての識別視標Taの三次元座標が算出されたか否かを判断し、座標未算出の識別視標Taが残っている場合はステップ406へ戻り、ステップ406〜409を繰り返して全ての識別視標Taの三次元座標を算出する。その後、ステップ411においてバンドル調整手段17により上述したバンドル調整の演算を行い、全ての識別視標Taの三次元座標の高精度化を図ることにより各計測部位の三次元座標を求める。
【００１６】
特許文献４の画像計測法によれば、複数の画像Igから撮像機10の位置・姿勢を算出することができるので、撮影位置・姿勢を別途測定する手間を省くことができ、計測作業の省力化及び時間短縮を図ることができる。ただし、識別視標Taの近傍に計測座標系となる基準視標群Tbを固定する必要があるため、基準視標群Tbを写し込むことができない場合は識別視標Taの三次元座標の計測ができなくなる問題点がある。上述したようにシールド掘進機の切羽付近には様々な器具等が配置されており、セグメントリングの形状計測を撮影上の障害物が存在する見通しの悪い条件下で行わなければならない場合も多い。基準視標群Tbを適宜移動させながら計測作業を進めることも可能であるが、計測作業の手間が増える。見通しの悪い条件下でのセグメントリングの形状計測には、基準視標群Tbを写し込む必要のない画像計測法を開発することが有効である。
【００１７】
そこで本発明の目的は、撮影上の障害物が存在する条件下でも適用できる筒状リングの歪み形状計測方法及びプログラムを提供することにある。
【課題を解決するための手段】
【００１８】
図１の実施例及び図２の流れ図を参照するに、本発明による筒状リングの歪み形状計測方法は、所定内周形状及び内周長さＧ（例えば内周長さＧ＝2πR）の筒状リングＡの歪み形状を計測する方法において、リング中心軸Ｃの直交断面Ｆと交差するリング内周面の交線Ｌ上に所定角度位置θkの計測視標Tkの群を取り付け、リング内側の可動撮像機10により異なる位置・姿勢から計測視標Tkが写り込む画像Igの群を撮影し、各画像Ig内の計測視標Tkの像の二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）とから各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）を検出し、各画像Ig内の計測視標Tkの像の二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）と各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）とに基づくバンドル調整により各計測視標Tkの三次元座標（X,Y,Z）を算出し、各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱ（図４参照）と所定内周長さＧとの縮尺Ｓに基づき補正した各計測視標Tkの補正三次元座標（Xs,Ys,Zs）により筒状リングＡの歪み形状を計測してなるものである。
【００１９】
好ましくは、図１に示すように、筒状リングＡの内周面上に所定角度位置θkの計測視標Tkと共に任意位置の補助視標Tsを取り付け、各画像Ig内の計測視標Tkの像の二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）とから各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）を検出すると共にその撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）と各画像Ig内の補助視標像Tsの二次元座標（x',y'）とから各補助視標Tsの位置座標（X',Y',Z'）を検出し、各画像Ig内の視標Tk、Tsの像の二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）と各補助座標Tsの検出位置座標（X',Y',Z'）と各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）とに基づくバンドル調整により計測視標Tk及び補助視標Tsの各々の三次元座標（X,Y,Z）を算出する。
【００２０】
また図１（Ｂ）のブロック図及び図２の流れ図を参照するに、本発明による筒状リングの歪み形状計測プログラムは、所定内周形状及び内周長さＧ（例えば内周長さＧ＝2πR）の筒状リングＡの歪み形状を計測するためコンピュータ13を、リング中心軸Ｃの直交断面Ｆと交差するリング内周面の交線Ｌ上に取り付けた所定角度位置θkの計測視標Tkの群をリング内側の異なる位置・姿勢から撮影した画像Igの群を入力する入力手段16、筒状リングの所定内周形状及び内周長さＧと各計測視標Tkの所定角度位置座標（X',Y',Z'）とを記憶する記憶手段15、各画像Ig内の計測視標Tkの像の二次元座標（x',y'）を検出し且つ検出二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）とから各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）を検出する検出手段17、各画像Ig内の計測視標Tkの像の二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）と各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）とに基づくバンドル調整により各計測視標Tkの三次元座標（X,Y,Z）を算出する座標算出手段21、各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱ（図４参照）と所定内周長さＧとの縮尺Ｓに基づき各計測視標Tkの三次元座標（X,Y,Z）を補正する座標補正手段22、並びに各計測視標Tkの補正三次元座標（Xs,Ys,Zs）により筒状リングＡの歪み形状を出力する出力手段25として機能させるものである。
【００２１】
好ましくは、図１に示すように所定内周長さＧの大径筒状リングＡの内側に小径筒状リングＢを同心状に組み立てる場合に、入力手段16により大径リングＡの内周面の中心軸CAと直交する交線LA上及び小径リングＢの内周面の中心軸CBと直交する交線LB上にそれぞれ取り付けた計測視標Tkの群を撮影した画像Igの群を入力し、座標補正手段22により大径リングＡの交線LA上の各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱと所定内周長さＧとの縮尺Ｓに基づき両リングＡ、Ｂの内周面の各計測視標Tkの三次元座標（X,Y,Z）を補正し、両リングＡ、Ｂの内周面の各計測視標Tkの補正三次元座標（Xs,Ys,Zs）により大径及び小径リングＡ、Ｂの間のクリアランスｅを計測するクリアランス計測手段24を設け、出力手段25により両リングＡ、Ｂの間のクリアランスｅを出力する。
【００２２】
更に好ましくは、図３に示すように、入力手段16により大径リングＡの内周面のリング中心軸CA方向に所定距離ｄだけ隔てた一対の交線LA1、LA2上及び小径リングＢの内周面のリング中心軸CB方向に所定距離ｄだけ隔てた一対の交線LB1、LB2上にそれぞれ取り付けた計測視標Tkの群を撮影した画像Igの群を入力し、記憶手段15に両リングＡ、Ｂの内周面の交線対間（LA1、LA2）、（LB1、LB2）の所定距離ｄを記憶し、各リングＡ、Ｂの内周面の交線対（LA1、LA2）、（LB1、LB2）上の各計測視標Tkの補正三次元座標（Xs,Ys,Zs）と交線対（LA1、LA2）、（LB1、LB2）の間の所定距離ｄとから両リングＡ、Ｂの中心軸CA、CBの交差角度δを計測する角度計測手段23を設け、クリアランス計測手段24により両リングＡ、Ｂの中心軸CA、CBの交差角度δに基づき両リングＡ、Ｂの間の中心軸CA、CB方向のクリアランスｅを計測する。
【発明の効果】
【００２３】
本発明による筒状リングの歪み形状計測方法及びプログラムは、所定内周形状及び内周長さの筒状リングＡの内周面の交線Ｌ上に所定角度位置θkの計測視標Tkの群を取り付け、異なる位置・姿勢から撮影した複数の画像Ig内の計測視標Tkの像の二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）とに基づき各画像Igの撮影位置・姿勢を検出し、その撮影位置・姿勢と各計測視標Tkの像の二次元座標（x',y'）と各計測視標Tkの所定角度位置座標（X',Y',Z'）とに基づくバンドル調整により各計測視標Tkの三次元座標（X,Y,Z）を算出し、各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱと所定内周長さとの縮尺Ｓに基づき各計測視標Tkの三次元座標（X,Y,Z）を補正し、各計測視標Tkの補正三次元座標（Xs,Ys,Zs）により筒状リングＡの歪み形状を計測するので、次の顕著な効果を奏する。
【００２４】
（イ）各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱが筒状リングＡの内周長さと一致するように補正するので、基準となる単位長さを示す視標を写し込まなくても筒状リングＡの歪み形状を定量的に計測できる。
（ロ）単位長さを示す視標を写し込む必要がないので、障害物等により写し込めない視標が存在する場合でも筒状リングＡの歪み形状を計測できる。
（ハ）筒状リングＡの内周長さは内周形状が変形しても維持されるので、内周形状の変形に拘らず筒状リングＡの歪み形状を精度よく計測できる。
（ニ）筒状リングＡに予め計測視標Tkを取り付けておけば、複数の画像Igを撮影するだけで歪み形状を継続的に計測することができ、計測作業の大幅な省力化・迅速化が実現できると共に他の作業と並行に進めることも可能である。
【００２５】
（ホ）従来の計測法では困難であったリング中央点の位置を求めることができるので、筒状リングの詳細なはらみ出しの向きや量を計測することができる。
（ヘ）大径筒状リングＡの内側に同心状に組み立てる小径筒状リングＢにも計測視標Tkを取り付ければ、同じ画像Igから両リングＡ、Ｂの精確な歪み形状を同時に計測できる。
（ト）また、大径筒状リングＡと小径筒状リングＢの歪み形状を同時に計測することにより、両リングＡ、Ｂの間の精確なクリアランスｅを計測できる。
（チ）シールド掘削機１のスキンプレートの歪み形状や、その内側に組み立てるセグメントリングの歪み形状の計測への有効利用が期待できる。
【発明を実施するための最良の形態】
【００２６】
図１（Ａ）は、この場合シールド掘削機１のスキンプレート２である筒状リングＡに本発明を適用した実施例を示す。図示例のスキンプレート２は内周形状が円形の筒状リングＡであり、その内径Ｒにより筒状リングＡの内周長さＧ＝2πRが定まる。ただし、本発明の適用対象の筒状リングＡは内周形状が円形のものに限定されず、内周形状が所定であって内周長さが一定の筒状リングＡに広く適用可能である。例えば内周形状が楕円又は長円形状のスキンプレート２やトンネル壁面等に本発明を適用することができ、筒状リングＡの断面形状の一部分に直線が含まれていてもよい。また、筒状リングの全周にわたる内周長さＧが分からない場合でも、その一部分の内周長さＧが分かれば、その内周部分に本発明を適用してその部分の歪み形状を計測することも可能である。同図（Ｂ）は本発明の実施に使用する装置の一例のブロック図を示し、図２は本発明における処理の流れ図の一例を示す。以下、図２の流れ図を参照して本発明による筒状リングＡの形状計測方法を説明する。
【００２７】
先ずステップ101において、図１（Ａ）に示すように、筒状リングＡのリング中心軸CAの直交断面FAと交差するリング内周面の交線LA上の複数の所定角度位置θkにそれぞれ計測視標Tkを取り付ける（ステップ101のａ）。計測視標Tkとして適当な反射板又はカラーマーキング等を用いることができるが、特許文献４の場合と同様に各視標Tkを自動識別可能とするため、相互識別可能な識別視標Taを用いることができる。相互識別可能な識別視標Taの一例は、図９（Ｂ）に示すドット分布型又は同図（Ｃ）に示す共心型の視標であり、ドット分布型の識別視標Taは500〜600の識別コードを容易に作ることができる。例えば同図（Ｂ）の識別視標Taにおいて、周辺部の３つの反射素材40O、40A、40Bは座標軸ＯＡＢを定める座標軸用反射素材であり、残りの３つの反射素材41a、41b、41cはＩＤコードを表すＩＤ用反射素材であり、３つのＩＤ用反射素材41a、41b、41cの配置により複数の識別コードを表す。
【００２８】
図示例では、筒状リングＡの内周面の交線LAをｎ等分した各角度位置θk（例えば48分割した場合は7.5度ピッチの角度位置）にそれぞれ計測視標Tkを取り付けている。ただし各計測視標Tkの取り付け角度位置θkのピッチは一定でなくてもよい。筒状リングＡの形状計測精度を高めるためには、その内周面の交線LA上にできるだけ多数（例えば48〜360個程度）の計測視標Tkを取り付けることが望ましく、計測視標Tkの数を増やすことで後述する縮尺Ｓに基づく補正の精度を高めることができる。例えば筒状リングＡをスキンプレート２とした場合、シールド掘進機１の製造工場等においてそのスキンプレート２の内周面の所定交線Ｌ上に沿って所定角度位置θkの計測視標Tkを予め取り付けておけば、計測時の視標取り付け作業の手間を省くことができる。一度取り付けた計測視標Tkは、筒状リングＡの同じ部位の形状を計測する際に繰り返し利用できる。
【００２９】
またステップ101において、所定角度位置θkの計測視標Tkと共に、筒状リングＡの内周面上の任意位置に補助視標Tsを取り付けることができる（ステップ101のｂ）。補助視標Tsを用いて筒状リングＡの内周面上の視標の数を増やすことにより、後述するバンドル調整の演算精度向上が期待できる。補助視標Tsとして、計測視標Tkと同様の反射板又は相互識別可能な識別視標Taを用いることができる。また、補助視標Tsは筒状リングＡの内周面上又は内側の計測を必要する所要部位に取り付けることができ、演算精度の向上と共に所要部位の三次元位置の計測に利用することもできる。この補助視標Tsの取り付け位置は、後述するように各画像Ig内の補助視標像Tsの二次元座標（x',y'）から検出する。
【００３０】
ステップ102において、筒状リングＡの所定内周形状と所定内周長さＧと各計測視標Tkの取り付け角度位置θkとを、コンピュータ13（図１（Ｂ）参照）の記憶手段15に記憶する。この内周長さＧ（＝2πR）と取り付け角度位置θkとから、例えば内周形状が円形の筒状リングＡ上の各計測視標Tkの角度位置座標（X',Y',Z'）を、リング中心軸CAと直交断面FAとの交点OAを原点とする座標系のXY平面上の座標（Rcosθk,Rsinθk,0）として算出することができる。また筒状リングＡの内周形状が円形以外であっても、その内周形状と内周長さＧと取り付け角度位置θkとが分かれば、同様に交点OAを原点とする座標系における計測視標Tkの角度位置座標（X',Y',Z'）を算出できる。ただし、座標系の原点は交点OAに限定されるものではない。この角度位置座標（X',Y',Z'）は、後述するバンドル調整において各計測視標Tkの近似三次元座標として利用する。
【００３１】
次にステップ103において、リング内側の可動撮像機10により、異なる位置・姿勢から視標Tk、Tsが写り込んだ複数の画像Igを撮影する。撮像機10の一例はCCDカメラ等のデジタルカメラであるが、従来の光学フィルム式カメラを用いることも可能であり、その場合は図１の撮像機10と入力手段16との間にアナログ画像をデジタルデータに変換するスキャナー等を設ける。撮像機10にはストロボやフラッシュ等の光源11を含めることができ（図９（Ａ）参照）、光源11を用いて視標Tk、Tsの撮影が可能な限り撮像機10の絞りを絞ることにより、視標Tk、Tsのみが浮き上がる画像Igを撮影することができる。望ましくは全ての画像Igを１台の撮像機10で撮影し、撮像機10の内部構造に起因する未知数（内部標定要素）の増加を防ぐ。また形状計測の精度を向上するためには、できるだけ異なる位置・姿勢から撮影した画像Igを用いることが望ましい。図示例のように大断面のスキンプレート２を内側から撮影する場合、リング中心軸CAより低い下方部や坑口側からの画像Igに偏りがちであるため、足場等の上方部や切羽側からの画像Igを加えて精度向上を試みる。ステップ104において、撮像機10を移動させながら撮影した視標Tk、Tsの複数の画像Igをコンピュータ13の入力手段16（図１（Ｂ）参照）に入力する。
【００３２】
ステップ105〜111は、画像Igを入力したコンピュータ13におけるプログラム処理を示す。先ずステップ105において、各画像Igを検出手段17（図１（Ｂ））に入力し、二次元座標検出手段18により各画像Ig内の計測視標Tk及び補助視標Tsの像の二次元座標（x',y'）を検出する（ステップ105のａ）。視標像の二次元座標の検出方法として従来技術に属する適当な方法を選択できるが、例えば適当な閾値による画像Igの２値化処理によって現れる各視標像の重心を二次元座標として検出する。この重心を用いた二次元座標の検出は、一部に切欠けのある視標像（切欠け視標像）を用いて重心を検出すると二次元座標の誤差が大きくなるため、画像Ig上の切欠けのない視標像（無切欠け視標像）のみを用いて行うことが望ましい。ただし、後述する座標抽出手段を用いることにより、切欠け視標像から誤差の少ない二次元座標を検出することも可能である。
【００３３】
ステップ105では、計測視標Tk及び補助視標Tsとして識別視標Ta（図９（Ｂ）及び（Ｃ）参照）を用いた場合に、二次元座標検出手段18によって各視標Tk、Tsの二次元座標（x',y'）の検出と同時にＩＤを検出することもできる。例えば同図（Ｂ）のドット分布型の識別視標Taを用いた場合は、２値化処理により現れた座標軸用反射素材40O、40A、40Bの視標像に基づき識別視標Ta上の二次元座標軸を決定し、座標軸用反射素材40O（又は40A、40B）の重心を用いて識別視標Taの二次元座標を検出し、ＩＤ用反射素材41a、41b、41cの配置に基づき識別視標TaのＩＤを読み取る。ただし、各視標Tk、視標Tsの識別は、画像Igや二次元座標の検出時の２値化画像等をディスプレイ等の出力装置14に出力し、作業者がディスプレイの画面を見ながら行ってもよい。
【００３４】
またステップ105において、二次元座標検出手段18により検出した各視標Tk、視標Tsの像の二次元座標（x',y'）を撮像位置・姿勢検出手段19に入力し、各計測視標Tkの像の検出二次元座標（x',y'）と所定角度位置座標（X',Y',Z'）＝（Rcosθk,Rsinθk,0）とから、各画像Igの撮影位置・姿勢を検出する（ステップ105のｂ）。撮像位置・姿勢検出手段19は、例えば６以上の計測視標Tkが写し込まれた画像Igを選択し、その画像Igに写し込まれた計測視標Tkの二次元座標（x',y'）と角度位置座標（X',Y',Z'）とを共線条件式(1)及び(2)（又は(3)及び(4)）へ代入することにより、その共線条件式の未知数としてその画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）を算出する。また図９のステップ405及び409と同様に、３以上の計測視標Tkが写し込まれた一対の画像Igを選択し、その画像Igの対に写し込まれた計測視標Tkの二次元座標（x',y'）と角度位置座標（X',Y',Z'）とを共線条件式(1)及び(2)（又は(3)及び(4)）へ代入することにより、それら一対の画像Igの撮影位置・姿勢を算出することもできる。この算出サイクルを各画像Igについて繰り返すことにより、全ての画像Igの撮影位置・姿勢を検出できる。
【００３５】
更にステップ105において、撮像位置・姿勢検出手段19により、各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）と各画像Ig内の補助視標像Tsの二次元座標（x',y'）とから、各画像Ig内の取り付け位置未定の補助視標Tsの三次元位置座標（X',Y',Z'）を検出することができる（ステップ105のｃ）。具体的には、各画像Igの共線条件式(1)及び(2)にその画像Ig内の取り付け位置未定の補助視標像Tsの二次元座標（x',y'）を代入することにより、計測視標Tkと同じ座標系における補助視標像Tsの三次元位置座標（X',Y',Z'）を算出する。この算出サイクルを全ての画像Igについて繰り返すことにより、全ての画像Ig内の取り付け位置未定の補助視標像Tsの位置座標（X',Y',Z'）を検出できる。
【００３６】
次いでステップ106において、座標算出手段21により、ステップ105において検出した各計測視標Tk及び補助視標Tsの位置座標（X',Y',Z'）に含まれる誤差をバンドル調整により最小化する。具体的には、各計測視標Tk及び補助座標Tsの三次元座標（X,Y,Z）をその近似値である位置座標（X',Y',Z'）に補正量を加えた値（X'＋ΔX,Y'＋ΔY,Z'＋ΔZ）とし、各画像Igの撮影位置・姿勢の近似値に補正量を加えた値（X'0＋ΔX0,Y'0＋ΔY0,Z'0＋ΔZ0）及び（ω'＋Δω,φ'＋Δφ,κ'＋Δκ）と、各視標像の二次元座標の検出値に誤差を加えた値（x'＋Δx,y'＋Δy）とを複数の共線条件式(1)及び(2)（又は(3)及び(4)）に代入し、テーラー展開により複数の線形方程式を作る。その複数の線形方程式から、各視標Tk、Tsの三次元座標の補正量（ΔX,ΔY,ΔZ）と、各画像Igの撮影位置・姿勢の補正量（ΔX0,ΔY0,ΔZ0）及び（Δω,Δφ,Δκ）と、各視標像の二次元座標の誤差（Δx,Δy）とを最小にする解を最小二乗法により求め、求めた最小補正量及び誤差を近似値及び検出値に加えたものを新たな近似値及び検出値として採用し、再度補正量及び誤差を最小にする解を求める。この過程を補正量及び誤差が十分に小さくなるまで繰り返す逐次近似解法により、各視標Tk、Tsの三次元座標（X,Y,Z）の収束解を求める。
【００３７】
更にステップ107において、座標算出手段21で算出した各計測視標Tk及び補助視標Tsの三次元座標（X,Y,Z）の座標系の単位長さを座標補正手段22により補正する。ステップ106で算出した各視標Tk、Tsの三次元座標（X,Y,Z）は、特許文献４における基準視標群Tbのような座標軸の単位長さの画像Ig中の基準がないので、座標系の単位長さが調整されていない。座標補正手段22は、筒状リングＡの内周長さが断面形状の変形にも拘らず維持されていることに基づき、各三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱが所定内周長さＧと一致するような縮尺Ｓを求め、その縮尺Ｓによって各三次元座標（X,Y,Z）を補正する。
【００３８】
ステップ107では先ず、図４（Ａ）に示すように、隣接する計測視標Tk、Tk+1の間の直線間隔（弦の長さ）ΔT＝Tk+1−Tkの総和として、各三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱ（＝ΣΔT）を計算する（ステップ107のａ）。ただし、図示例のように筒状リングＡの内周形状が円形等の曲面である場合は、計測視標Tk、Tk+1の間の直線距離ではなく、内周面に沿った弧状間隔（弧の長さ）を計算する必要がある。図５は、内周形状が円形の筒状リングＡに取り付けた計測視標Tk、Tk+1の角度ピッチθ（＝θk+1−θk）と、その計測視標Tk、Tk+1の間の弧状間隔に対する直線間隔の間隔比α（＝弦の長さ／弧の長さ）との関係を示すグラフである（図４（Ｂ）も参照）。例えば、筒状リングＡの内周面の交線LAを48分割した場合は、角度ピッチθ＝7.5度となるので、間隔比α＝（2(1−cosθ)）1/2／θ＝0.99929となる。同図は内周形状が円形の場合の関係を示すグラフであるが、円形以外の場合も同様に角度ピッチθと間隔比αとの関係式を求めることができる。従って、同図のような関係式をコンピュータ13の記憶しておけば、座標補正手段22において隣接する計測視標Tk、Tk+1の角度ピッチθに応じた間隔比αを求めて計測視標Tk、Tk+1の間の弧状間隔（＝ΔT／α）を計算し、その弧状間隔の総和として各計測視標Tkを結ぶ多角形の周囲長さＱ（＝Σ（ΔT／α））を計算することができる。望ましくは、計測視標Tkの取り付け角度ピッチθをできるだけ小さく（例えば１〜４度程度）することにより、間隔比αを１に近づけて周囲長さＱの精度を高める。
【００３９】
またステップ107において、各計測視標Tkを結ぶ多角形の周囲長さＱと筒状リングＡの所定内周長さＧとの縮尺Ｓを計算し（ステップ107のｂ）、その縮尺Ｓに基づき各計測視標Tk及び補助視標Tsの三次元座標（X,Y,Z）を補正する。例えば筒状リングＡの内周形状が円形である場合は、各計測視標Tkを結ぶ多角形の周囲長さＱ（＝Σ（ΔT／α））と筒状リングＡの内周長さＧ（＝2πR）との縮尺Ｓ（＝2πR・α／ΣΔT）を求め、その縮尺Ｓを各視標Tk、Tsの各三次元座標（X,Y,Z）に乗算することで補正三次元座標（Xs,Ys,Zs）を計算する。この補正により各三次元座標（X,Y,Z）の座標系の単位長さを調整し、筒状リングＡの実際のスケールと一致させることができる。筒状リングＡの内周形状が円形以外であっても、この補正方法により各視標Tkの補正三次元座標（Xs,Ys,Zs）を求めることができる。なお、ステップ107では計測視標Tkだけでなく各補助視標Tsの三次元座標（X,Y,Z）も同時に補正できる。
【００４０】
ステップ107において筒状リングＡの実際のスケールと一致させた各計測視標Tk、補助視標Tsの補正三次元座標（Xs,Ys,Zs）が求まるので、その補正三次元座標（Xs,Ys,Zs）により筒状リングＡの歪み形状を計測することができる。例えば、ステップ107で求めた各計測視標Tkの補正三次元座標（Xs,Ys,Zs）を出力手段25経由でディスプレイ・プリンタ等の出力装置14に表示又は印刷し、筒状リングＡの歪み形状を確認することができる（ステップ109）。図６は、各計測視標Tkの補正三次元座標（Xs,Ys,Zs）により図１のスキンプレート２（筒状リングＡ）の歪み形状を計測した結果を示す。同図から分かるように本発明によれば、筒状リングＡの極めて精確な歪み形状を計測することができる。本発明者の実験によれば、本発明によって半径6.5ｍ程度のスキンプレート２における0.5mm程度のはらみ出しや歪みの向き及び量を把握することが可能である。
【００４１】
ステップ108は、図３に示すようにリング中心軸CA方向に隔てた筒状リングＡの内周面の複数の交線LA1、LA2、……上にそれぞれ計測視標Tkを取り付けた場合に、各交線LA1、LA2、……毎にステップ107を繰り返し、各交線LA上の計測視標Tkの三次元座標（X,Y,Z）をそれぞれ筒状リングＡの実際のスケールと一致するように補正して補正三次元座標（Xs,Ys,Zs）を求める処理を示す。各交線LA上の計測視標Tkの三次元座標（X,Y,Z）は、ステップ105〜106のバンドル調整処理により同時に同時に算出することができる。複数の交線LA上の三次元座標（X,Y,Z）についてそれぞれステップ107の補正処理を施すことにより、リング中心軸CA方向に隔てた異なる断面における筒状リングＡの精確な歪み形状を計測することができ、後述するように筒状リングＡのリング中心軸CA方向の歪み（傾斜）を計測することもできる。
【００４２】
ステップ110〜111は、筒状リングＡの歪み形状の経時的変位を計測する処理を示す。ステップ110においてステップ103〜109による各視標Tk、Tsの補正三次元座標（Xs,Ys,Zs）を前回歪み形状として記憶したのち、ステップ103へ戻って所定時間経過後の画像Igの群を撮影し、その画像Igの群に対してステップ103〜109を繰り返すことにより所定時間経過後の筒状リングＡの今回歪み形状を計測する。今回歪み形状と前回歪み形状との経時的偏差を算出することとにより、筒状リングＡの歪み形状の経時的変位を計測できる。各視標Tk、Tsの補正三次元座標（Xs,Ys,Zs）の原点の位置ずれが想定される場合は、今回歪み形状と前回歪み形状と原点を一致させる処理を行ったのち、両歪み形状の偏差を算出すればよい。例えば、図１の実施例ではスキンプレート２の交線LA上の最下点視標Tkの位置ずれは少ないと想定されるので、例えば今回歪み形状及び前回歪み形状の最下点視標Tkの補正三次元座標（Xs,Ys,Zs）を原点とした上で、両歪み形状の偏差を算出することができる。ただし、原点の位置ずれの補正処理はこの方法に限定されるものではない。
【００４３】
本発明によれば、画像Ig中に座標軸の基準となるような視標を写し込む必要がないので、撮影上の障害物が存在する見通しの悪い条件下でも筒状リングＡの精確な歪み形状を計測することができる。また、筒状リングＡに予め計測視標Tkを取り付けておけば、複数の画像Igを撮影するだけで歪み形状を継続的に計測することができるので、計測作業の大幅な省力化・迅速化が実現できる。更に、従来の計測法では困難であったリング中央点の位置を求めることができるので、筒状リングＡの詳細なはらみ出し向きや量も計測できる。
【００４４】
こうして本発明の目的である「撮影上の障害物が存在する条件下でも適用できる筒状リングの歪み形状計測方法及びプログラム」の提供を達成することができる。
【実施例１】
【００４５】
図１（Ａ）の実施例では、シールド掘削機１のスキンプレート２の内側に同心状に組み立てるセグメントリング７にも計測視標Tkを取り付け、大径筒状リングＡであるスキンプレート２の歪み形状と、小径筒状リングＢであるセグメントリング７の歪み形状とを同時に計測している。すなわち、図２のステップ101において、大径筒状リングＡの内周面だけでなく、小径筒状リングＢのリング中心軸CBの直交断面FBと交差するリング内周面の交線LB上の複数の所定角度位置θk（例えば48分割した場合は7.5度ピッチの角度位置）にもそれぞれ計測視標Tkを取り付ける。また、筒状リングＢの内周面上の任意位置に補助視標Tsを取り付けることができる。例えば、図示例のように複数のセグメントピース８により構成されるセグメントリング７の場合は、各セグメントピース８の接合部位の近傍に補助視標Tsを取り付けることができる。
【００４６】
図２のステップ102〜106において、大径筒状リングＡの視標Tk、Tsと小径筒状リングＢの視標Tk、Tsとが共に写り込む画像Igの群を撮影し、その画像Igの群に対して上述した撮影位置・姿勢の検出処理（ステップ105）、及びバンドル調整による三次元座標算出処理（ステップ106）を施すことにより、両リングＡ、Ｂ上の各視標Tk、Tsの三次元座標（X,Y,Z）を同時に算出する。またステップ107において、大径筒状リングＡの各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱと所定内周長さＧとの縮尺Ｓに基づき、両リングＡ、Ｂの各視標Tk、Tsの三次元座標（X,Y,Z）をそれぞれ実際のスケールと一致するように補正して補正三次元座標（Xs,Ys,Zs）を求める。なお、小径筒状リングＢの内周長さGBが一定である場合は、ステップ107において、大径筒状リングＡの縮尺Ｓに基づく補正に代えて、小径筒状リングＢの各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さQBと所定内周長さGBとの縮尺SBを計算し、その縮尺SBに基づき小径筒状リングＢの各視標Tk、Tsの三次元座標（X,Y,Z）を補正してもよい。
【００４７】
図２のステップ107において大径筒状リングＡ及び小径筒状リングＢの各計測視標Tkの補正三次元座標（Xs,Ys,Zs）が求まるので、その補正三次元座標（Xs,Ys,Zs）により両リングＡ、Ｂの歪み形状を同時に計測することができる。また、ステップ107において補助視標Tsの補正三次元座標（Xs,Ys,Zs）も求まるので、その補正三次元座標（Xs,Ys,Zs）により小径筒状リングＢの内周面上の所要部位の三次元位置、例えばセグメントピース８の接合部位の三次元位置を計測してピース相互間の目違い等の検出に利用することができる。
【００４８】
更に、コンピュータ13に大径筒状リングＡと小径筒状リングＢとの間のクリアランスｅを計測するクリアランス計測手段24（図１（Ｂ）参照）を設ければ、図２のステップ109において両リングＡ、Ｂの各計測視標Tkの補正三次元座標（Xs,Ys,Zs）から両リングＡ、Ｂ間の任意角度向きのクリアランスを計測することができる。例えば、大径筒状リングＡの歪み形状と小径筒状リングＢの歪み形状とを同じ平面上に表わし（図６参照）、その平面上における大径筒状リングＡと小径筒状リングＢとの間のクリアランスを計測する。また、例えばセグメントリング７である小径筒状リングＢの厚さをコンピュータ13の記憶手段15に記憶しておけば、クリアランス計測手段24により小径筒状リングＢの内周面の補正三次元座標（Xs,Ys,Zs）から外周面の歪み形状を算出し、スキンプレート２の内周面とセグメントリング７の外周面と間のテールクリアランスを計測することもできる。図６は、セグメントリング７の厚さを考慮したスキンプレート２との間のテールクリアランスを表したものである。
【実施例２】
【００４９】
図３の実施例は、大径筒状リングＡの内側に小径筒状リングＢを同心状に組み立てる場合に、両リングＡ、Ｂのリング中心軸CA、CB方向の傾斜を考慮してクリアランスを計測する実施例を示す。この場合は、図２のステップ101において、大径リングＡの内周面のリング中心軸CA方向に所定距離ｄだけ隔てた直交断面FA1、FA2と交差する一対の交線LA1、LA2上にそれぞれ所定角度位置θkの計測視標Tkを取り付けると共に、小径リングＢの内周面のリング中心軸CB方向に所定距離ｄだけ隔てた直交断面FB1、FB2と交差する一対の交線LB1、LB2上にもそれぞれ所定角度位置θkの計測視標Tkを取り付ける。またステップ102において、両リングＡ、Ｂの所定内周形状と所定内周長さＧと各計測視標Tkの取り付け角度位置θkと共に、両リングＡ、Ｂの内周面の交線対（LA1、LA2）、（LB1、LB2）の間の所定距離（間隔）ｄを記憶する。
【００５０】
図２のステップ102〜106において、大径筒状リングＡの計測視標Tkと小径筒状リングＢの計測視標Tkとが共に写り込む画像Igの群を撮影し、その画像Igの群により両リングＡ、Ｂの交線対（LA1、LA2）、（LB1、LB2）上の各計測視標Tkの三次元座標（X,Y,Z）を同時に算出する。またステップ107〜108において、大径筒状リングＡの各計測視標Tkの三次元座標（X,Y,Z）を結ぶ多角形の周囲長さＱと所定内周長さＧとの縮尺Ｓに基づき、両リングＡ、Ｂの各計測視標Tkの三次元座標（X,Y,Z）を補正して補正三次元座標（Xs,Ys,Zs）を求める。ステップ108において一対の交線LA1、LA2毎に求まる縮尺Ｓが相違する場合も考えられるが、その相違は極めて小さいと考えられるので、何れかの縮尺Ｓ又はその平均値等を用いればよい。
【００５１】
図２のステップ109において、コンピュータ13に設けた角度計測手段23（図１（Ｂ）参照）により、大径筒状リングＡの交線対（LA1、LA2）上の各計測視標Tkの補正三次元座標（Xs,Ys,Zs）とその交線対（LA1、LA2）の間隔ｄとから大径筒状リングＡの中心軸CAの傾斜を求め、小径筒状リングＢの交線対（LB1、LB2）上の各計測視標Tkの補正三次元座標（Xs,Ys,Zs）とその交線対（LB1、LB2）の間隔ｄとから小径状リングＡの中心軸CAの傾斜を求め、両リングＡ、Ｂの中心軸CA、CBの交差角度δを計測する。この交差角度δを用いれば、両リングＡ、Ｂの交線LA、LB1の歪み形状から交線LA、LB1以外の歪み形状を推定することができるので、両リングＡ、Ｂの間の中心軸CA、CB方向の任意位置におけるクリアランスｅ、すなわち中心軸CA、CB方向のクリアランスｅの変位を計測することができる。また、小径筒状リングＢの厚さをコンピュータ13の記憶手段15に記憶しておけば、小径筒状リングＢの内周面の歪み形状から外周面の歪み形状を算出することができ、例えばスキンプレート２の内周面とセグメントリング７の外周面と間の中心軸Ｃ方向のテールクリアランスを計測することもできる。
【実施例３】
【００５２】
本発明は、画像Ig中に座標軸の基準となるような視標を必要としないので撮影上の障害物が存在する見通しの悪い条件下でも筒状リングＡ、Ｂの歪み形状を計測することが可能であるが、撮影上の障害物が多い場合は、多くの計測視標Tkの像に切欠けが発生して計測精度が低下するおそれがある。本発明の歪み形状の計測では、二次元座標が精確に検出できない切欠け視標像を除外することが望ましいが（ステップ105）、バンドル調整及び補正処理（図２のステップ106〜107）の演算精度を高めるためには計測視標Tkの数を増やす必要がある。図１（Ｂ）の実施例では、切欠け視標像から誤差の少ない二次元座標を検出する座標抽出手段20を設け、計測視標Tkの数を増やすことによって歪み形状の計測精度の向上を図っている。
【００５３】
座標抽出手段20を用いる場合は、図２のステップ102において記憶手段15に各計測視標Tkの所定形状と取り付け向きとを記憶し、ステップ105において撮像位置・姿勢検出手段19により各画像Ig内の無切欠け視標像の二次元座標（x',y'）と各視標Tkの角度位置座標（X',Y',Z'）とから各画像Igの撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）を検出したのち、その撮影位置（X'0,Y'0,Z'0）及び姿勢（ω',φ',κ'）と、各視標Tkの角度位置角度位置座標（X',Y',Z'）と、各計測視標Tkの所定形状及び所定取り付け向きとから、座標抽出手段20により切欠け視標の仮想画像It（以下、テンプレートItということがある）を作成する。テンプレートItは、筒状リングＡの各計測視標Tkの像を大体の位置に大体の形状として表している。そのテンプレートItと撮影画像Igとの重ね合わせにより切欠け視標像の二次元座標（x',y'）を抽出する。
【００５４】
図７は、座標抽出手段20による切欠け視標像の二次元座標（x',y'）の抽出方法の一例を示す。座標抽出手段20は、撮影画像Ig（又は撮像位置・姿勢検出手段19により作成された２値化画像。図面では実画像と表している。）とテンプレートItとを重ね合わせ、撮影画像Ig内の切欠け視標像がテンプレートItと重なるようにテンプレートItの位置を調整し、テンプレートItの調整位置から切欠け視標像の二次元座標（x',y'）を抽出する。
【００５５】
図７（Ａ）は、テンプレートItを撮影画像Ig内の切欠け視標像の外縁位置に一致させる方法である。とくに、撮影画像Ig内の切欠け視標像における切欠け存在部位が把握できる場合に有効な方法である。例えば、撮影画像Ig内の視標像の右下部分に切欠けが存在する場合に、切欠け視標像の上縁部及び左縁部がテンプレートItと重なるようにテンプレートItの位置を調整し、調整後のテンプレートItの重心位置から切欠け視標像の二次元座標（x',y'）を抽出する。簡便な方法であるが、障害物の位置等を把握しておく必要があるため自動化は難しい。
【００５６】
これに対し図７（Ｂ）は、テンプレートItの輪郭が撮影画像Ig内の切欠け視標像の輪郭と最も一致するようにテンプレートItの位置を調整する方法である。撮影画像Ig内の切欠け視標像の周辺でテンプレートItを１画素ずつ移動させながら、両者の輪郭形状が最も合う位置を探し出す。つまり、輪郭を構成する画素境界の重複する数（辺）が最も多い場合に相当する。輪郭形状が最も合う位置に調整したテンプレートItの重心位置から切欠け視標像の二次元座標（x',y'）を抽出する。処理にやや時間がかかるが、障害物の位置等が把握できなくても適用できるので、自動化処理が可能である。
【００５７】
座標抽出手段20で切欠け視標像の二次元座標（x',y'）を抽出することにより、計測視標Tkの数を増やしてバンドル調整及び補正処理の演算精度を高めることができる。具体的には、図２のステップ106において撮影画像Ig内の無切欠け視標像の二次元座標（x',y'）と切欠け視標像の抽出二次元座標（x',y'）とに基づいてバンドル調整により各視標Tkの三次元座標（X,Y,Z）を算出し、ステップ107において座標補正手段22により補正三次元座標（Xs,Ys,Zs）を求める。撮影条件にもよるが、座標抽出手段20により計測視標Tkの数を10％程度増やすことができれば、バンドル調整の演算精度を理論上約５％向上させることができる。
【図面の簡単な説明】
【００５８】
【図１】本発明による歪み形状計測方法の一実施例の説明図である。
【図２】本発明による歪み形状計測方法の一実施例の流れ図である。
【図３】本発明による歪み形状計測方法の他の実施例の説明図である。
【図４】本発明における各視標の三次元座標の補正方法の説明図である。
【図５】各視標の三次元座標の補正に用いるパラメタの説明図である。
【図６】本発明による計測結果の一例を示す説明図である。
【図７】画像中の切欠け視標像の二次元座標を抽出する方法の説明図である。
【図８】従来のバンドル調整法を用いた画像計測方法の流れ図の一例である。
【図９】図８の画像計測方法に用いるプログラム及び視標の説明図である。
【図１０】共線条件式の説明図である。
【図１１】バンドル調整法の説明図である。
【図１２】従来のセグメントリング形状の画像計測方法の説明図である。
【符号の説明】
【００５９】
１…シールド掘進機      ２…大径筒状リング（スキンプレート）
３…シールドジャッキ    ４…カッタヘッド
５…カッタチャンバ      ６…エレクタ
７……小径筒状リング（セグメントリング）
８…セグメントピース
10…可動撮像器          11…光源
12…移動手段            13…コンピュータ
14…出力装置            15…記憶手段
16…入力手段            17…検出手段
18…二次元座標検出手段  19…撮像位置・姿勢検出手段
20…二次元座標抽出手段  21…座標算出手段
22…座標補正手段        23…角度計測手段
24…クリアランス計測手段
25…出力手段
30…計測対象            31…視標像座標検出手段
32…撮影位置・姿勢検出手段
33…バンドル調整手段    34…撮影機位置・姿勢制御手段
36…レーザ距離計
40…座標軸用反射素材    41…ＩＤ用反射素材
Ａ…（大径）筒状リング  Ｂ…（小径）筒状リング
Ｃ…中心軸              ｄ…断面間距離
ｅ…クリアランス        Ｆ…直交断面
Ｇ…内周長さ
Ig…撮影画像（入力画像）
It…仮想画像（テンプレート）
Ｌ…交線                Ｑ…多角形の周囲長さ
Ｒ…筒状リング内径      Ｓ…縮尺
Tk…計測視標            Ts…補助視標
Ta…補助視標            Tb…基準視標群
【図１】
2007147498.tif 000004
【図２】
2007147498.tif 000005
【図３】
2007147498.tif 000006
【図４】
2007147498.tif 000007
【図５】
2007147498.tif 000008
【図８】
2007147498.tif 000009
【図９】
2007147498.tif 000010
【図１０】
2007147498.tif 000011
【図１１】
2007147498.tif 000012
【図１２】
2007147498.tif 000013
【図６】
2007147498.tif 000014
【図７】
2007147498.tif 000015