基于通用成像模型

suesmeisf

摘要：系统地介绍有理函数模型的定义和解算方法，推导基于有理函数模型的3维定位模型，并分别用画幅式相机、线阵列CCD
传感器影像进行了试验。试验结果表明有理函数模型的拟合误差可以忽略，3维定位能够达到严格成像模型同量级的精度，并能
8 K$ W2 [( o1 p* O& v; o替代严格成像模型完成摄影测量处理，利用不含物理意义的有理函数系数有效地实现传感器成像信息的隐藏。
) B3 A" Z; ^! a% Y关键词：有理函数模型；严格成像模型；画幅式；线阵列CCD；3维定位
% W$ T8 l, B$ `& v1 e# u一
、引 言
1 W1 a% t& Y! a传感器成像模型一般可分为两类：严格成像模
  {; c, e* ]! Q! [& m型和通用成像模型。摄影测量处理常采用严格成像
模型，但其前提是已知传感器物理构造及成像方式
# p  ~3 u: K2 h' T3 X7 q! R等参数。但是，为了保护敏感的设计参数和技术秘
密不被扩散，一些高性能传感器的镜头构造、成像方
& H7 v( g: j' v1 }2 o0 t* G式、卫星轨道等信息并未被公开，如美国空间成像公
5 A  t( z# d0 O3 z+ n司(Space Imaging)的IKONOS卫星等。为此，一些
厂家如Space Imaging和政府机构如美国国防部影
& T. l8 g! ~' x" q# e: j: c像与测绘局(NIMA)等，采用通用成像模型——有
理函数模型(Rational Function Model，RFM，下同)
6 ?" h- E2 Y3 J& P# H) n4 C替代传感器严格成像模型，达到影像的商业化目的。
近年来，IKONOS等高分辨率卫星影像的广泛
应用推动了人们对有理函数模型的全方位研究。有
学者将有理函数模型分为“正解”和“反解”两种形
式_3 ；也有人将有理函数模型的建立分为“基于地
% \( r+ ?2 U9 F' D形”和“独立于地形”两种方式_l’2 ；关于利用有理函
数模型实现立体影像重叠区地面点的3维定位理
+ \" j, {5 e. x! f* \8 n9 T# M论，文献[3]提出“基于反解”的定位算法，文献[2]在
此基础上提出“基于正解”算法，但其前提是左右影
像的有理函数模型具有相同的坐标标准化参数。
本文先介绍“独立于地形”即已建立传感器严格
成像模型情况下，“正解”有理函数模型的解算方法；
' T1 C1 b/ w' a; a& c9 P重点推导了不同标准化参数情况下，有理函数模型
0 f" c& O* u! ]3 f/ e  g' @* b的3维交会算法；最后通过画幅式和线阵列CCD影
7 N0 |: o% b1 D0 h* _- H像的试验，验证和分析有理函数模型拟合严格成像
- F/ t/ G; o+ ]模型的精度以及3维定位精度。
二、有理函数模型
1．定义
有理函数模型将像点坐标．s ， 分别表示为
* ?* o3 c) `8 [+ V/ [以相应地面坐标 ，y ，z 为自变量的多项式的
5 w( L: V4 S: n) u1 a" M9 y商，多项式最高阶数一般为3次，即
O 一 ! ! ! 2
* Q. w: x8 K& `& O1 \! ]* M* On — P2(X ，Y ，Zn)
6 x4 }5 [4 x5 z, W# zr 一 ! ! ! 2
—
! {9 ]8 [7 W5 Q6 a2 J- EP4(X ，Y ，Zn)
& v5 V* b& Q% g# g# b2 ]式中的(S ，L )和(X ，Y ，Z )是像坐标(S，L)和
地面坐标( ，y，z)经平移和缩放处理后的标准化
  [$ X6 [) s. \7 D* W7 P坐标，即
S =(S—So)／s
8 X- t4 i# k# X, ^/ r1 Y2 ]L =(L—L0)／L
$ Q2 U9 l3 |6 @, c/ i/ q  F# p/ iX，r=(X—X0)／
/ |# j0 o  [( a/ f; dY =(Y—Y0)／Y
Zn=(Z—Z0)／
(2)
" h$ E4 m* ~6 }6 l# y0 O# d6 W其中，(．s0， 0， 0，y0，Z0)为标准化平移参数，(．s ，
。， 。， y。， )为标准化比例参数，对于特定的某一
9 x% |" x, P& Z2 u# _幅影像而言它们是常量。
2．有理函数系数(赋ioIlal Function Coeffi．
cients-RFCs)的确定
文献[1]的试验表明，“基于地形”即通过实际量
收稿日期：2002．11．O7
& _, a( I! a2 x; U! t作者简介：刘军(1978一)，男，河南信阳人，硕士生，主要从事数字摄影测量方面的研究。
一 T～
维普资讯 http://www.cqvip.com
& @) a( ^+ G9 T- _) l, k2003年第4期 测绘通报
+ ]9 f) ]* Q' ^9 w& d5 d测的一定数量的地面控制点计算出的RFCs受控制
8 r+ i! b6 Q2 U0 e4 }0 u1 e9 j9 t点数量、空间分布的影响很大，精度不稳定，不能替
, X# K. u* z4 A( F! q) ^& H代严格成像模型完成摄影测量处理。通常解算
6 H0 ~  |8 ^, ~+ Y1 I( ]. eRFCs是在建立传感器严格成像模型后，利用严格
模型生成密集均匀的控制格网，以格网点作为控制
. ]2 {4 h" Y8 A" |6 Z' ?3 i" ]点，按照最d~-<-乘原理计算RFCs。此时有理函数模
" M+ G  O, A7 i3 U" k5 @+ d1 j型是对严格成像模型的拟合，其目的是方便用户使
( ~& @% N4 r6 D( _1 Y% L8 ?# l7 h用及传感器信息保密，这对处理成像关系特别复杂、
难以建立严格成像模型的高分辨率传感器影像(如
2 g: A4 F% _$ k4 UIKONOS等)具有更重要的意义，因为只需厂家按
$ g/ R: z. }* [, `" l照传感器成像方式建立复杂的严格成像模型，并在
7 \' t- o+ e8 T5 ~此基础上解算出RFCs，用户直接利用厂家提供的形
式简单的有理函数模型完成摄影测量处理。
7 w0 @. W; g+ E4 _4 @" z依据传感器严格成像模型，可按以下步骤构建
控制格网：
1．将影像均匀分成10X 10的网格即121个格
! t/ g. S% D! d网点，格网点像坐标作为控制点的像坐标(S，L)；
. [8 m& n; f$ a1 d( b( p) o4 h* W1 I* @2．概略确定出影像覆盖区域的地形起伏范围
(Z |n—Z )，将[Z i ，Z ]均匀分为6层，每层具
6 I3 q! @- t8 d( F/ T& f! d有相同的高程Z，并都有上述121个均匀分布的格
, m* e* O0 r  ~% ^5 L& k网点，由此产生726个空间分布均匀的格网点，格网
点的像坐标和高程z已知；
" c0 p% \+ Q- Q4 K2 c3．根据每个控制格网点的像坐标(S，L)和高
: `* S8 ]" S$ c1 p程z，利用严格成像模型计算出相应的地面坐标．】f，
y，这样得到726个格网点的全部坐标。
4 T1 _9 v! k, [5 g* _" l, r/ e3．利用有理函数模型解算地面点大地坐标
& H' f) m$ q& _, J# i4 N, w(1)坐标改正数的计算
综合式(1)、式(2)，将像点原始坐标(S，L)表
  k- _; l$ C. o示为
S=F(X，Y，Z)+So 1 ．
L：C(X，y，z)+Loj ‘3)
( Y: Y( ^' K2 G& P2 m# q' A1 Z! P式中
$ f- s7 _- ?- W5 q! v眦 z)_s。紫
G c 若
将式(3)按照泰勒公式展开到一次项，得到地面坐标
改正数的观测方程
s =
5 p0 O2 j, t5 _' ?1 A, X( TOF△z
$ E5 N) }' G  [- @& k+
2 [' S4 S4 g8 D- Q0 q, ]3 ?OF△ y
  I. ?1 N  {$ J+
5 V) b. |3 @- q( i, ~" }OF△
一(s一§)1 ，、
L =
OG△z
& r. Z' w  c% d, x- k+
: m' w& C: i% T( EOGA Y
+
OG△．】f
一( 一￡)J 一
9 e# v1 O+ A1 L( I$ s4 o9 z3 j& A根据立体像对同名像点坐标(Sf，Lf)，(S ，L )，可
列出以下4个误差方程
Sf
S
r
& ?" V: g- ~% [1 {$ cL
L
9 }- |3 p) l  ^5 B5 MOFI／OZ
2 P/ Z. c2 Q9 e; C# VOF ／OZ
' }/ O; n! f/ ]$ aOGI／OZ
+ G1 i. P6 `+ C, J+ p: sOG ／OZ
OFl／OY
5 q0 r' W- Z9 H. ~; O3 E$ t5 \6 _OF ／0Y
; R& D5 d+ m% }OGl／OY
3 v6 v2 m( d! t2 @OG ／0Y
' r0 ?# y8 u% M! b+ kSz—S z
; P0 E; ^( ?0 j3 cS 一S
/ O- F! K; |/ {: L: J2 Pz 一 z
- d$ w  {4 V% H' |+ U5 @5 \， 一 ，
" g: \' s; U3 a: W) a7 H% POFI／OX
OF，／OX
OGz／OX
% j4 b: L% v% ?  aOG ／OX
V =AA—Z (5)
则地面坐标改正数的最小二乘解为
A：[△Z △Y ZXX]T：(ATA)一 (AT Z) (6)
; m9 b5 ^' _% l: }, r(2)两种确定地面坐标初始值(x(O)，y(O)，
Z(o))的方法
. w2 n5 k. r3 n# y方法一：将左右影像有理函数模型的地面坐标。
, t3 Q4 ?% e0 x5 n/ g: b# p标准化平移参数的平均值作为初始值，即
X(。)=(Xol+Xo )／21
Y(0)：(Yoz+Yo )／2}
* b, a% m. \: M7 u2 C8 H4 z) g6 PZ(。)=(Zoz+Zo，)／2 J
方法二：利用有理函数模型的一次项计算地面
2 {3 }  Z1 N' y7 V8 ]4 l" ]坐标初始值。由
s L】 一L】 1 + b鲁l Z n+ b 2 Y n+ b簧3 Xn +1。 o。O
： { + 。
& |2 |1 L" {$ z; F3 Z6 x& p根据同名像点坐标(sz， )，(s ， ，)，可以列出4个
方程，按照最小二乘原理计算出相应的地面坐标(．】f，
2 f5 y9 F0 R7 A) `y，z)，将其作为迭代的初始值(．】f(∞，y(m，z(o )。
6 X6 k1 {5 I9 X- E) f0 m0 E三、实验及结果分析
: P) @1 ]5 `( G1 ]$ a* f" }1．试验数据
6 L# z7 x3 ]+ M数据I(航空画幅式像片)：深圳地区航空画幅
- t$ v7 P- w4 o2 Q* J; L( H式立体像对，立体覆盖区高程范围在0 400 m之
间。用LH Systems公司的SOCET SET数字摄影
3 Z. R- x, @, {: @测量软件系统解算出左右像片的外方位元素，建立
严格成像模型；用ATE软件模块以2O m的间距生
成规格网DEM，大小为8O列×111行。
数据II(线阵列CCD扫描影像)：某地区的
! l5 ?5 x8 N" f. F& |. mSPOT卫星立体影像，在影像覆盖区所在的
1：50 000地形图上量取89个明显标志点，并在左右
影像上量测相应的像点坐标。用其中l7个点作为
地面控制点，通过空间后方交会计算出两幅影像的
* C. `3 ~% I( [0 f9 P8 U/ A外方位元素，之后根据标志点的左右影像坐标(S ，
; _/ G9 E) S) VLt)，(S ，L )，前方交会出相应的地面坐标。通过
一
1●●●●● ●J
$ w! ^% Q" D2 {/ TZ y
△ △ △
—．．．．．．．．．．．．．．．．．
; V. q  J$ I+ U5 ^L
1 w8 r' k! W0 F# B维普资讯 http://www.cqvip.com
12 测绘通报 2003年第4期
2 R) h; H! l9 K- Q+ I7 Y0 l比较标志点的前方交会坐标和地图量测坐标知，该
& b* ^7 z; q' H) y% K: {& ^& l像对严格成像模型定位的均方差分别为RMSX=
25．72 m ，RM SY =21．83 m ，RM SZ =23．75 m 。
2．试验内容
# W, x% L0 ^5 l4 T2 F" }: y9 m(1)有理函数模型建立
8 y$ m8 o% D& ?. I4 t% k从数据I、数据II中任取一幅影像，按照上述方
法建立控制格网，利用格网点坐标解算出RFCs。类
似于控制格网的建立，将地形起伏范围均分为7层
+ o) c9 W+ C& y* Q* F( t: `7 h并去掉最高和最低两层，每层以2倍于控制格网的
; o8 G1 R, o/ @3 c" H密度，共生成2 205个检查点。将检查点的地面坐
9 d2 }  O% ?: P# C: u标代人解算出的有理函数模型中，计算出新的像点
- D  Y8 }$ M1 l. ~) e; A6 i坐标，将其与原始像点坐标比较，计算均方差和最大
误差，以此为依据分析有理函数模型拟合严格成像
. n. X) o' a* W& J模型的精度。
' D/ L8 C& z, E0 \. ~(2)地面点3维定位
! b; `) V6 u- C2 N) u' E: b% r先分别建立数据I、数据II立体像对左右影像
$ _( L, A3 `0 @6 W  z. {/ F8 W' Y% t) t  L的有理函数模型。按照严格成像模型，计算出数据
I中每个DEM 格网点的左右影像坐标( z，yt)，
% E- @) q# p5 h2 Z( ，y )；之后依据有理函数模型的3维定位原理，
计算同名像点( z，yt)，( ，y )对应的地面坐标，将
结果与该DEM 格网点的原始坐标比较，分析画幅
式影像有理函数模型的定位精度。同理，依据有理
函数模型定位原理，计算数据II中各标志点的同名
像点(s z，Lf)，(s ，L )对应的地面坐标，将其与该
标志点的前方交会结果比较，分析CCD线阵列影像
的有理函数模型定位精度。
) V  f6 U4 z6 q  ^% Y6 G3 J3．结果分析
2 i0 b, F* I. m" V. S+ D/ D为充分验证有理函数模型拟合严格成像模型的
精度情况，试验中考虑了有理函数模型的不同阶数
9 A) I9 c. E! e  v6 U和不同分母组合，表1、表2分别统计了数据I和数
4 j- _/ w2 g4 Z9 u# z. v1 b" c据II的试验结果。
$ T1 ~0 E4 v& C- y表1 画幅式像片的R脚拟合统计 mm
. A% ^3 ]8 E/ R3 d1 o从表1中可以看出：
1．带有分母多项式的一阶有理函数模型具有
2 B% Z! W) F: m  ]! ^最高的拟合精度，这是因为画幅式像片的严格成像
模型(共线方程)本身就是一阶有理函数模型；
2．带有分母多项式的有理函数模型(p2≠p4
% S7 |4 H- |' N6 w$ z及p2=p4)拟合精度较普通多项式(p2=p4—1)有
) Y  s. E. d7 B# S7 s6 D明显提高，用普通多项式取代严格成像模型的精度
$ j& ]: f- H( z+ [% M$ D) h& W) f; q没有保证。
从表2中可以得出以下结论：
# R0 |- @, w  @* _: {1．3阶有理函数模型拟合线阵列CCD影像严
7 _3 X  B" b. ~格成像模型可获得最高精度，1阶有理函数模型则
9 ~3 X4 n9 p( p+ q) D相对较差；
& X& L$ o$ f' a* c$ W5 z  M2．带有分母多项式的有理函数模型的拟合精
度较普通多项式有明显提高。
基于有理函数模型的摄影测量定位试验只考虑
% F0 x* W% h" W" W0 i1 f了2阶、3阶有理函数模型(带有分母多项式)，结果
维普资讯 http://www.cqvip.com
' r$ K9 U5 {+ K+ p4 v2003年第4期 测绘通报 13
6 ~" ~0 h" u) |/ z' ]$ R如表3、表4所示。其中，表3统计的是数据I中各
DEM格网点的左右像点，按照有理函数模型3维定
2 D' D7 o2 ]) h5 V  W% a位交会出的地面坐标与实际坐标之间的误差；表4
反映的是数据II中各标志点的左右像点，按照有理
; P- c: _- G' t函数模型3维定位交会出的地面坐标与其前方交会
, |' }  v8 V6 D1 T& s7 s8 n7 r& w结果之间的误差。
  T% G! t9 \! E8 d" e' H# w& J9 A表3 画幅式立体像对有理函数模型的3维定位精度
