一种用于克服有理函数模型病态性的改进的奇异值修正方法技术

本发明专利技术有效地克服了遥感影像有理函数模型成像过程中的病态性问题，可以应用于遥感影像成像领域。有理函数模型理论的提出已经有若干年，其系数解算过程中的病态性问题一直得不到有效的解决，本发明专利技术基于对系数矩阵奇异值的修改，提高小奇异值部分的稳定性，从而保证了解算过程的稳定。该方法可以用于多种遥感影像的正射纠正。

【技术实现步骤摘要】

本发有效地克服了遥感影像有理函数模型成像过程中的病态性问题，可以应用于遥感影像成像领域。有理函数模型理论的提出已经有若干年，其系数解算过程中的病态性 问题一直得不到有效的解决，本专利技术基于对系数矩阵奇异值的修改，提高小奇异值部分的 稳定性，从而保证了解算过程的稳定。
技术介绍
建立遥感影像与被测物体之间的成像模型，还原其消除干扰后的数据，使影像数 据和被测物体一一对应，是遥感数据处理的一部分，也是数据分析和应用的前提。传统的严 格物理成像模型根据传感器成像过程建立，模型中各分量物理意义明确，模型精度高，应用 广泛。但是，严格传感器模型强烈依赖于传感器，不同的传感器使用不同的传感器模型。随 着，科技的进步，遥感卫星日渐增多，这给传感器模型的建立和更新带来了很大的困难。此 外，就动态传感器而言，它们的方位参数是时刻变化的，对于这类传感器，严格物理模型也 有其局限性。再则，对于一些没有公开星历参数的卫星，严格物理模型就根本无法建立。 有理函数成像模型克服了以上严格物理成像模型的诸多不足。但是，运用传统的 最小二乘迭代算法解算有理函数模型时，由于方程病态性的存在，使得解算过程不稳定，解 的误差较大甚至错误。有理函数模型最小迭代解法大致过程如下 有理函数模型基本方程式" <formula>formula see original document page 3</formula> + fl15ZY"2 +a16Hf2 +fl17ZJ +a18rJ +a19ZJ (1-2) 构建误差方程<formula>formula see original document page 3</formula> (1-5) 5 = (i z y ;r …z3)(i & ' j = (fl0 a, …a19 " 62 … 619)T ￡> = (i z r义…r3 jt3)(i《《=(c0 ct … c19 《 J2利用(1-4)化简(1-1)可得(1-4)J_ Z I I 555 ￥j_工工工 D D D 5转化为矩阵形式为1 z, 1 z，1 z 《-A 《-"JT3r73一7 ''V-《7 .v-二5K一(1-6) 最后，利用最小二乘迭代法解算(1-6)。但是由于系数矩阵存在病态性，导致求得 得解不稳定。
技术实现思路
本专利技术立足于有理函数模型病态性的机理，通过增大其系数矩阵奇异值方阵中较 小的奇异值，克服其不稳定的部分，从而达到克服其病态性的目的。 本专利技术方法基于如下理论系数矩阵进行奇异值分解以后，相对大的奇异值表示 模型参数中比较可靠和稳定的部分，而相对较小的奇异值表示模型中不可靠或不稳定的部 分。本方法旨在保留较大的奇异值不变，增大较小的奇异值。目的是提高模型中不可靠和 不稳定的部分的稳定性，以此来减弱方程的病态性。具体做法为设t为截断奇异值法保留 的最小奇异值的门限，q为对应小于t的奇异值个数，用以下公式进行奇异值修正 a k' = a k， k《m_q(1-7) 奇异值，t-其中，m为奇异值的总个数，a a——为修改前的第k个奇异值，ak'为相应修改后的1000相应的修改后的奇异值矩阵为D' = diag(a / ， a 2' ... a k'(1-8)附图说明 图1是使用本专利技术方法的具体实例处理过程示意图； 图2是整个过程中的迭代求解部分的过程示意图。 具体实施过程 下文结合说明书附图1，以Landsat-5卫星遥感影像数据为例，对本专利技术的具体实施方式作详细说明。本专利技术所述方法包含但不限于所举实例。本方法使用于多种不同的卫星遥感影像的有理函数校正模型。 步骤1 :读取已经标准化后的数据； 步骤2 :利用已经读入的数据构建有理函数模型的系数矩阵； 步骤3 :求解初始值，初始值的精度直接影响后面迭代求解的收敛速度，第一次求 解最好结合奇异值分解，利用增大奇异值的方法求解； 步骤4 :设定循环条件，这里的循环条件是连续两次解的差的中误差，当中误差大 于设定阈值时进行循环迭代求解。另外，要设定最大循环次数，以防不收敛而陷入死循环； 步骤5 :利用最后的解构建有理函数模型； 步骤6 :利用检查数据进行模型的精度验证。权利要求的主要步骤是步骤(1)读取已经标准化后的数据；步骤(2)利用已经读入的数据构建有理函数模型的系数矩阵；步骤(3)求解初始值，初始值的精度直接影响后面迭代求解的收敛速度，第一次求解最好结合奇异值分解，利用增大奇异值的方法求解；根据步骤(4)设定循环条件，这里的循环条件是连续两次解的差的中误差，当中误差大于设定阈值时进行循环迭代求解。另外，要设定最大循环次数，以防不收敛而陷入死循环；步骤(5)利用最后的解构建有理函数模型；步骤(6)利用检查数据进行模型的精度验证。2. 根据权利要求1所示方法步骤(2)，其特征是需要构建联合解答系数矩阵，即对如 (1)式和(2)式所示模型的参数进行联合求解。<formula>formula see original document page 2</formula>+ a15ZZ2 + (316汉2 + a17Z3 + a18y3 + o19X3 (2)3. 根据权利要求l所示，方法步骤(3)，其特征是在进行迭代求解前，初值计算中使用增大奇异值算法。4. 根据权利要求l所示，方法步骤(4)，其特征是在迭代过程中，每一次都使用增大奇异值算法进行。5. 根据权利要求1所示，方法步骤(3)和(4)中涉及的增大奇异值算法，其特征在于 对系数矩阵经奇异值分解后的奇异值矩阵进行如公式(3)的修改<formula>formula see original document page 2</formula>其中，ak为修改前的第k个奇异值，ak'为相应修改后的奇异值。全文摘要本专利技术有效地克服了遥感影像有理函数模型成像过程中的病态性问题，可以应用于遥感影像成像领域。有理函数模型理论的提出已经有若干年，其系数解算过程中的病态性问题一直得不到有效的解决，本专利技术基于对系数矩阵奇异值的修改，提高小奇异值部分的稳定性，从而保证了解算过程的稳定。该方法可以用于多种遥感影像的正射纠正。文档编号G06F17/16GK101719118SQ20091022358公开日2010年6月2日 申请日期2009年11月24日 优先权日2009年11月24日专利技术者焦伟利, 陈立波 申请人:中国科学院对地观测与数字地球科学中心本文档来自技高网...

【技术保护点】
一种用于克服有理函数模型病态性的改进的奇异值修正方法的主要步骤是：步骤（１）：读取已经标准化后的数据；步骤（２）：利用已经读入的数据构建有理函数模型的系数矩阵；步骤（３）：求解初始值，初始值的精度直接影响后面迭代求解的收敛速度，第一次求解最好结合奇异值分解，利用增大奇异值的方法求解；根据步骤（４）：设定循环条件，这里的循环条件是连续两次解的差的中误差，当中误差大于设定阈值时进行循环迭代求解。另外，要设定最大循环次数，以防不收敛而陷入死循环；步骤（５）：利用最后的解构建有理函数模型；步骤（６）：利用检查数据进行模型的精度验证。

【技术特征摘要】

【专利技术属性】
技术研发人员：陈立波，焦伟利，
申请(专利权)人：中国科学院对地观测与数字地球科学中心，
类型：发明
国别省市：11[中国|北京]