专利名称：卫星影像二/三级产品的有理函数成像模型生成方法
技术领域：
本发明属于摄影成像技术领域，特别涉及卫星影像二/三级产品的有理函数成像 模型生成方法。
背景技术：
为了满足不同用户对卫星遥感影像不同的需求，本领域有专家提出对卫星 影像进行几何产品分级。但目前卫星的严密成像几何模型均是针对卫星的一级 产品而言，对于卫星的其他级别的几何产品，并没有对其严密成像几何模型的 研究，因此卫星其它级产品未能得到实际有效应用。但是，随着各种遥感卫星 的不断发展，遥感卫星几何产品的应用水平不断的提高，对卫星其它级产品进 行进一步的精纠正以便应用成为必要。
传统的卫星几何产品的进一步几何纠正方式一般利用参考影像选取适量的 控制点，进行多项式纠正。该种纠正方式对高程差引起的投影差表现并不理想， 已不能满足用户的应用精度需求，特别是新型传感器的发展，侧视、前后视线 阵传感器的出现，以及在高程起伏大的地区，高差引起的投影差较大。
由于严密成像几何模型具有较高的定位精度，因此一直是摄影测量学的首 选，卫星遥感影像处理常常要利用严密成像几何模型。然而， 一些高性能的传 感器系统虽然实现了商业化，但是由于卫星平台和有效载荷等关键技术参数是 航天大国的技术机密，高分辨率传感器的核心信息和卫星轨道参数并未公开，
如IKONOS等卫星，无法利用严格成像几何模型进行处理；同时传感器成像方
7式的多样化也对摄影测量软件提出了模型通用化的要求，这样在软件框架设计 时就不必为现有的或将来可能会出现的各种类型的传感器一一确定其严格几何 模型的形式，而可以统一采用通用模型进行处理，从而大大降低了程序设计的 复杂性，更易于软件升级和维护，特别在同时处理多源卫星遥感影像数据时， 通用成像几何模型更能显示出其优势。
常见的通用成像几何模型有以下几种多项式、直接线性变换、仿射变换、 有理函数成像模型(RPC)等。RPC模型是卫星遥感影像的通用成像几何模型 其中的一种，它适用于各类传感器包括最新的航空和航天传感器模型。基于
RPC的成像模型并不要求了解传感器的实际特性和成像过程，是一种能获得和
严格成像模型近似一致精度的形式简单的广义成像模型。
RPC模型具体定义关系式如下
，M/附,CP，Z，//) —Z)e"、(尸，丄,//)
其中(U)为正则化的影像坐标，(P，丄,//)为正则化的地面点坐标。 M/mJ尸,丄,//)、 Z)e"丄(尸,A丑)、Ww附,(尸,4//)、 "e"、(尸,丄,/0为三次多项式
+ <7|()//2 + nJPi// + fli2Z3 +fll3ZP2 +o,4i//2十a,5丄2p + a,6户3 + (,7尸//2
(户,i， = 62丄+ 63P + 64// + 65iP + 66Zi/ + 67 W/ + 68￡2 + 69尸2
+ 6l0//2 +Z>1,P￡// + 612i3 +613丄尸2 +614i//2 +6l5i275 + 616P3十Z^P//2
M/附'(户,1,//) = c, + c2i + c3/ + c4// + c5iP + c6Z// + c7P/f + c8i2 + c,P2
+ c1()//2 +c,1JPZ7/ + c,2i3 +<:|3丄尸2 +cMi//2 +￡715丄2尸+。6户3 +c17/>//2Z)e"、.(P,￡,//)=《+ d2￡ + ^尸+ +《iP + c/6ZJf + <i7JW/ + fi s丄2 + t^P2 其中，三次多项式的系数"p…,fl2。， A,…Ao， q,…,C20， 4,…,《。是RPC文件
中提供的模型参数，^和《通常为1，"分别为影像列数值和行数值。
所谓正则化，是一项处理无限大、发散以及一些不合理表示式的方法，其
方法透过引入一项辅助性的概念——正则化因子。在RPC模型中运用这种方法
对地面点和影像点的坐标进行处理。 正则化地面点坐标定义公式为
_ z 4r—SCALE
W _ //e妙,-服/GOT — OFF _ ~"服G7/71 —SCALE~~
其中，正贝U化因子丄v4r —OFF 、丄^r —SC4丄5、丄CWG —OFF、丄CWG —SC4Z五、 、 —SC4ZE是RPC文件中包含的地面点坐标正则化模型参
数。Latitude表示经度、Longitude表示纬度、Height表示高程，此三项即可代表 地面点的空间坐标。注某点沿铅垂线方向到大地水准面的距离，通常称为绝 对高程或海拔，本技术领域简称高程。
正则化影像坐标定义公式为
y _ 1/we-丄/JV五—OFF 一 i//V￡_SC/LLE
其中，正则化因子&4MP — (9FF、 5^MP_SC^L￡、 ￡/7V￡_(9FF、 —SC4丄五也
是RPC文件中包含的影像坐标正则化模型参数，sample代表影像列坐标，其数值即为"line代表影像行坐标，其数值即为/。
基于RPC模型具有模拟精度高，通用性好，应用方便，计算量小等等优点， 所以研究卫星各级产品有理函数成像模型有着重要的实用价值。
发明内容
本发明目的在于解决现有技术不足，提出了卫星影像的二级、三级产品的 几何模型这样概念和建立这些几何模型的方法。解决了目前卫星影像各级几何 产品进一步精纠正中无法应用其几何模型进行纠正这一卫星遥感影像应用中的 瓶颈问题。
本发明的技术方案如下
一种卫星影像二级产品的有理函数成像模型生成方法，所述二级产品是将 卫星拍摄所得一级影像按照卫星的地面分辨率投影在地球椭球面上得到的几何 产品，投影关系采用投影变换公式体现；所述一级影像是对遥感卫星拍摄图像 辐射纠正，但不进行几何处理得到的原始图像产品；二级产品提供的卫星产品 影像为二级影像；包括以下步骤，
步骤1.1，由某点的大地经纬度坐标(/^2,A。n2,/0通过一级影像严密数学模型反变 换公式或一级影像有理函数成像模型反变换公式，计算该点在一级影像上的像
素坐标(Xp为,/2);
步骤1.2，调整像素坐标(、,乂，/0中的高程^为/w。，通过一级影像严密数学模 型正变换公式或一级影像有理函数成像模型正变换公式，将像素坐标(x,^,U 变换为大地坐标(A。", Z)则， re);
步骤1.3，根据投影变换公式将大地坐标(D一D^A。,。)变化为投影面坐标步骤1.4，根据二级影像的影像信息进行坐标转换，由投影面坐标^自力"。—计算
二级影像像素坐标(A,;g;
步骤1.5，根据步骤1.1、 1.2、 1.3、 1.4所得结果构建物方均匀分布的控制点格 网，包括以下步骤，
步骤1.5.1，根据一级影像的覆盖范围，求出二级影像的覆盖范围；
步骤1.5.2，根据二级影像的覆盖范围，以给定的经纬度间隔，建立平面格
网；
步骤1.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高 程，根据这些高程确定二级影像的覆盖范围内的最小和最大的高程，并根据给 定的高程分层数建立立体格网；
步骤1.5.4，根据步骤l.l、 1.2、 1.3、 1.4建立的流程得到二级严密几何模型
的反算公式(^，h)=:r2-、A。,，z^,;o ，其中r2-'表示由大地经纬度坐标(/)/。,,^。 ,/0反算
到二级影像面坐标^，h)的转换关系；根据该反算公式算出步骤5.3所得立体格 网中各格网点所对应的像点坐标，从而建立物方均匀分布的控制点格网；
步骤1.6，根据步骤1.5建立的控制点格网求解二级产品的有理函数成像模型参
、W,
数；
步骤1.7，根据步骤1.6所得二级产品的有理函数成像模型参数，生成二级产品 的有理函数成像模型。
一种卫星影像三级产品的有理函数成像模型生成方法，所述三级产品是将 卫星拍摄所得一级影像经过控制点进行系统几何校正，并以按照卫星的地面分辨率投影在地球椭球面上得到的几何产品，投影关系采用投影变换公式体现； 所述一级影像是对遥感卫星拍摄图像辐射纠正，但不进行几何处理得到的原始 图像产品；三级产品提供的卫星产品影像为三级影像；包括以下步骤， 步骤2.1，由某点的大地经纬度坐标(A。。,i^3，/z)通过一级影像严密数学模型反变 换公式或一级影像有理函数成像模型反变换公式，计算该点在一级影像上的像
素坐标(Wi，力)；
步骤2.2，根据一级影像上的控制点，由像素坐标(x,，乂,/0通过仿射变换计算该 像点在纠正后对应的像面坐标(s讓;^,/^^);调整高程/z为/w，通过一级影像
严密数学模型正变换公式或一级影像有理函数成像模型正变换公式，将像素坐 标(sawp/e，//"e，/z)变换为大地坐标(D'to3，Z)'/。 3,;
步骤2.3，根据投影变换公式将大地坐标("。,3,"。 3,^变化为投影面坐标
e加，Z)w()/v/ );
步骤2.4，根据三级影像的影像信息进行坐标转换，由投影面坐标("自力"。一计算 三级影像像素坐标(x;,力)；
步骤2.5，根据步骤2.1、 2.2、 2.3、 2.4所得结果构建物方均匀分布的控制点格 网，包括以下步骤，
步骤2.5.1,根据一级影像的覆盖范围，求出三级影像的覆盖范围；
步骤2.5.2，根据三级影像的覆盖范围，以给定的经纬度间隔，建立平面格
网；
步骤2.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高 程，根据这些高程确定三级影像的覆盖范围内的最小和最大的高程，并根据给 定的高程分层数建立立体格网；步骤2.5.4，根据步骤2.1、 2.2、 2.3、 2.4所得结果建立三级严密几何模型的 反算公式(x3,^卜7rcD,。,z^,;o，其中7T表示由大地经纬度坐标反算到三级影像 面坐标转换关系；根据该反算公式算出步骤2.5.3所得立体格网中各格网点所对 应的像点坐标，从而建立物方均匀分布的控制点格网；
步骤2.6，根据步骤2.5建立的控制点格网求解三级产品的有理函数成像模型参 数；
步骤2.7，根据步骤2.6所得三级产品的有理函数成像模型参数，生成三级产品 的有理函数成像模型。
另一种卫星影像三级产品的有理函数成像模型生成方法，所述三级产品是 将卫星拍摄所得一级影像经过控制点进行系统几何校正，并以按照卫星的地面 分辨率投影在地球椭球面上得到的几何产品，投影关系采用投影变换公式体现； 所述一级影像是对遥感卫星拍摄图像辐射纠正，但不进行几何处理得到的原始 图像产品；三级产品提供的卫星产品影像为三级影像；包括以下步骤， 步骤3.1，由某点的大地经纬度坐标(D,。,3,/),。"a)通过二级影像严密数学模型反变 换公式或二级影像有理函数成像模型反变换公式，计算该点在二级影像上的像
素坐标(《，^，W;所述二级影像是二级产品提供的卫星产品影像，二级产品是 将卫星拍摄所得一级影像按照卫星的地面分辨率投影在地球椭球面上得到的几
何产品；
步骤3.2，根据二级影像上的控制点，由像素坐标";，^，W通过仿射变换计算该 像点在纠正后对应的像面坐标""A力)；调整高程A为t^通过二级影像严密 数学模型正变换公式或二级影像有理函数成像模型正变换公式，将像素坐标 (A，J^，W变换为大地坐标(D'w，D'—U ;步骤3.3，根据投影变换公式将大地坐标(D^,D;3，U变化为投影面坐标
步骤3.4，根据三级影像的影像信息进行坐标转换，由投影面坐标""一)计算 三级影像像素坐标(^,;g;
步骤3.5，根据步骤3.1、 3.2、 3.3、 3.4所得结果构建物方均匀分布的控制点格 网，包括以下步骤，
步骤3.5.1,根据二级影像的覆盖范围，求出三级影像的覆盖范围；
步骤3.5.2，根据三级影像的覆盖范围，以给定的经纬度间隔，建立平面格
网；
步骤3.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高 程，根据这些高程确定三级影像的覆盖范围内的最小和最大的高程，并根据给 定的高程分层数建立立体格网；
步骤3.5.4，根据步骤3.1、 3.2、 3.3、 3.4所得结果建立三级严密几何模型的
反算公式(13,少3) = 7^(/^,^。 ,/0 ，其中 r表示由大地经纬度坐标反算到三级影像
面坐标转换关系；根据该反算公式算出步骤3.5.3所得立体格网中各格网点所对 应的像点坐标，从而建立物方均匀分布的控制点格网；
步骤3.6，根据步骤3.5建立的控制点格网求解三级产品的有理函数成像模型参 数；
步骤3.7，根据步骤3.6所得三级产品的有理函数成像模型参数，生成三级产品 的有理函数成像模型。
本发明生成了卫星二级、三级产品的有理函数成像模型，提供通用几何模型对各级影像的严密模型进行模拟，使各模型能方便应用。所提供技术方按不 仅解决了现有严密模型只针对一级产品而产生的二级、三级产品无法用模型精 纠正的问题，而且使用有理函数成像模型使得各级产品商业化应用中优势明显。 本发明提供的方案将可在卫星遥感影像产品几何应用方面得到广泛应用，特别 是二级、三级的带几何模型的精纠正的应用。
图1是本发明实施例一流程图； 图2是本发明实施例二流程图； 图3是本发明实施例三流程图。
具体实施例方式
本发明针对的卫星产品分级情况为 一级影像是对遥感卫星拍摄图像辐射 纠正，但不进行几何处理得到的原始图像产品，也叫一级产品。二级产品是将 卫星拍摄所得一级影像按照卫星的地面分辨率投影在地球椭球面上得到的几何 产品，二级产品提供的卫星产品影像即为二级影像。三级产品是将卫星拍摄所 得一级影像经过控制点进行系统几何校正，并以按照卫星的地面分辨率投影在 地球椭球面上得到的几何产品，三级产品提供的卫星产品影像即为三级影像。
系统几何校正时一般用5或6个控制点。
本发明中，卫星二级、三级产品的有理函数成像模型，都是基于大地经纬 度坐标—前级影像上的像素坐标一调整高程后再求取大地经纬度坐标一投影面
坐标—本级影像像素坐标，然后建立物方均匀分布的控制点格网，求取出RPC 参数。因此各级产品的有理函数成像模型生成技术方案具有单一性，在本发明中要求保护。
具体来说，分别是利用二级几何模型的反变换建立物方均匀分布的控制点 格网，建立二级几何产品的RPC模型；利用有一级生成的三级产品的反变换建 立物方均匀分布的控制点格网，建立三级几何产品的RPC模型；利用二级生成 的三级产品的反变换建立物方均匀分布的控制点格网，建立三级几何产品的 RPC模型。具体实施时可采用计算机程序完成，以下结合实施例和附图1 3提 供详细说明本发明技术方案。
(1)实施例一提供卫星影像二级产品的有理函数成像模型生成方法，参见附图 h
步骤1.1，由某点的大地经纬度坐标(Z^,2,A。"2，/z)通过一级影像严密数学模型反变
换公式或一级影像有理函数成像模型反变换公式，计算该点在一级影像上的像
素坐标(X,乂，A)。
大地经纬度坐标的含义在大地坐标系中，某点大地经度，就是通过该点 的子午面与起始子午面所构成的二面角。由起始子午面起算，向东为正，叫东 经，向西为负，叫西经。这点的法线与赤道面的夹角，叫做大地纬度，由赤道 面起算，向北为正，叫北纬，向南为负，叫南纬。由大地经度和大地纬度所表 示的地面点坐标就成为经纬度坐标。经纬度坐标加上高程即可标注地面点的坐 标。
一级影像严密数学模型是指通过卫星成像的原理，地面上的一点通过卫星 传感器在卫星一级影像上形成像面上一点这两点的位置关系，这一位置关系可 用公式"力="(^"'，^'"，^>表示，即为给定一个大地经讳度坐标就可以唯一确定其在卫星一级影像面上的像素坐标，实现本步骤的计算过程。本发明实施例 中，因为按现有技术所得卫星一级产品的RPC模型的模拟精度可以达到足够高
度，可以用一级产品的RPC模型取代一级影像严密数学模型。即本步骤的计算
过程转换为基于一级影像有理函数成像模型反变换公式，只需直接读入相应的
RPC参数，把参数直接带入RPC公式 "e"JP，Z，i/)中进行求解。
步骤1.2，调整像素坐标(^，乂，/7)中的高程A为t、e，通过一级影像严密数学模
型正变换公式或一级影像有理函数成像模型正变换公式，将像素坐标(^^A。J
变换为大地坐标(Dm,D嗣,U 。
因为按现有技术所得卫星一级产品的RPC模型的模拟精度可以达到足够高 度，本步骤也可以采用一级影像有理函数成像模型正变换公式代替一级影像严
密数学模型正变换公式， 一级影像有理函数成像模型正变换公式可表示为 (D/a'，A"") = 7Xx，;a; 。实施例中计算过程是读入RPC参数，然后对RPC方程基于
3尸
P， L进行线性化: 可得其误差方程
式中，
5 =
3丄
尸=^*AP + ~^ 'M 3丄
3尸3丄
3尸y气
7 = &_/, 『
乂 r> 乂
、AZ乂
W为权矩阵,根据以上误差方程进行迭代求解。计算过程中^、s为二级影像生成中用的统一高 程， 一般情况下为0。 F,°、《分别是指未知数X、 y取初始值时的K、 F/直，AP、
M分别是指P、 L的改正数。
步骤1.3，根据投影变换公式将大地坐标(A。"，A。w,U变化为投影面坐标
这个过程是一个投影变换的过程，根据生成二级产品时将一级影像按照卫 星的地面分辨率投影在地球椭球面的投影方式，即投影变换公式，把大地经纬 度坐标转换成投影面(即地球椭球面)下的投影坐标。二级影像的投影方式可
以在二级产品的影像信息中读出。投影变换公式可表示为=。
步骤1.4，根据二级影像的影像信息进行坐标转换，由投影面坐标^自，""一)计算 二级影像像素坐标(A,h)。
由于影像中是以一定的地面分辨率进行排列的，所以二级影像的像面坐标
与其投影坐标间也具有线性的关系，该关系可以表达为下面的形式
^p^4)"其中，"A《，A均可从二级影像的影像信息中获取。a、 b为 二级影像左上点的大地坐标，《,^为X、 Y方向的分辨率。由于生成二级影像
时，记录了"，6，《，A，因此可以直接利用这些几何参数进行坐标转换。根据该公 式即可通过投影坐标求出像面坐标。
步骤1.5，根据步骤l.l、 1.2、 1.3、 1.4所得结果构建物方均匀分布的控制点格 网，包括以下步骤，
步骤1.5.1根据一级影像的覆盖范围，求出二级影像的覆盖范围。由于一级影像的覆盖范围和二级影像的覆盖范围一致，所以，在实施例求 解过程中，根据一级影像的四个角点的影像坐标求出一级影像的覆盖范围，设 为二级影像的覆盖范围。
步骤1.5.2，根据二级影像的覆盖范围，以给定的经纬度间隔，建立平面格
网；具体实施时，可以一点的经纬度间隔，建立平面格网。格网数一般大于30*30 即可。
步骤1.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高 程，根据这些高程确定二级影像的覆盖范围内的最小和最大的高程，并根据给 定的高程分层数建立立体格网。
步骤1.5.4，根据步骤l.l、 1.2、 1.3、 1.4建立的流程得到二级严密几何模型
的反算公式(^,少2)=7;-'(￡>,。,,1)/。 ,/0，其中7;—1表示由大地经纬度坐标(^。,,"/。",/0反算
到二级影像面坐标(X2,力)的转换关系；根据该反算公式算出步骤5.3所得立体格 网中各格网点所对应的像点坐标，从而建立物方均匀分布的控制点格网；
根据步骤l.l、 1.2、 1.3、 1.4建立的流程得到地理坐标到像面坐标的关系，
从而可以建立(x2, h) = r2-1 ("to, a。"， /o 。
步骤1.6，根据步骤1.5建立的控制点格网求解二级产品的有理函数成像模 型参数。
RPC参数的解求一般有多种方法，本发明实施例提供的一种方法如下
—"e"s,(尸，丄，//)
根据RPC表达式_ d^( ,"巧列出误差方程「=勤"，『式中，
3尺 3FV apv 3尺
Sat,
3尺
3c, 3F,
,3F少
3",改.3c,. 3c/
(/ = 0，19;j=l,19)
v 乂
/ =
-C 一F0
x = (a, 6c, 《) (/ = 0,19;_/ = 1,19)
W为权矩阵。",~ c,《即为背景技术中所述的RPC系数。
利用最小二乘原理x = ^7'5)—^^由此可计算出RPC模型参数。
若实施时由于参数的相关性，使得求解过程中方程出现病态， 一般的解决 方式是岭参数估计法。
步骤1.7，根据步骤1.6所得二级产品的有理函数成像模型参数，生成二级产品 的有理函数成像模型。
由于步骤1.6获得了二级影像的RPC参数，即可参照一级影像有理函数成
像模型正变化公式和反变换公式，建立起二级影像的有理函数成像模型，可以 表示为(x,力=T2—1 (Dto， Z)/。", /z)禾口 (Dto， ￡>/。") = r2 (x,少,/z)。
(2)实施例二提供卫星影像三级产品的有理函数成像模型生成方法之一，参见 附图2:
步骤2.1，由某点的大地经纬度坐标(i^3,Z)^,/0通过一级影像严密数学模型反变 换公式或一级影像有理函数成像模型反变换公式，计算该点在一级影像上的像 素坐标Op"，/ )。
该步骤在实施例一中已经有详细的描述，在此不再重复叙述。由于各级产品的有理函数成像模型生成方法中大部分步骤完全相应，故在下面步骤描述中 相应部分将不再具体描述。
步骤2.2，根据一级影像上的控制点，(由像素坐标(^,A，/0通过仿射变换计算 该像点在纠正后对应的像面坐标(3"附/^，//朋，/0;调整高程A为V^，通过一级影 像严密数学模型正变换公式或一级影像有理函数成像模型正变换公式，将像素
坐标(s awp/e, /,'"e, / )变换为大地坐标(D'to3,化。 3， /^ve);
所谓控制点，是指大地坐标已知，并知道其在相应影像上的影像坐标的点。 因为三级产品的生成需要将一级影像经过控制点进行系统几何校正，因此本步 骤需要根据控制点对像素坐标"，兀,/2 )进行仿射变换。实施例二根据仿射变换
关系式P/。+X'^^+^fe ，利用最小二乘原理计算影像的仿射变换参数
乂，,，《，e。，e'，、 x、 y分别表示控制点在一级影像上的影像坐标，""""一、分 别表示控制点大地坐标经过基于RPC模型反算后在影像上的影像坐标。仿射变 换是平移、旋转、縮放等运算的组合，仿射变换关系式的系数即为仿射变换参 数。
步骤2.3，根据投影变换公式将大地坐标(D',。,3，D',— U变化为投影面坐标
这个过程是一个投影变换的过程，根据三级影像的投影方式，把大地经讳 度坐标转换成投影面下的投影坐标。三级影像的投影方式可以在影像信息中读 出。
步骤2.4，根据三级影像的影像信息进行坐标转换，由投影面坐标(A'。、'，"""—计算 三级影像像素坐标",A)。由于影像中是以一定的地面分辨率进行排列的，所以三级影像的像面坐标 与其投影坐标间也具有线性的关系，该关系可以表达为下面的形式
x = (Z) — a)/血
"(Z)^-^/漆其中，"》，《，A均可从三级影像的影像信息中获取。根据该公
式即可通过投影坐标求出像面坐标。
步骤2.5，根据步骤2.1、 2.2、 2.3、 2.4所得结果构建物方均匀分布的控制点格 网，包括以下步骤，
步骤2.5.1，根据一级影像的覆盖范围，求出三级影像的覆盖范围；
步骤2.5.2，根据三级影像的覆盖范围，以给定的经纬度间隔，建立平面格
网；
步骤2.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高 程，根据这些高程确定三级影像的覆盖范围内的最小和最大的高程，并根据给 定的高程分层数建立立体格网；
步骤2.5.4，根据步骤2.1、 2.2、 2.3、 2.4所得结果建立三级严密几何模型的 反算公式(；")=^(/),。,，/),。 ,/0 ，其中7T表示由大地经纬度坐标反算到三级影像
面坐标转换关系；根据该反算公式算出步骤2.5.3所得立体格网中各格网点所对 应的像点坐标，从而建立物方均匀分布的控制点格网；
步骤2.6，根据步骤2.5建立的控制点格网求解三级产品的有理函数成像模型参 数；
步骤2.7，根据步骤2.6所得三级产品的有理函数成像模型参数，生成三级产品 的有理函数成像模型，可以表示为(x,少hr3"(A。,，Z)/。",/2)和(Z),。,，A。"hr3(x,3a)。(3)实施例三提供卫星影像三级产品的有理函数成像模型生成方法之二，参见 附图3:
步骤3.1，由某点的大地经纬度坐标(A。,3,A—Z0通过二级影像严密数学模型反变 换公式或二级影像有理函数成像模型反变换公式，计算该点在二级影像上的像 素坐标^，^，/0;所述二级影像是二级产品提供的卫星产品影像，二级产品是 将卫星拍摄所得一级影像按照卫星的地面分辨率投影在地球椭球面上得到的几
何产品；
由于精确度足够的关系，可以采用二级影像有理函数成像模型反变换公式 代替二级影像严密数学模型反变换公式，二级影像有理函数成像模型反变换公 式可以由实施例一所得结果二级产品的有理函数成像模型获取。
步骤3.2，根据二级影像上控制点，由像素坐标";，乂，W通过仿射变换计算该像 点在纠正后对应的像面坐标""力力)；调整高程A为、,t,，通过二级影像严密数 学模型正变换公式，将像素坐标"，^，W变换为大地坐标(^'3，^""3，U。
根据仿射变换关系式&=义"^+&，利用最小二乘原理计算影像的仿射变
换参数/0，,，/^^^2， — ^分别表示控制点在一级影像上的影像坐标，~ ^分
别表示控制点大地坐标经过基于RPC模型反算后在影像上的影像坐标；
通过二级影像严密数学模型正变换公式，将^2，^，^^)变换为大地坐标
(2A'3，ZA"3，U ，同样也可以采用二级影像有理函数成像模型正变换公式代替。
可以由实施例一所得结果，即二级产品的有理函数成像模型，获取二级影像有 理函数成像模型正变换公式。
步骤3.4，根据三级影像的影像信息进行坐标转换，由投影面坐标("。。、'，""一)计算三级影像像素坐标(13,;;3);
步骤3.5，根据步骤3.1、 3.2、 3.3、 3.4所得结果构建物方均匀分布的控制点格 网，包括以下步骤，
步骤3.5.1，根据二级影像的覆盖范围，求出三级影像的覆盖范围；
步骤3.5.2，根据三级影像的覆盖范围，以给定的经纬度间隔，建立平面格
网；
步骤3.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高 程，根据这些高程确定三级影像的覆盖范围内的最小和最大的高程，并根据给 定的高程分层数建立立体格网；
步骤3.5.4，根据步骤3.1、 3.2、 3.3、 3.4所得结果建立三级严密几何模型的 反算公式(x^^7r'(仏,d,一)，其中7T'表示由大地经纬度坐标反算到三级影像 面坐标转换关系；根据该反算公式算出步骤3.5.3所得立体格网中各格网点所对 应的像点坐标，从而建立物方均匀分布的控制点格网；
步骤3.6，根据步骤3.5建立的控制点格网求解三级产品的有理函数成像模型参 数；
步骤3.7，根据步骤3.6所得三级产品的有理函数成像模型参数，生成三级产品
的有理函数成像模型，可以表示为(x,少)=r3-1 (Aw,/o和(d/。,, a。") = r3 (x, y， w 。 实施例三中若采用二级影像严密数学模型反变换公式和二级影像严密数学
模型正变换公式，可以直接用步骤l.l、 1.2、 1.3、 1.4建立的流程得到二级严密 几何模型的反算公式(^;;2) = 72-'(A。,，Z)h/0作为二级影像严密数学模型反变换公 式。二级影像严密数学模型正变换公式 通过如下流程得到步骤a，由二级影像上任意一点影像坐标^,72,/0根据二级影像的影像信息变换 成投影面坐标^，D""一 ；
歩骤b，根据投影变换公式，将投影面坐标^^'力"。一转化为经纬度坐标 Z);2)，设置高程为/^、c ;
步骤c，根据一级影像严密数学模型反变换公式或一级影像有理函数成像模型反
变换公式，将经纬度坐标(D',a,2，Z)',。n2A,)转换为该点对应的一级影像上的影像坐
标(Wi);
步骤d，通过一级影像上的影像坐标(x,，力)加上高程A用一级影像严密数学模型 正变换公式或一级影像有理函数成像模型正变换公式生成对应的大地坐标
通过步骤a、 b、 c、 d建立的流程得到二级严密几何模型的正算公式 (A。,，LW) = r2(x2,3;2,/z)，其中72表示由大地经纬度坐标正算到二级影像面坐标转换
关系。该正算公式即可作为二级影像严密数学模型正变换公式。
权利要求
1.一种卫星影像二级产品的有理函数成像模型生成方法，所述二级产品是将卫星拍摄所得一级影像按照卫星的地面分辨率投影在地球椭球面上得到的几何产品，投影关系采用投影变换公式体现；所述一级影像是对遥感卫星拍摄图像辐射纠正，但不进行几何处理得到的原始图像产品；二级产品提供的卫星产品影像为二级影像；其特征是包括以下步骤，步骤1.1，由某点的大地经纬度坐标(Dlat2，Dlon2，h)通过一级影像严密数学模型反变换公式或一级影像有理函数成像模型反变换公式，计算该点在一级影像上的像素坐标(x1，y1，h)；步骤1.2，调整像素坐标(x1，y1，h)中的高程h为hbase，通过一级影像严密数学模型正变换公式或一级影像有理函数成像模型正变换公式，将像素坐标(x1，y1，hbase)变换为大地坐标(Dlat1，Dlon1，hbase)；步骤1.3，根据投影变换公式将大地坐标(Dlat1，Dlon1，hbase)变化为投影面坐标(Deast，Dnorth)；步骤1.4，根据二级影像的影像信息进行坐标转换，由投影面坐标(Deast，Dnorth)计算二级影像像素坐标(x2，y2)；步骤1.5，根据步骤1.1、1.2、1.3、1.4所得结果构建物方均匀分布的控制点格网，包括以下步骤，步骤1.5.1，根据一级影像的覆盖范围，求出二级影像的覆盖范围；步骤1.5.2，根据二级影像的覆盖范围，以给定的经纬度间隔，建立平面格网；步骤1.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高程，根据这些高程确定二级影像的覆盖范围内的最小和最大的高程，并根据给定的高程分层数建立立体格网；步骤1.5.4，根据步骤1.1、1.2、1.3、1.4建立的流程得到二级严密几何模型的反算公式(x2，y2)＝T2-1(Dlat，Dlon，h)，其中T2-1表示由大地经纬度坐标(Dlat，Dlon，h)反算到二级影像面坐标(x2，y2)的转换关系；根据该反算公式算出步骤5.3所得立体格网中各格网点所对应的像点坐标，从而建立物方均匀分布的控制点格网；步骤1.6，根据步骤1.5建立的控制点格网求解二级产品的有理函数成像模型参数；步骤1.7，根据步骤1.6所得二级产品的有理函数成像模型参数，生成二级产品的有理函数成像模型。
2. —种卫星影像三级产品的有理函数成像模型生成方法，所述三级产品是将卫星拍摄所得一级影像经过控制点进行系统几何校正，并以按照卫星的地面分辨率投影在地球椭球面上得到的几何产品，投影关系采用投影变换公式体现；所述一级影像是对遥感卫星拍摄图像辐射纠正，但不进行几何处理得到的原始图像产品；三级产品提供的卫星产品影像为三级影像；其特征是包括以下步骤，步骤2.1，由某点的大地经纬度坐标(A。,3,A,a)通过一级影像严密数学模型反变换公式或一级影像有理函数成像模型反变换公式，计算该点在一级影像上的像素坐标(W"力)；步骤2.2，根据一级影像上的控制点，由像素坐标(Xp乂,/2)通过仿射变换计算该像点在纠正后对应的像面坐标(s"m/^，/,'朋,/0;调整高程/z为/z^，通过一级影像严密数学模型正变换公式或一级影像有理函数成像模型正变换公式，将像素坐标(5"附/^,//加,/2)变换为大地坐标(D'/加3,";3,U ;歩骤2.3，根据投影变换公式将大地坐标(D;3,Z)V3,U变化为投影面坐标步骤2.4，根据三级影像的影像信息进行坐标转换，由投影面坐标(A，""")计算三级影像像素坐标(、，^);步骤2.5，根据步骤2.K 2.2、 2.3、 2.4所得结果构建物方均匀分布的控制点格网，包括以下步骤，步骤2.5.1，根据一级影像的覆盖范围，求出三级影像的覆盖范围；步骤2.5.2，根据三级影像的覆盖范围，以给定的经纬度间隔，建立平面格网；步骤2.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高程，根据这些高程确定三级影像的覆盖范围内的最小和最大的高程，并根据给定的高程分层数建立立体格网；步骤2.5.4，根据步骤2.1、 2.2、 2.3、 2.4所得结果建立三级严密几何模型的反算公式&，少3) = 7;-'(ZVA")，其中7T表示由大地经纬度坐标反算到三级影像面坐标转换关系；根据该反算公式算出步骤2.5.3所得立体格网中各格网点所对应的像点坐标，从而建立物方均匀分布的控制点格网；步骤2.6，根据步骤2.5建立的控制点格网求解三级产品的有理函数成像模型参数；步骤2.7,根据步骤2.6所得三级产品的有理函数成像模型参数，生成三级产品的有理函数成像模型。
3. —种卫星影像三级产品的有理函数成像模型生成方法，所述三级产品是将卫星拍摄所得一级影像经过控制点进行系统几何校正，并以按照卫星的地面分辨率投影在地球椭球面上得到的几何产品，投影关系采用投影变换公式体现；所述一级影像是对遥感卫星拍摄图像辐射纠正，但不进行几何处理得到的原始图像产品；三级产品提供的卫星产品影像为三级影像；其特征是包括以下步骤，步骤3.1，由某点的大地经纬度坐标(A。,3,A—Z0通过二级影像严密数学模型反变换公式或二级影像有理函数成像模型反变换公式，计算该点在二级影像上的像素坐标^，^，W;所述二级影像是二级产品提供的卫星产品影像，二级产品是将卫星拍摄所得一级影像按照卫星的地面分辨率投影在地球椭球面上得到的几何产品；步骤3.2，根据二级影像上的控制点，由像素坐标";，^，W通过仿射变换计算该像点在纠正后对应的像面坐标""h，W;调整高程A为l^通过二级影像严密数学模型正变换公式或二级影像有理函数成像模型正变换公式，将像素坐标(A , A , W变换为大地坐标("'紹,"。"3， U ;步骤3.3，根据投影变换公式将大地坐标(D"。,3,D'^,U变化为投影面坐标步骤3.4，根据三级影像的影像信息进行坐标转换，由投影面坐标(0^力"。一计算三级影像像素坐标(&jg;步骤3.5，根据步骤3.1、 3.2、 3.3、 3.4所得结果构建物方均匀分布的控制点格网，包括以下步骤，步骤3.5.1，根据二级影像的覆盖范围，求出三级影像的覆盖范围；步骤3.5.2，根据三级影像的覆盖范围，以给定的经纬度间隔，建立平面格网；步骤3.5.3，根据平面格网中的格网点在全球DEM上内插出各格网点的高程，根据这些高程确定三级影像的覆盖范围内的最小和最大的高程，并根据给定的高程分层数建立立体格网；步骤3.5.4，根据步骤3.1、 3.2、 3.3、 3.4所得结果建立三级严密几何模型的反算公式(x3,h)-7T'(ZVDs/O ，其中7T'表示由大地经纬度坐标反算到三级影像面坐标转换关系；根据该反算公式算出步骤3.5.3所得立体格网中各格网点所对应的像点坐标，从而建立物方均匀分布的控制点格网；步骤3.6，根据步骤3.5建立的控制点格网求解三级产品的有理函数成像模型参数； 步骤3.7，根据步骤3.6所得三级产品的有理函数成像模型参数，生成三级产品的有理函数成像模型。
全文摘要
本发明卫星影像二/三级产品的有理函数成像模型生成方法利用二级几何模型的反变换建立控制点格网，建立二级几何产品的RPC模型；利用有一级生成的三级产品的反变换建立控制点格网，建立三级几何产品的RPC模型；利用二级生成的三级产品的反变换建立控制点格网，建立三级几何产品的有理函数成像模型。本发明使二级、三级产品在进一步的应用中更具优势。
文档编号G01S7/48GK101598786SQ200910062910
公开日2009年12月9日 申请日期2009年6月30日 优先权日2009年6月30日
发明者斌 刘, 过 张, 江万寿 申请人:武汉大学