第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 組の対光線が交われば、相互標定が完成する。 ステレオモデル A 像点a2 像点a1 地面 Z Y X O2 O1 f 標定をする の場合には、接続 複数のモデル 変換する 変換で写真座標に ヘルマートやアフィン を入力 基準点の絶対座標 による PhotoShop 画像処理ソフト デジタル図化機 写真座標観測は、 被写体の3D座標パスポイントから 絶対標定 相互標定 内部標定 座標の測定 基準点の画像 の取り込み カメラからデータ による撮影 デジカメ シグナル設置 基準点測量 図4.1.1.1 ステレオ写真測量の原理 図 4.1.1.2 ステレオ計測の流れ
図4.1.1.3 に示す地上座標系ＵＶＷにおける点O1, O2, p1, p2を含む空間上の平面方程式を構成し、その不定方程 式から次に示す平面の条件が得られる。 | Uo1 Vo1 Wo1 1 Uo2 Vo2 Wo2 1 U1 U2 V1 V2 W1 W2 1 1 | = 0 ...(4.1.1.1) この行列式を共面条件式という。ここに，上に示す行列式において第 1,2,3,4 行はそれぞれ点O1, O2, p1, p2の UVW座標値である。（数学的証明：coplanarity.pdf） 4.1.2 相互標定 relative orientation 撮影したステレオ画像による、左右の対光線を空間に投影すると、画像の撮影位置及びカメラの回転角は不明 のためそれらの対光線は交わらない。そのときある投影面において左右の対光線との交点間の食い違いの長さベ クトルのｘ、ｙ成分において、それぞれ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 個選んで解けばよいことになる。 投影中心を結ぶ相互標定や左側投影中心を固定する相互標定は、いずれもその条件つきの相互標定といえる。 相互標定要素は左右のカメラにおいてカメラの座標軸ｘｙｚのまわりの回転角としてそれぞれω, φ, κおよびカメ ラ投影中心の移動量bx, by, bz（bx＝1としてよい）を設定する。左カメラでは添え文字１、右カメラでは添え文字 ２とすると、相互標定要素は次のとおり書ける。 回転角 (κ1, φ1, ω1, φ2, ω2) 移動量（by1, bz1, by2, bz2） 図4.1.1.3 共面条件 図 4.1.1.4 地上写真ステレオ計測
相互標定要素の組み合わせは、次の３つのグループから合計 5 個を選定して行う。
① シフトグループ (κ1, κ2, by）から 2 個（byは 2 個選んで解くことができない。by は 1 個のみ）
② スケールグループ（φ1, φ2, bz）から 2 個(bzは 2 個選べない。1 個のみ) ③ オメガグループ （ω1, ω2）から1 個 以下には、回転の順序を入れ替えた２種類の回転行列(M)の要素を求めておく。なお、各軸まわりの回転角Mi は次のとおりとする。 Mφ= ( cosφ 0 sinφ 0 1 0 −sinφ 0 cosφ ) , Mω= ( 1 0 0 0 cosω −sinω 0 sinω cosω ) , Mκ= ( cosκ −sinκ 0 sinκ cosκ 0 0 0 1 ) ...（4.1.2.1） (M = MφMωMκ順： Stereoplanigraph C8 の場合) Mφωκ= (
cosφ ∙ cosκ + sinφ ∙ sinω ∙ sinκ −cosφ ∙ sinκ + sinφ ∙ sinω ∙ cosκ sinφ ∙ cosω
cosω ∙ sinκ cosω ∙ cosκ −sinω
−sinφ ∙ cosκ + cosφ ∙ sinω ∙ sinκ sinφ ∙ sinκ + cosφ ∙ sinω ∙ cosκ cosφ ∙ cosω
) ...(4.1.2.2) (M = MωMφMκ順： WildA7,A8 の場合)
Mωφκ= (
cosφ ∙ cosκ −cosφ ∙ sinκ sinφ
cosω ∙ sinκ + sinω ∙ sinφ ∙ cosκ cosω ∙ cosκ − sinω ∙ sinφ ∙ sinκ −sinω ∙ cosφ sinω ∙ sinκ − cosω ∙ sinφ ∙ cosκ sinω ∙ cosκ + cosω ∙ sinφ ∙ sinκ cosω ∙ cosφ
) ...(4.1.2.3)
4.1.3 投影中心を結ぶ相互標定 relative orientation of intersecting projection centers (M = MφMωMκの場合) 図 4.1.3.1 に示すように投影中心O1, O2を結んだ軸をＸ軸とすれば、使用する標定要素はκ1, κ2, φ1, φ2, ω2の5 個 となる（ω2の代わりにω1を選んでもよい）。この場合、共面条件(4.1.1)においてVo1= Vo2, Wo1= Wo2とおけるか ら、共面条件式は次のように書ける。 f = | Vo1 Wo1 1 V1− Vo1 W1− Wo1 0 V2− Vo1 W1− Wo1 0 | = |YY1 Z1 2 Z2| = 0 ...(4.1.3.1) ここで、ＸＹＺ系はUVW系に平行であり、ＵＶＷ系の座標原点を投影中心O1に平行移動させた座標系を表し ている。 回転行列 M はM＝MφMωMκの場合と仮定してモデル座標系ＸＹＺとカメラ座標系(x，y，-c)との関係を表して おく。 左カメラについては、次のように表すことができる。 (XY11 Z1 ) = Mφ1∙ Mκ1( x1 y1 −c ) ...(4.1.3.2) または
X1= (cosφ1∙ cosκ1)x1+ (−cosφ1∙sinκ1)y1+ (sinφ1)(−c)
Y1= (sinκ1)x1+ (cosκ1)y1
Z1= (−sinφ1∙ cosκ1)x1+ (sinφ1∙ sinκ1)y1+ (cosφ1)(−c)
} ...(4.1.3.3) 右カメラに関しては、次のようになる。 (XY22 Z2 ) = Mφ2∙ Mω2∙ Mκ2( x2 y2 −c ) + ( 1 0 0 ) ...(4.1.3.4) または
X2= (cosφ2∙ cosκ2+ sinφ2∙ sinω2∙ sinκ2)x2+ (−cosφ2∙ sinκ2+ sinφ2∙ sinω2∙ cosκ2)y2+ (sinφ2∙ cosω2)(−c) + 1
Y2= (cosω2∙ sinκ2)x2+ (cosω2∙ cosκ2)y2+ (−sinω2)(−c)
Z2= (−sinφ2∙ cosκ2+ cosφ2∙ sinω2∙ sinκ2)x2+ (sinφ2∙ sinκ2+ cosφ2∙ sinω2∙ cosκ2)y2+ (cosφ2∙ cosω2)(−c)
} ...(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 Z̃1 Y ̃2 Z̃2| = a1∙ δκ1+ a2∙ δφ1+ a3∙ δκ2+ a4∙ δφ2+ a5∙ δω2 ...(4.1.3.7) ここで、 Ỹ1, Ỹ2, Z̃1, Z̃2：Y1, Y2, Z1, Z2の近似値、
a
i：関数fの微係数、 mijL: 左カメラの回転行列の要素、 mijR：右カメラの回転行列の要素 である。 そうして、式(4.1.3.7)の微係数を求めると、上の式は次のように得られる。 （M＝MωMφMκの場合) 左カメラに関して、 (XY11 Z1 ) = Mφ1∙ Mκ1( x1 y1 −c ) ...(4.1.3.9) f = |Ỹ1 Z̃1 Y ̃2 Z̃2| + | Y ̃1tanκ1− y1 sinκ1 Y ̃1sinφ1 Y ̃2 Z̃2 | δκ1+ | X1 0 0 Y2| δφ1+ | Y ̃1 Z̃1 m22R x2− m21Ry2 m32R x2− m31R y2 | δκ2 + |0 X̃2− 1 Y ̃1 0 | δφ2+ | Y ̃1 c cosω2− Y ̃2tanω2 Z̃1 Ỹ2cosφ2 | δω2 …(4.1.3.8) ステレオモデルA 像点a2 像点a1 地面 Z Y X O2 O1 f 図4.1.3.1 投影中心を結ぶ相互標定（共面条件）または
X1= (cosφ1∙ cosκ1)x1+ (−cosφ1∙ sinκ1)y1+ (sinφ1)(−c)
Y1= (sinκ1)x1+ (cosκ1)y1
Z1= (−sinφ1∙ cosκ1)x1+ (sinφ1∙ sinκ1)y1+ (cosφ1)(−c)
} …(4.1.3.10) 右カメラに関して、 (XY22 Z2 ) = Mω2∙ Mφ2∙ Mκ2( x2 y2 −c ) + ( 1 0 0 ) ...(4.1.3.11) または
X2= (cosω2∙ cosκ2)x2+ (−cosφ2∙ sinκ2)y2+ (sinφ2)(−c) + 1
Y2= (cosω2∙ sinκ2+ sinω2∙ sinφ2∙ cosκ2)x2+ (cosω2∙ cosκ2− sinω2∙ sinφ2∙ sinκ2)y2+ (−sinω2∙ cosφ2)(−c)
Z2= (sinω2∙ sinκ2− cosω2∙ sinφ2∙ cosκ2)x2+ (sinω2cosκ2+ cosω2∙ sinφ2∙ sinκ2)y2+ (cosω2∙ cosφ2)(−c)
} …(4.1.3.12) したがって、線形化した共面条件式は次のように表されます。 f = |Ỹ1 Z̃1 Y ̃2 Z̃2| + |
−y1+ Ỹ1cosκ1 Ỹ1sinφ1
Y ̃2 Z̃2 sinκ1 | δκ1+ | X ̃1 0 0 Ỹ2 | δφ1+ | Y ̃1 m22R x 2− m21Ry2 Z̃1 m32R x2− m31Ry2 | δκ2+ (X̃2− 1) | Y ̃1 sinω2 Z̃1 −cosω2 | δφ2+ | Y ̃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 に示すように左カメラを鉛直に固定し、右カメラの標定要素だけで相互標定する式を導いてみる。 このような標定法を接続法という。 式 4.1.1 に示す共面条件において接続法の場合には、次のように整理することができる。 f = | Uo1 Vo1 Wo1 1 Uo2− Uo1 Vo2−Vo1 Wo2−Wo1 0 U1− Uo1 U2− Uo1 V1−Vo1 V2−Vo1 W1−Wo1 W2−Wo1 0 0 | = 0 または f = |BXX1 BYY1 BZ1Z X2 Y2 Z2 | = 0 あるいは f 像点a2 X O1 ステレオモデル A 像点a1 地面 Z Y O1 f 図4.1.4.1 接続法
f = |bx1x by1y −cbz X2 Y2 Z2 | = 0 …(4.1.4.1) ここで、 1, by, bz：実際の基線成分（BX, BY, BZ）に対応する写真基線成分(bx＝1とおく)、 x1, y1, −c(= −f):左カメラ座標系（左カメラは鉛直に固定するから、ＸＹＺ系に一致する）、 X2, Y2, Z2:ＸＹＺ系に平行な右カメラの座標系 です。 まず、左右のカメラ座標は次のように表すことができる。 （M＝MφMωMκの場合) X1= x1 Y1= y1 Z1= −c } ...(4.1.4.2)
X2 = (cosφ2∙ cosκ2+ sinφ2∙ sinω2∙ sinκ2)x2+ (−cosφ2∙ sinκ2+ sinφ2∙ sinω2∙ cosκ2)y2+ (sinφ2∙ cosω2)(−c) + 1
Y2= (cosω2∙ sinκ2)x2+ (cosω2∙ cosκ2)y2+ (−sinω2)(−c)
Z2= (−sinφ2∙ cosκ2+ cosφ2∙ sinω2∙ sinκ2)x2+ (sinφ2∙ sinκ2+ cosφ2∙ sinω2∙ cosκ2)y2+ (cosφ2∙ cosω2)(−c)
} ...(4.1.4.3) この場合の線形化した相互標定式は次のように表すことができる。 f = | 1 x1 by y1 bz −c X ̃2 Ỹ2 Z̃2| o − |X̃y1 −c 2 Z̃2| δby+ | x1 y1 X ̃2 Ỹ2| δbz+ | 1 by bz x1 y1 −c m12R x2− m11Ry2 m22Rx2− m21R y2 m32R x2− m31Ry2 | δκ2 + | 1 by bz x1 y1 −c Z̃2 0 −(X̃2− 1) | δφ2+ | 1 by bz x1 y1 −c Y ̃2sinφ2 c cosω2− Y ̃2tanω2 Ỹ2cosφ2| δω2 ...(4.1.4.3) （M＝MωMφMκの場合) 左カメラに関して、 XY11= y= x11 Z1= −c } ...(4.1.4.4) 右カメラに関して、
X2= (cosφ2∙ cosκ2)x2+ (−cosφ2∙ sinκ2)y2+ (sinφ2)(−c) + 1
Y2= (cosω2∙ sinκ2+ sinω2∙ sinφ2∙ cosκ2)x2+ (cosω2∙ cosκ2− sinω2∙ sinφ2∙ sinκ2)y2+ (−sinω2∙ cosφ2)(−c)
Z2= (sinω2∙ sinκ2− cosω2∙ sinφ2∙ cosκ2)x2+ (sinω2∙ cosκ2+ cosω2∙ sinφ2∙ sinκ2)y2+ (cosω2∙ cosφ2)(−c)
} …(4.1.4.5) この場合の線形化した相互標定式は、次のようになる。 f = | 1 x1 by y1 bz −c X ̃2 Ỹ2 Z̃2 | o − |X̃y1 −c 2 Z̃2| δby+ | x1 y1 X ̃2 Ỹ2| δbz+ | 1 by bz x1 y1 −c m12Rx2− m11Ry2 m22R x2− m21Ry2 m32R x2− m31R y2 | δκ2+ | 1 by bz −x1 −y1 c c cosφ2+ (X
̃2− 1)tanφ2 −(X̃2− 1)sinω2 (X̃2− 1)cosω2| δφ2+ |
1 by bz
x1 y1 −c
0 −Z̃2 Ỹ2
| δω2 …(4.1.4.6)
4.2 モデル座標の計算 model coordinates
4.2.1 投影中心を結ぶ場合 intersecting projection centers
カメラ座標系(xi, yi, −c)との関係は次のように表すことができる。 (XYii Zi ) = λiL∙ ML( x1 y1 −c ) ...(4.2.1.1) (XYii Zi ) = λiR∙ MR( x2 y2 −c ) + ( 1 0 0 ) ...(4.2.1.2) ここで、回転行列(M)は式(5.5.3)、または式(5.5.4)に示すものと同様である。 上の2 式を等しいとおくと、次のようになる。 λiL∙ M L( x1 y1 −c ) = λiR∙ MR( x2 y2 −c ) + ( 1 0 0 ) ...(4.2.1.3) この３つの式を解けば、左右の写真像点の縮尺係数(λiL, λiR)は次式で表せる。
λ
i L=
ym2 xm1∙zm2−xm2∙zm1λ
iR=
ym1 xm1∙zm2−xm2∙zm1}
...(4.2.1.4) したがって、縮尺係数が求まれば、式4.2.1.1 および式 4.2.1.2 からモデル座標を算出できる。
4.2.2 左カメラを固定する場合 fixed left-hand projection center
この場合ML= I (単位行列)とおけるので、式 4.2.1.1 および式 4.2.1.2 は次のように書ける。 (XYii Zi ) = λiL( x1 y1 −c ) ...(4.2.2.1) ( Xi Yi Zi ) = λiR∙ MR( x2 y2 −c ) + ( bx by bz ) ...(4.2.2.2) これらの式が等しいとおいて、それぞれ縮尺係数を求めると次のようになる。 λiL( x1 y1 −c ) = λiR∙ MR( x2 y2 −c ) + ( bx by bz ) ...(4.2.2.3) そして ただし、通常はbx= 1とおいた。 4.3 絶対標定 absolute orientation 絶対標定は、4.2 項の計算で求めた平均モデル座標値(Xm, Ym, Zm)を使用して、観測対象の絶対座標値を求める 方法である。 モデル座標系(X, Y, Z)と地上座標系(E, N, H)との関係は、次式で表すことができる。 (NE H ) = S ∙ M ( X Y Z ) + ( Xo Yo Zo ) ...(4.3.1)
λ
i L=
bx∙zm2−bz∙xm2 x1zm2+c∙xm2λ
iR=
−bx∙c−bz∙x1 x1∙zm2+c∙xm2}
...(4.2.2.4)ここで、 Ｓ：各点の縮率、 (Xo, Yo, Zo)：投影中心座標 モデルの回転行列(M)はM = MΦ∙ MΩ∙ MΚの回転順とすれば、次の要素で表される。 MΦ= ( cosΦ 0 sinΦ 0 1 0 −sinΦ 0 cosΦ ) 、MΩ= ( 1 0 0 0 cosΩ −sinΩ 0 sinΩ cosΩ ) 、MΚ= ( cosΚ −sinΚ 0 sinΚ cosΚ 0 0 0 1 ) …(4.3.2)
M = ( cosΦ ∙ cosΚ + sinΦ ∙ sinΩ ∙ sinΚcosΩ ∙ sinΚ −cosΦ ∙ sinΚ + sinΦ ∙ sinΩ ∙ cosΚcosΩ ∙ cosΚ sinΦ ∙ cosΩ−sinΩ −sinΦ ∙ cosΚ + cosΦ ∙ sinΩ ∙ sinΚ sinΦ ∙ sinΚ + cosΦ ∙ sinΩ ∙ cosΚ cosΦ ∙ cosΩ
) ...(4.3.3)
そこで、式4.4.1 を展開すると、次式になる。
fE= S[(cosΦ ∙ cosΚ + sinΦ ∙ sinΩ ∙ sinΚ)X + (−cosΦ ∙ sinΚ + sinΦ ∙ sinΩ ∙ cosΚ)Y + (sinΦ ∙ cosΩ)Z] + Xo− E = 0
fN= S[(cosΩ ∙ sinΚ)X + (cosΩ ∙ cosΚ)Y + (−sinΩ)Z] + Yo− N = 0
fH= S[(−sinΦ ∙ cosΚ + cosΦ ∙ sinΩ ∙ sinΚ)X + (sinΦ ∙ sinΚ + cosΦ ∙ sinΩ ∙ cosΚ)Y + (cosΦ ∙ cosΩ)Z] + Zo− H = 0
} ...(4.3.4) 4.3.1 近似解法 approximate solution 式4.4.4 を用い、平面座標(X,Y)と高さ(H)を分離して回転角(Φ, Ω, Κ)と投影中心の位置、すなわちモデルの移動 量(Xo, Yo, Zo)を解いてみる。 (平面位置) 式(4.4.4)においてΦ=Ω=0、Ｚ＝0 とおけば、次式が得られる。 E = S(cosΚ)X + S(−sinΚ)Y + Xo N = S(sinΚ)X + S(cosΚ)Y + Yo} ...(4.3.1.1)
この式でa = S ∙ cosΚ, b = S ∙ sinΚ, Xo= c, Yo= dとおけば、次のような平面座標の観測方程式が得られる。
(E N) = ( X −Y 1 0 Y X 0 1) ( a b c d ) ...(4.3.1.2) または行列記号で書けば F=C・P …(4.3.1.3) この式に最小 2 乗法を適用し、正規方程式(CTf) = (CTC)Pを求め、これをP = (CTC)−1(CTf)により解けば、各標 定要素は次のように計算することができる。 a2+ b2= (S ∙ cosΚ)2+ (S ∙ sinΚ)2 ∴ S = √a2+ b2 ...(4.3.1.4) ∴ tanΚ =b a ...(4.3.1.5) ∴ Xo= c ...(4.3.1.6) ∴ Yo= d ...(4.3.1.7) (高さ) 式4.3.4 の第３式において、cosΚ=cosΦ=cosΩ=1,sinΩsinΚ=sinΦ=sinΚ=0 とおけば、高さは次のように表せる。 H = S[(−sinΦ)X + (sinΩ)Y + Z] + Zo …….(4.3.1.7)
この式において、A = −S ∙ sinΦ, B = S ∙ sinΩ, C = Zoとおけば、次のように高さの観測方程式が得られる。 (H − S ∙ Z) = (X Y 1) (AB C ) ….(4.3.1.8) または C = F ∙ Q ...(4.3.1.9) 平面位置と同様にして、この正規方程式は(FTh) = (FTF)Qから求められ、その解はQ = (FTF)−1(FTh) によって 解ける。Ｑから次のように標定要素が算出できる。 ∴ Φ = arcsin (−A S) ...(4.3.1.10) ∴ Ω = arcsin (B S) ...(4.3.1.11) ∴ Zo= C ...(4.3.1.12) 4.3.2 同時調整 simultaneous solution 式 4.3.4 の各標定要素を次のように近似値と補正値とで表しておく。 (4.3.2.1) { S = S̃ + δS Φ = Φ̃ + δΦ Ω = Ω̃ + δΩ Κ = Κ̃ + δΚ Xo= X̃o+ δXo Yo= Ỹo+ δYo Zo= Z̃o+ δZo そうして、補正値方程式（観測方程式、残差方程式）を書けば次のようになる。 ( vx vy vz ) = ( E N H ) − S ∙ M (X ̃ Y ̃ Z̃ ) − ( Xo Yo Zo ) + ( a1 a2 a3 a4 a5 a6 a7 b1 b2 b3 b4 b5 b6 b7 c1 c2 c3 c4 c5 c6 c7 ) ( δS δΦ δΩ δΚ δXo δYo δZo) ...(4.3.2.2) ここで、ai, bi, ciは各変数（標定要素）の微係数である。 絶対標定において、E,N,H は観測値でないので、次のように条件方程式で表して、上式 4.3.2.2 の係数をこ の式の偏微分によって解く。 ( fE fN fH ) = S ∙ M ( X Y Z ) + ( Xo Yo Zo ) − ( E N H ) …（4.3.2.3） (4.3.2.4) { a1=∂f∂SE=E-XSo a2=∂f∂ΦE=H-Zo a3=∂f∂ΩE=(N－Yo)sinΦ a4=∂f∂ΚE=S(m12X-𝑚11N) a5=∂X∂fE o=1 a6=∂Y∂fE o=0 a7=∂Z∂fE o=0
(4.3.2.5) { b1= ∂fN ∂S = N−Yo S b2= ∂fN ∂Φ= 0 b3= ∂fN
∂Ω = −S ∙ sinΩ (X ∙ sinΚ + Y ∙ cosΚ + Z tanΩ)
b4= ∂fN
∂Κ = S ∙ cosΩ(X ∙ cosΚ − Y ∙ sinΚ)
b5= ∂fN ∂Xo= 0 b6= ∂fN ∂Yo= 1 b7= ∂fN ∂Zo= 0 (4.3.2.6) { c1= ∂fH ∂S = H−Zo S c2= ∂fH ∂Φ= −(E − Xo) c3= ∂fH ∂Ω = (N − Yo)cosΦ c4= ∂fH ∂Κ = S(m32X − m31Y) c5= ∂fH ∂Xo= 0 c6= ∂fH ∂Yo= 0 c7= ∂fH ∂Zo= 1 式4.3.2.3 から正規方程式をつくり、順次補正値を解き、絶対標定に必要な標定要素を求めていく。 (NE H ) = S ∙ M ( X Y Z ) + ( Xo Yo Zo ) …（4.3.2.7）
4.4 計測全体の流れについて Summary of digital photogrammetry
本HP の目的は、写真測量による計測法を知ることであり、またデジタル一眼レフカメラを用いて簡単に被写 体の三次元座標を求めることができるという、その計算処理について述べてきた。その計測法の全体の流れにつ いて以下にまとめておく。 （カメラのキャリブレーション） 1） 非計測用カメラ(一眼レフ)を使用することを前提にすると、カメラキャリブレーション場を設置する。検定場 に配置する基準点の精度は1/10,000 以上とする。撮影距離 H=10ｍとすると、位置誤差は 1ｍｍ程度、角度 で言うと1/10,000×206,265”＝20”読み以上の精度の TS を使用する。H=5m では 5000ｍｍ/10,000＝0.5ｍ ｍの位置誤差、角度ではやはり20”の TS を用いることができる。ただし、TS での測距精度は 2ｍｍより大 きくなるので、これを考慮しなければならない。 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.20ｍｍ×4.65ｍｍとなる。 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.
