第 4 期，总第 2010 年 12 月 期 日 87 15 国 土 资 源 遥 感 REM OTE SENSING FOR LAND ＆ RESOURCES ， 2010 ， 2010 No． 4 Dec． 基于 RPC 模型的线阵卫星影像核线排列 及其几何关系重建 张 过1，潘红播1，江万寿1，秦绪文2 ( 1． 武汉大学测绘遥感信息工程国家重点实验室，武汉 430079 ; 2． 中国地质调查局，北京 100037 ) SPOT 5 HRS、IRS P5、GeoEye、IKONOS 摘要: 针对核线影像无几何模型的问题，提出基于有理多项式系数( 投影轨迹法制作线阵推扫式卫星核线影像及其几何模型的重建方法 像对和 像元之内; 利用重建几何模型和原始影像几何模型进行前方交会计算所得到的地面点平面中误差在 程中误差在 关键词: 投影轨迹法; 核线几何模型; 中图分类号: TP 75 文献标识码: A 模型; 线阵列影像 文章编号: 1001 － 070X 之内，验证了本文方法的可靠性 Rational Polynomial Coefficient 04 － 0001 － 05 0． 3 m 利用 2010 RPC 。 。 ( ) 同轨立体像对重建的核线影像几何模型上下视差的中误差在 SPOT 5 HRG、CBERS 2 － 03 ， RPC 1 m ) 模型，利用 异轨立体 个 之内，高 0． 2 0 引言 。 。 CCD 利用星载线阵 推扫式传感器获取遥感影 像，实现立体量测是遥感应用的一个重要方面 从 立体影像中提取三维信息，一个重要的约束条件就 是核线约束，在匹配检索同名点时，只用沿着核线方 向进行一维搜索即可，这不仅可以减少匹配时间，而 推扫式影 且可以提高匹配的可靠性 像多投影中心的成像几何特点，导致适用于单中心 投影的框幅式核线影像的制作方法不适用于星载线 阵列推扫式影像 但线阵 CCD 。 1 ］ 2 。 ］基于 Dowman 张祖勋［ 对于不同的传感器模型，其对应的核线模型是 不同的 提出将左扫描线作 为左核线的思路，认为右核线近似为多项式，并利用 若干同名点坐标，计算出对应的核线方程，但该方法 仅适用于异轨立体，对同轨立体不再适用［ 在基 于外方位元素呈线性变化的假设下，大量研究表明 ］，但并没有给出具体的核线 核线方程为双曲线［ 重采样方式 的共线方 ［ 6 程模型得到了核线方程为双曲线，但该方程无法应 用于核线影像重采样 当假设线推扫式传感器姿态 和速度不变时，利用平行投影法获取的核线对为直 线，但该方法需要地面控制点以计算平行投影法的 系数［ 速度 、 快等特点的通用传感器模型，被广泛应用于线阵列 模型是一种具有精度高 Natarjan ］基于 。Kim 由于 Orun RPC 和 ， ］ 8 。 。 。 3 － 5 7 11 10 ］ 9 。 。 RPC DEM ］利用 叶新魁等［ 张永军等［ 光学卫星中［ 模型提出基 于零高程点的星载推扫式立体影像核线确定方法， 获得了卫星影像投影轨迹核线的双曲线，但方程复 杂，无法进行核线重采样 ］利用投影轨 迹法，获取核线方程后进行分段拟合获取核线方程， 辅助，且在进行核线影 但该方法重采样时需要 像反算时需要核线影像上的同名点坐标及在原始影 像上相应点的坐标 ］利用投影轨迹法获 取核线在基准面上的方向，沿着核线方向进行投影 以获取近似核线 以上方法主要存在以下几点问 获取的核线方程无法应用于核线重采样，或 题: ① 辅助，而制作核线影像的主要目的是为 者需要 了方便将二维匹配简化为一维匹配，进而方便获取 核 以上所有方法检查核 ③ 像素精度的人工选点方 ，这是与制作核线影像的目的相矛盾的; DEM 线影像正变换［ 线影像的精度均采用 式，其方法不精确 ］过程复杂; 胡芬等［ DEM 0． 5 。 。 ② 12 9 针对上述问题，本文采用投影轨迹法获取核线， 用直线拟合方法计算核线方程，建立了原始影像与 在此基础上，进行核线重采 核线影像的对应关系 样，同时利用 模型重建核线影像的几何模型， 且提出利用核线模型是否存在上下视差作为检查核 线模 型 精 度 指 标 的 方 法，并 利 用 RPC 。 SPOT 5 HRG、 。 异轨立体数据， SPOT 5 HRS、GeoEye、 同轨立体数据进行试验，验证本文 CBERS 2 － 03 IKONOS、IRS P5 所提出方法的可行性 。 收稿日期: 2010 － 03 － 12 基金项目: 国 家 科 技 支 撑 项 目 ( 编 号: ; 修订日期: 2010 － 05 － 06 2007AA120203、2009AA12Z145 2011BAB01B01 ) 和对地观测技术国家测绘局重点实验室经费资助项目( 编号: 40801178 、 ) 国 家 自 然 科 学基 金 项 目 ( 编 号: ) 、 国 家 ) 2009 － 17 项 目 ( 编 号: 863 。 

·2· 国 土 资 源 遥 感 年 2010 1 投影轨迹法核线 q ［ 6 Kim ］提出的投影轨迹法是将高程差异引起的 点位轨迹定义为核线，左影像上一点 以其几何模 型投影到不同高程面上，得到对应的地面点，记为 ，并将这些地面点按右影像几何模型投影 Q1 、Q2 、Q3 到右影像上，记为 所 所示，若是核线对存在，则右 对应的右核线，如图 核线上的点所对应的左核线应当相同 利用右影像 几何模型反解法［ 点投影到 不同高程面，利用左影像几何模型将这些地面点投 影到左影像上，即可获取同名核线对 用直线拟合 核线上的离散点，获得核线方程 ，获取了左影像上点 ］将右影像上 q1 、q2 、q3 q1 、q2 、q3 。 。 1 q 9 。 图 1 投影轨迹法示意图 Fig． 1 Sketch for projection track method 13 ］ 按框幅式影像的核线制作流程，在生成核线影 像之前需要进行相对定向［ 对于线阵列传感器 。 影像来说，可以利用基于 模型的像面仿射变换 进行无控制点定向处理，以保证同名光线对对相 交［ 。 由核线定义可知，核线方向与摄影基线的方向 RPC ］ 9 相近似 线 阵 列 影 像 通 常有 异 轨 立 体 ( 。 Across － 。 ，采用侧视的方式) 和同轨立体( ， track Along － track 采用前后视方式) 两种立体观测，不同成像方式的 卫星影像，其核线在影像上的分布状态不同，前者的 核线方向与扫描方向大致相同，后者的核线方向则 为获取核线影像，需要计算 近似垂直于扫描方向 整幅影像的核线方程 点位置影 而核线方程受 沿近似垂直于核线方向前进: 对于异轨 响，故使 ，应当沿着垂直于扫描方向前进; 对于同轨立 立体 体，则应当沿着扫描方向前进 为使拟合的核线尽 可能准确， 当核线对方程 q 与影像区域没有交点的时候，核线方程可以舍去 这样就获得了所有核线对的方程，即 点应当位于影像中间 。 。 。 。 q q q { i ai x + bi y + ci = 0 ( ， … ， n ) i = 1 ( ) 1 式中， = 0 a' i x + b' i y + c' i a、b、c 为核线方程的编号; 为直线方程的系数; 为核线上点 x、y 为核线方程的数目 的坐标; 在重采样时，本文采用了文献［ 的方法逐行进行，再沿着 斜像片上获取影像 q 方向进行线性内插，获得整个区域的核线影像 ” 。 n ］ 13 “ 。 直接在倾 前进 2 核线影像几何模型重建 。 i。 核线影像上的点与原始影像上的点是一一对应 的，而原始影像的几何模型建立了原始影像上点与地 面点的对应关系 根据以上关系可建立核线影像上 点与地面点的对应关系，重建核线影像的几何模型 。 其中，建立核线影像与原始影像对应关系是确定原始 影像点所在核线方程的序号 当核线影像上一点变 换到原始影像时， 可以从其坐标值获得，但是对于现 i 有的摄影测量系统来说，无法直接利用核线方程进行 核线影像前方交会 模型作为通用 传感器模型的特点，重建核线影像的几何模型，使核 线影像可以直接应用于三维信息提取 2． 1 像素点与对应地面点解析关系的确定 本文将利用 像的核线方程，因此对于核线影像上一点 由于核线影像上每一行对应于所求解的整幅影 ， x0 ，但核线的序号为 值即可确定核线的序号 ) 的 RPC M0 。 。 ( i 的 y M0 值进行取整运算，如图 2 y y0 整数，这就需要对 所示 。 图 2 异轨核线影像与原始影像上同名点对应关系 ( 左: 核线影像; 右: 原始影像) Fig． 2 Corresponding relation of tie － points in across － track stereo images and original images 对于异轨立体数据，核线方程是近似于扫描方 值与原始影像上对应点 即可确 ， y0 值相等［ ) 的 ］ 向的，那么 ) 的 x0 ( 13 x0 利用核线方程及 ) 的坐标为 x0 M0 x ， y 的近似值 ( M 定 x y 。 ， y x ( y'。M x = x0 { y = y' + y0 － int y0 而对于同轨立体数据来说，核线方程是近似垂 ) ( ) 2 ( 

第 4 期 张 过，等: 基于 RPC 模型的线阵卫星影像核线排列及其几何关系重建 ·3· 直于扫 描 方 向 的，那 么， x0 值相等，因此利用 与 原 始 影 像 上 对 应 点 值即可确定 的 M0 x 数［ 9 ］，具体步骤如下: ( 1 ) 根据核线影像的覆盖范围，利用美国地质 调查局提供的全球 分辨率 1 km ( DEM Global 30 － ) ，计算该区域的 。 。 Arc － second Digital Elevation Model 最大 最小椭球高 、 ) 在核线影像平面上，以一定的间隔建立规 ( 2 利用核线方程，计算出这些格网点所对应 则格网 原始影像上的坐标，再利用原始影像的 模型将 这些点投影到以一定间隔分层的高程面上，这样即 同时，加 建立核线影像上点与地面点的对应关系 密控制格网和高程面，建立独立的检查点 RPC 。 。 ) 利用控制点坐标计算核线影像与地面坐标 的归一化参数，将控制点和检查点归一化 。 ］，获取核 模型参数［ 9 ) 利用岭估计法求解 ( 线影像的 模型参数，并用检查点进行精度检查 RPC RPC 。 ( 3 4 3 实验结果与分析 SPOT 5 HRS、GeoEye、IRS P5、IKONOS 采用 轨数据， 试验，各个影像的参数如表 SPOT 5 HRG、CBERS 2 － 03 所示 。 1 同 异轨数据进行 M M ， y ( x 的 ) 的 值 y x 。 利用式( 2 ) 即建立了核线影像与原始影像的对 由于原始影像的几何模型参数是已知的， 应关系 即可建立核线影像上一点与地面点的对应关系，即 核线影像的正变换过程，其流程如图 所示 。 3 。 图 3 核线影像的正变换 Fig． 3 The forward transformation of epipolar images 2 ， y0 首先，根据公式( ) 变换为原始影像上对应的点( ( ， y x0 据原始影像的几何模型，将原始影像上点( 到某高程面上，获取其对应地面点坐标( 2． 2 RPC 参数的求解 ) 将核线影像上的像素坐标 ) ; 然后，根 ) 投影 ， y ， h x ， lon lat 。 x ) 采用与地形无关的求解方式求解 模型参 RPC 表 1 实验数据的详细说明 数据源 立体覆盖 区域中心 SPOT 5 HRS IRS P 5 IKONOS GeoEye SPOT 5 HRG CBERS 2 － 03 107°27'19． 758″E 34°24'3． 683″N 116°31'44． 802″E 40°24'7． 869″N 147°16'30． 720″E 40°51'52． 560″S 147°15'31． 500″E 42°51'38． 520″S 116°23'23． 615″E 39°55'28． 692″N 116°27'49． 844″E 39°58'1． 159″N Tab． 1 Detailed specification of datasets in experiments 左影像大小 像素 像素 右影像大小 像素 像素 最小 高程 最大 高程 立体模式 / × / × / m / m 12 000 × 12 000 12 000 × 12 000 384 2 452 12 000 × 12 000 12 000 × 12 000 － 14 1315 12 122 × 13 148 12 122 × 13 148 － 30 1 284 26 928 × 31 668 26 960 × 31 668 － 30 1 284 12 002 × 12 002 12 000 × 12 002 － 27 10 002 × 10 000 10 002 × 10 000 － 28 234 111 同轨 同轨 同轨 同轨 异轨 异轨 影像分 辨率 / m 定向［ ］精 8 像素 度 / 10 0． 160 2． 5 0． 182 1 0． 154 0． 5 0． 146 5 3 0． 115 0． 235 地貌 山地 山地 丘陵 丘陵 平原 平原 其 中， 为 宝 鸡 地 区 影 像， 和 2 － 03 SPOT 5 HRS 为北京市区影像， SPOT 5 HRG 为澳大利亚 CBERS 和 GeoEye 为北京市 丘陵和山地，具有代 、 定向精度为利用立体像对上若干分布均匀的 地形包括平原 地区影像， Hobart IKONOS 怀柔地区影像 表性 P 5 。 。 连接点，采用基于 RPC 控制点定向的中误差［ 15 模型的像面仿射变换进行无 ， 16 ］ 本文采用分母不同的三阶 。 RPC 模型求解核线 个像素， 参数［ 影像 高程分层为 RPC 9 ］，控制点格网间隔为 层，其求解精度如表 200 所示 2 。 15 表 2 核线 RPC 模型控制点和检查点精度 传感器 核 线 影 像 左 SPOT 5 HRS 右 Tab． 2 Accuracy of control points and check points for epipolar RPC model 控制点 检查点 X 最大 中误差 最大 Y 平面 中误差 最大 中误差 最大 X 中误差 最大 Y 中误差 最大 中误差 ( pixels) 平面 － 2． 11E － 02 9． 47E － 05 － 5． 98E － 02 2． 10E － 02 6． 00E － 02 2． 10E － 02 － 8． 72E － 02 4． 18E － 04 4． 10E － 02 1． 35E － 02 8． 93E － 02 1． 35E － 02 － 8． 72E － 02 4． 19E － 04 4． 10E － 02 1． 35E － 02 8． 93E － 02 1． 35E － 02 7． 20E － 04 5． 03E － 06 － 6． 04E － 02 2． 26E － 02 6． 04E － 02 2． 26E － 02 

·4· 国 土 资 源 遥 感 年 2010 ( 续表) 平面 IKONOS 传感器 核 线 影 像 左 IRS P5 右 左 右 左 右 左 SPOT 5 HRG 右 左 CBERS 2 － 03 右 由表 GeoEye 控制点 检查点 X 最大 中误差 最大 Y 平面 中误差 最大 中误差 最大 X 中误差 最大 Y 中误差 最大 中误差 5． 09E － 05 1． 17E － 05 1． 17E － 03 2． 51E － 04 1． 17E － 03 2． 51E － 04 4． 22E － 05 1． 08E － 05 9． 36E － 04 2． 33E － 04 9． 36E － 04 2． 33E － 04 － 1． 33E － 03 1． 89E － 04 － 2． 78E － 03 3． 03E － 04 2． 80E － 03 3． 57E － 04 － 1． 06E － 03 1． 74E － 04 － 2． 77E － 03 2． 81E － 04 2． 81E － 03 3． 31E － 04 － 1． 65E － 05 4． 55E － 06 － 2． 01E － 05 7． 43E － 06 2． 46E － 05 8． 71E － 06 1． 23E － 05 3． 85E － 06 － 1． 70E － 05 7． 03E － 06 2． 09E － 05 8． 01E － 06 － 1． 95E － 05 5． 09E － 06 － 5． 47E － 06 1． 95E － 06 2． 02E － 05 5． 45E － 06 － 1． 45E － 05 4． 29E － 06 － 4． 62E － 06 1． 84E － 06 1． 52E － 05 4． 67E － 06 － 1． 86E － 05 3． 51E － 06 － 3． 66E － 02 1． 51E － 02 3． 66E － 02 1． 51E － 02 － 1． 54E － 05 3． 16E － 06 － 4． 96E － 02 1． 52E － 02 4． 96E － 02 1． 52E － 02 4． 92E － 05 5． 15E － 06 4． 28E － 02 1． 49E － 02 4． 28E － 02 1． 49E － 02 4． 18E － 05 4． 71E － 06 － 3． 38E － 02 1． 39E － 02 3． 38E － 02 1． 39E － 02 2． 99E － 03 1． 95E － 05 1． 29E － 03 4． 70E － 04 2． 99E － 03 4． 71E － 04 － 1． 53E － 02 7． 97E － 05 － 4． 05E － 03 4． 44E － 04 1． 58E － 02 4． 51E － 04 9． 53E － 05 1． 98E － 05 7． 29E － 05 1． 90E － 05 9． 72E － 05 2． 74E － 05 7． 41E － 05 1． 85E － 05 5． 87E － 05 1． 79E － 05 7． 61E － 05 2． 57E － 05 － 2． 08E － 02 1． 09E － 04 － 3． 86E － 03 2． 09E － 05 2． 12E － 02 1． 11E － 04 － 2． 12E － 03 1． 20E － 05 － 3． 98E － 04 5． 77E － 06 2． 15E － 03 1． 33E － 05 2 4． 68E － 06 4． 14E － 06 4． 98E － 04 － 6． 66E － 05 4． 93E － 04 可以看出，在采用分母不相等且在三阶 模型参数求解的检 2． 26 × 10 － 2 像素，控制点平面精度 RPC 多项式的情况下，核线影像 查点平面精度为 为 0． 05 个像素 2． 10 × 10 － 2像素，均小于 理论上，核线模型是不存在上下视差的，但是无 控制点定向和直线拟合过程会引入误差 本文将在 左影像范围按一定间隔布设格网，并将其投影到不 同高程面上，获取其对应的地面点，然后将这些地面 。 。 6． 25E － 06 2． 24E － 04 2． 22E － 04 3． 33E － 06 4． 36E － 06 － 3． 16E － 05 5． 49E － 06 点投影到右影像上去，此时按照核线影像定义右影 值相等，将该两 像上点的 值相减获得上下视差，作为评价重建核线影像几 个 何模型的准确性指标 值应当与左影像上点的 y y y ( 在左影像空间的 200 × 15 素选取一个点，在高程方向分为 几何模型精度如表 最大值在 像素之内，中误差优于 所示 。 3 方向和 本实验格网布设方案为 方向每 200 × 像 y 层) ，获得的核线 核线几何模型上下视差 200 15 x 。 0． 5 表 3 核线模型上下视差 Tab． 3 Vertical parallax of epipolar model 像素 。 0． 2 ( pixels) 传感器 最大差值 中误差 SPOT 5 HRS IRS P5 － 5． 05E － 02 － 2． 50E － 02 1． 19E － 02 4． 57E － 03 IKONOS 6． 63E － 02 1． 96E － 02 GeoEye SPOT 5 HRG CBERS 2 － 03 － 3． 02E － 01 8． 62E － 02 4． 19E － 01 1． 39E － 01 2． 88E － 01 1． 20E － 01 14 ］ 9 。 RPC 进行三维信息提取［ 为验证核线模型的准确性，本文利用生成的核 线影像及 通过核线影 个均匀分布的同名点进行前方交会，获取同 像上 名点的地面坐标，将这些同名点利用核线方程投影 到原始影像，利用原始影像的 模型前方交会， 获取其投影点的地面坐标，将两次计算的地面点坐 标相减来评价核线影像几何模型的实际精度，结果 如表 所示 RPC 4 。 表 4 核线影像前方交会与原始影像前方交会坐标差 Tab． 4 Difference between epipolar images forward intersection and original images forward intersection ( m) 数据源 高程 纬度 经度 最大差值 中误差 最大差值 中误差 最大差值 中误差 SPOT 5 HRS 0． 103 0． 052 － 0． 499 0． 248 0． 003 0． 002 IRS P5 0． 176 0． 11 － 0． 883 0． 617 IKONOS 0． 399 GeoEye 0． 563 0． 234 0． 311 － 1． 421 0． 837 － 1． 204 0． 678 0． 166 0． 479 0． 586 0． 12 0． 24 0． 289 4 结论 ( 1 ) 采用不同地形 核线影像几何模型，控制点中误差最大为 素，最大误差为 不同星载立体影像重建的 像 像素; 检查点中误差最大为 0． 021 、 0． 089 像素，最大误差为 0． 026 模型重建几何精度不损失 。 0． 089 像素，核线影像 RPC ( 2 ) 利用投影轨迹法获取的星载推扫式影像核 像素之内，中误 像素，验证本文的核线影像和核线模型 线几何模型，其上下视差最大在 差小于 制作方法不受地形影响 0． 2 0． 5 ( 3 地面坐标差值在平 面 上 小 于 ) 核线影像前方交会和原始影像前方交会的 ，高 程 差 值 小 于 ，说明利用本文方法重建的核线几何模型可 1 m 0． 3 m 以直接用于三维信息提取而不损失精度 。 。 SPOT 5 HRG CBERS 2 － 03 0． 08 0． 035 0． 043 0． 027 0． 035 0． 025 0． 446 0． 268 － 0． 106 0． 061 0． 136 0． 081 参考文献: ［ ］ 张祖勋，周月琴 1 ． 用拟合法进行 可以看出，计算的地面点在平面上坐标差的中 误差在 1 m 之内，高程差在 0． 3 m 之内 。 ］ ［ J ． 武汉测绘科技大学学报， ) : 1989 三线阵 20 － 25． 2 卫星影像的模拟研究 ［ ］ 江万寿，张剑清，张祖勋 2 ． ］ ［ J ． 武汉大学学报: 信息科学版， 2002 ， 27 ( 4 ) : 414 － 419． SPOT ， 14 ( CCD 影像的近似核线排列 

第 4 期 张 过，等: 基于 RPC 模型的线阵卫星影像核线排列及其几何关系重建 ·5· ［ ］ 3 Gupta R ， Hartley R I． Linear Pushbroom Cameras Pattern Analysis ［ ］ 张 过 9 ． ］ 缺少控制点的高分辨率卫星遥感影像几何纠正［ D ． and Machine Intelligence ［ ］ J ysis and Machine Intelligence ［ ］ 4 Habib A F ， Morgan M ， Jeong S ， ． IEEE Transactions on Pattern Anal- ， 1997 ， et al． Analysis of Epipolar Geom- 963 － 975． ) : 19 9 ( ［ ］ 5 ［ ］ 6 etry in Linear Array Scanner Scees ． The Photogrammetric Re- ［ ］ J ， 2005 ， ( ) : 20 cord 27 － 47． 109 ， Park W． Epipolarity Analysis for Linear Pushb- Lee Hae － Yeoun ］ ［ ， J ． International Symposium on Remote Sensing room Imagery : ， 17 2001 593 － 598． Kim T． A Study on the Epipolarity of Linear Pushbroom Images ］ ［ ， J 66 ． Photogrammetric Engineering and Remote Sensing ( ) : 2000 ， 8 961 － 966． 武汉: 武汉大学， 2005． ［ ］ 叶新魁，文贡坚，王继阳，等 10 基于零高程点的卫星影像核线 ． ］ 确定方法［ J ． 测绘信息与工程， 2009 ( ， 34 2 ) : 28 － 31． ［ 11 ］ 张永军，丁亚洲 ． ］ 影像的生成［ J ． 基于有理多项式系数的线阵卫星近似核线 ) : 武汉大学学报: 信息科学版， ， ( 2009 34 9 1068 － 1071． ［ ］ 胡 芬，王 密，李德仁，等 12 基于投影基准面的线阵推扫式卫星 ． 立体像对近 似 核 线 影 像 生 成 方 法［ ］ ． J ( ) : 数字摄影测量学［ ］ ． M ． 5 428 － 436 ］ 张祖勋，张剑清 ［ 13 社， 2001 : 117 － 118． 测 绘 学 报， 2009 ， 38 武汉: 武汉大学出版 ［ ］ 7 Morgan M ， Kim K O ， Soo J ， et al ． Epipolar Resampling of Space ［ ］ 14 － borne Linear Array Scanner Scenes Using Parallel Projection ［ ］ J ( ． Photogrammetric Engineering ＆ Remote Sensing ) : ， 72 2006 ， 11 1255 － 1263． ［ ］ 8 Ono T． Epipolar Resampling of High Resolution Satellite Imagery ［ ］ C I /3 and IV /4 on Sen- ∥Joint Workshop of ISPRS WG I /1 ， 1999． sors and Mapping from Space． Hanover ， Grodecki J． ［ ］ C / / ASPRS Annual Conference ， St Louis ， 2001． IKONOS Stereo Feature Extraction － RPC Approach ［ ］ 15 Grodecki J ， Dial G． Block Adjustment of High － resolution Satellite Images Described by Rational Polynomials ， Engineering and Remote Sensing 2003 ， 69 ． Photogrammetric ［ ］ J : 59 － 68． ［ 16 ］ 祝小勇，张 过，秦绪文 ) : 土资源遥感， ( 2009 2 国产光学卫星影像 ． ］ 制作［ J ． 国 RPC 32 － 34． Epipolar Resampling and Epipolar Geometry Reconstruction of Linear Array Scanner Scenes Based on RPC Model ZHANG Guo1， PAN Hong － bo1， JIANG Wan － shou1， QIN Xu － wen2 ( 1． State Key Laboratory of Information Engineering in Surveying ， 2． China Geological Survey Beijing 100037 China ， ) ， Mapping and Remote Sensing Wuhan 430079 ， ， China ; : Abstract In view of the nonexistence of geometric model for epipolar image ， this paper ， based on the RPC model ， puts forward a practical way of building geometric model from epipolar image generated by the projection track method． With the support of across － track stereo image data of SPOT 5 HRG CBERS 2 － 03 and along － track ， stereo image data of SPOT 5 HRS P5 IKONOS this paper indicates that the RMS of vertical parallax of the epipolar image is within 0． 2 pixel that the planar RMS calculated from forward intersection by the epipolar ， ， ， GeoEye ， ， ， model and the original image model is within 1m and that the altitude RMS is within 0． 3 m． Projection track method Epipolar geometric model RPC model Linear array scanner scenes ; ; ; : Key words 第一作者简介: 张 过( ) ，男，副教授，主要从事航天摄影测量研究工作 1976 － 。 ( 责任编辑: 刁淑娟)