摄影测量基本公式

1) 空间直角坐标变换

用像片上的点坐标解求相应地面点坐标时，在建立起各种空间直角坐标系的基础上，需要在两两不同的坐标系中进行坐标变换，由高等数学可知，空间直角坐标的变换是正交变换，一个坐标系按三个角元素顺次旋转即可变换为另一个同原点的坐标系。

设像点a在像空间坐标系中的坐标为(x，y，-f)，而在像空间辅助坐标系中的坐标为(X，Y，Z)，两者之间的正交变换关系可以用下式表示：

Xxa1YRyb1Zfc1a2b2c2a3xy b3c3fR是一个33的正交矩阵，由九个方向余弦组成。

1) 2) 像点与相应地面点之间的坐标关系

如图（8-7）O-XYZ为选定的地面右手坐标系，地面

点A与投影中心S在此坐标系中的坐标分别为X, Y, Z和XS, YS, ZS；A在像片的构像a，在像空间坐标系与像空间辅助坐标系中的坐标分别为(x, y, -f)和(X’, Y’, Z’)，其中像空间辅助坐标系s-xyz与地面坐标系O-XYZ的相应坐标轴分别平行，由于s, a, A三点位于一直线上，由图中各相似三角形关系可得

X'Y'Z'1

XXsYYsZZsXXsX'以矩阵表示为YYsY' ZZsZ'有由～式可得像点的像空间坐标系与像空间辅助坐标系的关系为

XXsa1YYbs1ZZsc1上式的逆变换为

a2b2c2a3xy b3c3f1

xXXsa111yR1YYbs1fZZsc1整理可得

a2b2c2a3XXsYY b3sc3ZZsxfyfa1(XAXS)b1(YAYS)c1(ZAZS)a3(XAXS)b3(YAYS)c3(ZAZS)

a2(XAXS)b2(YAYS)c2(ZAZS)a3(XAXS)b3(YAYS)c3(ZAZS)上式即为像点与物点中心投影构像方程。由于像点、投影中心和地面点三点共线，故又称共线条件方程式，此式是摄影测量中重要基本公式之一，应用十分广泛。

2) 空间后方交会

利用地面控制点来解算像片的外方位元素的工作称为空间后方交会。可以利用共线条件方程进行空间后方交会。

xfyfa1(XAXS)b1(YAYS)c1(ZAZS)a3(XAXS)b3(YAYS)c3(ZAZS)

a2(XAXS)b2(YAYS)c2(ZAZS)a3(XAXS)b3(YAYS)c3(ZAZS)其中，X, Y, Z为地面控制点的地面测量坐标或地面摄影测量坐标，可用普通测量的方法得到，或通过坐标变换取得，x, y为地面控制点的相应像点坐标，可量测得到。f为摄影主距，已知。未知数为六个外方位元素。由于共线条件方程的严密关系式是非线性函数不便于计算机迭代计算，因此要将共线条件方程线性化，一般用台劳级数展开取至小一次项，

xxxxxxx(x)XYZΔωκSSSXSYSZSωκ yyyyyyXSYSZSΔκy(y)XYZωκSSS其中(x), (y)为函数的近似值，可用未知数的初始值XS,YS,ZS,,,代入共线条件方程得到。XS,, YS,, ZS,, , , 为外方位元素初始值的改正数，解算过程需反复迭代

运算。

竖直摄影下的空间后交会公式：

000000 2

fxx2xyx(x)XZf(1)ySS2HHff fyxyy2y(y)YSZSf(12)xHHff这就是后方交会的实用公式，计算时至少需三个地面控制点，每个控制列两个方程，从

而解算出6个外方位元素，若有多余控制点时，把像点坐标作为观测值，可把控制点坐标作为真值，列误差方程，按最小二乘法平差计算。误差方程式为：

fxx2xyVXZf(1)ylxSSx2HHff 2fyxyyVYZf(1)xlyySS2HHffa1(XXs)lx(x)xfx式中，

a(XX)sly(y)yf1y写成矩阵形式：V=AX-L 相应法方程式：ATAX=ATL 解：X=(ATA)-1ATL 3) 空间前方交会：

如图8-8，过左右摄站点分别做像空

间辅助坐标系S1-X1Y1Z1和S2-X2Y2Z2，使它们的坐标轴与地面上选定的一空间坐标系D-XYZ的相应坐标轴平行。地面

任意点A和左右投影中心S1与S2在地面上的空间直角坐标系中的坐标分别为XA，YA，ZA、XS1，Y S1，Z S1和X S2，Y S2，Z S2。像点a1和a2为A点在左右像片上的构像，在各自像空间辅助坐标系中的坐标为X1，Y1，Z1和X2，Y2，Z2，A在左片像空间辅助坐标系S1-X1Y1Z1中的坐标为XA-XS，YA-YS，ZA-ZS，从而有

S1AXAXS1YAYS1ZAZS1N1SaX1Y1Z111 S2AXAXS2YAYS2ZAZS2N2SaXYZ22222

3

X1x1其中：YX21R1y1x2YRy2Z1f22f

Z2XAXS1X1XS2X2XS1Bx得YNNX2AYS11Y1YS22Y2YS1ByN2Y2 ZAZS1Z1ZS2Z2ZS1BzZ2由上式得 Bx = N1X1-N2X2 (1) By = N1Y1-N2Y2 (2) Bz = N1Z1-N2Z2

(3)

取(1)、(3)式或(1)、(2)式，解得

上两式就是用立体像对确定地面点坐标得计算公式，也叫空间前方交会公式。 4