TH 1A 数据处理GXL解决方案.pdfVIP

基于PCI GXL 的TH1A 数据正射纠正 第一次全国地理国情普查时点校正采用了大量的天绘一号（TH01 ）1A 级别 的卫星数据，该数据具有幅宽大、更新快、提供及时等优势，但由于该级别数据 提供的RPC 存在较大系统误差，常规纠正方法无法满足项目精度要求，本文采用 GXL 软件，通过大量控制点反算RPC，较好的解决了该问题。 天绘一号卫星影像产品分为初级、中级、高级3 类。1A 级、2 级为初级产品， 1B 级、3A 级为中级产品，正射影像（DOM）、数字高程模型（DEM）为高级产 品。按有效载荷类型分为5 米三线阵影像、2 米全色影像和10 米多光谱影像。 表1 天绘一号卫星产品列表 级别 描述 备注 1A 级 原始数据经过相对辐射校正后得到的卫星影像 TIF+RPC/TIF+XML 产品，包括影像数据、元数据和浏览图。 2 级 1A 级数据产品利用系统参数经过几何校正后得 到的卫星影像产品，包括影像数据、元数据和浏 览图 1B 级 1A 级数据产品经过摄影测量处理后得到的卫星 TIF+RPC+XML 影像产品，包括影像数据、元数据、RPC 参数和 浏览图 3A 级 1A 级数据产品利用系统参数和地面控制点经过 几何校正后得到的卫星影像产品，包括影像数 据、元数据和浏览图 本次时点校正采用的天绘数据以1A 级为主，1B 级为辅。1B 级数据精度高， 处理比较容易，采用RPC+少量控制点即可达到精度要求。而1A 级数据由于RPC 没有经过摄影测量处理，存在较大的畸变误差，无论采用多少控制点均不能满足 精度要求，为解决该问题，用户普遍采用几何纠正和局部纠正的方法进行处理。 采用几何校正的方法对于平坦地形是可行的，但对于山区，有时地形起伏的 影响可能比系统畸变还要大，而几何校正不考虑DEM，即使采用高阶拟合也很 难满足整体精度要求。采用局部纠正的方法虽然可以解决部分问题，但无疑大大 增加了处理的繁琐程度。以上两种方法都不能有效解决TH1A 的纠正问题，因此， 即使没批次任务中仅有几景1A 数据，也会大大增加用户的工作量，甚至影响整 批成果的提交。 从TH1A 和1B 数据的特点我们可以看出，两种的级别的数据影像文件完全 一致，唯一不同的是RPC 文件，如果我们可以通过软件计算一个与1B 级数据精 度相当的RPC，则可以彻底解决1A 级数据的纠正问题。而GXL 则恰恰具备这种 能力。  可靠的控制点采集算法( ＦＦＴＰ) 能否快速采集到足够的、准确的、分布均匀的控制点是计算RPC 的关键。 地理国情监测底图普遍采用高分辨率卫星影像甚至航片，而时点校正采用TH 数据分辨率为2 米且图面质量普遍较差，人工选点效率低，常规软件的基于 灰度相关的自动采点很难采集到足够的控制点，因此具备计算RPC 的能力往 往也得不到好的结果。GXL 除了具备灰度相关的控制点自动采集算法外，还 提供了基于频率的控制点自动匹配算法（FFTP），该算法具有以下特点： 1. 提取速度快，是常规控制点采集算法效率的5 倍以上。 2. 对不同传感器、不同分辨率、不同坐标系统等均提取到可靠的控制点。 3. 对噪声和光照变化不敏感 4. 精度可达1/1000 像素 5. 可采集成千上万的GCP 基于FFTP 的以上特点，通过仔细调整系统参数，在几乎所有的情况下 都可以采集到足够的控制点，经过GXL 系统对控制点进行优化，基本可以得 到满足计算RPC 所需的控制点。  ２．根据控制点计算ＲＰＣ文件 目前，卫星成像模型主要是严格物理模型和有理函数模型，严格物理模型基 于共线方程，描述了影像和地面之间的严格几何成像关系，其参数具有明确的物 理意义，但是，严格物理模型与卫星传感器相关，形式复杂，缺乏通用性，需要 专业软件的才能支持。而有理函数模型是一种通用的广义成像模型，形式简单， 参数没有物理意义，精度与严格物理模型相近，商业软件基本都支持。因此，鉴 于