第 4 章ステレオ写真測量 stereo photogrammetry 4.1 ステレオ写真測量とは写真測量による計測は 対象物が平面の場合には単写真で計測することができるが 対象物が立体の場合にはその対象物を 2 か所以上から左右の画像が 60% 以上の重複度で重なるように撮影して行うステレオ写真測量の原理を利用する 計測対象物によっては 2 枚一組以上のステレオ画像を用いて三次元計測を行う 被写体の三次元座標は 写真面より投影中心をとおり実空間で互いに交わった光束の交点を 空間前方交会法を用いて求めるものである 通常 地図作成 (1/500~1/25,000) は ステレオ写真測量法が用いられている ステレオ写真は 撮影方法に制約があり 第 5 章 (5.1 フィルム航空カメラ 5.2DMC 航空カメラ 5.3 近接一眼レフ平行撮影 5.4 近接マルチ回転マッチング ) で説明する撮影計画を立てておく必要がある 撮影計画は 測量目的に見合った良質な写真 ( 画像 ) を撮影するために行い 精度 効率 及び経済の観点に立って作業を進めるために必要である 4.1.1 共面条件とは coplanarity condition 撮影時には撮影した左右の写真の像による光線は空間で交会している その写真を撮影時と同じ位置に置けば 左右の光線を再現して空間で交会させれば三次元計測が行える 左右の対光線を交会させる標定を 相互標定といい 解析写真測量では 共面条件 を用いる 共面条件とは 図 4.1.1.1 に示すように撮影したカメラにおいて 投影中心 O1 O2 像点 a1 a2 の 4 点が同一平面 ( 図 4.1.1.3 のハッチング面 ) に存在すれば 必ず光線 O1a1 と O2a2 は点 A で交わるという条件である 5 組の対光線が交われば 相互標定が完成する Z Y 基準点測量 シグナル設置 O1 f 像点 a1 像点 a2 O2 X デジカメ による撮影 カメラからデータの取り込み 基準点の画像座標の測定 内部標定 相互標定 写真座標観測は デジタル図化機 画像処理ソフト PhotoShop による ヘルマートやアフィン変換で写真座標に変換する 基準点の絶対座標を入力 A ステレオモデル 地面 絶対標定 パスポイントから被写体の3D 座標 複数のモデル の場合には 接続標定をする 図 4.1.1.1 ステレオ写真測量の原理 図 4.1.1.2 ステレオ計測の流れ
図 4.1.1.3 共面条件 図 4.1.1.4 地上写真ステレオ計測 図 4.1.1.3 に示す地上座標系 UVW における点 O 1, O 2, p 1, p 2 を含む空間上の平面方程式を構成し その不定方程 式から次に示す平面の条件が得られる U o1 V o1 W o1 1 U o2 V o2 W o2 1 = 0...(4.1.1.1) U 1 V 1 W 1 1 U 2 V 2 W 2 1 この行列式を共面条件式という ここに, 上に示す行列式において第 1,2,3,4 行はそれぞれ点 O 1, O 2, p 1, p 2 の UVW 座標値である ( 数学的証明 :coplanarity.pdf) 4.1.2 相互標定 relative orientation 撮影したステレオ画像による 左右の対光線を空間に投影すると 画像の撮影位置及びカメラの回転角は不明のためそれらの対光線は交わらない そのときある投影面において左右の対光線との交点間の食い違いの長さベクトルのx y 成分において それぞれ x 視差 ( 横視差 ;x-parallax) y 視差 ( 縦視差 ;y-parallax) という 相互標定は 左右のカメラの回転角及び / 又は移動量の標定要素を用い 左右のステレオ画像の縦視差を消去して行くことである 5 個の対光線を交会させるということは 5 個の標定要素を解けば ステレオ画像の相互標定が行ったことになる このことは解析法における 共面条件が成立すれば縦視差はなくなることに対応している 図 4.1.1.3 のハッチングの部分を 核理論において核面 (epipolar plane) といい 核面と写真面との交線を核線 (epipolar line) になる 空間上の点は左右の画像にある核線上に存在する 相互標定は 空間立体幾何学で説明できる 相互標定のための未知量である標定要素について次のように定義する 1 つのカメラにおいて 空間上そのカメラが移動したときの移動量 b を XYZ 座標軸のベクトル成分に分けて bx, by, bz の要素にする 次に カメラが回転 ( 空間での傾斜 ) したとき Z 軸周りの回転角 κ Y 軸周りの回転角 φ 及び X 軸周りの回転角 ωの成分を表す したがって 1 つのカメラの未知量は bx を除き以上の 5 個となる ステレオ画像では左右 2 台のカメラを使用するので 合計 10 個の未知量を持つが 相互標定では左右の対光線を空間上で交会させるので その 10 個の標定要素から条件付で 5 個選んで解けばよいことになる 投影中心を結ぶ相互標定や左側投影中心を固定する相互標定は いずれもその条件つきの相互標定といえる 相互標定要素は左右のカメラにおいてカメラの座標軸 xyzのまわりの回転角としてそれぞれω, φ, κおよびカメラ投影中心の移動量 bx, by, bz(bx=1としてよい ) を設定する 左カメラでは添え文字 1 右カメラでは添え文字 2とすると 相互標定要素は次のとおり書ける 回転角 (κ 1, φ 1, ω 1, φ 2, ω 2 ) 移動量 (b, bz 1, by 2, bz 2 )
相互標定要素の組み合わせは 次の 3 つのグループから合計 5 個を選定して行う 1 シフトグループ (κ 1, κ 2, by) から 2 個 (by は 2 個選んで解くことができない by は 1 個のみ ) 2 スケールグループ (φ 1, φ 2, bz) から 2 個 (bz は 2 個選べない 1 個のみ ) 3 オメガグループ (ω 1, ω 2 ) から 1 個 以下には 回転の順序を入れ替えた 2 種類の回転行列 (M) の要素を求めておく なお 各軸まわりの回転角 M i は次のとおりとする cosφ 0 sinφ 1 0 0 cosκ sinκ 0 M φ = ( 0 1 0 ), M ω = ( 0 cosω sinω), M κ = ( sinκ cosκ 0)...(4.1.2.1) sinφ 0 cosφ 0 sinω cosω 0 0 1 (M = M φ M ω M κ 順 : Stereoplanigraph C8 の場合 ) cosφ cosκ + sinφ sinω sinκ osφ sinκ + sinφ sinω cosκ sinφ cosω M φωκ = ( cosω sinκ cosω cosκ sinω )...(4.1.2.2) sinφ cosκ + cosφ sinω sinκ sinφ sinκ + cosφ sinω cosκ cosφ cosω (M = M ω M φ M κ 順 : WildA7,A8 の場合 ) cosφ cosκ osφ sinκ sinφ M ωφκ = ( cosω sinκ + sinω sinφ cosκ cosω cosκ sinω sinφ sinκ sinω cosφ)...(4.1.2.3) sinω sinκ cosω sinφ cosκ sinω cosκ + cosω sinφ sinκ cosω cosφ 4.1.3 投影中心を結ぶ相互標定 relative orientation of intersecting projection centers (M = M φ M ω M κ の場合 ) 図 4.1.3.1 に示すように投影中心 O 1, O 2 を結んだ軸を X 軸とすれば 使用する標定要素は κ 1, κ 2, φ 1, φ 2, ω 2 の 5 個 となる (ω 2 の代わりに ω 1 を選んでもよい ) この場合 共面条件 (4.1.1) において V o1 = V o2, W o1 = W o2 とおけるか ら 共面条件式は次のように書ける V o1 W o1 1 f = V 1 V o1 W 1 W o1 0 = Y 1 Z 1 = 0...(4.1.3.1) Y V 2 V o1 W 1 W o1 0 2 Z 2 ここで XYZ 系は UVW 系に平行であり UVW 系の座標原点を投影中心 O 1 に平行移動させた座標系を表し ている 回転行列 M は M=M φ M ω M κ の場合と仮定してモデル座標系 XYZ とカメラ座標系 (x,y,-c) との関係を表して おく 左カメラについては 次のように表すことができる X 1 Y 1 ( ) = M φ1 M κ1 ( )...(4.1.3.2) Z 1 または X 1 = (cosφ 1 cosκ 1 ) + (osφ 1 sinκ 1 ) + (sinφ 1 )() Y 1 = (sinκ 1 ) + (cosκ 1 ) }...(4.1.3.3) Z 1 = ( sinφ 1 cosκ 1 ) + (sinφ 1 sinκ 1 ) + (cosφ 1 )() 右カメラに関しては 次のようになる X 2 x 2 1 ( Y 2 ) = M φ2 M ω2 M κ2 ( y 2 ) + ( 0)...(4.1.3.4) Z 2 0 または
X 2 = (cosφ 2 cosκ 2 + sinφ 2 sinω 2 sinκ 2 )x 2 + (osφ 2 sinκ 2 + sinφ 2 sinω 2 cosκ 2 )y 2 + (sinφ 2 cosω 2 )() + 1 Y 2 = (cosω 2 sinκ 2 )x 2 + (cosω 2 cosκ 2 )y 2 + ( sinω 2 )() } Z 2 = ( sinφ 2 cosκ 2 + cosφ 2 sinω 2 sinκ 2 )x 2 + (sinφ 2 sinκ 2 + cosφ 2 sinω 2 cosκ 2 )y 2 + (cosφ 2 cosω 2 )()...(4.1.3.5) 5 個の標定要素を最小 2 乗法により求めるために 次のように未知量を近似値と補正量とによって表しておく κ 1 = κ 1 + δκ 1 φ 1 = φ 1 + δφ 1 κ 2 = κ 2 + δκ 2 φ 2 = φ 2 + δφ 2 ω 2 = ω 2 + δω 2 }...(4.1.3.6) したがって 共面条件式 (5.5.5) を上に示す 5 つの要素に関して線形化した式は 次の形で表すことができる f = Ỹ 1 Ỹ 2 ここで Z 1 = a 1 δκ 1 + a 2 δφ 1 + a 3 δκ 2 + a 4 δφ 2 + a 5 δω 2...(4.1.3.7) Z 2 Ỹ 1, Ỹ 2, Z 1, Z 2:Y 1, Y 2, Z 1, Z 2 の近似値 a i : 関数 f の微係数 m L ij : 左カメラの回転行列の要素 m ij R : 右カメラの回転行列の要素 である そうして 式 (4.1.3.7) の微係数を求めると 上の式は次のように得られる f = Ỹ 1 Z 1 + Ỹ 1tanκ 1 Ỹ sinκ 1 sinφ 1 1 δκ 1 + X 1 0 δφ Ỹ 2 Z 2 Ỹ 2 Z 0 Y 1 + 2 2 + 0 X 2 1 δφ 2 + Ỹ c 1 Ỹ cosω 2 tanω 2 2 δω 2 Ỹ 1 0 Z 1 Ỹ 2 cosφ 2 m 22 Ỹ 1 Z 1 R x 2 m R 21 y 2 m R 32 x 2 m R δκ 2 31 y 2 Z Y O1 f 像点 a1 像点 a2 O2 X A ステレオモデル 地面 図 4.1.3.1 投影中心を結ぶ相互標定 ( 共面条件 ) (M=M ω M φ M κ の場合 ) 左カメラに関して X 1 Y 1 ( ) = M φ1 M κ1 ( )...(4.1.3.9) Z 1
または X 1 = (cosφ 1 cosκ 1 ) + (osφ 1 sinκ 1 ) + (sinφ 1 )() Y 1 = (sinκ 1 ) + (cosκ 1 ) } (4.1.3.10) Z 1 = ( sinφ 1 cosκ 1 ) + (sinφ 1 sinκ 1 ) + (cosφ 1 )() 右カメラに関して X 2 x 2 1 ( Y 2 ) = M ω2 M φ2 M κ2 ( y 2 ) + ( 0)...(4.1.3.11) Z 2 0 または X 2 = (cosω 2 cosκ 2 )x 2 + (osφ 2 sinκ 2 )y 2 + (sinφ 2 )() + 1 Y 2 = (cosω 2 sinκ 2 + sinω 2 sinφ 2 cosκ 2 )x 2 + (cosω 2 cosκ 2 sinω 2 sinφ 2 sinκ 2 )y 2 + ( sinω 2 cosφ 2 )()} Z 2 = (sinω 2 sinκ 2 cosω 2 sinφ 2 cosκ 2 )x 2 + (sinω 2 cosκ 2 + cosω 2 sinφ 2 sinκ 2 )y 2 + (cosω 2 cosφ 2 )() したがって 線形化した共面条件式は次のように表されます f = Ỹ 1 Ỹ 2 + Ỹ 1 cosκ 1 Ỹ 1 sinφ 1 Z 1 + Z 2 Ỹ 2 Z 2 δκ 1 + X 1 0 sinκ 1 δφ 1 + Ỹ 1 m R 22 x 2 m R 21 y 2 0 Ỹ 2 Z 1 m R 32 x 2 m 31 R y 2 δκ 2 + (4.1.3.12) (X 2 1) Ỹ 1 sinω 2 Z 1 osω 2 δφ 2 + Ỹ 1 Z 1 Z 2 Ỹ 2 δω 2 (4.1.3.13) 4.1.4 左カメラを固定する相互標定 relative orientation of fixed left-hand projection center 図 4.1.4.1 に示すように左カメラを鉛直に固定し 右カメラの標定要素だけで相互標定する式を導いてみる このような標定法を接続法という Z Y O1 O1 像点 a2 f X f 像点 a1 A ステレオモデル 地面 図 4.1.4.1 接続法 式 4.1.1 に示す共面条件において接続法の場合には 次のように整理することができる または あるいは U o1 V o1 W o1 1 U f = o2 U o1 V o2 V o1 W o2 W o1 0 = 0 U 1 U o1 V 1 V o1 W 1 W o1 0 U 2 U o1 V 2 V o1 W 2 W o1 0 B X B Y B Z f = X 1 X 2 Y 1 Y 2 Z 1 = 0 Z 2
ここで です b x b y b z f = = 0 (4.1.4.1) X 2 Y 2 Z 2 1, by, bz: 実際の基線成分 (B X, B Y, B Z ) に対応する写真基線成分 (bx=1 とおく ),, (= f): 左カメラ座標系 ( 左カメラは鉛直に固定するから XYZ 系に一致する ) X 2, Y 2, Z 2 :XYZ 系に平行な右カメラの座標系 まず 左右のカメラ座標は次のように表すことができる (M=M φ M ω M κ の場合 ) X 1 = Y 1 = }...(4.1.4.2) Z 1 = X 2 = (cosφ 2 cosκ 2 + sinφ 2 sinω 2 sinκ 2 )x 2 + (osφ 2 sinκ 2 + sinφ 2 sinω 2 cosκ 2 )y 2 + (sinφ 2 cosω 2 )() + 1 Y 2 = (cosω 2 sinκ 2 )x 2 + (cosω 2 cosκ 2 )y 2 + ( sinω 2 )() } Z 2 = ( sinφ 2 cosκ 2 + cosφ 2 sinω 2 sinκ 2 )x 2 + (sinφ 2 sinκ 2 + cosφ 2 sinω 2 cosκ 2 )y 2 + (cosφ 2 cosω 2 )() この場合の線形化した相互標定式は次のように表すことができる...(4.1.4.3) 1 f = b y b z X 2 Ỹ 2 Z 2 o X 2 δb y + δb Z 2 X 2 Ỹ z + 2 δκ 2 R x 2 m R 11 y 2 m R 22 x 2 m R 21 y 2 m R 32 x 2 m R 31 y 2 m 12 + δφ 2 + δω c 2...(4.1.4.3) Z 2 0 (X 2 1) Ỹ 2 sinφ 2 Ỹ cosω 2 tanω 2 Ỹ 2 cosφ 2 2 (M=M ω M φ M κ の場合 ) 左カメラに関して X 1 = Y 1 = }...(4.1.4.4) Z 1 = 右カメラに関して X 2 = (cosφ 2 cosκ 2 )x 2 + (osφ 2 sinκ 2 )y 2 + (sinφ 2 )() + 1 Y 2 = (cosω 2 sinκ 2 + sinω 2 sinφ 2 cosκ 2 )x 2 + (cosω 2 cosκ 2 sinω 2 sinφ 2 sinκ 2 )y 2 + ( sinω 2 cosφ 2 )()} Z 2 = (sinω 2 sinκ 2 cosω 2 sinφ 2 cosκ 2 )x 2 + (sinω 2 cosκ 2 + cosω 2 sinφ 2 sinκ 2 )y 2 + (cosω 2 cosφ 2 )() この場合の線形化した相互標定式は 次のようになる (4.1.4.5) 1 f = b y b z X 2 Ỹ 2 Z 2 o X 2 δb y + δb Z 2 X 2 Ỹ z + 2 δκ 2 + R x 2 m R 11 y 2 m R 22 x 2 m R 21 y 2 m R 32 x 2 m R 31 y 2 m 12 c δφ c 2 + δω 2 (4.1.4.6) + cosφ (X 2 1)tanφ 2 (X 2 1)sinω 2 (X 2 1)cosω 2 2 0 Z 2 Ỹ 2 4.2 モデル座標の計算 model coordinates 4.2.1 投影中心を結ぶ場合 intersecting projection centers 左カメラの投影中心 O 1 を座標原点 (0,0,0) とする右手座標系 XYZ をモデル座標系に選ぶと この系と左右の
カメラ座標系 (x i, y i, ) との関係は次のように表すことができる X i ( Y i ) = λ L i M L ( )...(4.2.1.1) Z i X i x 2 1 ( Y i ) = λ R i M R ( y 2 ) + ( 0)...(4.2.1.2) Z i 0 ここで 回転行列 (M) は式 (5.5.3) または式(5.5.4) に示すものと同様である 上の 2 式を等しいとおくと 次のようになる x 2 1 λ L i M L ( ) = λ R i M R ( y 2 ) + ( 0)...(4.2.1.3) 0 この3つの式を解けば 左右の写真像点の縮尺係数 (λ L i, λ R i ) は次式で表せる λ i L = λ i R = y m2 x m1 z m2 x m2 z m1 y m1 }...(4.2.1.4) x m1 z m2 x m2 z m1 したがって 縮尺係数が求まれば 式 4.2.1.1 および式 4.2.1.2 からモデル座標を算出できる 4.2.2 左カメラを固定する場合 fixed left-hand projection center この場合 M L = I ( 単位行列 ) とおけるので 式 4.2.1.1 および式 4.2.1.2 は次のように書ける X i Y i ( ) = λ i Z i X i ( Y i ) = λ i Z i L ( R M R ( )...(4.2.2.1) x 2 b x y 2 ) + ( b y )...(4.2.2.2) b z これらの式が等しいとおいて それぞれ縮尺係数を求めると次のようになる そして x 2 y 2 λ L i ( ) = λ R i M R ( ただし 通常は b x = 1 とおいた b x ) + ( b y λ L i = b x z m2 b z x m2 z m2 +c x m2 λ R i = b }...(4.2.2.4) x c b z z m2 +c x m2 b z )...(4.2.2.3) 4.3 絶対標定 absolute orientation 絶対標定は 4.2 項の計算で求めた平均モデル座標値 (X m, Y m, Z m ) を使用して 観測対象の絶対座標値を求める方法である モデル座標系 (X, Y, Z) と地上座標系 (E, N, H) との関係は 次式で表すことができる E X X o ( N) = S M ( Y) + ( Y o )...(4.3.1) H Z Z o
ここで S: 各点の縮率 (X o, Y o, Z o ): 投影中心座標 モデルの回転行列 (M) はM = M Φ M Ω M Κ の回転順とすれば 次の要素で表される cosφ 0 sinφ 1 0 0 cosκ sinκ 0 M Φ = ( 0 1 0 ) M Ω = ( 0 cosω sinω) M Κ = ( sinκ cosκ 0) (4.3.2) sinφ 0 cosφ 0 sinω cosω 0 0 1 cosφ cosκ + sinφ sinω sinκ osφ sinκ + sinφ sinω cosκ sinφ cosω M = ( cosω sinκ cosω cosκ sinω )...(4.3.3) sinφ cosκ + cosφ sinω sinκ sinφ sinκ + cosφ sinω cosκ cosφ cosω そこで 式 4.4.1 を展開すると 次式になる f E = S[(cosΦ cosκ + sinφ sinω sinκ)x + (osφ sinκ + sinφ sinω cosκ)y + (sinφ cosω)z] + X o E = 0 f N = S[(cosΩ sinκ)x + (cosω cosκ)y + ( sinω)z] + Y o N = 0 f H = S[( sinφ cosκ + cosφ sinω sinκ)x + (sinφ sinκ + cosφ sinω cosκ)y + (cosφ cosω)z] + Z o H = 0...(4.3.4) } 4.3.1 近似解法 approximate solution 式 4.4.4 を用い 平面座標 (X,Y) と高さ (H) を分離して回転角 (Φ, Ω, Κ) と投影中心の位置 すなわちモデルの移動 量 (X o, Y o, Z o ) を解いてみる ( 平面位置 ) 式 (4.4.4) において Φ=Ω=0 Z=0 とおけば 次式が得られる E = S(cosΚ)X + S( sinκ)y + X o N = S(sinΚ)X + S(cosΚ)Y + Y o }...(4.3.1.1) この式で a = S cosκ, b = S sinκ, X o = c, Y o = d とおけば 次のような平面座標の観測方程式が得られる a ( E Y 1 0 ) = (X N Y X 0 1 ) ( b c )...(4.3.1.2) d または行列記号で書けば F=C P (4.3.1.3) この式に最小 2 乗法を適用し 正規方程式 (C T f) = (C T C)Pを求め これをP = (C T C) 1 (C T f) により解けば 各標定要素は次のように計算することができる a 2 + b 2 = (S cosκ) 2 + (S sinκ) 2 S = a 2 + b 2...(4.3.1.4) tanκ = b a...(4.3.1.5) X o = c...(4.3.1.6) Y o = d...(4.3.1.7) ( 高さ ) 式 4.3.4 の第 3 式において cosκ=cosφ=cosω=1,sinωsinκ=sinφ=sinκ=0 とおけば 高さは次のように表せる H = S[( sinφ)x + (sinω)y + Z] + Z o.(4.3.1.7)
この式において A = S sinφ, B = S sinω, C = Z o とおけば 次のように高さの観測方程式が得られる A (H S Z) = (X Y 1) ( B).(4.3.1.8) C または C = F Q...(4.3.1.9) 平面位置と同様にして この正規方程式は (F T h) = (F T F)Q から求められ その解は Q = (F T F) 1 (F T h) によって 解ける Q から次のように標定要素が算出できる Φ = arcsin ( A S )...(4.3.1.10) Ω = arcsin ( B S )...(4.3.1.11) Z o = C...(4.3.1.12) 4.3.2 同時調整 simultaneous solution 式 4.3.4 の各標定要素を次のように近似値と補正値とで表しておく (4.3.2.1) S = S + δs Φ = Φ + δφ Ω = Ω + δω Κ = Κ + δκ X o = X o + δx o Y o = Ỹ o + δy o { Z o = Z o + δz o そうして 補正値方程式 ( 観測方程式 残差方程式 ) を書けば次のようになる δs δφ v x E X X o a 1 a 2 a 3 a 4 a 5 a 6 a 7 δω ( v y ) = ( N) S M ( Ỹ) ( Y o ) + ( b 1 b 2 b 3 b 4 b 5 b 6 b 7 ) δκ v z H Z Z o c 1 c 2 c 3 c 4 c 5 c 6 c 7 δx o δy o ( δz o ) ここで a i, b i, c i は各変数 ( 標定要素 ) の微係数である...(4.3.2.2) 絶対標定において E,N,H は観測値でないので 次のように条件方程式で表して 上式 4.3.2.2 の係数をこ の式の偏微分によって解く f E X ( f N ) = S M ( Y) + ( f H Z X o Y o E ) ( N) (4.3.2.3) Z o H (4.3.2.4) { a 1 = f E = E-X o S S a 2 = f E Φ =H-Z o a 3 = f E Ω =(N-Y o)sinφ a 4 = f E Κ =S(m 12X-m 11 N) a 5 = f E X o =1 a 6 = f E Y o =0 a 7 = f E Z o =0
(4.3.2.5) { b 3 = f N Ω b 1 = f N = N Y o S S b 2 = f N Φ = 0 = S sinω (X sinκ + Y cosκ + Z b 4 = f N = S cosω(x cosκ Y sinκ) Κ b 5 = f N X o = 0 b 6 = f N Y o = 1 b 7 = f N Z o = 0 tanω ) (4.3.2.6) { c 1 = f H = H Z o S S c 2 = f H Φ = (E X o) c 3 = f H Ω = (N Y o)cosφ c 4 = f H Κ = S(m 32X m 31 Y) c 5 = f H X o = 0 c 6 = f H Y o = 0 c 7 = f H Z o = 1 式 4.3.2.3 から正規方程式をつくり 順次補正値を解き 絶対標定に必要な標定要素を求めていく X o E X ( N) = S M ( Y) + ( Y o ) H Z Z o (4.3.2.7) 4.4 計測全体の流れについて Summary of digital photogrammetry 本 HP の目的は 写真測量による計測法を知ることであり またデジタル一眼レフカメラを用いて簡単に被写体の三次元座標を求めることができるという その計算処理について述べてきた その計測法の全体の流れについて以下にまとめておく ( カメラのキャリブレーション ) 1) 非計測用カメラ ( 一眼レフ ) を使用することを前提にすると カメラキャリブレーション場を設置する 検定場に配置する基準点の精度は 1/10,000 以上とする 撮影距離 H=10mとすると 位置誤差は 1mm 程度 角度で言うと 1/10,000 206,265 =20 読み以上の精度の TS を使用する H=5m では 5000mm/10,000=0.5m mの位置誤差 角度ではやはり 20 の TS を用いることができる ただし TS での測距精度は 2mmより大きくなるので これを考慮しなければならない 2) カメラの画面の大きさ 概略焦点距離 (f) と撮影距離 (H) との関係を考慮し 基準点位置を決め TS 等で正確に検定場の XYZ 座標 ( 数学座標系右手系 ) を算出する 3) その基準点位置には TS の反射板を貼り付け これを撮影の際のシグナル ( 標識 ) とする 4) 検定場の基準点数は 8 点以上 20 点ぐらいとする 5) 使用カメラにより検定場の基準点 ( ターゲット 標識 ) を撮影する 6) 基準点の画像座標値をコンパレータにより計測し XYZ 座標値を用いて 共線条件式 (DLT) により焦点距離 主点位置 カメラレンズ収差を算出する ( 一眼レフにより計測対象物の撮影 )
7) 対象物周辺に任意の座標系 ( 平面直角座標系 ) または GPS/TS 測量機によりで測地座標系を設定し 5 点以上の基準点を設ける 8) 基準点には サイズ 背景色を考慮して 必ず標識板 ( シグナル ターゲット マーク ) を取り付ける 9) デジタル一眼レフカメラを用いて 異なる場所 (3 カメラ点以上 ) から計測対象物を撮影する このとき 基線比は 1/10<B/H<1/2 となるように撮影する 10) 撮影した画像をパソコンに取り込む ( デジタル写真測量 ) 11)PhotoShop などのアプリケーションを用い 写真画像上のマーク及び求めたい点を計測し 画像座標値 (x,y) を記録する 12)2.4 項に示したヘルマート変換式 又はアフィン変換式を用いて Abobe 社の PhotoShop などで計測した画像座標値を写真座標値に変換する この場合 写真画像座標系は CCD サイズを実寸単位で求めた大きさである たとえば 3648 2736 で 1 画素当たり 1.7 ミクロンの場合 実寸は 6.20mm 4.65mmとなる 13) 写真座標値を用い 相互標定を行う 14) モデル座標の計算を行う 15) モデル座標値と絶対座標値を用い 絶対標定を行う 16) ベクトル計測はステレオ観測により実施する ベクトル ( 等高線 ) より DEM を作成するか ステレオモデルの自動相関により DEM を構築する 17) 元画像と DEM によりオルソ画像を作成する ( 文章作成 1997 年 修正更新追加 2009 年 5 月 11 月 2010 年 3 月 22 日作成者小林和夫 ) 参考文献小林和夫 上野郁夫 ポケコンプログラムによる測量計算法 山海堂 1994 年 2 月 Kobayashi, K., and C. Mori, Relations between the coefficients in the projective transformation equations and the orientation elements of a photograph, PE & RS, Vol.63, No.9, September 1997, pp. 1151-1127.