・其
他－
一种基于相位一致性相关的多源遥感影像配准方法
范登科，潘
励，叶沅鑫
（武汉大学遥感信息工程学院，武汉４３００７９）
摘要：针对几何校正过程中多源遥感影像同名点匹配率低的问题，提出一种基于相位一致性相关的遥感影像配准方法。该方法首先使用多尺度Ｈａｒｒｉｓ提取出不受高斯平滑影响、住置稳定的角点，而后以金字塔分层映射为搜索策略，在参考影像上预测可能包含同名点的子区域，通过计算该子区域与待匹配点领域的相位一致性，同时引入相关系数作为相似性度量获取同名点对，最终实现不同光谱影像间的配准。实验表明，谊方法的误匹配率较低。精
度和稳定性高于传统原始影像灰度相关及同类配准方法，适用于不同传感器或不同光谱通道间影像的匹配。
关键词：影像配准：多尺度卷积；相位一致性；相关系数
中图分类号：ＴＰ３９１
文献标识码：Ａ
文章编号：１００４—２９５４（２０１２）０２一０１１８—０６
ＡＮｅｗ
Ｍｕｌｔｉ－ｓｏｕｒｃｅＲｅｍｏｔｅＳｅｎｓｉｎｇＩｍａｇｅＲｅｇｉｓｔｒａｔｉｏｎＭｅｔｈｏｄ
Ｂａｓｅｄ
ｏｎ
ＣｏｒｒｅｌａｔｉｏｎｏｆＰｈａｓｅＣｏｎｇｒｕｅｎｃｙ
ＦａｎＤｅｎｇｋｅ，ＰａｎＬｉ，ＹｅＹｕａｎｘｉｎ
（Ｓｃｈｏｏｌ
ｏｆＲｅｍｏｔｅ
ＳｅｎｓｉｎｇａｎｄＩｎｆｏｒｍａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ，ＷｕｈａｎＵｎｉｖｅｒｓｉｔｙ，Ｗｕｈａｎ４３００７９）
Ａｂｓｔｒａｃｔ：Ｔｈｉｓｐａｐｅｒ
ｐｒｏｖｉｄｅｓ
ａｎｅｗ
ｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ
ｍｅｔｈｏｄ
ｔｉｅ
ｂａｓｅｄ
ｏｎ
ｔｈｅ
ｃｏｒｒｅｌａｔｉｏｎｏｆｐｈａｓｅ
ｃｏｎｇｒｕｅｎｃｙｔｏｓｏｌｖｅｔｈｅｐｒｏｂｌｅｍｔｈａｔｔｈｅｉｍａｇｅｓｉｓｒａｔｈｅｒｌｏｗｄｕｒｉｎｇｇｅｏｍｅｔｒｉｃｗｉｔｈｏｕｔ
ｍａｔｃｈｉｎｇ
ｒａｔｉｏｏｆｔｈｅ
ｐｏｉｎｔｓｏｆｍｕｌｔｉ－ｓｏｕｒｃｅｒｅｍｏｔｅｓｅｎｓｉｎｇ
ｃｏｒｎｅｒｓ
ｃｏｒｒｅｃｔｉｏｎｃｏｕｒｓｅ．Ｉｎｔｈｉｓ
ａｒｅ
ｍｅｔｈｏｄ，ｔｈｅ
ｏｆｓｔａｂｌｅｐｏｓｉｔｉｏｎｓ
ｉｎｆｌｕｅｎｃｅｏｆ
ｓｔｒｕｃｔｕｒｅ
Ｇａｕｓｓｉａｎ
ａｓ
ｓｍｏｏｔｈｉｎｇｆｉｒｓｔｌｙ
ｆｉｌｔｅｒｅｄ
ｂｙｍｕｌｔｉ—ｓｃａｌｅＨａｒｒｉｓ，ａｎｄｔｈｅｎ
ｌａｙｅｒｅｄ
ｐｙｒａｍｉｄ
ｔｈｅｔｉｅ
ｉｓｕｔｉｌｉｚｅｄ
ｔｈｅｓｅａｒｃｈｉｎｇｓｔｒａｔｅｇｙｔｏｐｒｅｄｉｃｔｔｈｅｓｕｂｒｅｇｉｏｎｏｆｒｅｆｅｒｅｎｃｅｉｍａｇｅｉｎｗｈｉｃｈ
ｐｏｉｎｔｍａｙｂｅｌｏｃａｔｅｓ．Ｓｏ，ｂｙｃａｌｃｕｌａｔｉｎｇｔｈｅｐｈａｓｅｃｏｎｇｒｕｅｎｃｙｂｅｔｗｅｅｎｔｈｉｓｓｕｂｒｅｇｉｏｎａｎｄｔｈａｔｏｆ
ｐｏｉｎｔｔｏｂｅｍａｔｃｈｅｄ，ａｎｄｂｙｉｎｔｒｏｄｕｃｉｎｇｔｈｅｃｏｒｒｅｌａｔｉｏｎ
ｃｏｅｆｆｉｃｉｅｎｔ
ａｓ
ｔｈｅｓｉｍｉｌａｒｉｔｙｍｅａｓｕｒｅｍｅｎｔ
ｃａｎ
ｔｏ
ｏｂｔａｉｎ
ｔｈｅｔｉｅ
ｐｏｉｎｔ—ｐａｉｒｓ，ｔｈｅｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ
ｂｅｔｗｅｅｎ
ｄｉｆｆｅｒｅｎｔｓｐｅｃｔｒｕｍｓ
ｂｅ
ａｃｈｉｅｖｅｄ．Ｅｘｐｅｒｉｍｅｎｔ
ｃｏｎｔｒａｓｔ
ｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｉｓｍｅｔｈｏｄｇｏｔｌｏｗｍｉｓｍａｔｃｈｉｎｇｒａｔｉｏｂｕｔｈｉｇｈａｃｃｕｒａｃｙａｎｄｈｉｇｈｓｔａｂｉｌｉｔｙｉｎ
ｔｏ
ｔｈｅｔｒａｄｉｔｉｏｎａｌｇｒａｙｓｅａｌｅｃｏｒｒｅｌａｔｉｏｎｍｅｔｈｏｄｏｆｏｒｉｇｉｎａｌｉｍａｇｅａｎｄｏｔｈｅｒｓｉｍｉｌａｒｍｅｔｈｏｄｓ，ａｎｄｃｏｕｌｄｂｅ
ｏｒ
ｓｕｉｔａｂｌｅｆｏｒｔｈｅｉｍａｇｅｍａｔｃｈｉｎｇｂｅｔｗｅｅｎｄｉｆｆｅｒｅｎｔｓｅｎｓｏｒｓｄｉｆｆｅｒｅｎｔｓｐｅｃｔｒａｌｃｈａｎｎｅｌｓ．
Ｋｅｙｗｏｒｄｓ：ｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ；ｍｕｌｔｉ—ｓｃａｌｅｃｏｎｖｏｌｕｔｉｏｎ；ｐｈａｓｅｃｏｎｇｒｕｅｎｃｙ；ｃｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔ
１
概述
图像配准作为一项基础性的影像几何处理技术，
行性和研究价值的配准策略和方法，也因此产生了多种将配准方法进行分类的方式，包括基于配准影像数据源类型的分类、基于影像间几何畸变类型的分类等。最常见的是基于匹配时所使用影像信息类型的分类方式…，该方式将配准分为特征匹配和区域匹配两类。
特征匹配是当前配准技术研究的热点方向，它通过计算影像的梯度、矩等强度或对象的几何形态信息来描述局部范围内的显著性特征，然后选择适用于该特征匹配的相似性测度建立同名点的对应关系，最终完成对几何畸变的改正。常见的用于遥感影像配准的点特征有：Ｈａｒｒｉｓ、Ｆｏｒｓｔｅｒｎｅｒ、Ｍｏｒｏｖａｃ。线形或边缘特征有：Ｃａｎｎｙ、Ｓｏｂｅｌ，面状或区域特征有Ｈｕ不变矩、形状上下文。ＳＩＦＴ¨。是近几年图像配准研究领域兴起的一种具有尺度和旋转不变性描述的特征算子，它以其多
在医学、计算机视觉等学科领域中发挥着重要作用。尽管它在遥感中的研究和应用起步较晚，但是由于遥感影像在采集环境、观测角、地形起伏等多方面因素作用下所体现出的复杂性远远大于其他图像，再加上其本身海量及多源异构的特性，对配准技术在该领域的研究和应用提出新的挑战。目前，国内外众多学者针对遥感影像匹配问题的不同方面，提出了很多具有可
收稿日期：２０１１一ｌＯ一２０
基金项目：国家９７３项目支持。高分辨宰遥感数据精处理和空间信息智
能转化的理论与方法１（２０１２ＣＢ７１９９０５）
作者简介：范登科（１９¨一）．男．博士研究生，Ｅ．ｍｍｌ：ｆｄｋｈｂ＠１６３．㈣。
１１８
铁通话准设计ＲＡ／ＬＷＡｌｒ
ＳＴＡＮＤＡＲＤＤＥＳＩＧＮ
２０１２（２）
范登科．潘励，叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
・其
７
他・
尺度、稳定的点特征表达特性赢得许多研究学者的青睐．并由此衍生出许多基于ＳＩＦＴ的配准方法．如ＰＣＡ—ＳＩＦＴ［３］、ＣＳｌＦＴ”１、ＳＲ—ＳＩＦＴ“等。然而，由于大多数特征算子都是以图像梯度的强度和方向为基础的，且很少顾及区域性信息，造成多源异构遥感影像配准的同名点匹配率和纠正精度不能满足实际应用生产的需要，因此特征匹配方法具有应＿ｌ书局限性和问题针对性。
区域匹配作为一种早期发展成熟的、方便易用的图像配准方法，被广泛应用于各类商业化图像处理软件中。它通过统计运算的方式对特定区域内的图像信息进行描述，并借助相关关系模型度量彼此问的相似程度，相对于特征匹配需要稳定且充足的特征作为配准前提条件来说，区域匹配具有很强的通用性。目前研究和应用最多的区域匹配方法有相关系数和互信息，它们分别是基于影像原始灰度和信息熵进行计算和匹配的。Ｇａｏ使用相关系数作为匹配策略”１，借助影像金字塔和二次多项式实现对ＬａｎｄｓａｔＧｅｏＣｏｖｅｒ数据子像素精度的自动配准；Ｐｅｔｅｒ”ｏ同样以相关系数为核心匹配手段．完成对ＨｙＭａｐ、ＡＩＲＳＡＲ等２４组同源和非同源影像对的几何校正；Ｓａｈｉｌ”１通过迭代计算ＴｅｒｒａＳＡＲ—Ｘ和ｌｋｏｎｏｓ影像间的互信息，研究分析了城镇区域异构图像配准的精度。由于非同源异构影像问在原始灰度上通常是非线性变化的，且信息熵的计算量比较大，导致区域匹配方法的推广应用受到一定的制约。
本文综合特征匹配和区域匹配的优势，提出一种适用于非同源异构遥感影像的配准方法，该方法以不同滤波尺度下位置稳定不变的角点作为控制点．借助金字塔分层映射搜索策略．通过计算领域和搜索区域内的相位一致性，引入相关系数作为相似性测度实现同名点的匹配。并最终给出利用可见光、近红外及ＳＡＲ影像对该方法进行配准实验的验证和分析。２配准方法２．１控制点提取
为了保证提取足够多且位置相对稳定的控制点进行匹配和模型参数估计。需要使用通用性强、可靠性高的角点特征提取算子在待配准的输入影像上执行计算和检测。当前，Ｈａｒｒｉｓ”１和ＦＡＳＴ０１０１被证明具有较高的角点提取重复率，适用于各类光谱组成结构的影像。本文选择Ｈａｒｒｉｓ算子在输入影像上提取特征控制点，Ｈａｒｒｉｓ通过局部图像信息的自相关方式，在各个方向都计算了灰度变化．其角点提取的原理是：首先使用标准偏离差为口的高斯滤波模板与原始灰度影像１进行卷积运算。得到平滑后的影像７，然后计算７的一阶梯度矩阵Ｍ．公式为
妖道耘准设计ＲＡ儿ＷＡｙ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
Ｍ：『玖”）
【，，（ｘ，Ｙ），，（ｘ，Ｙ）
ｒ（Ⅵ）７，（埘’１（１）
Ｊ
咏ｚ，Ｙ）
其中，。（ｚ．Ｙ）和，。（ｚ．Ｙ）分镕Ⅱ代表平滑影像，中点（ｚ，Ｙ）处沿ｚ和Ｙ方向的梯度，通过解算矩阵肼的２个特征值Ａ．和Ａ：，使用判别函数Ｒ得到对于点（，，Ｙ）的特征判断，Ｒ的计算公式为
Ｒ＝＾，Ａ，一ｋ（＾．＋Ａ２）４
（２）
其中ｋ为介于０．０４—０．０６的常数，如果２个特征值都比较大，说明在该点处沿任意方向运动会引起判别函数Ｒ的急剧变化，判断为具有强响应的角点，否则为弱边缘或平滑域点．在使用阈值ｒ对该值进行局部非最大抑制后，可以获得稳定的特征点。
由于使用不同标准偏差口的高斯卷积模板会使平滑处理后的影像，变模糊．在特定滤波器尺度下检测的角点可能会在其他尺度的滤波影像中发生位置偏移甚至消失。为了提取不受高斯滤波尺度影响的稳定角点。需要使用临近尺度的高斯模板重新提取角点进行检测和判断。假设角点（‰，Ｙ。）是使用尺度（标准偏差值）为∥。的高斯滤波模板所提取的Ｈａｒｒｉｓ特征点，构造尺度为盯ｌ＝矿ｏ－Ｓｘｃｒｏ和盯２＝盯ｏ＋占ｘ盯。的高斯滤波模板对原始影像，进行卷积运算，提取的角点分别为（ｇ。。Ｙ．）和（ｘ：，Ｙ：），当提取的角点满足条件
Ｄ（（ｘｏ，Ｙｏ），（ｘｌ，ＹＩ））≤２＆ｒ
（３）
并且
Ｄ（（ｘｏ．Ｙｏ），（ｘ２，Ｙ２））≤２＆ｒ
（４）
时，认为（‰，Ｙ。）是一个不受卷积尺度影响的稳定角点。选为控制点，否则被排除。公式（３）、（４）中Ｄ代表欧式距离。即以角点（ｘ．，Ｙ．）和（ｚ：，Ｙ：）相对于（ｘ。，Ｙ。）的位置偏移作为特征点的判定依据。图１展示了分别使用原始Ｈａｒｒｉｓ算子和多尺度高斯滤波后Ｈａｒｒｉｓ提取的特征点分布情况，可以看出，在加人了公式（３）、（４）的判别后，位于道路边缘和开阔地中心的点受平滑尺度变化的影响，位置变化显著，被判断为不稳定的特征点排除．保留的点为稳定的不受平滑滤波干扰的特
征点。
２．２同名点匹配
由于Ｈａｒｒｉｓ特征算子提取的控制点为梯度变化显
咏准帕整仉为Ｉ５的商扣
随扳滤波后攫取结果
囤ｔ
６ｘ嘞为ｏ．５的临近Ｋ堰高抑髓扳滤波聍筛选挺墩络垠
Ｈａｒｒｉｓ算子提取特征点分布示意
２０１２（黝
・其他・
范登科，潘励，叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
著的角点，因此在其领域内影像的纹理信息比较丰富，
灰度变化显著，相关系数在小范围内容易形成变化剧烈
∑∑Ｗｓ（ｘ，ｙ）ＬＡ。（ｚ，ｙ）Ａ４，。（ｘ，ｙ）一ＬＪ
生生生——ｉ—－————————一（５）口
的极值，可以选择作为同名点的匹配测度。然而，非同Ｒ
、
源异构遥感影像间的光谱组成存在明显差异，在灰度强㈨０＝１
∑∑Ａ¨（ｚ，ｙ）＋８
度上呈现非线性变化，导致相关系数不能稳定描述相似式中，ＰＣ（ｚ，Ｙ）为点（ｚ，Ｙ）处的相位一致性，０为滤波度高的同名区域。针对这一问题，本文引入相位一致性方向，在区间［０，１ｒ］之间以固定步长变化，ｎ为滤波器
对影像进行转换，得到相关性程度高的特征图像。
尺度ｓ的最大值，矾（ｘ，Ｙ）为频率扩展的权重系数，
相位一致性是一种使用傅里叶谐波分量描述信号Ａ。（ｘ，Ｙ）代表使用沿日方向尺度为ｓ的ＬｏｇＧａｂｏｒ滤局部强度的特征，它利用特定方向口不同尺度ｓ的Ｌｏｇ波器进行分解后谐波分量的振幅值，ＬＪ表示仅取正值
Ｃａｂｏｔ滤波器，将信号分解为频率域下与日和ｓ相对应的运算，负值全部赋为０，Ｔ为噪声阈值，ｓ是一个避免
的傅里叶谐波分量．通过加权叠合不同尺度的分量，获除零的常数，ａ４，。（ｚ，ｙ）可由分解后谐波分量的相位得具有原始信号强弱特性描述的响应。Ｍｏｒｒｏｎｅ证明角计算获得，公式为
了相位一致性在描述图像特征信息上与人类神经视觉Ａ４，¨（ｚ，ｙ）＝ＣＯＳ（咖Ⅲ（ｚ，Ｙ）一币（ｘ，Ｙ））一
具有相一致的敏感度…ｏ。Ｗｏｎｇ指出。相位一致性特Ｉｓｉｎ（币。（ｚ，Ｙ）一咖（＾，Ｙ））Ｉ
（６）
征具有抵抗图像对比度和亮度变化的能力，并利用其其中咖（ｘ，ｙ）代表沿口方向尺度ｓ从１到ｎ之间的所构造的最小和最大距特征，分别实现了遥感影像特征有相位小的加权平均值。使用具有不同光谱组成结点的提取和同名点的匹配”…。应用于图像二维离散构、不同传感方式采集获得到的同名区域影像进行实信号，相位一致性的计算公式为
验分析。分别计算其相位一致性得到的结果如图２ＰＣ（ｘ．Ｙ）＝
所示。
（劬第３授段反射率彻第３渡段相位一致性
（ｃ）第５渡段反射率（ｄ）第５渡段相位一致性
囤２
Ｌａｎｄｓａｔ５
ＴＭ影像殛其相位一致性特征围像
从图２可以看出，边缘和角点处的相位一致性明分辨率，缩小图幅范围。首先，在输入影像的顶层金字显大于平坦开阔地区，通常在角点处取得极值。同时．塔图像上提取特征点．并搜索参考影像对应尺度的金相位一致性特征影像反映的并不是二值化的边缘，而字塔，找出相关系数大的窗口中心作为匹配点．使用简代表的是不同尺度和方向频率域卷积运算下傅里叶谐单的二维平面转换关系模拟初始的几何变形，根据最波分量的累积值。因此，以特征点为中心的窗口发生小二乘准则删除残差大的误匹配点对；而后。利用不同位置改变时，不会引起相关系数的剧烈变化，且由于角尺度金字塔图像的映射关系，在下一层更高分辨率尺点处相位一致性的显著特征，使得极值比较容易获得。度的金字塔图像上获取同名区域，执行与上述一致的经计算对应图像的相关系数，图２中（ａ）与（ｃ）之间的精确匹配，直到达到底层即原始分辨率影像为止。两相关系数为一０．０４５，而（ｂ）与（ｄ）达到０．６０５，可见采层金字塔结构的同名区域映射关系及几何变换模型估用相位一致性特征度量影像间的相关性明显优于原始计如图３所示。
灰度。
２．３模型参数估计与重采样
遥感影像在采集过程中一般表现出大视角、宽视利用匹配获取的同名点对．依据最小二乘原理估场的特性，成像的图幅范围大，使全局匹配的效率降计几何畸变模型的参数，剔除残差大的同名点对以保低，计算复杂度增大，需要使用特定的搜索策略优化同证模型估计精度在可接受的范围内，最终使用特定的名点匹配过程。基于金字塔结构的同名区域预测与分插值方法实现对输入影像的重采样处理，输出在地理层索引映射是目前最受欢迎的、应用于遥感影像匹配位置上与参照影像配准的校正影像。在这一过程中，中的搜索策略．方法是以一定比例逐级降低图像空间
全局性几何畸变模型的选择显得十分重要。
１２０
铁道嚣准设计ＲＡＩＬＷＡＹＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
２０１２（２）
范登科．潘励．叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
．其
他・
最大相关的窗口，以其中心作为匹配同名点；（４）＿｝｝Ｊ简单二维平面变换模型估计畸变改正关系，预测下一级
，
金字塔图像的同名区域；（５）在下级金字塔图像的子区域内迭代执行步骤（２）一（４），直到完成原始分辨率影像上同名点的匹配；（６）根据待配准影像的类型确定全局几何变换模型，采用最小二乘原理迭代解算和
圈３基于金字塔结构的遥豫影像同名区域索引和搜索策略
修正模型精度；（７）重采样输出配准后的结果影像。图４展示了配准策略的完整流程。
由于传感器姿态和侧视角的差异，采集获得的多源异构遥感影像在全局几何畸变上表现出不同形式。中分辨率光学影像的成像方式为推扫式，其传感器所处轨道高度及线阵扫描方式决定了影像几乎不受侧视角和地形起伏的影响．在全局影像上仅存在平移、旋转和尺度畸变，因此可以使用仿射模型描述待配准光学影像问的几何变形关系。仿射变换模型公式为
ｆ３１：＾『００８１『。１＋［ＣｌＬｙＪＬｓｉｎ０一ｃ０８０ＪＬ训
Ｌ
８８１“８
（７）
图４配准方法流程
ｒＪ
式中，（ｚ，ｙ）和（Ｘ，Ｙ）分别为对应同名点在输入影像和参考影像的图像坐标，＾代表尺度缩放因子，０表示影像的相对旋转角，（ｃ，ｒ）是图像在二维平面的相对平移量。因此仅需要两对同名点即可以求解该模型参数。而同分辨率ＳＡＲ影像的采集使用主动传感方式，虽然轨道高度避免了地形起伏的影响，但是其成像图幅范围较小．并且存在一定的侧视角，因此对于光学影像与ＳＡＲ影像以及不同时相ＳＡＲ影像问的几何畸变，通常采用投影变形模型来描述，一次项投影模型公式为
ｇｘ＋ｈｙ＋１
ｇｘ＋ｈｙ＋ｌ
３实验与分析
３．１实验背景和数据
实验以ＷｉｎｄｏｗｓＸＰ系统为平台，在Ｖｉｓｕａｌ法．程序运行的硬件环境为Ｉｎｔｅｒ
２
Ｃ＋＋
２００８集成开发环境中设计开发程序实现上述配准方
Ｃｏｒｅ２ＣＰＵ２．６６ＧＨｚ，
ＧＢ内存。实验数据选用３组遥感影像，每一组中输
ｘ：旦生譬２掣．ｙ：坐±丝（８）
入影像和参考影像分别采集自异掏光谱范围或不同传感器。前２组数据均由光学遥感影像组成，采集的区域为湖北省大东湖和梁子湖水域．其中第ｌ组数据来自于同一传感器，为了体现空间位置上的差异性，手工对输入影像旋转ｌｏ。，第２组数据来自不同传感器，存在微小旋转及１０～２０像素的平移变形。第３组数据为光学影像与ＳＡＲ影像之间的配准，为了正确计算相关系数．对ＳＡＲ影像进行重采样．确保其与光学影像的分辨率一致。３组实验数据的详细信息如表ｌ
所示。
式中包括。一＾总计８个模型参数，需要４对同名点进行求解，分母部分描述了成像平面的倾斜变形。在最小二乘原理基础上对模型参数进行估计后，使用参考影像中的控制点坐标反算输人影像控制点的模型坐标，并与真实坐标求差后计算均方根误差ＲＭＳＥ，ＲＭＳＥ值的大小代表了该同名点对偏离模型估计的程度．删除ＲＭＳＥ值最大的多个点对，重新估计模型参数，迭代执行这一过程直到同名点对集合中最大的ＲＭＳＥ值小于预先设定的阈值。以获得满足纠正精度要求的模型参数。最终，双线性插值方法被用于重采样过程输出配准结果。２．４总体配准策赔
本文提出的基于相位一致性相关的多源遥感影像配准方法可概括为以下７个步骤：（１）根据影像大小。设定特定的分辨率比例，分别构建输人和参考影像的多级金字塔图像；（２）在输入影像的顶层金字塔上使用多尺度Ｈａｒｒｉｓ提取控制点．并计算其领域内窗口的相位一致性；（３）计算参照影像顶层金字塔的相位一致性特征图像，使用相关系数作为匹配测度全局搜索
帙遵标准艟汁ＲＡ儿ＷＡＹ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
在进行影像配准运算处理过程中．需要对相关参数进行初始化．具体项目及其初始值为：金字塔各级比例Ｉ：３，金字塔级数３，原始分辨率下搜索区域大小２００像素，相关系数阈值０．７５．相关系数计算窗口大小１３像素．ＲＭＳＥ阈值Ｉ．５。
３．２配准精度分析
使用以上３组数据对本文配准方法进行实验验证，并对该方法的适用性和精度进行评价。表２列出了经实验获得的配准评价结果。
由表２可以看出．在配准精度方面，第３组数据的效果最好，最大残整和总体均方根误差指标都控制在
１２１
２０１２（引
・算他・
范登科，潘励，叶沅蠢一一种基于相位一致性相关的多源遥感影像配准方法
表１实验数据信息概览
ｌＬａｎｄｓａ＊５ＴＭ
ｂａｎｄ３／０．６３－０
６９
ｐｍ
３０ｍ＋ｌＯ。３５００ｘ３５００２００７—０４—１９
Ｌａｎｄｓａｔ５ＴＭｂａｎｄ５／１５５—１７５ｌ＾ｍ３０ｍ３０００ｘ３ＯＯＯ２００７—０４—１９２ＨＪ—ｌＡＣＣＤＩｂ衄ｄ４／０．７６—０．９０“ｍ
３０ｍ３２５０ｘ３２５０２０１１一Ｏｌ一１６Ｌａｎｄｓａｔ５ＴＭｂａｎｄ５／１５５—１
７５＊ｍ
３０ｍ３０００ｘ３０００２００７－０４一１９３
Ｌａｎｄｓａｔ５ＴＭ
ｂａｎｄ
１２３／可见光
３０ｍ
Ｉ７００ｘ１６５０２００６—０９一１０
ＡＬＯＳｐａｉｓａ
ｂａｎｄＩ／ＳＡＲ
３０ｍ
１８００ｘ１７００２００５—１１—１０
表２配准方法的实验评价结果
范围内，且在对比度高的局部区域内相对集中。符合相
位一致性描述亮度和对比度差异的不变特征。
综合分析３组实验数据的配准结果，本文方法更适用于对非同源异构的两景光学影像进行配准。对于光学与ＳＡＲ影像的配准有一定局限性和区域选择性，
１个像素以内，然而同名点对的误匹配率却达到了同时将相关系数作为匹配测度引入到该方法中。对于
９５％以上，仅有５对正确匹配的同名点用于模型参数全局仿射和投影变形有一定的抵抗能力。估计，且它们分布集中于图５（ｅ）、（ｆ）中左下角的岛屿３．３对比验证
海岸线上，占影像中主体的陆地范围内并没有正确匹为了进一步验证本文方法的配准效果．引入传统基配的同名点对，这是因为大陆中地物类型丰富而海岛于灰度的相关系数（ＮＣＣ）、尺度不变特征（ＳＩＦＴ’）、改进内单一．ＳＡＲ影像之于前者表现出更多的噪声．而相尺度不变特征（ＳＲ—ＳＩＦＴ）３种方法分别对３组实验数据位一致性的计算结果对噪声十分敏感．造成陆地中相进行配准实验，统计各项实验结果的误匹配率和均方根位一致性的相关程度并不是十分显著；反之，海岛在中误差（ＲＭＳＥ），得到的对比情况如表３所示。
ＳＡＲ影像中纹理结构简单，海岸线轮廓特征明显，有表３
３种常用方法实验结果对比（“一”表示配准失败）
利于相位一致性稳定描述局部图像信息。
配准
实验１
实验２实验３
方法
误匹配率／％ＲＭＳＥ谩匹配率／％ＲＭＳＥ误匹配卓／％ＲＭＳＥ
ＮＣＣ
５８
３８
０８７１４４．１７０．９４３圈鬻糕一一ＳＩｒｒ６６７３０９９３５５．９７０９８４
一
一ＳＲ—ＳＩＦＴ
５５
３６
０９７７
４１．８４
０９９６——
一
实验过程中设定ＲＭＳＥ阈值为１．０像素，以保证几何校正模型参数估计的准确性和一致性。当３种方法达到相同精度水平时，结果显示实验２的误匹配率最低，原鬻翟甄
因是两景近红外影像间灰度相关性程度最高．且与特征表述相关的梯度强度和方向变化较为稳定．与此同时，ＮＣＣ和ＳＲ—ＳＩＦＴ方法的误匹配率明显低于本文方法。相反，３种方法对实验ｌ数据进行配准的误匹配率高于或接近本文方法，可见，对于光谱组成差异越显著的影像来｛¨Ｊ？’｝＿＿｛；＿ｖ
－。‘１ｒ。、，“…
＿・－
“≮
说，本文方法的适用性更强，效果更好。３种方法不能完图５
３组实验藏据影像及同名点分布情况
成对实验３数据配准的事实更进一步验证了这一结论。前２组实验数据的配准结果在精度和误匹配率上因此，与３种配准方法进行对比分析得出的结论是，基于达到～致水平。尽管第２组数据的总体均方根中误差相位一致性相关的配准方法可以应对传统方法所不能解略微偏高，但其匹配得到的同名点对数目是第ｌ组的决的异构遥感影像间难以配准的问题。３倍，是因为第２组数据的影像对均由近红外波段构４结语
成，相对于第ｌ组的可见光一近红外组合，相关程度更高。然而，在匹配获得大量同名点对的同时．发生误匹本文提出一种基于相位一致性相关的多源遥感影
配的点对数目也相应提升，且从２组数据对比情况来像配准方法，该方法以多级图像金字塔结构为索引和
看，这一比例水平并不受光谱组成结构的影响．因此要搜索策略，首先使用多尺度高斯滤波判别和遘选稳定迸一步降低误匹配率，仍需要对匹配方法和策略进行的Ｈａｒｒｉｓ特征点．在控制点领域和相关系数极值搜索改进。另一方面，由于第２组影像采集自不同时问．地区域内，引人相位一致性特征进行相似性度量，在降低物变化也会造成同名点的误匹配。从图５（ａ）～（ｄ）误匹配概率的同时获得更多精确匹配的同名点对。经
中可以看出，正确匹配的同名点均匀分布于整个图幅
实验分析和对比验证．该方法适用于中分辨率多源异
１２２
铁道捧准设计ＲＡＩＬＷＡＹ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
２Ｄｊ２（２）
范登科，潘励，叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
．其他・
构光学一光学、光学一ＳＡＲ遥感影像间的配准，精度满ｂａｓｅｄ
ｏｎ
ＳＩＦＴ［Ｊ］ＥｌｅｃｔｒｏｎｉｃｓＬｅｔｔｅｎ．２００８．４４（２）：１０７一１０８
足实际生产过程中影像几何校正处理的需要，且高于［６］ＧａｏＦ．Ｍａ蚰ｋ
Ｊ，Ｗｏｌｆｅ
Ｒ
Ｅ
Ａｕｔｏｍａｔｅｄｒｅｇｉｓｔｒａｔｉｏｎ
ａｍｌｏｒｅｈｏｎｃｔｉｆｉｃａｔｉｏｎ
Ｌａｎｄｓａｔａｎｄ
Ｌａｎｄｓａｌ—ｌｉｋｅ
ｄａｔａ
其他基于区域匹配的配准方法。下一步的研究工作将ｐａｃｋａｇｅ
ｆｏｒ
ｐｒｏｃｅｓｓｉｎｇ【Ｊ］．Ｊｏｕｒｎａｌ
ｏｆＡｐｐｌｉｅｄＲｅｍｏｔｅ
Ｓｅｎｓｉｎｇ．２００９．３（１）
专注于提高运算效率，改进算法以降低误匹配率，并适［７］ＢｕｎｔｉｎｇＰ．Ｌ“ｍ∞Ｆ．Ｌｕｃａｓ
Ｒ
Ａ
ｍｕｌｔｉ－ｒｅｓｏｌｕｔｉｏｎ㈣‘ｂａｓｅｄ
用于不同分辨率、地形起伏等复杂条件下的多源遥感ｔｅｃｈｎｉｑｕｅｆｏｒａｕｔｏｍａｔｉｃｍｕｌｔｉ－ｍｏｄａｌｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ［Ｊ］Ｉｍａｇｅ
ａｎｄ
影像配准。同时希望本文方法能为图像配准领域的学Ｖｉｓｉｏｎ
Ｃｏｍｐｕｔｉｎｇ．２０１０．２８（８）：１２０３—１２１９．
者提供研究基础和借鉴。［８］Ｓｕｒｉ
Ｓ．Ｒｅｉｎａｒｔｚ
Ｐ．Ｍｕｔｕａｌ－ｉｎｆｅ．ｎａｔｉｏｎ－ｂｍｄ
ｒｅｇｉｓｔｒａｔｉｏｎ
ｏｆ
参考文献：
ＴｅｎａＳＡＲ．Ｘａｎｄ
Ｉｋｏｎｏｓ
ｉｍａｇｅｒｙｉｎ
ｕｒｂａｎ
ｌｉｔｅｒ［Ｊ］．ＩＥＥＥ
ＣＩ】ｚｉｔｏｖ＾Ｂ．Ｆｌ㈣Ｊ
Ｔｒａｎｓａｃｔｉｏｎｓ
ｏｎ
Ｇｅｏｓｃｉｅｎｃｅ
ａｎｄ
Ｒｅｍｏｔｅ
Ｓｅｎｓｉｎｇ．２０１０．４８（２）：
Ｉｍａｇｅ”ｇｌｓｘｒｍｉｏｎ
ｍｅｔｈＭ‘：Ａｓｕｒｖｅｙ【Ｊ］｜ｍａｒｅ
９３９—９４９
ａｎｄＶｉｓｌｏｎ
Ｃｏｍｐｕｔｉｎｇ．２００３．２１（Ｉ１）：９７７—１０００
［９］Ｈａｒｒｉｓ
ｃ．ＳｔｅｐｈｅｎｓＭＫ
Ａｃｏｍｂｉｎｅｄｃｏｍｅｒ
ｓｏｄ
ｅｄｇｅ
ｄｅｔｅｃｔｏｒ［Ｃ］．
［２］Ｌｏｗｓ
ＤＧ
Ｄｉｓｔｉｎｃｔｉｖｅｉｍａｇｅｆｅａｔｕ亿ｓ
ｆｒｏｍ
ｓｃａｌｅ－ｉｎｖａｒｉ肌ｔｋｅｙｐｎｉｎｔｓ
Ｐｒｏｃ
（Ｊ】ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌ
ＡｌｖｅｙＶｉｓｉｏｎＣｏｎｆ．．１９８８：１４７—１５２．ｏｆＣｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎ，２００４．６０（２）：９１一１１０
ＰＣＡ—ＳＩＦＴ：Ａ㈣ｄｉｓｔｉｎｃｔｉｖｅ
［１０］ＲｏｅｔｅｎＥ．ＤｒｕｍｍｏｎｄＴ
Ｍａｃｈｉｎｅ［３］Ｋｅ
ｌｅａｒｎｉｎｇ
ｆｏｒ
ｈｉｇｈ—ｓｐｅｅｄ
ｃｏｒｎｅｒ
Ｙ．ＳｕｋｔｈａｎｋａｒＲ
ｒｅｐ＂ｓｅｎｔａｔｉｏｎ
ｄｅｔｅｃｔｉｏｎ［Ｃ］ｐｒｏｃｅｅｄｉｎｇｓ
ｏｆｔｈｅＥｕｒｏｐｅａｎＣｏｎｆｅｒｅｎｃｅ
ｏｎ
Ｃｏｍｐｕｔｅｒ
ｆｏｒｌｏｃａｌｉｍａｇｅ
ｄｅ虻ｄ阱ｏｎ［Ｃ］Ｐｒｏｃｅｅｄｉｎｇｓ
ｏｆｔｈｅＩＥＥＥ
Ｃｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎ．２００６：４３０—４４３
ＳｏｃｉｅｔｙＣｏｎｆｅｒｅｎｃｅｏｎ
Ｃｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎ蛐ｄ
ＰａｔｔｅｍＲ∞ｏｇｎｉｔｉｏｎ．
［ＩＩ］Ｍｏ肿ｎｅ
Ｍｃ．ＯｗｅｎｓＲＡ．Ｆｅａｔｕｒｅｄｅｔｅｃｔｉｏｎｆｒｏｍｌｏｃａｌｅｎｅｒｇｙ［Ｊ］
２００４：５０６—５１３
［４］Ａｂｄｅｌ・Ｈａｋｉｍ
ＡＥ．ＦａｒａｇＡＡ
ＣＳＩＦＴ：Ａ
ＳＩＦＴ
ＰａｔｔｅｒｎＲｅｃｏｇｎｉｔｉｏｎＬｅｔｔｅｒｓ．１９８７．６（５）：３０３—３１３．
ｄｅ虻ｒｌｐｔｏｒｗｉｔｈｃｏｌｏｒ
［１２］Ｗｏｎｇ
Ａ．Ｃｌａｕｓｉ
Ｄ
Ａ
ＡＩＳＩＲ：Ａｕｔｏｍａｔｅｄｉｎｔｅｒ．柏ｎｍｒ／ｉｎｔｅｒ—ｂ蛆ｄ
ｉｎｖａｒｉａｎｔｃｈ盯ａｅｔｅｒｉｓｔｉｅｓ［Ｃ］Ｐｒｏｃｅｅｄｉｎｇｓ
ｏｆ
ｔｈｅ
ＩＥＥＥ
Ｃｏｍｐｕｔｅｒ
ｓａｔｅｌｌｉｔｅｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎｕｓｉｎｇｒｏｂｕｓｔｃｏｍｐｌｅｘｗａｖｅｌｅｔｆｅａｔｕｒｅｒｅｐｒｏ－
Ｓｏｃｉｅｔｙ
Ｃｏｎｆｅｒｅｎｃｅ
ｏｎ
Ｃｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎａｎｄ
ＰａｔｔｅｍＲｅｃｏｇｎｉｔｉｏｎ．
盹ｎｔａｔｉｏｎａ［Ｊ］Ｐａｔｔｅｒｎ
ＲｅｃｏｇｎｉｔｉｏｎＬｅｔｔｅｍ，２０１０．３１（１０）：１１６０一
２００６：｜９７９一１９８３
１１６７．
［５】Ｙｉｚ．ＺｈｉｇｕｏＣ，ＹａｎｇＸ
Ｍｕｌｔｉ—ｓｐｅｃｔｒａ］ｍｍｏｔｅｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ
（上接第１１２页）
翻车棚的端头设置翻车机监控室．便于操作人员对翻的取送车，棚的跨度需大于２条重调机走行轨道间距
车机作业时的监控。控制室内还设有配电室、夹轮器（２５．４１１１），建筑模数取整后，棚的跨度取２７
ｍ。
液压站等房同。翻车机卸车系统工艺平面如图９
翻车机控制室设置在２条重车线的中间，在靠近
所示。
１８０ＤＤ
２５００［
６ｌ瑚６０００６０００
６０００
６０∞
甲镌盎’谶迸芈方向
洲州
，●Ｉ————』
｜．．｛
Ｌ亡三Ｉｌ圭些
￡一位城料枷…。笠Ｌ—Ｌ丑墅磬
《习ｉ蕊。弘嗣矧
一．ｕ翌ｌ
＂
旺匠ｉ口：Ⅱ旺
口ｊＩＩｎ…————■一■．－……“
Ｚ『央轮器罗呷ｉ
晨
ｊ．｜孑｝—＝Ｆ—孔
｝’
’
。
。
’
’
’
’
’／。
’
。怍革≠±手ｒ＋一｝＝二ｊ州Ｉ仃Ｉ警≤必ｌ錾Ｉ…｝
一～一
塑塑。
…一■＿
川Ｌ
安全止Ｈ群
圈９翻车机卸车系统工艺平面布置（单位：ｍｍｌ
５结论
参考文献：
［Ｉ］张磊．马明札．燃料远行与检謦［Ｍ］．北京：中国电力出版
翻车机卸车工艺作为铁路货场一种新型的卸车工杜．２００６．
［２］高志刚．新“Ｃ”型折退式双车翻车机在定洲发电公司＝期工程的
艺，具有生产能力大、自动化程度高、对车辆损伤小和应用［Ｊ】神毕科技．２０１０，８（Ｉ）
环境污染小等特点，应推广翻车机在大宗散装物料铁［３］＊撬军．蛱路货车翻车机使用中存在的问题及解挟方击［Ｊ］铁道
路货场的应用。进行翻车机卸车工艺设计时，应优先车辆．２００９．４７（９）．
采用折返式布置；如采用２台翻车机，则应采用２台翻
［４】
张凡华鼻车机设鲁选捌分析［Ｊ］华电技术．２００８．３０（６）．
车机对称。相背的布置形式。
【５］石小东．幢路散货运■不同卸车方式比较［Ｊ］．内簟古煤炭烃脐．
２００８（３）
侠嬉檬准设计ＲＡＩＬＷＡＹ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
２０１２（２）
１２３
一种基于相位一致性相关的多源遥感影像配准方法
作者：作者单位：刊名：英文刊名：年，卷(期)：
范登科， 潘励， 叶沅鑫， Fan Dengke， Pan Li， Ye Yuanxin武汉大学遥感信息工程学院,武汉,430079铁道标准设计
Railway Standard Design2012(2)
1. Zitová B;Flusser J Image registration methods:A survey[外文期刊] 2003(11)2. Harris C;Stephens M K A combined corner and edge detector 1988
3. Suri S;Reinartz P Mutual-information-based registration of TerraSAR-X and Ikonos imagery in urban area 2010(02)4. Bunting P;Labrosse F;Lucas R A multi-resolution area-based technique for automatic multi-modal image registration2010(08)
5. Gao F;Masek J;Wolfe R E Automated registration and orthorectification package for Landsat and Landsat-like dataprocessing 2009(01)
6. Yi Z;Zhiguo C;Yang X Multi-spectral remote image registration based on SIFT[外文期刊] 2008(02)7. Abdel-Hakim A E;Farag A A CSIFT:A SIFT descriptor with color invariant characteristics 20068. Ke Y;Sukthankar R PCA-SIFT:A more distinctive representation for local image descriptor 20049. Lowe D G Distinctive image features from scale-invariant keypoints[外文期刊] 2004(02)
10. Wong A;Clausi D A AISIR:Automated inter-sensor/inter-band satellite image registration using robust complex waveletfeature representations 2010(10)
11. Morrone M C;Owens R A Feature detection from local energy 1987(05)12. Rosten E;Drummond T Machine learning for high-speed corner detection 2006
本文链接：http://d.g.wanfangdata.com.cn/Periodical_tdbzsj201202034.aspx
・其
他－
一种基于相位一致性相关的多源遥感影像配准方法
范登科，潘
励，叶沅鑫
（武汉大学遥感信息工程学院，武汉４３００７９）
摘要：针对几何校正过程中多源遥感影像同名点匹配率低的问题，提出一种基于相位一致性相关的遥感影像配准方法。该方法首先使用多尺度Ｈａｒｒｉｓ提取出不受高斯平滑影响、住置稳定的角点，而后以金字塔分层映射为搜索策略，在参考影像上预测可能包含同名点的子区域，通过计算该子区域与待匹配点领域的相位一致性，同时引入相关系数作为相似性度量获取同名点对，最终实现不同光谱影像间的配准。实验表明，谊方法的误匹配率较低。精
度和稳定性高于传统原始影像灰度相关及同类配准方法，适用于不同传感器或不同光谱通道间影像的匹配。
关键词：影像配准：多尺度卷积；相位一致性；相关系数
中图分类号：ＴＰ３９１
文献标识码：Ａ
文章编号：１００４—２９５４（２０１２）０２一０１１８—０６
ＡＮｅｗ
Ｍｕｌｔｉ－ｓｏｕｒｃｅＲｅｍｏｔｅＳｅｎｓｉｎｇＩｍａｇｅＲｅｇｉｓｔｒａｔｉｏｎＭｅｔｈｏｄ
Ｂａｓｅｄ
ｏｎ
ＣｏｒｒｅｌａｔｉｏｎｏｆＰｈａｓｅＣｏｎｇｒｕｅｎｃｙ
ＦａｎＤｅｎｇｋｅ，ＰａｎＬｉ，ＹｅＹｕａｎｘｉｎ
（Ｓｃｈｏｏｌ
ｏｆＲｅｍｏｔｅ
ＳｅｎｓｉｎｇａｎｄＩｎｆｏｒｍａｔｉｏｎＥｎｇｉｎｅｅｒｉｎｇ，ＷｕｈａｎＵｎｉｖｅｒｓｉｔｙ，Ｗｕｈａｎ４３００７９）
Ａｂｓｔｒａｃｔ：Ｔｈｉｓｐａｐｅｒ
ｐｒｏｖｉｄｅｓ
ａｎｅｗ
ｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ
ｍｅｔｈｏｄ
ｔｉｅ
ｂａｓｅｄ
ｏｎ
ｔｈｅ
ｃｏｒｒｅｌａｔｉｏｎｏｆｐｈａｓｅ
ｃｏｎｇｒｕｅｎｃｙｔｏｓｏｌｖｅｔｈｅｐｒｏｂｌｅｍｔｈａｔｔｈｅｉｍａｇｅｓｉｓｒａｔｈｅｒｌｏｗｄｕｒｉｎｇｇｅｏｍｅｔｒｉｃｗｉｔｈｏｕｔ
ｍａｔｃｈｉｎｇ
ｒａｔｉｏｏｆｔｈｅ
ｐｏｉｎｔｓｏｆｍｕｌｔｉ－ｓｏｕｒｃｅｒｅｍｏｔｅｓｅｎｓｉｎｇ
ｃｏｒｎｅｒｓ
ｃｏｒｒｅｃｔｉｏｎｃｏｕｒｓｅ．Ｉｎｔｈｉｓ
ａｒｅ
ｍｅｔｈｏｄ，ｔｈｅ
ｏｆｓｔａｂｌｅｐｏｓｉｔｉｏｎｓ
ｉｎｆｌｕｅｎｃｅｏｆ
ｓｔｒｕｃｔｕｒｅ
Ｇａｕｓｓｉａｎ
ａｓ
ｓｍｏｏｔｈｉｎｇｆｉｒｓｔｌｙ
ｆｉｌｔｅｒｅｄ
ｂｙｍｕｌｔｉ—ｓｃａｌｅＨａｒｒｉｓ，ａｎｄｔｈｅｎ
ｌａｙｅｒｅｄ
ｐｙｒａｍｉｄ
ｔｈｅｔｉｅ
ｉｓｕｔｉｌｉｚｅｄ
ｔｈｅｓｅａｒｃｈｉｎｇｓｔｒａｔｅｇｙｔｏｐｒｅｄｉｃｔｔｈｅｓｕｂｒｅｇｉｏｎｏｆｒｅｆｅｒｅｎｃｅｉｍａｇｅｉｎｗｈｉｃｈ
ｐｏｉｎｔｍａｙｂｅｌｏｃａｔｅｓ．Ｓｏ，ｂｙｃａｌｃｕｌａｔｉｎｇｔｈｅｐｈａｓｅｃｏｎｇｒｕｅｎｃｙｂｅｔｗｅｅｎｔｈｉｓｓｕｂｒｅｇｉｏｎａｎｄｔｈａｔｏｆ
ｐｏｉｎｔｔｏｂｅｍａｔｃｈｅｄ，ａｎｄｂｙｉｎｔｒｏｄｕｃｉｎｇｔｈｅｃｏｒｒｅｌａｔｉｏｎ
ｃｏｅｆｆｉｃｉｅｎｔ
ａｓ
ｔｈｅｓｉｍｉｌａｒｉｔｙｍｅａｓｕｒｅｍｅｎｔ
ｃａｎ
ｔｏ
ｏｂｔａｉｎ
ｔｈｅｔｉｅ
ｐｏｉｎｔ—ｐａｉｒｓ，ｔｈｅｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ
ｂｅｔｗｅｅｎ
ｄｉｆｆｅｒｅｎｔｓｐｅｃｔｒｕｍｓ
ｂｅ
ａｃｈｉｅｖｅｄ．Ｅｘｐｅｒｉｍｅｎｔ
ｃｏｎｔｒａｓｔ
ｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｔｈｉｓｍｅｔｈｏｄｇｏｔｌｏｗｍｉｓｍａｔｃｈｉｎｇｒａｔｉｏｂｕｔｈｉｇｈａｃｃｕｒａｃｙａｎｄｈｉｇｈｓｔａｂｉｌｉｔｙｉｎ
ｔｏ
ｔｈｅｔｒａｄｉｔｉｏｎａｌｇｒａｙｓｅａｌｅｃｏｒｒｅｌａｔｉｏｎｍｅｔｈｏｄｏｆｏｒｉｇｉｎａｌｉｍａｇｅａｎｄｏｔｈｅｒｓｉｍｉｌａｒｍｅｔｈｏｄｓ，ａｎｄｃｏｕｌｄｂｅ
ｏｒ
ｓｕｉｔａｂｌｅｆｏｒｔｈｅｉｍａｇｅｍａｔｃｈｉｎｇｂｅｔｗｅｅｎｄｉｆｆｅｒｅｎｔｓｅｎｓｏｒｓｄｉｆｆｅｒｅｎｔｓｐｅｃｔｒａｌｃｈａｎｎｅｌｓ．
Ｋｅｙｗｏｒｄｓ：ｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ；ｍｕｌｔｉ—ｓｃａｌｅｃｏｎｖｏｌｕｔｉｏｎ；ｐｈａｓｅｃｏｎｇｒｕｅｎｃｙ；ｃｏｒｒｅｌａｔｉｏｎｃｏｅｆｆｉｃｉｅｎｔ
１
概述
图像配准作为一项基础性的影像几何处理技术，
行性和研究价值的配准策略和方法，也因此产生了多种将配准方法进行分类的方式，包括基于配准影像数据源类型的分类、基于影像间几何畸变类型的分类等。最常见的是基于匹配时所使用影像信息类型的分类方式…，该方式将配准分为特征匹配和区域匹配两类。
特征匹配是当前配准技术研究的热点方向，它通过计算影像的梯度、矩等强度或对象的几何形态信息来描述局部范围内的显著性特征，然后选择适用于该特征匹配的相似性测度建立同名点的对应关系，最终完成对几何畸变的改正。常见的用于遥感影像配准的点特征有：Ｈａｒｒｉｓ、Ｆｏｒｓｔｅｒｎｅｒ、Ｍｏｒｏｖａｃ。线形或边缘特征有：Ｃａｎｎｙ、Ｓｏｂｅｌ，面状或区域特征有Ｈｕ不变矩、形状上下文。ＳＩＦＴ¨。是近几年图像配准研究领域兴起的一种具有尺度和旋转不变性描述的特征算子，它以其多
在医学、计算机视觉等学科领域中发挥着重要作用。尽管它在遥感中的研究和应用起步较晚，但是由于遥感影像在采集环境、观测角、地形起伏等多方面因素作用下所体现出的复杂性远远大于其他图像，再加上其本身海量及多源异构的特性，对配准技术在该领域的研究和应用提出新的挑战。目前，国内外众多学者针对遥感影像匹配问题的不同方面，提出了很多具有可
收稿日期：２０１１一ｌＯ一２０
基金项目：国家９７３项目支持。高分辨宰遥感数据精处理和空间信息智
能转化的理论与方法１（２０１２ＣＢ７１９９０５）
作者简介：范登科（１９¨一）．男．博士研究生，Ｅ．ｍｍｌ：ｆｄｋｈｂ＠１６３．㈣。
１１８
铁通话准设计ＲＡ／ＬＷＡｌｒ
ＳＴＡＮＤＡＲＤＤＥＳＩＧＮ
２０１２（２）
范登科．潘励，叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
・其
７
他・
尺度、稳定的点特征表达特性赢得许多研究学者的青睐．并由此衍生出许多基于ＳＩＦＴ的配准方法．如ＰＣＡ—ＳＩＦＴ［３］、ＣＳｌＦＴ”１、ＳＲ—ＳＩＦＴ“等。然而，由于大多数特征算子都是以图像梯度的强度和方向为基础的，且很少顾及区域性信息，造成多源异构遥感影像配准的同名点匹配率和纠正精度不能满足实际应用生产的需要，因此特征匹配方法具有应＿ｌ书局限性和问题针对性。
区域匹配作为一种早期发展成熟的、方便易用的图像配准方法，被广泛应用于各类商业化图像处理软件中。它通过统计运算的方式对特定区域内的图像信息进行描述，并借助相关关系模型度量彼此问的相似程度，相对于特征匹配需要稳定且充足的特征作为配准前提条件来说，区域匹配具有很强的通用性。目前研究和应用最多的区域匹配方法有相关系数和互信息，它们分别是基于影像原始灰度和信息熵进行计算和匹配的。Ｇａｏ使用相关系数作为匹配策略”１，借助影像金字塔和二次多项式实现对ＬａｎｄｓａｔＧｅｏＣｏｖｅｒ数据子像素精度的自动配准；Ｐｅｔｅｒ”ｏ同样以相关系数为核心匹配手段．完成对ＨｙＭａｐ、ＡＩＲＳＡＲ等２４组同源和非同源影像对的几何校正；Ｓａｈｉｌ”１通过迭代计算ＴｅｒｒａＳＡＲ—Ｘ和ｌｋｏｎｏｓ影像间的互信息，研究分析了城镇区域异构图像配准的精度。由于非同源异构影像问在原始灰度上通常是非线性变化的，且信息熵的计算量比较大，导致区域匹配方法的推广应用受到一定的制约。
本文综合特征匹配和区域匹配的优势，提出一种适用于非同源异构遥感影像的配准方法，该方法以不同滤波尺度下位置稳定不变的角点作为控制点．借助金字塔分层映射搜索策略．通过计算领域和搜索区域内的相位一致性，引入相关系数作为相似性测度实现同名点的匹配。并最终给出利用可见光、近红外及ＳＡＲ影像对该方法进行配准实验的验证和分析。２配准方法２．１控制点提取
为了保证提取足够多且位置相对稳定的控制点进行匹配和模型参数估计。需要使用通用性强、可靠性高的角点特征提取算子在待配准的输入影像上执行计算和检测。当前，Ｈａｒｒｉｓ”１和ＦＡＳＴ０１０１被证明具有较高的角点提取重复率，适用于各类光谱组成结构的影像。本文选择Ｈａｒｒｉｓ算子在输入影像上提取特征控制点，Ｈａｒｒｉｓ通过局部图像信息的自相关方式，在各个方向都计算了灰度变化．其角点提取的原理是：首先使用标准偏离差为口的高斯滤波模板与原始灰度影像１进行卷积运算。得到平滑后的影像７，然后计算７的一阶梯度矩阵Ｍ．公式为
妖道耘准设计ＲＡ儿ＷＡｙ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
Ｍ：『玖”）
【，，（ｘ，Ｙ），，（ｘ，Ｙ）
ｒ（Ⅵ）７，（埘’１（１）
Ｊ
咏ｚ，Ｙ）
其中，。（ｚ．Ｙ）和，。（ｚ．Ｙ）分镕Ⅱ代表平滑影像，中点（ｚ，Ｙ）处沿ｚ和Ｙ方向的梯度，通过解算矩阵肼的２个特征值Ａ．和Ａ：，使用判别函数Ｒ得到对于点（，，Ｙ）的特征判断，Ｒ的计算公式为
Ｒ＝＾，Ａ，一ｋ（＾．＋Ａ２）４
（２）
其中ｋ为介于０．０４—０．０６的常数，如果２个特征值都比较大，说明在该点处沿任意方向运动会引起判别函数Ｒ的急剧变化，判断为具有强响应的角点，否则为弱边缘或平滑域点．在使用阈值ｒ对该值进行局部非最大抑制后，可以获得稳定的特征点。
由于使用不同标准偏差口的高斯卷积模板会使平滑处理后的影像，变模糊．在特定滤波器尺度下检测的角点可能会在其他尺度的滤波影像中发生位置偏移甚至消失。为了提取不受高斯滤波尺度影响的稳定角点。需要使用临近尺度的高斯模板重新提取角点进行检测和判断。假设角点（‰，Ｙ。）是使用尺度（标准偏差值）为∥。的高斯滤波模板所提取的Ｈａｒｒｉｓ特征点，构造尺度为盯ｌ＝矿ｏ－Ｓｘｃｒｏ和盯２＝盯ｏ＋占ｘ盯。的高斯滤波模板对原始影像，进行卷积运算，提取的角点分别为（ｇ。。Ｙ．）和（ｘ：，Ｙ：），当提取的角点满足条件
Ｄ（（ｘｏ，Ｙｏ），（ｘｌ，ＹＩ））≤２＆ｒ
（３）
并且
Ｄ（（ｘｏ．Ｙｏ），（ｘ２，Ｙ２））≤２＆ｒ
（４）
时，认为（‰，Ｙ。）是一个不受卷积尺度影响的稳定角点。选为控制点，否则被排除。公式（３）、（４）中Ｄ代表欧式距离。即以角点（ｘ．，Ｙ．）和（ｚ：，Ｙ：）相对于（ｘ。，Ｙ。）的位置偏移作为特征点的判定依据。图１展示了分别使用原始Ｈａｒｒｉｓ算子和多尺度高斯滤波后Ｈａｒｒｉｓ提取的特征点分布情况，可以看出，在加人了公式（３）、（４）的判别后，位于道路边缘和开阔地中心的点受平滑尺度变化的影响，位置变化显著，被判断为不稳定的特征点排除．保留的点为稳定的不受平滑滤波干扰的特
征点。
２．２同名点匹配
由于Ｈａｒｒｉｓ特征算子提取的控制点为梯度变化显
咏准帕整仉为Ｉ５的商扣
随扳滤波后攫取结果
囤ｔ
６ｘ嘞为ｏ．５的临近Ｋ堰高抑髓扳滤波聍筛选挺墩络垠
Ｈａｒｒｉｓ算子提取特征点分布示意
２０１２（黝
・其他・
范登科，潘励，叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
著的角点，因此在其领域内影像的纹理信息比较丰富，
灰度变化显著，相关系数在小范围内容易形成变化剧烈
∑∑Ｗｓ（ｘ，ｙ）ＬＡ。（ｚ，ｙ）Ａ４，。（ｘ，ｙ）一ＬＪ
生生生——ｉ—－————————一（５）口
的极值，可以选择作为同名点的匹配测度。然而，非同Ｒ
、
源异构遥感影像间的光谱组成存在明显差异，在灰度强㈨０＝１
∑∑Ａ¨（ｚ，ｙ）＋８
度上呈现非线性变化，导致相关系数不能稳定描述相似式中，ＰＣ（ｚ，Ｙ）为点（ｚ，Ｙ）处的相位一致性，０为滤波度高的同名区域。针对这一问题，本文引入相位一致性方向，在区间［０，１ｒ］之间以固定步长变化，ｎ为滤波器
对影像进行转换，得到相关性程度高的特征图像。
尺度ｓ的最大值，矾（ｘ，Ｙ）为频率扩展的权重系数，
相位一致性是一种使用傅里叶谐波分量描述信号Ａ。（ｘ，Ｙ）代表使用沿日方向尺度为ｓ的ＬｏｇＧａｂｏｒ滤局部强度的特征，它利用特定方向口不同尺度ｓ的Ｌｏｇ波器进行分解后谐波分量的振幅值，ＬＪ表示仅取正值
Ｃａｂｏｔ滤波器，将信号分解为频率域下与日和ｓ相对应的运算，负值全部赋为０，Ｔ为噪声阈值，ｓ是一个避免
的傅里叶谐波分量．通过加权叠合不同尺度的分量，获除零的常数，ａ４，。（ｚ，ｙ）可由分解后谐波分量的相位得具有原始信号强弱特性描述的响应。Ｍｏｒｒｏｎｅ证明角计算获得，公式为
了相位一致性在描述图像特征信息上与人类神经视觉Ａ４，¨（ｚ，ｙ）＝ＣＯＳ（咖Ⅲ（ｚ，Ｙ）一币（ｘ，Ｙ））一
具有相一致的敏感度…ｏ。Ｗｏｎｇ指出。相位一致性特Ｉｓｉｎ（币。（ｚ，Ｙ）一咖（＾，Ｙ））Ｉ
（６）
征具有抵抗图像对比度和亮度变化的能力，并利用其其中咖（ｘ，ｙ）代表沿口方向尺度ｓ从１到ｎ之间的所构造的最小和最大距特征，分别实现了遥感影像特征有相位小的加权平均值。使用具有不同光谱组成结点的提取和同名点的匹配”…。应用于图像二维离散构、不同传感方式采集获得到的同名区域影像进行实信号，相位一致性的计算公式为
验分析。分别计算其相位一致性得到的结果如图２ＰＣ（ｘ．Ｙ）＝
所示。
（劬第３授段反射率彻第３渡段相位一致性
（ｃ）第５渡段反射率（ｄ）第５渡段相位一致性
囤２
Ｌａｎｄｓａｔ５
ＴＭ影像殛其相位一致性特征围像
从图２可以看出，边缘和角点处的相位一致性明分辨率，缩小图幅范围。首先，在输入影像的顶层金字显大于平坦开阔地区，通常在角点处取得极值。同时．塔图像上提取特征点．并搜索参考影像对应尺度的金相位一致性特征影像反映的并不是二值化的边缘，而字塔，找出相关系数大的窗口中心作为匹配点．使用简代表的是不同尺度和方向频率域卷积运算下傅里叶谐单的二维平面转换关系模拟初始的几何变形，根据最波分量的累积值。因此，以特征点为中心的窗口发生小二乘准则删除残差大的误匹配点对；而后。利用不同位置改变时，不会引起相关系数的剧烈变化，且由于角尺度金字塔图像的映射关系，在下一层更高分辨率尺点处相位一致性的显著特征，使得极值比较容易获得。度的金字塔图像上获取同名区域，执行与上述一致的经计算对应图像的相关系数，图２中（ａ）与（ｃ）之间的精确匹配，直到达到底层即原始分辨率影像为止。两相关系数为一０．０４５，而（ｂ）与（ｄ）达到０．６０５，可见采层金字塔结构的同名区域映射关系及几何变换模型估用相位一致性特征度量影像间的相关性明显优于原始计如图３所示。
灰度。
２．３模型参数估计与重采样
遥感影像在采集过程中一般表现出大视角、宽视利用匹配获取的同名点对．依据最小二乘原理估场的特性，成像的图幅范围大，使全局匹配的效率降计几何畸变模型的参数，剔除残差大的同名点对以保低，计算复杂度增大，需要使用特定的搜索策略优化同证模型估计精度在可接受的范围内，最终使用特定的名点匹配过程。基于金字塔结构的同名区域预测与分插值方法实现对输入影像的重采样处理，输出在地理层索引映射是目前最受欢迎的、应用于遥感影像匹配位置上与参照影像配准的校正影像。在这一过程中，中的搜索策略．方法是以一定比例逐级降低图像空间
全局性几何畸变模型的选择显得十分重要。
１２０
铁道嚣准设计ＲＡＩＬＷＡＹＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
２０１２（２）
范登科．潘励．叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
．其
他・
最大相关的窗口，以其中心作为匹配同名点；（４）＿｝｝Ｊ简单二维平面变换模型估计畸变改正关系，预测下一级
，
金字塔图像的同名区域；（５）在下级金字塔图像的子区域内迭代执行步骤（２）一（４），直到完成原始分辨率影像上同名点的匹配；（６）根据待配准影像的类型确定全局几何变换模型，采用最小二乘原理迭代解算和
圈３基于金字塔结构的遥豫影像同名区域索引和搜索策略
修正模型精度；（７）重采样输出配准后的结果影像。图４展示了配准策略的完整流程。
由于传感器姿态和侧视角的差异，采集获得的多源异构遥感影像在全局几何畸变上表现出不同形式。中分辨率光学影像的成像方式为推扫式，其传感器所处轨道高度及线阵扫描方式决定了影像几乎不受侧视角和地形起伏的影响．在全局影像上仅存在平移、旋转和尺度畸变，因此可以使用仿射模型描述待配准光学影像问的几何变形关系。仿射变换模型公式为
ｆ３１：＾『００８１『。１＋［ＣｌＬｙＪＬｓｉｎ０一ｃ０８０ＪＬ训
Ｌ
８８１“８
（７）
图４配准方法流程
ｒＪ
式中，（ｚ，ｙ）和（Ｘ，Ｙ）分别为对应同名点在输入影像和参考影像的图像坐标，＾代表尺度缩放因子，０表示影像的相对旋转角，（ｃ，ｒ）是图像在二维平面的相对平移量。因此仅需要两对同名点即可以求解该模型参数。而同分辨率ＳＡＲ影像的采集使用主动传感方式，虽然轨道高度避免了地形起伏的影响，但是其成像图幅范围较小．并且存在一定的侧视角，因此对于光学影像与ＳＡＲ影像以及不同时相ＳＡＲ影像问的几何畸变，通常采用投影变形模型来描述，一次项投影模型公式为
ｇｘ＋ｈｙ＋１
ｇｘ＋ｈｙ＋ｌ
３实验与分析
３．１实验背景和数据
实验以ＷｉｎｄｏｗｓＸＰ系统为平台，在Ｖｉｓｕａｌ法．程序运行的硬件环境为Ｉｎｔｅｒ
２
Ｃ＋＋
２００８集成开发环境中设计开发程序实现上述配准方
Ｃｏｒｅ２ＣＰＵ２．６６ＧＨｚ，
ＧＢ内存。实验数据选用３组遥感影像，每一组中输
ｘ：旦生譬２掣．ｙ：坐±丝（８）
入影像和参考影像分别采集自异掏光谱范围或不同传感器。前２组数据均由光学遥感影像组成，采集的区域为湖北省大东湖和梁子湖水域．其中第ｌ组数据来自于同一传感器，为了体现空间位置上的差异性，手工对输入影像旋转ｌｏ。，第２组数据来自不同传感器，存在微小旋转及１０～２０像素的平移变形。第３组数据为光学影像与ＳＡＲ影像之间的配准，为了正确计算相关系数．对ＳＡＲ影像进行重采样．确保其与光学影像的分辨率一致。３组实验数据的详细信息如表ｌ
所示。
式中包括。一＾总计８个模型参数，需要４对同名点进行求解，分母部分描述了成像平面的倾斜变形。在最小二乘原理基础上对模型参数进行估计后，使用参考影像中的控制点坐标反算输人影像控制点的模型坐标，并与真实坐标求差后计算均方根误差ＲＭＳＥ，ＲＭＳＥ值的大小代表了该同名点对偏离模型估计的程度．删除ＲＭＳＥ值最大的多个点对，重新估计模型参数，迭代执行这一过程直到同名点对集合中最大的ＲＭＳＥ值小于预先设定的阈值。以获得满足纠正精度要求的模型参数。最终，双线性插值方法被用于重采样过程输出配准结果。２．４总体配准策赔
本文提出的基于相位一致性相关的多源遥感影像配准方法可概括为以下７个步骤：（１）根据影像大小。设定特定的分辨率比例，分别构建输人和参考影像的多级金字塔图像；（２）在输入影像的顶层金字塔上使用多尺度Ｈａｒｒｉｓ提取控制点．并计算其领域内窗口的相位一致性；（３）计算参照影像顶层金字塔的相位一致性特征图像，使用相关系数作为匹配测度全局搜索
帙遵标准艟汁ＲＡ儿ＷＡＹ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
在进行影像配准运算处理过程中．需要对相关参数进行初始化．具体项目及其初始值为：金字塔各级比例Ｉ：３，金字塔级数３，原始分辨率下搜索区域大小２００像素，相关系数阈值０．７５．相关系数计算窗口大小１３像素．ＲＭＳＥ阈值Ｉ．５。
３．２配准精度分析
使用以上３组数据对本文配准方法进行实验验证，并对该方法的适用性和精度进行评价。表２列出了经实验获得的配准评价结果。
由表２可以看出．在配准精度方面，第３组数据的效果最好，最大残整和总体均方根误差指标都控制在
１２１
２０１２（引
・算他・
范登科，潘励，叶沅蠢一一种基于相位一致性相关的多源遥感影像配准方法
表１实验数据信息概览
ｌＬａｎｄｓａ＊５ＴＭ
ｂａｎｄ３／０．６３－０
６９
ｐｍ
３０ｍ＋ｌＯ。３５００ｘ３５００２００７—０４—１９
Ｌａｎｄｓａｔ５ＴＭｂａｎｄ５／１５５—１７５ｌ＾ｍ３０ｍ３０００ｘ３ＯＯＯ２００７—０４—１９２ＨＪ—ｌＡＣＣＤＩｂ衄ｄ４／０．７６—０．９０“ｍ
３０ｍ３２５０ｘ３２５０２０１１一Ｏｌ一１６Ｌａｎｄｓａｔ５ＴＭｂａｎｄ５／１５５—１
７５＊ｍ
３０ｍ３０００ｘ３０００２００７－０４一１９３
Ｌａｎｄｓａｔ５ＴＭ
ｂａｎｄ
１２３／可见光
３０ｍ
Ｉ７００ｘ１６５０２００６—０９一１０
ＡＬＯＳｐａｉｓａ
ｂａｎｄＩ／ＳＡＲ
３０ｍ
１８００ｘ１７００２００５—１１—１０
表２配准方法的实验评价结果
范围内，且在对比度高的局部区域内相对集中。符合相
位一致性描述亮度和对比度差异的不变特征。
综合分析３组实验数据的配准结果，本文方法更适用于对非同源异构的两景光学影像进行配准。对于光学与ＳＡＲ影像的配准有一定局限性和区域选择性，
１个像素以内，然而同名点对的误匹配率却达到了同时将相关系数作为匹配测度引入到该方法中。对于
９５％以上，仅有５对正确匹配的同名点用于模型参数全局仿射和投影变形有一定的抵抗能力。估计，且它们分布集中于图５（ｅ）、（ｆ）中左下角的岛屿３．３对比验证
海岸线上，占影像中主体的陆地范围内并没有正确匹为了进一步验证本文方法的配准效果．引入传统基配的同名点对，这是因为大陆中地物类型丰富而海岛于灰度的相关系数（ＮＣＣ）、尺度不变特征（ＳＩＦＴ’）、改进内单一．ＳＡＲ影像之于前者表现出更多的噪声．而相尺度不变特征（ＳＲ—ＳＩＦＴ）３种方法分别对３组实验数据位一致性的计算结果对噪声十分敏感．造成陆地中相进行配准实验，统计各项实验结果的误匹配率和均方根位一致性的相关程度并不是十分显著；反之，海岛在中误差（ＲＭＳＥ），得到的对比情况如表３所示。
ＳＡＲ影像中纹理结构简单，海岸线轮廓特征明显，有表３
３种常用方法实验结果对比（“一”表示配准失败）
利于相位一致性稳定描述局部图像信息。
配准
实验１
实验２实验３
方法
误匹配率／％ＲＭＳＥ谩匹配率／％ＲＭＳＥ误匹配卓／％ＲＭＳＥ
ＮＣＣ
５８
３８
０８７１４４．１７０．９４３圈鬻糕一一ＳＩｒｒ６６７３０９９３５５．９７０９８４
一
一ＳＲ—ＳＩＦＴ
５５
３６
０９７７
４１．８４
０９９６——
一
实验过程中设定ＲＭＳＥ阈值为１．０像素，以保证几何校正模型参数估计的准确性和一致性。当３种方法达到相同精度水平时，结果显示实验２的误匹配率最低，原鬻翟甄
因是两景近红外影像间灰度相关性程度最高．且与特征表述相关的梯度强度和方向变化较为稳定．与此同时，ＮＣＣ和ＳＲ—ＳＩＦＴ方法的误匹配率明显低于本文方法。相反，３种方法对实验ｌ数据进行配准的误匹配率高于或接近本文方法，可见，对于光谱组成差异越显著的影像来｛¨Ｊ？’｝＿＿｛；＿ｖ
－。‘１ｒ。、，“…
＿・－
“≮
说，本文方法的适用性更强，效果更好。３种方法不能完图５
３组实验藏据影像及同名点分布情况
成对实验３数据配准的事实更进一步验证了这一结论。前２组实验数据的配准结果在精度和误匹配率上因此，与３种配准方法进行对比分析得出的结论是，基于达到～致水平。尽管第２组数据的总体均方根中误差相位一致性相关的配准方法可以应对传统方法所不能解略微偏高，但其匹配得到的同名点对数目是第ｌ组的决的异构遥感影像间难以配准的问题。３倍，是因为第２组数据的影像对均由近红外波段构４结语
成，相对于第ｌ组的可见光一近红外组合，相关程度更高。然而，在匹配获得大量同名点对的同时．发生误匹本文提出一种基于相位一致性相关的多源遥感影
配的点对数目也相应提升，且从２组数据对比情况来像配准方法，该方法以多级图像金字塔结构为索引和
看，这一比例水平并不受光谱组成结构的影响．因此要搜索策略，首先使用多尺度高斯滤波判别和遘选稳定迸一步降低误匹配率，仍需要对匹配方法和策略进行的Ｈａｒｒｉｓ特征点．在控制点领域和相关系数极值搜索改进。另一方面，由于第２组影像采集自不同时问．地区域内，引人相位一致性特征进行相似性度量，在降低物变化也会造成同名点的误匹配。从图５（ａ）～（ｄ）误匹配概率的同时获得更多精确匹配的同名点对。经
中可以看出，正确匹配的同名点均匀分布于整个图幅
实验分析和对比验证．该方法适用于中分辨率多源异
１２２
铁道捧准设计ＲＡＩＬＷＡＹ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
２Ｄｊ２（２）
范登科，潘励，叶沅鑫一一种基于相位一致性相关的多源遥感影像配准方法
．其他・
构光学一光学、光学一ＳＡＲ遥感影像间的配准，精度满ｂａｓｅｄ
ｏｎ
ＳＩＦＴ［Ｊ］ＥｌｅｃｔｒｏｎｉｃｓＬｅｔｔｅｎ．２００８．４４（２）：１０７一１０８
足实际生产过程中影像几何校正处理的需要，且高于［６］ＧａｏＦ．Ｍａ蚰ｋ
Ｊ，Ｗｏｌｆｅ
Ｒ
Ｅ
Ａｕｔｏｍａｔｅｄｒｅｇｉｓｔｒａｔｉｏｎ
ａｍｌｏｒｅｈｏｎｃｔｉｆｉｃａｔｉｏｎ
Ｌａｎｄｓａｔａｎｄ
Ｌａｎｄｓａｌ—ｌｉｋｅ
ｄａｔａ
其他基于区域匹配的配准方法。下一步的研究工作将ｐａｃｋａｇｅ
ｆｏｒ
ｐｒｏｃｅｓｓｉｎｇ【Ｊ］．Ｊｏｕｒｎａｌ
ｏｆＡｐｐｌｉｅｄＲｅｍｏｔｅ
Ｓｅｎｓｉｎｇ．２００９．３（１）
专注于提高运算效率，改进算法以降低误匹配率，并适［７］ＢｕｎｔｉｎｇＰ．Ｌ“ｍ∞Ｆ．Ｌｕｃａｓ
Ｒ
Ａ
ｍｕｌｔｉ－ｒｅｓｏｌｕｔｉｏｎ㈣‘ｂａｓｅｄ
用于不同分辨率、地形起伏等复杂条件下的多源遥感ｔｅｃｈｎｉｑｕｅｆｏｒａｕｔｏｍａｔｉｃｍｕｌｔｉ－ｍｏｄａｌｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ［Ｊ］Ｉｍａｇｅ
ａｎｄ
影像配准。同时希望本文方法能为图像配准领域的学Ｖｉｓｉｏｎ
Ｃｏｍｐｕｔｉｎｇ．２０１０．２８（８）：１２０３—１２１９．
者提供研究基础和借鉴。［８］Ｓｕｒｉ
Ｓ．Ｒｅｉｎａｒｔｚ
Ｐ．Ｍｕｔｕａｌ－ｉｎｆｅ．ｎａｔｉｏｎ－ｂｍｄ
ｒｅｇｉｓｔｒａｔｉｏｎ
ｏｆ
参考文献：
ＴｅｎａＳＡＲ．Ｘａｎｄ
Ｉｋｏｎｏｓ
ｉｍａｇｅｒｙｉｎ
ｕｒｂａｎ
ｌｉｔｅｒ［Ｊ］．ＩＥＥＥ
ＣＩ】ｚｉｔｏｖ＾Ｂ．Ｆｌ㈣Ｊ
Ｔｒａｎｓａｃｔｉｏｎｓ
ｏｎ
Ｇｅｏｓｃｉｅｎｃｅ
ａｎｄ
Ｒｅｍｏｔｅ
Ｓｅｎｓｉｎｇ．２０１０．４８（２）：
Ｉｍａｇｅ”ｇｌｓｘｒｍｉｏｎ
ｍｅｔｈＭ‘：Ａｓｕｒｖｅｙ【Ｊ］｜ｍａｒｅ
９３９—９４９
ａｎｄＶｉｓｌｏｎ
Ｃｏｍｐｕｔｉｎｇ．２００３．２１（Ｉ１）：９７７—１０００
［９］Ｈａｒｒｉｓ
ｃ．ＳｔｅｐｈｅｎｓＭＫ
Ａｃｏｍｂｉｎｅｄｃｏｍｅｒ
ｓｏｄ
ｅｄｇｅ
ｄｅｔｅｃｔｏｒ［Ｃ］．
［２］Ｌｏｗｓ
ＤＧ
Ｄｉｓｔｉｎｃｔｉｖｅｉｍａｇｅｆｅａｔｕ亿ｓ
ｆｒｏｍ
ｓｃａｌｅ－ｉｎｖａｒｉ肌ｔｋｅｙｐｎｉｎｔｓ
Ｐｒｏｃ
（Ｊ】ＩｎｔｅｒｎａｔｉｏｎａｌＪｏｕｒｎａｌ
ＡｌｖｅｙＶｉｓｉｏｎＣｏｎｆ．．１９８８：１４７—１５２．ｏｆＣｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎ，２００４．６０（２）：９１一１１０
ＰＣＡ—ＳＩＦＴ：Ａ㈣ｄｉｓｔｉｎｃｔｉｖｅ
［１０］ＲｏｅｔｅｎＥ．ＤｒｕｍｍｏｎｄＴ
Ｍａｃｈｉｎｅ［３］Ｋｅ
ｌｅａｒｎｉｎｇ
ｆｏｒ
ｈｉｇｈ—ｓｐｅｅｄ
ｃｏｒｎｅｒ
Ｙ．ＳｕｋｔｈａｎｋａｒＲ
ｒｅｐ＂ｓｅｎｔａｔｉｏｎ
ｄｅｔｅｃｔｉｏｎ［Ｃ］ｐｒｏｃｅｅｄｉｎｇｓ
ｏｆｔｈｅＥｕｒｏｐｅａｎＣｏｎｆｅｒｅｎｃｅ
ｏｎ
Ｃｏｍｐｕｔｅｒ
ｆｏｒｌｏｃａｌｉｍａｇｅ
ｄｅ虻ｄ阱ｏｎ［Ｃ］Ｐｒｏｃｅｅｄｉｎｇｓ
ｏｆｔｈｅＩＥＥＥ
Ｃｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎ．２００６：４３０—４４３
ＳｏｃｉｅｔｙＣｏｎｆｅｒｅｎｃｅｏｎ
Ｃｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎ蛐ｄ
ＰａｔｔｅｍＲ∞ｏｇｎｉｔｉｏｎ．
［ＩＩ］Ｍｏ肿ｎｅ
Ｍｃ．ＯｗｅｎｓＲＡ．Ｆｅａｔｕｒｅｄｅｔｅｃｔｉｏｎｆｒｏｍｌｏｃａｌｅｎｅｒｇｙ［Ｊ］
２００４：５０６—５１３
［４］Ａｂｄｅｌ・Ｈａｋｉｍ
ＡＥ．ＦａｒａｇＡＡ
ＣＳＩＦＴ：Ａ
ＳＩＦＴ
ＰａｔｔｅｒｎＲｅｃｏｇｎｉｔｉｏｎＬｅｔｔｅｒｓ．１９８７．６（５）：３０３—３１３．
ｄｅ虻ｒｌｐｔｏｒｗｉｔｈｃｏｌｏｒ
［１２］Ｗｏｎｇ
Ａ．Ｃｌａｕｓｉ
Ｄ
Ａ
ＡＩＳＩＲ：Ａｕｔｏｍａｔｅｄｉｎｔｅｒ．柏ｎｍｒ／ｉｎｔｅｒ—ｂ蛆ｄ
ｉｎｖａｒｉａｎｔｃｈ盯ａｅｔｅｒｉｓｔｉｅｓ［Ｃ］Ｐｒｏｃｅｅｄｉｎｇｓ
ｏｆ
ｔｈｅ
ＩＥＥＥ
Ｃｏｍｐｕｔｅｒ
ｓａｔｅｌｌｉｔｅｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎｕｓｉｎｇｒｏｂｕｓｔｃｏｍｐｌｅｘｗａｖｅｌｅｔｆｅａｔｕｒｅｒｅｐｒｏ－
Ｓｏｃｉｅｔｙ
Ｃｏｎｆｅｒｅｎｃｅ
ｏｎ
Ｃｏｍｐｕｔｅｒ
Ｖｉｓｉｏｎａｎｄ
ＰａｔｔｅｍＲｅｃｏｇｎｉｔｉｏｎ．
盹ｎｔａｔｉｏｎａ［Ｊ］Ｐａｔｔｅｒｎ
ＲｅｃｏｇｎｉｔｉｏｎＬｅｔｔｅｍ，２０１０．３１（１０）：１１６０一
２００６：｜９７９一１９８３
１１６７．
［５】Ｙｉｚ．ＺｈｉｇｕｏＣ，ＹａｎｇＸ
Ｍｕｌｔｉ—ｓｐｅｃｔｒａ］ｍｍｏｔｅｉｍａｇｅｒｅｇｉｓｔｒａｔｉｏｎ
（上接第１１２页）
翻车棚的端头设置翻车机监控室．便于操作人员对翻的取送车，棚的跨度需大于２条重调机走行轨道间距
车机作业时的监控。控制室内还设有配电室、夹轮器（２５．４１１１），建筑模数取整后，棚的跨度取２７
ｍ。
液压站等房同。翻车机卸车系统工艺平面如图９
翻车机控制室设置在２条重车线的中间，在靠近
所示。
１８０ＤＤ
２５００［
６ｌ瑚６０００６０００
６０００
６０∞
甲镌盎’谶迸芈方向
洲州
，●Ｉ————』
｜．．｛
Ｌ亡三Ｉｌ圭些
￡一位城料枷…。笠Ｌ—Ｌ丑墅磬
《习ｉ蕊。弘嗣矧
一．ｕ翌ｌ
＂
旺匠ｉ口：Ⅱ旺
口ｊＩＩｎ…————■一■．－……“
Ｚ『央轮器罗呷ｉ
晨
ｊ．｜孑｝—＝Ｆ—孔
｝’
’
。
。
’
’
’
’
’／。
’
。怍革≠±手ｒ＋一｝＝二ｊ州Ｉ仃Ｉ警≤必ｌ錾Ｉ…｝
一～一
塑塑。
…一■＿
川Ｌ
安全止Ｈ群
圈９翻车机卸车系统工艺平面布置（单位：ｍｍｌ
５结论
参考文献：
［Ｉ］张磊．马明札．燃料远行与检謦［Ｍ］．北京：中国电力出版
翻车机卸车工艺作为铁路货场一种新型的卸车工杜．２００６．
［２］高志刚．新“Ｃ”型折退式双车翻车机在定洲发电公司＝期工程的
艺，具有生产能力大、自动化程度高、对车辆损伤小和应用［Ｊ】神毕科技．２０１０，８（Ｉ）
环境污染小等特点，应推广翻车机在大宗散装物料铁［３］＊撬军．蛱路货车翻车机使用中存在的问题及解挟方击［Ｊ］铁道
路货场的应用。进行翻车机卸车工艺设计时，应优先车辆．２００９．４７（９）．
采用折返式布置；如采用２台翻车机，则应采用２台翻
［４】
张凡华鼻车机设鲁选捌分析［Ｊ］华电技术．２００８．３０（６）．
车机对称。相背的布置形式。
【５］石小东．幢路散货运■不同卸车方式比较［Ｊ］．内簟古煤炭烃脐．
２００８（３）
侠嬉檬准设计ＲＡＩＬＷＡＹ
ＳＴＡＮＤＡＲＤ
ＤＥＳＩＧＮ
２０１２（２）
１２３
一种基于相位一致性相关的多源遥感影像配准方法
作者：作者单位：刊名：英文刊名：年，卷(期)：
范登科， 潘励， 叶沅鑫， Fan Dengke， Pan Li， Ye Yuanxin武汉大学遥感信息工程学院,武汉,430079铁道标准设计
Railway Standard Design2012(2)
1. Zitová B;Flusser J Image registration methods:A survey[外文期刊] 2003(11)2. Harris C;Stephens M K A combined corner and edge detector 1988
3. Suri S;Reinartz P Mutual-information-based registration of TerraSAR-X and Ikonos imagery in urban area 2010(02)4. Bunting P;Labrosse F;Lucas R A multi-resolution area-based technique for automatic multi-modal image registration2010(08)
5. Gao F;Masek J;Wolfe R E Automated registration and orthorectification package for Landsat and Landsat-like dataprocessing 2009(01)
6. Yi Z;Zhiguo C;Yang X Multi-spectral remote image registration based on SIFT[外文期刊] 2008(02)7. Abdel-Hakim A E;Farag A A CSIFT:A SIFT descriptor with color invariant characteristics 20068. Ke Y;Sukthankar R PCA-SIFT:A more distinctive representation for local image descriptor 20049. Lowe D G Distinctive image features from scale-invariant keypoints[外文期刊] 2004(02)
10. Wong A;Clausi D A AISIR:Automated inter-sensor/inter-band satellite image registration using robust complex waveletfeature representations 2010(10)
11. Morrone M C;Owens R A Feature detection from local energy 1987(05)12. Rosten E;Drummond T Machine learning for high-speed corner detection 2006
本文链接：http://d.g.wanfangdata.com.cn/Periodical_tdbzsj201202034.aspx