第30卷第1期
2010年1月
文章编号:1001-9081(2010)01-0098-03
计算机应用
JournalofComputerApplications
Vol.30No.1
Jan.2010
几种图像相似性度量的匹配性能比较
陈卫兵
(南通职业大学电子工程系,江苏南通226007)
(ntcwb@mail.ntvc.edu.cn)
摘 要:针对景象匹配中匹配性能和匹配实时性会受相似性度量选择影响的问题,从常用的相似性度量(归一化
积相关、相位相关、均平方差和去均值均平方差)入手,对图像施加噪声和进行畸变(如图像旋转变化、图像比例变化、光照强度变化和云层遮挡等),通过相似性度量匹配性能的仿真试验,从匹配性、适应度和实时性等方面对各相似性度量进行比较,并对结果进行了归纳总结与证明。
关键词:图像匹配;相似性度量;匹配性能中图分类号:TP391 文献标志码:A
Comparisonofmatchingcapabilitiesinsiilarityents
CHENWei2(DepartmentofElectronicEngineering,NNantongJ226007,China)
Abstract:Imagematchingreal2eonhowtochoosesimilaritymeasurementmethodofimagetotopicin2depth,thepapercarriedoutaseriesofsimulation
experimentsonoffoursimilaritymeasurements(NProd,PC,MSDandequalizationMSD)throughonimagedistortions.Theimagedistortionsincludepictureorientation,imagescalechange,illuminationintensitychange,cloudscover,andsoon.Theexperimentalresultswerecomparedwiththosefromthematchingperformance,adaptivecapacityandreal2timecapability.Analysisandtheoreticaldemonstrationoftheresultswerealsogiven.
Keywords:imagematching;similaritymeasurement;matchingperformance
0 引言
景象匹配技术在飞行器制导定位等领域一直是人们研究的热门课题,具有广泛的应用前景。研究可靠性高、鲁棒性好的图像匹配算法是提高匹配制导性能的有效途径。按照
Brown的理论
[1]
度量,两幅图像越相似,则距离越小。常用的距离度量有。
1)平均绝对差度量MAD。D(u,v)=
m×n
m
n
∑∑|
j=1k=1m
n
X(j+u,k+v)-Y(j,k)|
(1)
,用于景象匹配各种匹配算法都是4个元素
2)平均平方差度量MSD。D(u,v)=
m×n
的不同选择的组合:特征空间、相似性度量、搜索空间、搜索算法。其中相似性度量是指衡量匹配图像特性之间相似性的尺度,在景象匹配中占有重要的位置。基于各种相似性度量的算法研究不在少数,但对相应的适用场合没做详细的分析介绍,如何根据参与匹配的实际情况,选择合适的相似性度量,对有效地抑制噪声、畸变干扰,提高匹配概率与匹配精度有重要的意义。
本文通过对几种常用的相似性度量的匹配性能进行仿真并加以分析证明,给出了几种普遍存在的因素对相似性度量的影响以及相似性度量的适应性和实时性,为实际图像匹配应用中选择何种相似性度量提供了理论依据。
∑∑(X(j+u,k+v)
j=1k=1m
n
2
-Y(j,k))
(2)
3)去均值平均绝对差度量。D(u,v)=
m×n
∑∑(X(j+u,k+v)
j=1k=1
-(3)
X (u,v)-Y(j,k)+Y )
2
其中:X (u,v)为基准子图的平均值,Y 为实时图的平均值。
4)Hausdorff距离度量。
Hausdorff距离度量是适合于二值图像(如边缘特征图
像)相似性的度量。给定有限的两个点集A和B,A和B之间的
Hausdorff距离定义为:
H(A,B)=max[h(A,B),h(B,A)]
(4)
h(A,B)称有向
1 常用的相似性度量
图像的相似性度量有好多种,但应用广泛,相对成熟的相似性度量可以归结为两大类:距离度量和相关度量。设矩阵X为基准图矢量,大小为M×N,Y代表实时图矢量,大小为m×
n,则几种常见的相似性度量可表示如下。
其中:h(A,B)
Hausdorff距离
=max-min‖a-b‖,
a∈A
b∈B
[2]
。
1.2 相关度量
相关度量是定义在某种相关意义下的度量[3]。
1)积相关度量。
a)归一化积相关度量(NormalizedProduct,NProd)。
1.1 距离度量
距离度量是相比较的两幅图像定义在某种距离意义下的
收稿日期:2009-06-25;修回日期:2009-08-04。
基金项目:江苏省自然科学基金资助项目(BK2008188);江苏省六大人才高峰项目(072E2013)。 作者简介:陈卫兵(1966-),男,江苏南通人,副教授,硕士,主要研究方向:人工智能。
第1期
R(u,v)=
m
n
陈卫兵:几种图像相似性度量的匹配性能比较99
均值为0,因此仿真时对均值不为0的噪声不作考虑。
X(j+u,k+v)Y(j,k)
∑∑
对截取图添加均值为0,方差σ为0.
032、0.072、0.149、
[∑∑X
j=1k=1
mn
2
(j+u,k+v)
][∑∑Y
j=1k=1
1/2
mn
2
(j,k)
]
1/2
(5)
0.200产生的模拟实时图如图2所示。
b)去均值归一化积相关度量。
m
n
∑∑(X(j+u,k+v)
R(u,v)=
mn
j=1k=1
-X (u,v))Y(j,k)-X (u,v))
2
[∑∑(X(j+u,k+v)
m
n
]
1/2
×
图2 添加噪声的仿真用实时图
2.3.2 图像旋转变化的仿真
(6)
景象匹配导航中,由于惯性仪表的测量误差,拍摄的实时图与基准图总存在一定的角度偏差,这种偏差可用对截取图施加旋转变化来仿真,但是不能直接从截取图中旋转得到,因为截取图直接旋转会使截取图变大,。采取的措施是:;
以此位置为中心,,;对旋转后的图采用上述方法所得具有旋转畸变的实时图如图3所示。
[∑∑(Y(j,k)
j=1k=1
-Y )
2
]
1/2
2)相位相关度量(PhaseCorrelation,PC)。
设gr(x,y)为参考图,gs(x,y)为实时图,Gr(u,v)、Gs(u,
v)分别为它们的Fourier变换,令它们的互功率谱为Gr3Gs,
3
则相位相关函数R为:
R=F
-1
G33
|Gr3Gs=F
-1
[e
φr-φs)j(
]7]。其中:F-1表示Fourier,r2由于相似性度量较多,进行一一仿真工作量过大,因此选用最具代表性的四种相似性度量(归一化积相关、相位相关、均平方差,去均值均平方差)进行仿真研究。
2.1 选用的基准图及截取图
图3 图像旋转的仿真用实时图
为了便于计算匹配性能(匹配概率、匹配误差),仿真时,选用从基准图截取一小块图像,添加噪声及施加畸变的方法产生仿真用实时图[5]。选用的基准图为分辨率1m的某地区可见光卫星照片,大小为412像素×419像素,如图1(a)所示
,从中截取的图像,大小为152像素×157像素,如图1(b)所示。
2.3.3 图像比例变化的仿真
景象匹配导航中,由于测高系统存在测量误差,拍摄的实时图经过比例校正后,仍然存在一定大小的比例误差,这种误差可采用对截取图施加缩放来仿真,具体措施为:记下截取图的大小及其中心位置;以此位置为中心
,重新在基准图上截取一块较大的图像,以保证缩小后实时图大小的区域内没有空白区域;对该截取图缩放;对缩放后的图像中心截取实时图大小的一块图像即为畸变的实时图。采用上述方法所得具有缩放畸变的实时图如图4所示。
图1 仿真所用的基准图像及截取图像
2.2 选用的噪声及畸变模型
图4 具有缩放畸变的仿真实时图
为了能够尽可能地反映基准图与实时图的差异,仿真选用式(8)所示的噪声及畸变模型。
Y(j,k)=αX(u0+j+δ,k)j,v0+k+δk)+n(j
(8)
2.3.4 光照强度变化的仿真
拍摄基准图与实时图的时间的不同,天气不同,都会造成光照强度的变化,
表现为在基准图与实时图灰度值的变化,即式(8)中α的变化,α取不同值所得的仿真用实时图如图5所示。
其中:X为基准图矢量,Y为实时图矢量,基准图与实时图差异除了具有位置平移(u0,v0)外,还存在随机噪声n(j,k),灰度畸变α,以及比例缩放、旋转等畸变引起位置的微小变化δj、δk等。
2.3 噪声及各种畸变的仿真2.3.1 图像噪声的仿真
图像的噪声是以式(8)中n(j,k)的值来体现的,仿真时,取n(j,k)为均值是0,方差为σ的随机噪声。其实,若n(j,k)不服从均值为0,实时图经过去均值预处理后,噪声就会变为
图5 光照强度变化的仿真用实时图
2.3.5 云层遮挡的仿真
景象匹配导航中,拍摄实时图时,出现云层遮挡的情况也
100 计算机应用
2.4 仿真结果
第30卷
不可避免,可采取从截取图中去掉一部分图像来仿真云层遮挡的情况。所得的仿真用实时图如图6所示
。
对截取图分别施加噪声及畸变用不同相似性度量进行匹
配的结果如表1所示。其中:匹配误差,即匹配到的位置(u1,v1)与截取位置(u0,v0)之间的距离
ε=(u0-u1)2+(v0-v1)2
度量值为匹配点的度量值,对于相关度量为相关曲面的峰值,对于距离度量为距离曲面的峪值。有下划线的结果表示明显的误匹配。
表1 实验结果
图6 云层遮挡的仿真用实时图
畸变
类型
匹配误差
畸变值
0.032
度量值
去均值
MSD0.000000.000000.[1**********]05.303.162302.236101.000001.000000.00000NProd0.967520.936620.897430.0.0.0.440.927603.162300.944660.924560.993840.993840.993840.937850.938870.93887
实时图图图图图图图图图图图图图图图
2(a)2(b)2(c)2(d)3(a)3(b)3(c)3(d)4)4(b4(c)5(a)5(b)5(c)6(a)6(b)6(c)
NProd0.000000.000000.000000.000001.000003.605605.385208.162100.000000.000000.000000.000000.00000相位相关
0.000000.000000.000000.000003.162304.472105.[1**********].770308.602301.000001.000001.000000.000000.000000.00000
MSD0.000000.000000.000000.000001.000003.605303.302.236101.000000.000000.00000相位相关
0.123320.094470.0680.0.390.044220.024630.017470.013680.015660.151130.151130.151130.092990.080740.08074
MSD(×107)1.555403.[***********].421302.701803.336602.513902.618603.809800.222510.293940.370110.175220.208730.20873
去均值MSD
(×107)1.555303.046705.019706.088601.846602.421102.701703.336502.506402.611003.784100.488770.935741.858901.565200.181570.18152
添加
噪声
0.0720.1490.2003.02°5.02°7.02°10.00°0.915
旋转
缩放1.1.0.70550.50600.37689%13%20%
光照变化云层遮挡
3 分析与证明
通过分析仿真结果得到如下结论。
1)添加噪声不会影响匹配性能。
从仿真结果可以看出,添加噪声后利用各种相似性度量匹配结果的误差都为0,即使噪声添加到使图像人眼都分辨不出来,匹配误差仍然为0。
mn
(X(j+u,k+v)-Y(j,k))2= D(u,v)=∑∑m×nj=1k=1
m×nj=1km×nj=1km×nj=1km×nj=1k
mmm
证明 基准子图和实时图只存在噪声,则畸变模型简化为如下所示:
(9)Y(j,k)=X(u0+j,v0+k)+n(j,k)
为了证明方便,设n(j,k)是服从均值为0,方差为δ的噪声,当n(j,k)的均值不为0,可以通过对实时图去均值处理,使n(j,k)是服从均值为0。
对于平方差度量:
∑∑(X(j+u,k
=1n
mn
2
+v)-(X(u0+j,v0+k)+n(j,k)))=
∑∑(X(j+u,k
=1n
+v)-X(u0+j,v0+k))
2
+
m×n
j=1k=1
∑∑n(j,k)
mn
2
-
∑∑(X(j+u,k
=1n
+v)-X(u0+j,v0+k))n(j,k)=
22
+v)-X(u0+j,v0+k))+δ-0(.)
∑∑(X(j+u,k
=1m
n
(10)
其中0(.)表示近似为0,因此D(u,v)与噪声无关。而对于归一化积相关有:
R(u,v)=
m
n
∑∑X(j+u,k
2
+v)[X(u0+j,v0+k)+n(j,k)]
m
n
[∑∑X
j=1k=1
(j+u,k+v)
m
n
][∑∑[X(u
j=1k=1
1/2
+j,v0+k)+n(j,k)]
2
]]
1/2
=
∑∑X(j+u,k
(j+u,k+v)
1/2
+v)X(u0+j,v0+k)+0(.)
m
n
[∑∑X
j=1k=1
mn
2
][∑∑[X(u
j=1k=1
+j,v0+k)+n(j,k)]
2
1/2
(11)
(下转第110页)
110 计算机应用第30卷
精确匹配,匹配的成功率得到了提高,由此,将更多的资源纳入为任务选择的范围内,使得系统在吞吐量、任务满足率和资源调度率上均有所提高。就失败率来看,供求基本平衡和供不应求两种情况下分别增长了0.35%和0.80%,这和其他参数的增长相比较,增长幅度甚微,是可承受的。它反映出用户在QoS参数的理解上的模糊性和不确定性,导致一般在定义中提出略高的要求,因此,软化几乎对QoS没有影响。
最大匹配算法从两个方面取得了进展:1)在满足QoS需求的基础上,适度扩大了资源选择范围,减少了因为QoS参数而导致的匹配问题,使得资源调度率和任务满足率提高;
2)确定了依照QoS选择资源的优先顺序,为任务选择了适当
证明,对QoS的软化处理提高了系统吞吐量和资源匹配成功率,最终实现了整个系统资源调度率的提高。参考文献:
[1] ALOISIOG,CAFAROM,EPICOCOI,etal.Resourceandservice
discoveryintheiGridinformationservice[EB/OL].[2009-06-15].http://www.gridlab.org/WorkPackages/wp210/Documents/igrid2iccsa.pdf.
[2] BRADLEYA,CURRANNK,PARRZG.Resourcediscoveryand
managementincomputationalGRIDenvironments[J].InternationalJournalofCommunicationSystems,2007,19(6):639-657.[3] CHENP,XUZ,ZHANGB.AsolutiontoQoScontrolandavaila2
bilitypromotionincomplexgridcomputing[C]//Proceedingsof12thIEEEInternationalConferenceonNetworks.NewYork:IEEE,2004:403-407.
[4] LICHUNLIN,LILAYUAN.UtilitymultipleQoSguaranteed
resourceoptimonmputing[C]//ProceedingstheIonalCoTheoryandApplica2tionson,ICoSociety,2009:165-169.R,D,GIDDYJ.Nimrod/G:Anarchitecture
ormanagementandschedulingsysteminaglobalcompu2tationalgrid[C]//Proceedingsofthe4thInternationalConference/ExhibitiononHighPerformanceComputingintheAsia2PacificRe2gion.NewYork:IEEE,2000:283-289.
[6] PLESTYSR,VILUTISG,SANDONAVICIUSD.Themeasurement
ofgridQoSparameters[C]//Proceedingsof29thInternationalCon2ferenceonInformationTechnologyInterfaces.NewYork:IEEE,2007:703-707.
的资源,而非当前网格中最好的资源,使系统工作在最优状态,由此进一步提高系统的可用性。与传统经典算法比较,最大匹配算法在任务一级取得了次优解,但在系统一级得到更优解。算法有利于克服网格中因自私性带来的缺点,在用户和系统效率之间取得折中,更有利于系统的稳定和健壮,尤其适用于资源供不应求状态下。
4 结语
,按照参数性质,。针对弹性参数,进行了软化处理,使其在不失定性含义的前提下,扩展了资源选择的范围,提高了任务和资源在QoS参数上的匹配率,由此提高系统的吞吐量,实现网格系统的调度优化。实验
(上接第100页)
可见式(11)的R(u,v)与式(6)具有相同位置的峰值点,即添加噪声也不影响NProd度量的匹配概率。
2)相关度量对于光照强弱的变化不敏感,而距离度量却比较敏感。
实时图和基准子图的相似程度可以直观地从两图的图像矢量在欧氏空间的几何关系得出来
。
关算法正确找到匹配位置。
4)对于旋转、缩放等畸变,几种相似性度量都比较敏感。
对实时图施加旋转、缩放畸变,各种相似性度量都不同程度地出现了匹配误差,甚至误匹配,可见旋转、缩放畸变是影响匹配概率的主要因素,因此要尽量消除旋转、缩放畸变后再进行匹配,才能取得较好的匹配性能。
4 结语
不同的相似性度量对于不同的畸变类型适应度不同,在进行景象匹配时,应根据主要的畸变类型选择适应度比较强的相似性度量。若对于各种畸变的影响综合衡量,相关度量优于距离度量,但相关度量的计算量大于距离度量,因此选择
图7 向量V1和V2在欧氏空间的关系图
相似性度量时还要考虑是否能达到实时性要求。参考文献:
[1] BROWNLG.Asurveyofimageregistrationtechniques[J].ACM
ComputerSurveys,1992,24(4):325-376.
[2] 赵锋伟,沈振康,李吉成.景象匹配技术研究[J].系统工程与电
如图7所示,若实时图Y(j,k)相当于向量V1,基准子图
X(u+j,v+k)相当于向量V2,则平方差度量D(u,v)相当于d,归一化积相关度量R(u,v)相当于θ。由此关系可以看出,
光照强弱只影响V1的模,从而影响d的变化,而不引起θ的变化,因此光照变化会引起距离度量的误匹配,而对相关度量没有影响。这与仿真结果相符。
3)相位相关度量对于云层具有较高的适应能力。
子技术,2002,24(12):111-114.
[3] YANGXG,MIAOD,CAOF.Studyonthematchingsimilarity
measuremethodforimagetargetrecognition[C]//The2ndInterna2tionalConferenceonFuzzySystemsandKnowledgeDiscovery.Ber2lin:Springer,2005:289-292.[4] ZITOVAB,FLUSSERJ.
Imageregistrationmethods:Asurvey
[J].ImageandVisionComputing,2007,21(11):977-1000.[5] 钟声伟,王玉珍,屈志毅.图像拼接中特征块的选取[J].计算
以上仿真都基于灰度的匹配算法,其中相位相关舍去了图像的能量信息,与图像的边缘密切相关,云层遮挡影响了实时图的灰度分布情况,对NProd、MSD和去均值MSD算法影响较大,但云层遮挡只是使实时图减小,其遮挡部分包含了很少的边缘信息,减小后的实时图包含的信息也足以使相位相
机工程,2002,28(9):165-167,232.
第30卷第1期
2010年1月
文章编号:1001-9081(2010)01-0098-03
计算机应用
JournalofComputerApplications
Vol.30No.1
Jan.2010
几种图像相似性度量的匹配性能比较
陈卫兵
(南通职业大学电子工程系,江苏南通226007)
(ntcwb@mail.ntvc.edu.cn)
摘 要:针对景象匹配中匹配性能和匹配实时性会受相似性度量选择影响的问题,从常用的相似性度量(归一化
积相关、相位相关、均平方差和去均值均平方差)入手,对图像施加噪声和进行畸变(如图像旋转变化、图像比例变化、光照强度变化和云层遮挡等),通过相似性度量匹配性能的仿真试验,从匹配性、适应度和实时性等方面对各相似性度量进行比较,并对结果进行了归纳总结与证明。
关键词:图像匹配;相似性度量;匹配性能中图分类号:TP391 文献标志码:A
Comparisonofmatchingcapabilitiesinsiilarityents
CHENWei2(DepartmentofElectronicEngineering,NNantongJ226007,China)
Abstract:Imagematchingreal2eonhowtochoosesimilaritymeasurementmethodofimagetotopicin2depth,thepapercarriedoutaseriesofsimulation
experimentsonoffoursimilaritymeasurements(NProd,PC,MSDandequalizationMSD)throughonimagedistortions.Theimagedistortionsincludepictureorientation,imagescalechange,illuminationintensitychange,cloudscover,andsoon.Theexperimentalresultswerecomparedwiththosefromthematchingperformance,adaptivecapacityandreal2timecapability.Analysisandtheoreticaldemonstrationoftheresultswerealsogiven.
Keywords:imagematching;similaritymeasurement;matchingperformance
0 引言
景象匹配技术在飞行器制导定位等领域一直是人们研究的热门课题,具有广泛的应用前景。研究可靠性高、鲁棒性好的图像匹配算法是提高匹配制导性能的有效途径。按照
Brown的理论
[1]
度量,两幅图像越相似,则距离越小。常用的距离度量有。
1)平均绝对差度量MAD。D(u,v)=
m×n
m
n
∑∑|
j=1k=1m
n
X(j+u,k+v)-Y(j,k)|
(1)
,用于景象匹配各种匹配算法都是4个元素
2)平均平方差度量MSD。D(u,v)=
m×n
的不同选择的组合:特征空间、相似性度量、搜索空间、搜索算法。其中相似性度量是指衡量匹配图像特性之间相似性的尺度,在景象匹配中占有重要的位置。基于各种相似性度量的算法研究不在少数,但对相应的适用场合没做详细的分析介绍,如何根据参与匹配的实际情况,选择合适的相似性度量,对有效地抑制噪声、畸变干扰,提高匹配概率与匹配精度有重要的意义。
本文通过对几种常用的相似性度量的匹配性能进行仿真并加以分析证明,给出了几种普遍存在的因素对相似性度量的影响以及相似性度量的适应性和实时性,为实际图像匹配应用中选择何种相似性度量提供了理论依据。
∑∑(X(j+u,k+v)
j=1k=1m
n
2
-Y(j,k))
(2)
3)去均值平均绝对差度量。D(u,v)=
m×n
∑∑(X(j+u,k+v)
j=1k=1
-(3)
X (u,v)-Y(j,k)+Y )
2
其中:X (u,v)为基准子图的平均值,Y 为实时图的平均值。
4)Hausdorff距离度量。
Hausdorff距离度量是适合于二值图像(如边缘特征图
像)相似性的度量。给定有限的两个点集A和B,A和B之间的
Hausdorff距离定义为:
H(A,B)=max[h(A,B),h(B,A)]
(4)
h(A,B)称有向
1 常用的相似性度量
图像的相似性度量有好多种,但应用广泛,相对成熟的相似性度量可以归结为两大类:距离度量和相关度量。设矩阵X为基准图矢量,大小为M×N,Y代表实时图矢量,大小为m×
n,则几种常见的相似性度量可表示如下。
其中:h(A,B)
Hausdorff距离
=max-min‖a-b‖,
a∈A
b∈B
[2]
。
1.2 相关度量
相关度量是定义在某种相关意义下的度量[3]。
1)积相关度量。
a)归一化积相关度量(NormalizedProduct,NProd)。
1.1 距离度量
距离度量是相比较的两幅图像定义在某种距离意义下的
收稿日期:2009-06-25;修回日期:2009-08-04。
基金项目:江苏省自然科学基金资助项目(BK2008188);江苏省六大人才高峰项目(072E2013)。 作者简介:陈卫兵(1966-),男,江苏南通人,副教授,硕士,主要研究方向:人工智能。
第1期
R(u,v)=
m
n
陈卫兵:几种图像相似性度量的匹配性能比较99
均值为0,因此仿真时对均值不为0的噪声不作考虑。
X(j+u,k+v)Y(j,k)
∑∑
对截取图添加均值为0,方差σ为0.
032、0.072、0.149、
[∑∑X
j=1k=1
mn
2
(j+u,k+v)
][∑∑Y
j=1k=1
1/2
mn
2
(j,k)
]
1/2
(5)
0.200产生的模拟实时图如图2所示。
b)去均值归一化积相关度量。
m
n
∑∑(X(j+u,k+v)
R(u,v)=
mn
j=1k=1
-X (u,v))Y(j,k)-X (u,v))
2
[∑∑(X(j+u,k+v)
m
n
]
1/2
×
图2 添加噪声的仿真用实时图
2.3.2 图像旋转变化的仿真
(6)
景象匹配导航中,由于惯性仪表的测量误差,拍摄的实时图与基准图总存在一定的角度偏差,这种偏差可用对截取图施加旋转变化来仿真,但是不能直接从截取图中旋转得到,因为截取图直接旋转会使截取图变大,。采取的措施是:;
以此位置为中心,,;对旋转后的图采用上述方法所得具有旋转畸变的实时图如图3所示。
[∑∑(Y(j,k)
j=1k=1
-Y )
2
]
1/2
2)相位相关度量(PhaseCorrelation,PC)。
设gr(x,y)为参考图,gs(x,y)为实时图,Gr(u,v)、Gs(u,
v)分别为它们的Fourier变换,令它们的互功率谱为Gr3Gs,
3
则相位相关函数R为:
R=F
-1
G33
|Gr3Gs=F
-1
[e
φr-φs)j(
]7]。其中:F-1表示Fourier,r2由于相似性度量较多,进行一一仿真工作量过大,因此选用最具代表性的四种相似性度量(归一化积相关、相位相关、均平方差,去均值均平方差)进行仿真研究。
2.1 选用的基准图及截取图
图3 图像旋转的仿真用实时图
为了便于计算匹配性能(匹配概率、匹配误差),仿真时,选用从基准图截取一小块图像,添加噪声及施加畸变的方法产生仿真用实时图[5]。选用的基准图为分辨率1m的某地区可见光卫星照片,大小为412像素×419像素,如图1(a)所示
,从中截取的图像,大小为152像素×157像素,如图1(b)所示。
2.3.3 图像比例变化的仿真
景象匹配导航中,由于测高系统存在测量误差,拍摄的实时图经过比例校正后,仍然存在一定大小的比例误差,这种误差可采用对截取图施加缩放来仿真,具体措施为:记下截取图的大小及其中心位置;以此位置为中心
,重新在基准图上截取一块较大的图像,以保证缩小后实时图大小的区域内没有空白区域;对该截取图缩放;对缩放后的图像中心截取实时图大小的一块图像即为畸变的实时图。采用上述方法所得具有缩放畸变的实时图如图4所示。
图1 仿真所用的基准图像及截取图像
2.2 选用的噪声及畸变模型
图4 具有缩放畸变的仿真实时图
为了能够尽可能地反映基准图与实时图的差异,仿真选用式(8)所示的噪声及畸变模型。
Y(j,k)=αX(u0+j+δ,k)j,v0+k+δk)+n(j
(8)
2.3.4 光照强度变化的仿真
拍摄基准图与实时图的时间的不同,天气不同,都会造成光照强度的变化,
表现为在基准图与实时图灰度值的变化,即式(8)中α的变化,α取不同值所得的仿真用实时图如图5所示。
其中:X为基准图矢量,Y为实时图矢量,基准图与实时图差异除了具有位置平移(u0,v0)外,还存在随机噪声n(j,k),灰度畸变α,以及比例缩放、旋转等畸变引起位置的微小变化δj、δk等。
2.3 噪声及各种畸变的仿真2.3.1 图像噪声的仿真
图像的噪声是以式(8)中n(j,k)的值来体现的,仿真时,取n(j,k)为均值是0,方差为σ的随机噪声。其实,若n(j,k)不服从均值为0,实时图经过去均值预处理后,噪声就会变为
图5 光照强度变化的仿真用实时图
2.3.5 云层遮挡的仿真
景象匹配导航中,拍摄实时图时,出现云层遮挡的情况也
100 计算机应用
2.4 仿真结果
第30卷
不可避免,可采取从截取图中去掉一部分图像来仿真云层遮挡的情况。所得的仿真用实时图如图6所示
。
对截取图分别施加噪声及畸变用不同相似性度量进行匹
配的结果如表1所示。其中:匹配误差,即匹配到的位置(u1,v1)与截取位置(u0,v0)之间的距离
ε=(u0-u1)2+(v0-v1)2
度量值为匹配点的度量值,对于相关度量为相关曲面的峰值,对于距离度量为距离曲面的峪值。有下划线的结果表示明显的误匹配。
表1 实验结果
图6 云层遮挡的仿真用实时图
畸变
类型
匹配误差
畸变值
0.032
度量值
去均值
MSD0.000000.000000.[1**********]05.303.162302.236101.000001.000000.00000NProd0.967520.936620.897430.0.0.0.440.927603.162300.944660.924560.993840.993840.993840.937850.938870.93887
实时图图图图图图图图图图图图图图图
2(a)2(b)2(c)2(d)3(a)3(b)3(c)3(d)4)4(b4(c)5(a)5(b)5(c)6(a)6(b)6(c)
NProd0.000000.000000.000000.000001.000003.605605.385208.162100.000000.000000.000000.000000.00000相位相关
0.000000.000000.000000.000003.162304.472105.[1**********].770308.602301.000001.000001.000000.000000.000000.00000
MSD0.000000.000000.000000.000001.000003.605303.302.236101.000000.000000.00000相位相关
0.123320.094470.0680.0.390.044220.024630.017470.013680.015660.151130.151130.151130.092990.080740.08074
MSD(×107)1.555403.[***********].421302.701803.336602.513902.618603.809800.222510.293940.370110.175220.208730.20873
去均值MSD
(×107)1.555303.046705.019706.088601.846602.421102.701703.336502.506402.611003.784100.488770.935741.858901.565200.181570.18152
添加
噪声
0.0720.1490.2003.02°5.02°7.02°10.00°0.915
旋转
缩放1.1.0.70550.50600.37689%13%20%
光照变化云层遮挡
3 分析与证明
通过分析仿真结果得到如下结论。
1)添加噪声不会影响匹配性能。
从仿真结果可以看出,添加噪声后利用各种相似性度量匹配结果的误差都为0,即使噪声添加到使图像人眼都分辨不出来,匹配误差仍然为0。
mn
(X(j+u,k+v)-Y(j,k))2= D(u,v)=∑∑m×nj=1k=1
m×nj=1km×nj=1km×nj=1km×nj=1k
mmm
证明 基准子图和实时图只存在噪声,则畸变模型简化为如下所示:
(9)Y(j,k)=X(u0+j,v0+k)+n(j,k)
为了证明方便,设n(j,k)是服从均值为0,方差为δ的噪声,当n(j,k)的均值不为0,可以通过对实时图去均值处理,使n(j,k)是服从均值为0。
对于平方差度量:
∑∑(X(j+u,k
=1n
mn
2
+v)-(X(u0+j,v0+k)+n(j,k)))=
∑∑(X(j+u,k
=1n
+v)-X(u0+j,v0+k))
2
+
m×n
j=1k=1
∑∑n(j,k)
mn
2
-
∑∑(X(j+u,k
=1n
+v)-X(u0+j,v0+k))n(j,k)=
22
+v)-X(u0+j,v0+k))+δ-0(.)
∑∑(X(j+u,k
=1m
n
(10)
其中0(.)表示近似为0,因此D(u,v)与噪声无关。而对于归一化积相关有:
R(u,v)=
m
n
∑∑X(j+u,k
2
+v)[X(u0+j,v0+k)+n(j,k)]
m
n
[∑∑X
j=1k=1
(j+u,k+v)
m
n
][∑∑[X(u
j=1k=1
1/2
+j,v0+k)+n(j,k)]
2
]]
1/2
=
∑∑X(j+u,k
(j+u,k+v)
1/2
+v)X(u0+j,v0+k)+0(.)
m
n
[∑∑X
j=1k=1
mn
2
][∑∑[X(u
j=1k=1
+j,v0+k)+n(j,k)]
2
1/2
(11)
(下转第110页)
110 计算机应用第30卷
精确匹配,匹配的成功率得到了提高,由此,将更多的资源纳入为任务选择的范围内,使得系统在吞吐量、任务满足率和资源调度率上均有所提高。就失败率来看,供求基本平衡和供不应求两种情况下分别增长了0.35%和0.80%,这和其他参数的增长相比较,增长幅度甚微,是可承受的。它反映出用户在QoS参数的理解上的模糊性和不确定性,导致一般在定义中提出略高的要求,因此,软化几乎对QoS没有影响。
最大匹配算法从两个方面取得了进展:1)在满足QoS需求的基础上,适度扩大了资源选择范围,减少了因为QoS参数而导致的匹配问题,使得资源调度率和任务满足率提高;
2)确定了依照QoS选择资源的优先顺序,为任务选择了适当
证明,对QoS的软化处理提高了系统吞吐量和资源匹配成功率,最终实现了整个系统资源调度率的提高。参考文献:
[1] ALOISIOG,CAFAROM,EPICOCOI,etal.Resourceandservice
discoveryintheiGridinformationservice[EB/OL].[2009-06-15].http://www.gridlab.org/WorkPackages/wp210/Documents/igrid2iccsa.pdf.
[2] BRADLEYA,CURRANNK,PARRZG.Resourcediscoveryand
managementincomputationalGRIDenvironments[J].InternationalJournalofCommunicationSystems,2007,19(6):639-657.[3] CHENP,XUZ,ZHANGB.AsolutiontoQoScontrolandavaila2
bilitypromotionincomplexgridcomputing[C]//Proceedingsof12thIEEEInternationalConferenceonNetworks.NewYork:IEEE,2004:403-407.
[4] LICHUNLIN,LILAYUAN.UtilitymultipleQoSguaranteed
resourceoptimonmputing[C]//ProceedingstheIonalCoTheoryandApplica2tionson,ICoSociety,2009:165-169.R,D,GIDDYJ.Nimrod/G:Anarchitecture
ormanagementandschedulingsysteminaglobalcompu2tationalgrid[C]//Proceedingsofthe4thInternationalConference/ExhibitiononHighPerformanceComputingintheAsia2PacificRe2gion.NewYork:IEEE,2000:283-289.
[6] PLESTYSR,VILUTISG,SANDONAVICIUSD.Themeasurement
ofgridQoSparameters[C]//Proceedingsof29thInternationalCon2ferenceonInformationTechnologyInterfaces.NewYork:IEEE,2007:703-707.
的资源,而非当前网格中最好的资源,使系统工作在最优状态,由此进一步提高系统的可用性。与传统经典算法比较,最大匹配算法在任务一级取得了次优解,但在系统一级得到更优解。算法有利于克服网格中因自私性带来的缺点,在用户和系统效率之间取得折中,更有利于系统的稳定和健壮,尤其适用于资源供不应求状态下。
4 结语
,按照参数性质,。针对弹性参数,进行了软化处理,使其在不失定性含义的前提下,扩展了资源选择的范围,提高了任务和资源在QoS参数上的匹配率,由此提高系统的吞吐量,实现网格系统的调度优化。实验
(上接第100页)
可见式(11)的R(u,v)与式(6)具有相同位置的峰值点,即添加噪声也不影响NProd度量的匹配概率。
2)相关度量对于光照强弱的变化不敏感,而距离度量却比较敏感。
实时图和基准子图的相似程度可以直观地从两图的图像矢量在欧氏空间的几何关系得出来
。
关算法正确找到匹配位置。
4)对于旋转、缩放等畸变,几种相似性度量都比较敏感。
对实时图施加旋转、缩放畸变,各种相似性度量都不同程度地出现了匹配误差,甚至误匹配,可见旋转、缩放畸变是影响匹配概率的主要因素,因此要尽量消除旋转、缩放畸变后再进行匹配,才能取得较好的匹配性能。
4 结语
不同的相似性度量对于不同的畸变类型适应度不同,在进行景象匹配时,应根据主要的畸变类型选择适应度比较强的相似性度量。若对于各种畸变的影响综合衡量,相关度量优于距离度量,但相关度量的计算量大于距离度量,因此选择
图7 向量V1和V2在欧氏空间的关系图
相似性度量时还要考虑是否能达到实时性要求。参考文献:
[1] BROWNLG.Asurveyofimageregistrationtechniques[J].ACM
ComputerSurveys,1992,24(4):325-376.
[2] 赵锋伟,沈振康,李吉成.景象匹配技术研究[J].系统工程与电
如图7所示,若实时图Y(j,k)相当于向量V1,基准子图
X(u+j,v+k)相当于向量V2,则平方差度量D(u,v)相当于d,归一化积相关度量R(u,v)相当于θ。由此关系可以看出,
光照强弱只影响V1的模,从而影响d的变化,而不引起θ的变化,因此光照变化会引起距离度量的误匹配,而对相关度量没有影响。这与仿真结果相符。
3)相位相关度量对于云层具有较高的适应能力。
子技术,2002,24(12):111-114.
[3] YANGXG,MIAOD,CAOF.Studyonthematchingsimilarity
measuremethodforimagetargetrecognition[C]//The2ndInterna2tionalConferenceonFuzzySystemsandKnowledgeDiscovery.Ber2lin:Springer,2005:289-292.[4] ZITOVAB,FLUSSERJ.
Imageregistrationmethods:Asurvey
[J].ImageandVisionComputing,2007,21(11):977-1000.[5] 钟声伟,王玉珍,屈志毅.图像拼接中特征块的选取[J].计算
以上仿真都基于灰度的匹配算法,其中相位相关舍去了图像的能量信息,与图像的边缘密切相关,云层遮挡影响了实时图的灰度分布情况,对NProd、MSD和去均值MSD算法影响较大,但云层遮挡只是使实时图减小,其遮挡部分包含了很少的边缘信息,减小后的实时图包含的信息也足以使相位相
机工程,2002,28(9):165-167,232.