基于RPC模型的线阵卫星影像核线排列及其几何关系重建
基于RPC模型的天绘一号卫星影像高精度立体定位
—
代入
（ １ ）中， 整理得 ：
Ｓ Ｃ Ａ Ｌ Ｅ・ Ｆ （ Ｐ， Ｌ 。 Ｈ 、 ＋Ｌ Ｉ Ｎ Ｅ ０ Ｆ Ｆ
ｆ ４１
—
＼ － Ｃ＝ ｓ Ａ ＭＰ ｓ Ｃ Ａ Ｌ Ｅ・ Ｇ （ Ｐ。 Ｌ。 Ｈ 、 ）＋ｓ Ａ ＭＰ ｏ Ｆ Ｆ
将上 式 的两个 方程 按 照 Ｔ ａ ｙ ｌ ｏ ｒ 公 式 展开 至一
Ｊ ｉ ａ ｎ Ｊ ｉ ａ ｎ ｆ ｅ ｎ ｇ 一，Ｌ ｉ Ｘｉ ｎ ｔ ａ ｏ 一，Ｈｕ Ｙａ ｎ ’
１． Ｘｉ ’ ａ ｎ Ｒｅ ｓ ｅ ａ ｒ ｃ ｈ Ｉ ｎ ｓ ｔ ｉ ｔ ｕ ｔ ｅ ｏ ｆ Ｓｕ ｒ ｖ ｅ ｙ ｉ ｎ ｇ ａ ｎｄ Ｍａ ｐｐ ｉ ｎｇ，Ｘｉ ’ ａ ｎ ７１ ０ ０５ ４， Ｃｈ ｉ ｎ ａ ２． Ｓ ｔ ａ ｔ ｅ Ｋｅ ｙ Ｌａ ｂｏ ｒ ａ ｔ ｏ ｒ ｙ ｏ ｆ Ｇｅ ｏ — Ｉ ｎ ｆ ｏ ｒ ｍａ ｔ ｉ ｏ ｎ Ｅｎ ｇ ｉ ｎｅ ｅ ｒ ｉ ｎ ｇ，Ｘｉ ’ ａ ｎ ７ １ ００ ５ ４ ，Ｃｈ ｉ ｎ ａ
ｄ，
～
ＨＥＩ ＧＨＴ ＳＣ ＡＬ Ｅ
１ Ｊｒ ｌ
Ｏ Ｌ ｏ ｎ
Ｏ ｃ ，
Ｏ Ｌａ ｔ
Ｏ ｃ ，
Ｏ Ｈｇｔ
Ｏ ｃ ７
ｃ—ｓ ＡＭ Ｐ
—
ｏＦＦ
ｓ ＡＭ Ｐ ｓＣＡＬＥ ｒ — ＬＩ ＮＥ ｏＦＦ
（ ３ ）
Ｖ＝
ｃ ｌ
，
Ｏ Ｌ ｏ ｎ
Ｏ Ｌａ ｔ
） 、 （ Ｃ ， ， ｒ ｒ ） ， 由式 （ ６ ）可得 如下误 差方 程 ：
Ｌｏ ｎ
： —
－
ＬＯＮＧ
—
一种基于 RPC 模型的 SAR 影像几何精纠正新方法
一种基于 RPC 模型的 SAR 影像几何精纠正新方法刘美琳【摘要】Synthetic Aperture Radar (SAR ) images have advantages of being usefulness all-weather and all-phase. Therefore,they have been widely used in terrain mapping and change detection.Rational Polynomial Coefficient (RPC)model is a novel method being used in SAR image geo-location,which is a highly accurate general model comparing with the physical model and can be used to take place of the physical model.This paper utilizes GCPs (Ground Controlling Points)to achieve RPC model refining and puts forward a novel method of SAR geocorrection based on the refined RPC Model.The test site is in Genhe district in Inner Mongolia province in China and the test images are TerraSAR-X and Radarsat-2 images.This paper concludes that:(1)RPC Model has a highly accuracy compared with R-D model,which is all most 10e-6 pixel.(2)Based on the refined RPC model,the geo-location accuracy of TerraSAR-X image is 0.01 pixel and the one of Rasarsat-2 image is 0.1 pixel. (3)The quality of orthoimages of geocorrection can be several meters based on refined RPC model.%RPC（Rational Polynomial Coefficient）是最新提出的应用于 SAR （Synthetic Aperture Radar）影像定位的模型算法，该模型通过对严格物理模型的高精度拟合，达到替代严格物理模型的目的。
基于RPC模型的线阵卫星影像核线排列及其几何关系重建
第4期，总第87期国土资源遥感No．4，20102010年12月15日REMOTE SENSING FOR LAND ＆RESOURCESDec．，2010基于RPC 模型的线阵卫星影像核线排列及其几何关系重建张过1，潘红播1，江万寿1，秦绪文2（1．武汉大学测绘遥感信息工程国家重点实验室，武汉430079；2．中国地质调查局，北京100037）摘要：针对核线影像无几何模型的问题，提出基于有理多项式系数（Rational Polynomial Coefficient ，RPC ）模型，利用投影轨迹法制作线阵推扫式卫星核线影像及其几何模型的重建方法。
利用SPOT 5HRG 、CBERS 2－03异轨立体像对和SPOT 5HRS 、IRS P5、GeoEye 、IKONOS 同轨立体像对重建的核线影像几何模型上下视差的中误差在0．2个像元之内；利用重建几何模型和原始影像几何模型进行前方交会计算所得到的地面点平面中误差在1m 之内，高程中误差在0．3m 之内，验证了本文方法的可靠性。
关键词：投影轨迹法；核线几何模型；RPC 模型；线阵列影像中图分类号：TP 75文献标识码：A文章编号：1001－070X （2010）04－0001－05收稿日期：2010－03－12；修订日期：2010－05－06基金项目：国家科技支撑项目（编号：2011BAB01B01）、国家自然科学基金项目（编号：40801178）、国家863项目（编号：2007AA120203、2009AA12Z145）和对地观测技术国家测绘局重点实验室经费资助项目（编号：2009－17）。
0引言利用星载线阵CCD 推扫式传感器获取遥感影像，实现立体量测是遥感应用的一个重要方面。
从立体影像中提取三维信息，一个重要的约束条件就是核线约束，在匹配检索同名点时，只用沿着核线方向进行一维搜索即可，这不仅可以减少匹配时间，而且可以提高匹配的可靠性。
但线阵CCD 推扫式影像多投影中心的成像几何特点，导致适用于单中心投影的框幅式核线影像的制作方法不适用于星载线阵列推扫式影像。
RPC模型在星载SAR与星载InSAR几何处理中的研究与应用
RPC 模型在星载SAR 与星载InSAR 几何处理中的研究与应用摘要：在卫星遥感领域中，星载合成孔径雷达（SAR）和星载干涉合成孔径雷达（InSAR）是两种被广泛应用的雷达技术。
这两种技术需要漂移校正、几何校正、相位解缠等步骤的处理。
本文针对这些处理步骤中的几何处理阶段，介绍了RPC 几何模型及其在星载SAR 和星载InSAR 中的研究和应用。
关键词：RPC 模型；星载SAR；星载InSAR；几何处理引言：合成孔径雷达（SAR）和干涉合成孔径雷达（InSAR）是卫星遥感技术中常用的雷达技术。
这两种技术能够获取高分辨率的地表图像和高精度的DEM 信息。
在SAR 和InSAR 数据处理过程中，需要进行漂移校正、几何校正、相位解缠等步骤。
其中，几何处理是前置步骤，是确保后续处理流程正确进行的必要步骤。
遥感图像的几何处理，即是将遥感数据中的像元坐标系转换为地理坐标系的过程，也就是将遥感图像像素的位置与实际地物位置相互关联。
RPC 模型是一种常见的几何变换模型，在星载SAR 和星载InSAR 几何处理中有着广泛的应用。
本文将介绍RPC 模型的原理、应用和研究进展。
一、RPC 几何模型的原理RPC（Rational Polynomial Coefficients）模型是一种基于多项式函数描述的几何变换模型，能够将像元坐标系和地理坐标系之间的映射关系进行描述。
RPC 模型基于多项式函数描述，通过计算多项式系数进行像素坐标到地理坐标的转换。
由于其功能强大，操作简单，因此RPC 模型被广泛应用于卫星遥感数据的几何处理中。
c^2 RPC 模型的本质是一个二维的多项式函数，其数学表达式如下：X = P_1 + P_2 * c + P_3 * r + P_4 * r * c + P_5 * r^2 + P_6 * Y = P_7 + P_8 * c + P_9 * r + P_10 * r * c + P_11 * r^2 + P_12* c^2式中，r 和c 分别为像元的行列号，(X,Y)为像素坐标对应地理坐标系下的坐标，P_i 为多项式系数。
重构RPC在提高卫星影像精度的应用
多项式描述为：
{ｘ＝ＮｕｍＬ（Ｘ，Ｙ，Ｚ）／ＤｅｎＬ（Ｘ，Ｙ，Ｚ） ｙ＝ＮｕｍＳ（Ｘ，Ｙ，Ｚ）／ＤｅｎＳ（Ｘ，Ｙ，Ｚ）
（１）
式中：（ｘ，ｙ）为正则化的原始卫星图像坐标；（Ｘ，Ｙ，Ｚ）为
正则化的地面点大地测量坐标；ＮｕｍＬ 和 ＮｕｍＳ 分别为求 解正则化图像横纵坐标 ｘ和 ｙ所用的分子多项式；ＤｅｎＬ 和 ＤｅｎＳ分别为求解正则化图像横纵坐标 ｘ和 ｙ所用的 分母多项式。
本文使用的 ＰＩＥｏｒｔｈｏｒ的重构 ＲＰＣ功能，通过控制点 平差，对存在问题的卫星影像 ＲＰＣ进行重新解算优化， 以解决同景卫星影像数据因 ＲＰＣ问题导致的精度达不 到要求的情况。
１ 原理概述
１．１ 几何校正原理 几何校正一般利用一系列数学模型来降低遥感影像
成像时产生 的 几 何 畸 变，获 得 更 高 精 度 的 遥 感 影 像。 基 于有理函数模型的几何精校正的重点是利用严密的几何 模型，通过计算分析，在影像数据中提取数量充足、均匀 分布的控制点，获得足够的参数用于定向分析，提高模拟 精度，而且应用更灵活。 １．２ 有理函数模型
高分辨率卫星影像通常使用严格轨道模型完成影像 的高精度几何定位，但卫星制造者基于保密原因，会提供 一种包含有理多 项 式 系 数 的 文 件，作 为 卫 星 影 像 的 几 何 校正处理的参数，代替以共线条件为严格成像几何模型， 并将这种使用 ＲＰＣ文件系数的模型称为有理函数模型。 有理函数模型相对传感器和卫星平台，没有直接关联，不 涉及传感器成像的物理意义，而是通过有理多项式，描述 像点和相应地面 点 之 间 的 几 何 关 系，具 有 优 良 的 拟 合 特 性，但影像自带的 ＲＰＣ参数有较大的误差，需要建立误 差补偿模型进行补偿。
上述试验结果证明，当遇到因为 原始 ＲＰＣ，导 致 定
一种基于RPC模型的高分卫星影像侧视角计算方法
因和影响的程度,从 原 理 上 论 证 侧 视 角 与 高 程 精 度
袁修孝
等人基 于 CCD 线 阵 列 侧 视 角 匀 速 变
[
４]
收稿日期:
２０１９
Ｇ
０３
Ｇ
０６
第一作者简介:尚大帅(
１９８８－ ),男,工程师,硕士 ．
的关系.这些研究 多 是 从 理 论 出 发,需 要 卫 星 平 台
参数、相机 参 数 以 及 成 像 条 件 参 数 等 的 支 持,计 算
r
dbyt
hep
i
xe
lc
oo
rd
i
na
t
eo
ft
hec
en
t
r
eo
fimageandt
hespe
c
i
f
i
ed
yi
e
l
eva
t
i
onva
l
ue,wh
i
chi
ss
imp
l
ei
np
r
i
nc
i
l
e,sma
l
li
nc
a
l
cu
l
a
t
i
on,e
a
s
oimp
l
emen
tandha
sadvan
t
age
si
n
p
yt
eng
LIZh
i
i
an
f
eng,
YIN Ya
uan
yong,WANGJ
j
(
Xi
′
anS
t
基于有理函数模型的线阵推扫式卫星近似核线影像的生成
论
成方法。该方法 以有理函数模 型描述推扫式传感器的成像 几何 关 系， 并采 用投 影轨迹 法建立线阵推扫式影像 的扩展
核 线模型 。试验表明该方法能生成子像素级的近似核线影像 ， 算法简单可靠。
关键词 ： 有理 函数模型 ； 核线； 投影轨迹法 ； 线阵推扫式成像
ｌ 引 言
Ｉ Ｋ Ｏ Ｎ Ｏ Ｓ 、 Ｑ ｕ ｉ ｃ ｋ Ｂ ｉ ｒ ｄ 、 Ａ Ｌ Ｏ Ｓ等高分 辨率 线阵推 扫式 卫 星的 出现 ， 使 得利 用 卫 星 立体 影像 实现 地 面 目标 的 高精 度定 位和 大 比例 尺 测 图成 为 可 能 。然 而 ， 线 阵推
２ ０ １ ３年 １ ０月 第５ 期 来自城市勘
测
０ｃ ｔ ． ２０１ ３ Ｎｏ ． ５
Ｕｒ ｂ ａ ｎ Ｇｅ ｏ ｔ ｅ ｃ ｈｎ ｉ ｃ ａ ｌ Ｉ ｎ ｖ ｅ ｓ ｔ ｉ ｇ ａ ｔ ｉ ｏ ｎ＆ Ｓ ｕｒ ｖ ｅ ｙ ｉ ｎ ｇ
文章编号 ： １ ６ ７ ２ — ８ ２ ６ ２ （ ２ ０ １ ３ ） ０ ５ — ６ Ｏ — ｏ ３
发， 经过左像 的投影 中心 Ｓ （ Ｘ ， ， Ｚ ） 成像 于左 像上 Ｐ 点， 则该成像光线上 每个点经过 Ｓ（ Ｘ ， ， ） 投影 到右
像上 ， 得到 的投 影 点轨 迹将 形成 一 条 曲线 ， 将这 条 曲线 定义 为像 点 Ｐ的核线 。类 似 于框 幅式 成像 的核线几 何 关系， 点 Ｐ在右像上 的同名点 Ｐ必定位于这条 “ 核 曲线 ” 上， 这就是线 阵推扫式影像核线 的几何定义 。
论 最 为严密 。
像 的物理 几何过 程是非 常重要 的 。随着高 分辨 率传 感 器 技术 的发 展 ， 高 分 辨 率商 业 卫 星 如 Ｉ Ｋ Ｏ Ｎ Ｏ Ｓ 、 Ｑ ｕ ｉ ｃ ｋ —
基于通用RPC文件的卫星影像正射纠正应用
基于通用RPC文件的卫星影像正射纠正应用摘要：利用航空或航天摄影像片制作正射影像图，是当前地理信息系统建设和遥感应用中非常重要的内容。
本文以Alos卫星影像为例，介绍在envi软件中生成通用RPC和正射影像制作过程，本方法通过实际生产，验证了其方法生产的正射影像精度可靠，解决了无RPC文件时对山区影像很好消除地形高差引起的变形问题。
关键词：Abstract: This paper take Alos satellite images for example, introduced in the envi software to generate a general-purpose RPC and orthophoto production process, and by actual production, verify its method of production of orthophoto accuracy to solve the RPC file on the mountain image very well eliminate the problem of deformation caused by the terrain elevation.中图分类号: P185.18文献标识码：A 文章编号：一．DOM是指消除各种影像变形后与地面物体形状相似，成一定比例并带有坐标的数字影像地图，也称数字正射影像图(digital orthor map)。
在测绘生产中，4D产品DLG、DRG、DEM、DOM是目前的主要内容，正射影像图DOM 做为其中一种新产品，以其直观、信息丰富等特点，成为测绘产中一重要部分。
在DOM产品指标中，像点精度、地面分辨率、色彩、反差是主要的内容，而DOM影像点坐标精确度，是一项很重要指标，如果不同像片的DOM坐标不一致，将无法做为影像图使用，也无法完成相互之间的拼接，也就无法制作合格的DOM产品。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第4期，总第87期国土资源遥感No．4，20102010年12月15日REMOTE SENSING FOR LAND ＆RESOURCESDec．，2010基于RPC 模型的线阵卫星影像核线排列及其几何关系重建张过1，潘红播1，江万寿1，秦绪文2（1．武汉大学测绘遥感信息工程国家重点实验室，武汉430079；2．中国地质调查局，北京100037）摘要：针对核线影像无几何模型的问题，提出基于有理多项式系数（Rational Polynomial Coefficient ，RPC ）模型，利用投影轨迹法制作线阵推扫式卫星核线影像及其几何模型的重建方法。
利用SPOT 5HRG 、CBERS 2－03异轨立体像对和SPOT 5HRS 、IRS P5、GeoEye 、IKONOS 同轨立体像对重建的核线影像几何模型上下视差的中误差在0．2个像元之内；利用重建几何模型和原始影像几何模型进行前方交会计算所得到的地面点平面中误差在1m 之内，高程中误差在0．3m 之内，验证了本文方法的可靠性。
关键词：投影轨迹法；核线几何模型；RPC 模型；线阵列影像中图分类号：TP 75文献标识码：A文章编号：1001－070X （2010）04－0001－05收稿日期：2010－03－12；修订日期：2010－05－06基金项目：国家科技支撑项目（编号：2011BAB01B01）、国家自然科学基金项目（编号：40801178）、国家863项目（编号：2007AA120203、2009AA12Z145）和对地观测技术国家测绘局重点实验室经费资助项目（编号：2009－17）。
0引言利用星载线阵CCD 推扫式传感器获取遥感影像，实现立体量测是遥感应用的一个重要方面。
从立体影像中提取三维信息，一个重要的约束条件就是核线约束，在匹配检索同名点时，只用沿着核线方向进行一维搜索即可，这不仅可以减少匹配时间，而且可以提高匹配的可靠性。
但线阵CCD 推扫式影像多投影中心的成像几何特点，导致适用于单中心投影的框幅式核线影像的制作方法不适用于星载线阵列推扫式影像。
对于不同的传感器模型，其对应的核线模型是不同的。
张祖勋［1］基于Dowman 提出将左扫描线作为左核线的思路，认为右核线近似为多项式，并利用若干同名点坐标，计算出对应的核线方程，但该方法仅适用于异轨立体，对同轨立体不再适用［2］。
在基于外方位元素呈线性变化的假设下，大量研究表明核线方程为双曲线［3－5］，但并没有给出具体的核线重采样方式。
Kim ［6］基于Orun 和Natarjan 的共线方程模型得到了核线方程为双曲线，但该方程无法应用于核线影像重采样。
当假设线推扫式传感器姿态和速度不变时，利用平行投影法获取的核线对为直线，但该方法需要地面控制点以计算平行投影法的系数［7，8］。
由于RPC 模型是一种具有精度高、速度快等特点的通用传感器模型，被广泛应用于线阵列光学卫星中［9］。
叶新魁等［10］利用RPC 模型提出基于零高程点的星载推扫式立体影像核线确定方法，获得了卫星影像投影轨迹核线的双曲线，但方程复杂，无法进行核线重采样。
张永军等［11］利用投影轨迹法，获取核线方程后进行分段拟合获取核线方程，但该方法重采样时需要DEM 辅助，且在进行核线影像反算时需要核线影像上的同名点坐标及在原始影像上相应点的坐标。
胡芬等［12］利用投影轨迹法获取核线在基准面上的方向，沿着核线方向进行投影以获取近似核线。
以上方法主要存在以下几点问题：①获取的核线方程无法应用于核线重采样，或者需要DEM 辅助，而制作核线影像的主要目的是为了方便将二维匹配简化为一维匹配，进而方便获取DEM ，这是与制作核线影像的目的相矛盾的；②核线影像正变换［9］过程复杂；③以上所有方法检查核线影像的精度均采用0．5像素精度的人工选点方式，其方法不精确。
针对上述问题，本文采用投影轨迹法获取核线，用直线拟合方法计算核线方程，建立了原始影像与核线影像的对应关系。
在此基础上，进行核线重采样，同时利用RPC 模型重建核线影像的几何模型，且提出利用核线模型是否存在上下视差作为检查核线模型精度指标的方法，并利用SPOT 5HRG 、CBERS 2－03异轨立体数据，SPOT 5HRS 、GeoEye 、IKONOS 、IRS P5同轨立体数据进行试验，验证本文所提出方法的可行性。
国土资源遥感2010年1投影轨迹法核线Kim ［6］提出的投影轨迹法是将高程差异引起的点位轨迹定义为核线，左影像上一点q 以其几何模型投影到不同高程面上，得到对应的地面点，记为Q 1、Q 2、Q 3，并将这些地面点按右影像几何模型投影到右影像上，记为q 1、q 2、q 3，获取了左影像上点q 所对应的右核线，如图1所示，若是核线对存在，则右核线上的点所对应的左核线应当相同。
利用右影像几何模型反解法［9］将右影像上q 1、q 2、q 3点投影到不同高程面，利用左影像几何模型将这些地面点投影到左影像上，即可获取同名核线对。
用直线拟合核线上的离散点，获得核线方程。
图1投影轨迹法示意图Fig．1Sketch for projection track method按框幅式影像的核线制作流程，在生成核线影像之前需要进行相对定向［13］。
对于线阵列传感器影像来说，可以利用基于RPC 模型的像面仿射变换进行无控制点定向处理，以保证同名光线对对相交［9］。
由核线定义可知，核线方向与摄影基线的方向相近似。
线阵列影像通常有异轨立体（Across －track ，采用侧视的方式）和同轨立体（Along －track ，采用前后视方式）两种立体观测，不同成像方式的卫星影像，其核线在影像上的分布状态不同，前者的核线方向与扫描方向大致相同，后者的核线方向则近似垂直于扫描方向。
为获取核线影像，需要计算整幅影像的核线方程。
而核线方程受q 点位置影响，故使q 沿近似垂直于核线方向前进：对于异轨立体q ，应当沿着垂直于扫描方向前进；对于同轨立体，则应当沿着扫描方向前进。
为使拟合的核线尽可能准确，q 点应当位于影像中间。
当核线对方程与影像区域没有交点的时候，核线方程可以舍去。
这样就获得了所有核线对的方程，即a i x +b i y +c i =0a'i x +b'i y +c'i {=0（i =1，…，n ）（1）式中，a 、b 、c 为直线方程的系数；x 、y 为核线上点的坐标；i 为核线方程的编号；n 为核线方程的数目。
在重采样时，本文采用了文献［13］“直接在倾斜像片上获取影像”的方法逐行进行，再沿着q 前进方向进行线性内插，获得整个区域的核线影像。
2核线影像几何模型重建核线影像上的点与原始影像上的点是一一对应的，而原始影像的几何模型建立了原始影像上点与地面点的对应关系。
根据以上关系可建立核线影像上点与地面点的对应关系，重建核线影像的几何模型。
其中，建立核线影像与原始影像对应关系是确定原始影像点所在核线方程的序号i 。
当核线影像上一点变换到原始影像时，i 可以从其坐标值获得，但是对于现有的摄影测量系统来说，无法直接利用核线方程进行核线影像前方交会。
本文将利用RPC 模型作为通用传感器模型的特点，重建核线影像的几何模型，使核线影像可以直接应用于三维信息提取。
2．1像素点与对应地面点解析关系的确定由于核线影像上每一行对应于所求解的整幅影像的核线方程，因此对于核线影像上一点M 0（x 0，y 0）的y 值即可确定核线的序号i ，但核线的序号为整数，这就需要对M 0的y 值进行取整运算，如图2所示。
图2异轨核线影像与原始影像上同名点对应关系（左：核线影像；右：原始影像）Fig．2Corresponding relation of tie －points inacross －track stereo images and original images对于异轨立体数据，核线方程是近似于扫描方向的，那么M 0（x 0，y 0）的x 0值与原始影像上对应点M （x ，y ）的x 值相等［13］。
利用核线方程及x 0即可确定y 的近似值y'。
M （x ，y ）的坐标为x =x 0y =y'+y 0－int （y 0{）（2）而对于同轨立体数据来说，核线方程是近似垂·2·第4期张过，等：基于RPC 模型的线阵卫星影像核线排列及其几何关系重建直于扫描方向的，那么，x 0与原始影像上对应点M （x ，y ）的y 值相等，因此利用M 0的x 值即可确定M 的x 值。
利用式（2）即建立了核线影像与原始影像的对应关系。
由于原始影像的几何模型参数是已知的，即可建立核线影像上一点与地面点的对应关系，即核线影像的正变换过程，其流程如图3所示。
图3核线影像的正变换Fig．3The forward transformation of epipolar images首先，根据公式（2）将核线影像上的像素坐标（x 0，y 0）变换为原始影像上对应的点（x ，y ）；然后，根据原始影像的几何模型，将原始影像上点（x ，y ）投影到某高程面上，获取其对应地面点坐标（lat ，lon ，h ）。
2．2RPC 参数的求解采用与地形无关的求解方式求解RPC 模型参数［9］，具体步骤如下：（1）根据核线影像的覆盖范围，利用美国地质调查局提供的全球1km 分辨率DEM （Global 30－Arc －second Digital Elevation Model ），计算该区域的最大、最小椭球高。
（2）在核线影像平面上，以一定的间隔建立规则格网。
利用核线方程，计算出这些格网点所对应原始影像上的坐标，再利用原始影像的RPC 模型将这些点投影到以一定间隔分层的高程面上，这样即建立核线影像上点与地面点的对应关系。
同时，加密控制格网和高程面，建立独立的检查点。
（3）利用控制点坐标计算核线影像与地面坐标的归一化参数，将控制点和检查点归一化。
（4）利用岭估计法求解RPC 模型参数［9］，获取核线影像的RPC 模型参数，并用检查点进行精度检查。
3实验结果与分析采用SPOT 5HRS 、GeoEye 、IRS P5、IKONOS 同轨数据，SPOT 5HRG 、CBERS 2－03异轨数据进行试验，各个影像的参数如表1所示。
表1实验数据的详细说明Tab．1Detailed specification of datasets in experiments数据源立体覆盖区域中心左影像大小/像素ˑ像素右影像大小/像素ˑ像素最小高程/m 最大高程/m 立体模式影像分辨率/m 定向［8］精度/像素地貌SPOT 5HRS 107ʎ27'19．758ᵡE 34ʎ24'3．683ᵡN 12000ˑ1200012000ˑ120003842452同轨100．160山地IRS P 5116ʎ31'44．802ᵡE 40ʎ24'7．869ᵡN 12000ˑ1200012000ˑ12000－141315同轨2．50．182山地IKONOS 147ʎ16'30．720ᵡE 40ʎ51'52．560ᵡS 12122ˑ1314812122ˑ13148－301284同轨10．154丘陵GeoEye 147ʎ15'31．500ᵡE 42ʎ51'38．520ᵡS 26928ˑ3166826960ˑ31668－301284同轨0．50．146丘陵SPOT 5HRG 116ʎ23'23．615ᵡE 39ʎ55'28．692ᵡN 12002ˑ1200212000ˑ12002－27234异轨50．115平原CBERS 2－03116ʎ27'49．844ᵡE 39ʎ58'1．159ᵡN10002ˑ1000010002ˑ10000－28111异轨30．235平原其中，SPOT 5HRS 为宝鸡地区影像，CBERS 2－03和SPOT 5HRG 为北京市区影像，GeoEye 和IKONOS 为澳大利亚Hobart 地区影像，P 5为北京市怀柔地区影像。