m
; Q, V  l0 d8 ~表4 线阵列CCD立体像对的有理函数模型3维
定位精度与“后交一前交”精度比较 m
9 ?4 C' ?$ u8 I: v9 M( ?# W, l通过对表3、表4的分析及试验过程中出现的
- e7 S: M- k9 z问题，有以下几点结论：
( o* `" K3 I2 L. [: f" {8 M8 o1．按照文中给出方法建立的有理函数模型具
有与严格成像模型相当的定位精度，可以替代严格
& Q  S1 `1 }& B! x3 A) @成像模型完成摄影测量处理；同时由于RFCs不含
有物理意义，实现了对传感器成像信息的隐藏；
2．对于画幅式影像，随着其有理函数模型阶数
的提高3维定位精度有所下降，这是因为低阶的有
理函数模型与画幅式影像的严格成像模型形式更相
近；而对于线阵列CCD行中心投影影像，由于其严
格成像模型比较复杂，高阶的有理模型定位精度相
$ y1 _! O# C9 Q2 y  N) d7 z对较高；
6 w- N  }2 X" z% |4 i8 }% {9 G3．为使3维定位的迭代过程收敛，计算出的
RFCs在检查点像坐标．s，L上的拟合误差必须相
当。
四、总 结
本文在总结近两年来国外对有理函数模型的研
究成果基础上，推导出坐标标准化参数不同时有理
3 Z5 \6 X7 r7 R5 o! c, v函数模型的3维定位和算法。文中的试验证明，适
! _$ o& x4 m& I6 V4 U2 u当方式构建的有理函数模型能替代传感器严格成像
模型完成摄影测量处理，同时不损失应有的精度，并
能有效地实现传感器镜头构造、成像方式、轨道参数
9 [: T9 e# F* q* e* C. L& B! N等机密信息的隐藏。有理函数模型的应用使我国高
8 M: L8 p$ E6 h5 X- b6 f! b9 R) J5 o分辨率遥感影像的商业化成为可能。
- U( k7 G% \: i5 Z) p. j参考文献：
[1] TAO C V，HU Yong．A Comprehensive Study of the
Rational Function Model for Photogrammetric Processing
- S. A7 M+ k& A: o* O[j]．PE＆RS，2001，67：1 347-1 357．
8 q% p3 x& D! v5 ~8 E$ c12] TAO C V，HU Yong．3D Reconstruction Algorithms
with the Rational Function Model and Their Applications
( ^1 x6 c: G# ]3 Qfor IKONOS Stereo Imagery lA]．Proceedings of“High
2 P* `0 c6 e4 x2 }6 RResolution Mapping from Space 2001”[C]．Hannover：
Joint W orkshop on ISPRS Working Groups I／2，I／5 and
I／7．2001．252-263．
13] YANG Xinghe．Accuracy of Rational Function Approximation
in Photogrammetry l A]．ASPRS 2000 Annual
Conference Proceedings[C]．Washington D．C．：Amer—
clan Society for Photogrammetry an d Remote Sensing
# {1 |4 r0 E6 s( C: w3 b/ |(ASPRS)，2000．
《中国地理信息产业机构指南》开始编篡
地理信息产业作为信息产业的重要组成部分，是融全球定位系统、遥感、地理信息系统、测量学、制图学、图形图像学、计
4 Z' E! q9 J! W+ _7 S- X7 f算机科学与技术等多学科为一体，集科研、教学、生产、销售与技术服务等多种功能于一身的产业。为促进地理信息产业界的
' }& b: ]1 a( g. F9 J交流与合作，全国测绘科技信息网决定编制《中国地理信息产业机构指南》。主要介绍地理信息产业相关单位，包括地质、地
% ?- N% a0 m& S+ B4 Y% @理、煤炭、石油、水利、水电、农业、林业、海洋、地籍与国土资源、城乡勘察与规划、建筑、交通、环境等部门的大地测量、工程测
量、摄影测量、地图制图、全球定位系统、航空航天遥感、地理信息系统、计算机软硬件、