专利名称：：基于相位一致特征的多光谱图像配准检测修正方法
技术领域：
：本发明属于图像预处理方法，具体涉及一种基于相位一致特征的多光谱遥感图像配准检测修正方法。
背景技术：
：多光谱卫星遥感图像是指卫星对同一地域拍摄的不同光谱下的图像，最简单的多光谱图像就是对同一地点拍摄的R、G、B三色图像。这些图像在使用过程中可能会将其再叠放在一起，此时由于图像间微小的位置差异可能导致叠放的图像出现"重影"。这种情况非常不利于图像的精确分析。常用的图像配准方法分为基于图像灰度和基于图像特征两大类。基于图像灰度的方法，直接利用图像的灰度信息，建立两幅图像之间相似性度量，不需要对图像进行特征提取，省去了特征匹配这个难点，由于用到了全部的图像信息，因此提高了配准的精度和鲁棒性，但同时也存在计算量大、速度较慢等缺点。基于图像特征的方法，通过提取参考图像和待配准图像的一些共同特征，建立特征之间的对应关系，从而提取控制点对，实现图像的自动配准。这类方法的优点是它提取了图像的显著特征，大大压縮了图像的信息量，使得计算量小、速度较快，而且对图像灰度的变化具有鲁棒性，但特征匹配的难度较高，易出现误匹配。与基于灰度的方法相比，基于特征的方法不是直接在图像的灰度层次上计算两块图像之间的相似性，而是在特征的层次上进行配准，因此它对图像灰度变化不太敏感，更适合于灰度相差较大的多光谱不同谱段图像之间的配准。根据其关键步骤特征选择和特征匹配方法的不同，具体的基于特征的图像自动配准方法是多种多样的，但这类方法的基本处理流程是一致的首先要对待配准图像进行特征提取，再利用提取的特征完成两幅图像特征之间的匹配，通过特征的匹配关系建立图像之间的配准映射变换。由于卫星多光谱遥感图像不同谱段间的差异复杂性，现有的基于特征的图像自动配准方法仍然不能同时稳定可靠地解决特征稳定提取与特征误匹配问题。
发明内容本发明的目的是针对现有技术缺陷，提供一种能精确对多光谱图像波段配准精度进行检测和修正，修正后图像可以很好匹配，且计算量不大的基于相位一致特征的多光谱遥感图像配准检测修正方法。本发明是这样实现的基于相位一致特征的多光谱遥感图像配准检测修正方法，包括下述步骤(1)指定参考波段图像和待评价波段图像在同一组多光谱图像中指定任意一幅作为参考波段图像，其余图像作为待评价波段图像；(2)在参考波段图像和待评价波段图像中进行相位一致特征点提取在步骤(2)处理后的参考波段图像和其中一幅需要修正的待评价波段图像中提取特征点，特征点的数量范围为50200个，在步骤(2)中根据相位一致性大小，由大到小的开始选取特征点，取满预设特征点的数量；(3)计算相关系数阵S首先对每个特征点均以该点为中心取像素窗，像素窗的大小为7X715X15，参考波段图像和待评价波段图像中像素窗的大小必须一致，对每个参考波段图像中的像素窗均计算与待评价波段图像中每个像素窗中的相关系数，该相关系数组成相关系数阵S;(4)计算参数T参数T为a6r=少001其中a、b、c、d、x、y均为待求参数，这几个参数按照下述方法计算得到。在步骤(3)中得到的S阵中取相关度最大的6个点，分别得到这6个相关系数对应的参考波段图像和待评价波段图像特征点，参考图像中的6个特征点的灰度值按照X";2Xa-3X"X"、6'的形式排列成A阵，待评价波段图像中的6个特征点的灰度少W>^3111凡.4兀5111值按照、51、61的形式排列成B阵，设定TA=B，求出T阵中的每个未、1、2义/"3~4>Vl>V2少r41111知参数；(5)计算空间一致阵C对每个待评价波段图像特征点的坐标均计算经过T阵变换后的坐标值，检验变换后的坐标值与参考波段图像中每个特征点坐标的欧氏距离，若距离小于等于预设值则判定为l，大于预设值则判定为O，这些值构成空间一致阵C(6)计算匹配矩阵M和新的参数T'令空间一致阵C与步骤(3)中计算得到的相关系数阵S点乘，得到的就是匹配矩阵M;用匹配矩阵M中的所有点按照步骤(5)的方法计算新的变换模型T'参数；(7)判定AT是否小于预设值令Aa=a-a，，Ab=b-b，，Ac=c-c，，Ad=d-d，，Ax=x-x，，Ay=y-y，，令Ar=^/a"2++Ac2+++Ay2判断AT是否小于0.0001，若AT小于0.0001则执行步骤(8)，否则执行步骤(5)计算参数T，此时用于计算参数T的矩阵是匹配矩阵M;(8)对图像进行校正用匹配成功的匹配特征点对待评价波段图像进行现行校正，校正方法选用现有的线性校正方法。如上所述的基于相位一致特征的多光谱遥感图像配准检测修正方法，其中，在步骤(2)前先在整幅图像中先均匀的设定小区域，特征点在这些设定的小区域内选取，所述的区域大小为100X100300X300的正方形像素区。如上所述的基于相位一致特征的遥感图像检测修正方法，其中，所述的均匀区域分布是5X5、10X10等。如上所述的一种基于相位一致特征的遥感图像检测修正方法，其中，在步骤步骤(3)选取特征点时在相位一致图中用曲面插值法精确定位特征点。如上所述的一种基于相位一致特征的遥感图像检测修正方法，其中，步骤(4)中的像素窗大小为9X9。本发明的显著效果是该方法不仅利用了相位一致特征对遥感图像亮度和对比度具有不变性，能捕捉图像的不变特征，同时适用于阶跃边缘、线边缘、屋脊型边缘及介于阶跃边缘和线边缘之间的边缘信息等多种边缘特征的检测，并且相位一致是一个无量纲的量，可以对大多数图像使用较为固定的阈值，相位一致特征检测算法对线形物体产生单线响应，抗噪能力强，检测到的特征不会因为噪声而出现波动，具有更稳定的特征定位精度的优点；而且采用特征相似和空间一致相结合匹配策略实现了多光谱图像波段配准高精度检测，根据检测结果修正了遥感图像地面数据处理系统中波段配准精度，使多光谱配准精度提高到了0.3像元。图1是在图像中设定小区域的示意图(5X5)。具体实施例方式基于相位一致特征的多光谱遥感图像配准检测修正方法，包括下述步骤(1)指定参考波段图像和待评价波段图像所述的参考波段图像和待评价波段图像都是卫星拍摄的多光谱图像，这些图像是同一地点不同波段的图像。指定其中的任意一幅图像作为参考波段图像，其余图像均作为待评价波段图像。最终将所有待评价波段图像指定过程不需要任何其它要求。[OOSS](2)相位一致特征点提取采用相位一致算法(phasecongre皿cy)对参考波段图像和待评价波段图像进行相位一致特征提取。在步骤(2)中特征点的数量范围为50200个，特征点从步骤(2)中得到的图像中相位一致性由大到小的开始选取，取满预设特征点的数量。参考波段图像和待评价波段图像中特征点的数量可以相同也可以不同，例如在参考波段图像中特征点数量为50个，待评价波段图像中特征点数量为200个。(3)计算相关系数阵S首先对每个特征点(包括参考波段图像和待评价波段图像中的特征点)均以该点为中心取像素窗，像素窗的大小为7X715X15，例如7X7、9X9、11X11、15X15。即使由于特征点位置过近而使像素窗部分重叠也不影响本方法的继续执行。像素窗中的值是步骤(2)中得到的相位值。参考波段图像和待评价波段图像中像素窗的大小必须一致。对每个参考波段图像中的像素窗均计算待评价波段图像中每个像素窗中的相关系数。例如在参考波段图像中取50个特征点，本步骤中的像素窗取7X7，那么在参考波段图像中会有50个6像素窗；类似的若在待评价波段图像中取200个特征点，那么在待评价波段图像中会有200个像素窗，本步骤计算参考波段图像中50个像素窗中的每个像素窗与待评价波段图像中200个像素窗中的每个像素窗的相关系数，该相关系数组成相关系数阵S。对本例中S应为50X200的矩阵，矩阵的行为参考波段图像中的像素窗的序号，矩阵的列为待评价波段图像中像素窗的序号，矩阵的值为相应像素窗之间的相关系数。(4)计算参数T参数T的形式为<formula>formulaseeoriginaldocumentpage7</formula>其中a、b、c、d、x、y均为待求参数。这几个参数按照下述方法计算得到。在步骤(3)中得到的S阵中取相关度最大的6个点(即矩阵数值最大的6个点)分别得到这6个相关系数对应的参考波段图像和待评价波段图像特征点，即参考图像中的特征点6个，待评价波段图像中的特征点6个。参考图像中的6个特征点的灰度值按照x,iir2A:r3xMx,5x,6'的形式排列成A阵，待评价波段图像中的6个特征点的灰度<formula>formulaseeoriginaldocumentpage7</formula>的形式排列成B阵，设定TA=B，由于A阵和B阵中所有的值均为已知，因此可以唯一的求出T阵中的每个；(5)计算空间一致阵C对每个待评价波段图像特征点的坐标均计算经过T阵变换后的坐标值，检验变换后的坐标值与参考波段图像中每个特征点坐标的欧氏距离，若距离小于等于预设值则判定为l，大于预设值则判定为O，这些值构成空间一致阵C。同样以参考波段图像中取50个特征点，待评价波段图像中取200个特征点为例，每个待评价波段图像中的特征点均有坐标，用[xsi，ysi]表示，其中i表示第i个特征点，本例中i取值范围为1200。按照公式T[xsi，ysi，l]'=[x，si，y，si，l]，计算变换后的坐标[x'si，y'si]，然后计算该坐标与参考波段图像中取50个特征点中每个特征点坐标[Xrk，yA]的距离，其中k表示第k个特征点，本例中k取值范围为150。若距离小于等于预设值则判定结果为1，否则判定结果为0。本例中所述的预设值是一维高斯函数方差取0.5时的值，当然也可以根据实际情况取其它数值。所得的判定结果组成空间一致阵C，该阵也应是50X200的矩阵。计算完成后C阵应呈现类似对角阵的形式，即其中部分值以对较阵的形式为l，其余为0或散点1值。(6)计算匹配矩阵M和新的参数T'令空间一致阵C与步骤(3)中计算得到的相关系数阵S点乘，得到的就是匹配矩阵M。用匹配矩阵M中的所有点按照步骤(4)的方法计算新的参数T'。本步骤使用的参数多于计算参数T'需要的数量，可以用拟合的方法计算T'中的每个参数。计算后"=[(7)判定AT是否小于预设值令<formula>formulaseeoriginaldocumentpage8</formula>判断AT是否小于0.0001，若AT小于0.0001则执行步骤(8)，否则执行步骤(4)计算参数T，只是此时用于计算参数T的矩阵不再是相关系数阵S，而是匹配矩阵M，其余计算过程均按照步骤(4)中的步骤计算。[OO56](8)对图像进行校正用匹配成功的匹配特征点对待评价波段图像进行现行校正，校正方法可以选用现有的线性校正方法。所述的匹配特征点是指在最后一次循环中步骤(7)中匹配矩阵M判定为非零的参考波段图像和待评价波段图像对应的特征点。所述的线性校正方法是指用这些匹配成功的特征点求公式y=kx+b中的参数k和参数b，其中x为待评价波段图像中的特征点，y为参考波段图像中的特征点，匹配成功的特征点往往多于计算参数k，b所需要的值，此时可以用拟合的方式求得所需参数。求得参数k，b后就可以用该参数对地面数据处理系统进行图像校正了。为了避免在步骤(2)选取特征点中特征点的选择不均匀，在整幅图像中先均匀的设定小区域，如5X5、10X10等。如附图1所示。特征点只能在这些设定的小区域内选取，后续计算时仅针对参考波段图像和待评价波段图像相对应的小区域进行计算，这样既可以保证在大幅卫星遥感图像上特征点均匀分布，避免特征点过于随机选取，又可以大大减少计算量。所述的区域范围是100X100300X300的正方形像素区，例如100X100，150X150，200X200，250X250，300X300都是可以选择的。为了进一步提高校正精度，在步骤(2)选取特征点时在相位一致图上采用曲面插值法定位特征点，提高了特征点的定位精度，使得基于相位一致特征的多光谱遥感图像配准检测修正方法的检测精度提高到0.001个像素，而上述方法的精度仅为1个像素。权利要求基于相位一致特征的多光谱遥感图像配准检测修正方法，包括下述步骤(1)指定参考波段图像和待评价波段图像在同一组多光谱图像中指定任意一幅作为参考波段图像，其余图像作为待评价波段图像；(2)在参考波段图像和待评价波段图像中提取相位一致特征点在步骤(2)处理后的参考波段图像和其中一幅需要修正的待评价波段图像中提取特征点，特征点的数量范围为50～200个，特征点从步骤(2)中得到的图像中相位一致性由大到小的开始选取，取满预设特征点的数量；(3)计算相关系数阵S首先对每个特征点均以该点为中心取像素窗，像素窗的大小为7×7～15×15，参考波段图像和待评价波段图像中像素窗的大小必须一致，对每个参考波段图像中的像素窗均计算待评价波段图像中每个像素窗中的相关系数，该相关系数组成相关系数阵S；(4)计算参数T参数T为<mrow><mi>T</mi><mo>=</mo><mfencedopen='['close=']'><mtable><mtr><mtd><mi>a</mi></mtd><mtd><mi>b</mi></mtd><mtd><mi>x</mi></mtd></mtr><mtr><mtd><mi>c</mi></mtd><mtd><mi>d</mi></mtd><mtd><mi>y</mi></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>其中a、b、c、d、x、y均为待求参数，这几个参数按照下述方法计算得到。在步骤(3)中得到的S阵中取相关度最大的6个点，分别得到这6个相关系数对应的参考波段图像和待评价波段图像特征点，参考图像中的6个特征点的灰度值按照<mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>4</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>5</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>r</mi><mn>6</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>4</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>5</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>r</mi><mn>6</mn></mrow></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced>的形式排列成A阵，待评价波段图像中的6个特征点的灰度值按照<mfencedopen='['close=']'><mtable><mtr><mtd><msub><mi>x</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>s</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>s</mi><mn>4</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>s</mi><mn>5</mn></mrow></msub></mtd><mtd><msub><mi>x</mi><mrow><mi>s</mi><mn>6</mn></mrow></msub></mtd></mtr><mtr><mtd><msub><mi>y</mi><mrow><mi>s</mi><mn>1</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>s</mi><mn>2</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>s</mi><mn>3</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>s</mi><mn>4</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>s</mi><mn>5</mn></mrow></msub></mtd><mtd><msub><mi>y</mi><mrow><mi>s</mi><mn>6</mn></mrow></msub></mtd></mtr><mtr><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced>的形式排列成B阵，设定TA＝B，求出T阵中的每个未知参数初值；(5)计算空间一致阵C对每个待评价波段图像特征点的坐标均计算经过T阵变换后的坐标值，检验变换后的坐标值与参考波段图像中每个特征点坐标的欧氏距离，若距离小于等于预设值则判定为1，大于预设值则判定为0，这些值构成空间一致阵C；(6)计算匹配矩阵M和新的变换模型T’参数令空间一致阵C与步骤(4)中计算得到的相关系数阵S点乘，得到的就是匹配矩阵M；用匹配矩阵M中的所有点按照步骤(4)的方法计算新的变换模型T’参数；(7)判定ΔT是否小于预设值令Δa＝a-a’，Δb＝b-b’，Δc＝c-c’，Δd＝d-d’，Δx＝x-x’，Δy＝y-y’，令<mrow><mi>ΔT</mi><mo>=</mo><msqrt><mrow><mi>Δ</mi><msup><mi>a</mi><mn>2</mn></msup></mrow><mo>+</mo><mi>Δ</mi><msup><mi>b</mi><mn>2</mn></msup><mo>+</mo><mi>Δ</mi><msup><mi>c</mi><mn>2</mn></msup><mo>+</mo><mi>Δ</mi><msup><mi>d</mi><mn>2</mn></msup><mo>+</mo><mi>Δ</mi><msup><mi>x</mi><mn>2</mn></msup><mo>+</mo><mi>Δ</mi><msup><mi>y</mi><mn>2</mn></msup></msqrt></mrow>判断ΔT是否小于0.0001，若ΔT小于0.0001则执行步骤(8)，否则执行步骤(4)计算参数T，此时用于计算参数T的矩阵是匹配矩阵M；(8)对图像进行校正用匹配成功的匹配特征点对待评价波段图像进行现行校正，校正方法选用现有的线性校正方法。2.如权利要求1所述的一种基于相位一致特征的多光谱遥感图像配准检测修正方法，其特征在于在步骤(2)前先在整幅图像中先进行均匀的提取检测区域，如5X5、10X10等，特征点在这些设定的小区域内选取。3.如权利要求2所述的一种基于相位一致特征的多光谱遥感图像配准检测修正方法，其特征在于所述的区域范围是IOOX100300X300的正方形像素区。4.如权利要求1所述的一种基于相位一致特征的多光谱遥感图像配准检测修正方法，其特征在于在步骤(2)选取特征点时用相位一致图曲面插值法精确定位特征点。5.如权利要求1所述的一种基于相位一致特征的遥感图像检测修正方法，其特征在于步骤(3)中的像素窗大小为9X9。全文摘要本发明属于一种基于相位一致特征的多光谱图像配准检测修正方法。一种无热化的红外光学镜头包括保护窗，其特征在于还包括同光轴顺次排列的入射镜、透镜、出射镜和衍射镜。使用本发明的效果是在不增加测温反馈系统及机械调焦等机构的情况下，只依靠光学器件本身的配合就能达到下述性能使用温度范围-40℃～60℃，焦距150mm，视场4.74°±1°，工作波段3.7μm～4.8μm，单调线性相对畸变≤1％，相对孔径≥1/2，后工作距23mm。尤其是-40℃～60℃大范围的温度应用范围是传统红外光学镜头在不增加测温反馈系统及机械调焦等机构的前提下不能达到的。文档编号G06T5/00GK101751673SQ20091026582公开日2010年6月23日申请日期2009年12月24日优先权日2009年12月24日发明者徐建艳,王洪海,陆书宁申请人:中国资源卫星应用中心