基于广义似然比的多航过sar相干变化检测方法
【专利摘要】本发明公开了一种基于广义似然比的多航过SAR相干变化检测方法，包括以下步骤：S1：选取多航过SAR图像对记为{f1,f2,…,fK}；S2：选取多航过SAR图像像素对；S3：对协方差矩阵进行最大似然估计；S4：进行似然比假设检验；S5：依次选取多航过SAR图像像素对，重复步骤S4，得到变化检测结果。本发明假设SAR成像区域在多幅图像采集期间发生变化和未发生变化时其对应的复像素对分别服从不同的圆对称复高斯分布，对假设中圆对称复高斯分布的协方差矩阵进行估计，然后确定检测统计量，并将该检测统计量和门限比较，检验上述两个假设的成立与否，也就是检测成像区域有无变化，可以实现微弱变化检测和变化过程的观测。
【专利说明】基于广义似然比的多航过SAR相干变化检测方法
【技术领域】
[0001]本发明属于合成孔径雷达(SAR)【技术领域】，特别涉及一种基于广义似然比的多航过SAR相干变化检测方法。
【背景技术】
[0002]变化检测技术可以广泛应用于监测森林植被、土壤水分、积雪量的变化；监测农作物的长势变化、土地覆盖的变化；监测各种灾害发生前后的变化，如地震区域的定位和灾害评估；监测海冰的运动、山岳冰川的移动以及滑坡运动；军事目标区域的动态监测、战场打击评估等。但是在云、雾等恶劣天气条件下难以通过光学图像获取变化信息。由于SAR是一种全天时、全天候的现代高分辨率微波遥感成像雷达，通过对获取的SAR图像进行变化检测，能够在第一时间为决策提供有力的支持。然而当目标变化微弱时，传统SAR图像变化检测方法检测不到变化。另外，受森林、海浪等微动杂波的影响，传统SAR图像变化检测方法使目标的变化淹没在周围背景的变化中，区分不出目标变化与背景变化的差异，使得变化检测性能急剧下降。
[0003]多航过SAR是目前微波遥感技术研究与应用的一个重要领域，它深入挖掘信息有效地提高了雷达对地物信息的获取能力，对多航过雷达遥感图像进行信息处理的研究具有重要的理论价值和广阔的应用前景。
[0004]近年来，国外很多学者开始探索利用多航过极化微波遥感获取同一地区的多幅图像，经过信号处理检测目标微弱变化。在文献“Leslie M Novak, Change detect1nfor mult1-polarizat1n, mult1-pass SAR, in Defense and Security.1nternat1nalSociety for Optics and Photonics, 2005，pp.234 - 246” 中提出了一种基于广义似然比的变化检测方法，该方法利用多航过图像的幅度信息提高了变化检测性能，但是此非相干变化检测方法通过场景的后向散射功率的变化来进行变化检测，没有利用相位信息，对SAR图像变化的灵敏度低。然而雷达系统是相干系统，回波的幅度和相位对目标的变化都非常敏感。单纯利用幅度信息进行变化检测的方法丢失了大量的有用的相位信息已经不能满足精确变化检测需求。在文献 “Mark Preiss and Nicholas JS Stacy, Coherent ChangeDetect1n !Theoretical Descript1n and Experimental Results, Tech.Rep, DTICDocument, 2006”中提出了一种相干变化检测方法，该方法同时利用图像的幅度和相位信息，能够检测图像微弱变化。但是该方法只利用两幅图像的信息，没有利用多幅图像的时相信息。
【发明内容】
[0005]本发明的目的在于克服两航过非相干变化检测方法中只利用幅度信息以及仅利用两航过SAR图像信息导致变化检测概率低或检测不到变化的问题，提供一种假设SAR成像区域在多幅图像采集期间发生变化和未发生变化时其对应的复像素对分别服从不同的圆对称复高斯分布，再确定检测统计量并与门限比较以检验上述两个假设的成立与否的基于广义似然比的多航过SAR相干变化检测方法。
[0006]本发明的目的是通过以下技术方案来实现的:基于广义似然比的多航过SAR相干变化检测方法，包括以下步骤:
[0007]S1:选取多航过SAR图像对:选取不同时间对同一地区多次观测获得的K幅SAR图像并进行配准后记为{f\，f2,…，fK}，K彡3 ;
[0008]S2:选取多航过SAR图像像素对，依次选取f\，f2,…，fK对应的复矩阵第m行第η
列的值，并记为向量，假设向量x_?CN(0，Γ)，其中Γ =E [XfflnXj]
为协方差矩阵，设H0表示目标区域未发生变化，得到协方差矩阵Γ 0,设H1表示目标区域发生变化，得到协方差矩阵F1 ;
[0009]S3:对协方差矩阵进行最大似然估计:分别对协方差矩阵Γ ^和Γ i进行最大似然估计；
[0010]S4:进行似然比假设检验；
[0011]S5:依次选取多航过SAR图像像素对，重复步骤S4，得到变化检测结果。
[0012]进一步地，所述的步骤S2的具体实现方法为:由于图像存储为矩阵形式，因此f1； f2,…，fK均为复图像，依次选取f\，f2,…，fK对应的复矩阵第m行第η列的值分别记为令复向量…,/!Γ, m = 1，2，…，M，n= 1，2，...，N，其中[]τ
表示转置运算，Μ、N为图像对应的复矩阵的大小；
[0013]设向量Xmn服从K维圆对称复高斯分布，即X--CN(0，Γ)，其中协方差矩阵Γ =E [XnmXnmaI，因此Xmn的概率密度函数表示为:
[0014]P(KJ = pJj^expC—xw/rx,f/)
[0015]其中，E [XmnXmn11]表示求Hh的均值，X*表示求Xmn的复共轭转置，I Γ I表示Γ的行列式，exp表示指数运算；
[0016]设H。表示目标区域未发生变化，此时对应的向量Xmn?CN(0，Γ。)，其中，
[0017]
彳Λ:σ,σ:6χρΗΦ12)…P,^σ,σ, CxpC^Dljr)]
P _ ACr1CT2 exp(jcl).0CT22…/λΑσ:σ人 exp(-」Φ2【)
Aq-....ρ
?λm.I
exp{j#!K) —exp{j.M)… σ/ J
[0018]其中，j表示复数单元，P ab ?1，Oab ?0° , a = I, 2,..., K, b = 1,2,…，K, a
<b ；
[0019]设H1表示目标区域发生变化，此时对应的向量Xmn?CN(0，Γ\)，其中，
[0020]
σ,2cxpH#；2)…/^擎火―expH<D?A—)
P = #12σ1σ；6χρ(]φ,12)σ；2…Ρ2?σ；σΑ； εχρ(-]Φ'2Α.)，
>iK<￥4exp(ja)lir)p'^cr,CrctCXiXjay2J…σΛ2观
[0021]p' ?~0，Φ' ab#0°，由于Γ。、未知，因此为复合假设。
[0022]进一步地，所述的步骤S3对协方差矩阵Γ ^和Γ i进行最大似然估计的具体方法为:
[0023]S31:对Γ ^进行最大似然估计，根据最大似然估计理论，对于假设Htl，根据多航过SAR图像对{f\，f2,…，fK}中明显未变化区域得到参数为Γ 0的复高斯随机向量Y的L个相互独立的观测Yi,每个观测Yi的概率密度函数为P (Yi Γ0), Γ0的似然函数Iik ( Γ 0)表示为:
[_ 級" !
[0025]令备Uk(T0) = O,求得协方差矩阵Γ 0的最大似然估计:
O
γ?ΙΠI\7 H
[0026]F0 = — / ,XX I
[0027]S32:对进行最大似然估计，对于假设H1，根据多航过SAR图像对{f\，f2,…，fK}中明显变化区域得到的参数为Γ I的复高斯随机向量Y的S个相互独立的观测Yi,每个观测1的概率密度函数为P(YiI Γ\)，的似然函数IiHr1)可以表示为:
[_ 勝(Γ|)=如(Xir1)=^jexpI —'
Ji/ I Jl 11
[0029]令最mj\、=Q，求出协方差矩阵r i的最大似然估计:
Λ I丄
Γλλολ? F 一 VVV;/
1_0030」11; * j * j ◎
Λ I I
[0031]进一步地，所述的步骤S4包括以下子步骤:
[0032]S41:确定检验统计量，根据步骤S2的假设和步骤S3的参数估计得到:
[0033]PiXmn I H(1) =」^exp〔-Χ?,/ Γ() 1X臓
Zjr0I VJ
? ? \
IΛ ι
r?W ! U \ —V HY
[0034]Ρ\ I ** I /— exp ?人腿置 I
? Γ, j VJ
[0035]取K幅多航过SAR图像中第m行第η列像素点的Q个相互独立的像素对XK?，...,X匕，根据似然比假设检验理论，令似然比:
_p(x；Mn,XL,-XL；H0) Qp(^mnIH0)
[0036]Z=..........——。..........:；..( =—TTT
ρ(Κ?’χ-5■.-xL;Ηι) ρ(χ,ιιπιΗι)
[0037]将PdI H?)、PiXm IH1)代入上式化简得，
[0038]eXp^Tr|r? ^Γ^￡χ；?,χ；?； jj
[0039]取对数并忽略常数项，得检测统计量:
If.' -S Λ..?] 卜)
[0040]Z = Trj^rfl -T]hTr|rdG|
[0041]其中TrfdG丨表示矩阵G的迹，也就是矩阵对角元素的和，
Α Λ "1 Λ "1 夺..η
rd=r0 -F1 ,G=mL.>
[0042]S42:进行假设检验，定义一个和用来检测的图像一样大小的矩阵，记为R，选择门限T，判断门限T与检测统计量Z的大小:当Z>T时，判定假设H1成立，变化检测结果就是该像素对应的区域发生变化，令R对应的像素值为255 ;否则判定假设Htl成立，变化检测结果就是该像素对应的区域未发生变化，令R对应的像素值为O。
[0043]进一步地，所述的步骤S5中得到变化检测结果的方法为:依次选取多航过SAR图像像素对，重复步骤S4，直到确定矩阵R中所有的像素值，R即为变化检测结果。
[0044]本发明的有益效果是:
[0045]本发明利用广义似然比假设检验理论，假设SAR成像区域在多幅图像采集期间发生变化和未发生变化时其对应的复像素对分别服从不同的圆对称复高斯分布，对假设中圆对称复高斯分布的协方差矩阵进行估计，然后确定检测统计量，并将该检测统计量和门限比较，检验上述两个假设的成立与否，也就是检测成像区域有无变化。
[0046]本发明的检测方法充分利用多航过SAR复图像丰富的幅度和相位信息，与现有的两航过相干变化检测方法和多航过非相干变化检测方法相比，可以实现微弱变化检测和变化过程的观测，在地球遥感和地质灾害监测等领域十分适用。
【专利附图】
【附图说明】
[0047]图1为本发明的检测方法的流程图；
[0048]图2为本发明实施例采用的某一地区的SAR图像；
[0049]图3为本发明实施例采用的某一地区的SAR图像f2 ；
[0050]图4为本发明实施例采用的某一地区的SAR图像f3 ；
[0051]图5为本发明实施例变化检测结果图像；
[0052]图6为本发明实施例不同航过相干变化检测的ROC曲线。
【具体实施方式】
[0053]下面结合附图进一步说明本发明的技术方案，但本发明所保护的内容不局限于以下所述。
[0054]为了方便描述本发明的内容，首先作以下解释:
[0055]1、复高斯分布
[0056]a)复高斯分布的定义
[0057]假设X和Y是k维实空间的随机向量，向量vect[X Y]是2k维正态随机向量。那么以X为实部、Y为虚部的复随机向量Z = X+Yj具有复高斯分布，记为Z?CN(l.!，r，C)。
[0058]μ = Ε[Ζ], Γ = Ε[(Ζ-μ ) (Zl)H]，C = Ε[(Ζ_μ) (Ζ_ μ )，]，
[0059]其中:j是复数单元，S卩一I开根，E[Z]表示求Z的均值，均值μ可以是任意k维的复矢量；Zh表示求Z的复共轭转置，协方差矩阵Γ必须是厄尔米特的和非负定的；Z’表示求Z矩阵转置，相关矩阵C是对称的。
[0060]b)圆对称复高斯分布
[0061]圆对称复高斯分布对应的参数为μ = 0，C = O。若k维随机向量Z = X+iY服从圆对称复高斯分布，通常写做Z?CN(0，Γ)，其概率密度函数为:
[0062]/(Z)=exp-,r,z
[0063]其中I Γ I表示Γ的行列式，exp表示指数运算，Γ-1表示Γ的逆。
[0064]2、复合假设
[0065]虚无与对立假设中不只包含一个母数值。在一般的假设检验中还存在一个未知的参数。
[0066]3、最大似然估计
[0067]假设得到参数α的随机变量y的m个观测。每个观测Ji有pdfp (Yi | α )，并且由于m个观测是独立的，联合pdfppla)称之为似然函数，可以表示为
m
[0068]H, I gO=EbU IcO
1-Ι
[0069]其中Π表示连乘运算。那么通过求似然函数的最大值可以确定参数α的最大似然估计?Μο计算时，这个最大值通过令似然函数的导数为零求得，也就是
[0070]= 0
[0071]由于对数函数是单调函数，所以常用一个更简单的过程来求似然函数对数的最大值(称为对数-似然)，也就是
Q—
[0072]—bip( v|
da
[0073]最大似然估计的前提是似然函数导数存在。
[0074]4、虚警概率、检测概率
[0075]本发明中虚警概率是指场景未发生变化但被检测为变化的概率。检测概率指场景发生变化同时该变化也被检测出的概率。
[0076]5、ROC 曲线
[0077]通过改变门限得到多组检测概率和虚警概率，以虚警概率为横坐标，检测概率为纵坐标，根据这些点连成的曲线称为ROC曲线。
[0078]本发明提供一种基于广义似然比的多航过SAR相干变化检测方法，如图1所示，包括以下步骤:
[0079]S1:选取多航过SAR图像对:选取不同时间对同一地区多次观测获得的K幅SAR图像并进行配准后记为{f\，f2,…，fK}，K彡3 ;
[0080]S2:选取多航过SAR图像像素对，依次选取f\，f2,…，fK对应的复矩阵第m行第η
列的值，并记为向量假设向量Xnn?CN(0，Γ)，其中Γ =E[XmnXnm11]
为协方差矩阵，设H0表示目标区域未发生变化，得到协方差矩阵Γ 0,设H1表示目标区域发生变化，得到协方差矩阵F1 ;
[0081]S3:对协方差矩阵进行最大似然估计:分别对协方差矩阵Γ ^和Γ i进行最大似然估计；
[0082]S4:进行似然比假设检验；
[0083]S5:依次选取多航过SAR图像像素对，重复步骤S4，得到变化检测结果。
[0084]进一步地，所述的步骤S2的具体实现方法为:由于图像存储为矩阵形式，因此f1； f2,…，fK均为复图像，依次选取f\，f2,…，fK对应的复矩阵第m行第η列的值分别记为
fLL.'.?，/I，令复向量 I =IfLfL-.-JiY，m= 1，2，…，Μ，η= 1，2，…，N，其中[]τ
表示转置运算，Μ、N为图像对应的复矩阵的大小；
[0085]设向量Xmn服从K维圆对称复高斯分布，即X--CN(0，Γ)，其中协方差矩阵Γ =E [XnmXnmaI，因此Xmn的概率密度函数表示为:
[0086]PQimn) = ^ αφ(-Χ?/ΓΧ?/)
[0087]其中，E [XmnXmn11]表示求U?11的均值，X*表示求L的复共轭转置，I Γ I表示Γ的行列式，exp表示指数运算；
[0088]设Htl表示目标区域未发生变化，此时对应的向量Xmn?CN(0，Γ^)，其中，
[0089]
Γ Cr,2P,2^,^2exp(-j<I)1;)…p;A.q?exp(—j<l)IA.)
r —…AA.A%exp(—j(U
夏 0 _....，
Iama*
P—xp(j_2[)…σ/^
[0090]其中，j表示复数单元，P ab ?1，Oab ?0° ,a= 1,2,..., K, b = 1,2,…，K, a
<b ；
[0091]设H1表示目标区域发生变化，此时对应的向量Xmn?CN(0，Γ\)，其中，
[0092]
σ.,'2Ρησ?σ2…Ρικσισκ ?Ρ(-]Φ'ιΚ)
Γ — AAA.exp(jtu'i2)σ22- ρ2κσ2σκ1Κ)
I I """..<.,
__/>]Κσ,σλ.cxp(jcI)',A) ρ乂σ,σλ..οφ(.?Φ2Α.)…σ'/冊
[0093]p' ?~0，Φ' ab#0°，由于Γ。、未知，因此为复合假设。
[0094]进一步地，所述的步骤S3对协方差矩阵Γ ^和Γ i进行最大似然估计的具体方法为:
[0095]S31:对Γ ^进行最大似然估计，根据最大似然估计理论，对于假设Htl，根据多航过SAR图像对{f\，f2,…，fK}中明显未变化区域得到参数为Γ 0的复高斯随机向量Y的L个相互独立的观测Yi,每个观测Yi的概率密度函数为P (Yi Γ0), Γ0的似然函数Iik ( Γ 0)表示为:
/
/-..1 扣丨'<,'、；
[0096]Μ(Γ(,) = Π p(Y j TΛ=.-rexp 1'
0 ,.1 0； ^|r0f- F
[0097]令iMlT(l) = 0，求得协方差矩阵Γ ^的最大似然估计:
[0098]r0 =^Xx-Y/'；
L.?::.ι
[0099]S32:对进行最大似然估计，对于假设H1，根据多航过SAR图像对{f\，f2,…，fK}中明显变化区域得到的参数为Γ I的复高斯随机向量Y的S个相互独立的观测Yi,每个观测1的概率密度函数为P(YiI Γ\)，的似然函数IiHr1)可以表示为:
S
Iy" I': γ
[0100]Iik(Tl) = Up(Xiri) =......」.一rexp ' ':'
' H " I^s |Γ I
[0101]令€^0 = 0，求出协方差矩阵^的最大似然估计:
[0102]T1β
^ ? I
[0103]进一步地，所述的步骤S4包括以下子步骤:
[0104]S41:确定检验统计量，根据步骤S2的假设和步骤S3的参数估计得到:
J?/s 4-、
[0105]ρ(\,? IH0)=——Γ...............exp Γ0 Xmil
I r01 VJ
IfA\
r I η--? i O \1 PWI Y H Y
[0106]PK^rm I r^l/ ~~ " Λ —“P 画 * I ^tImi
TTA丨[丨 VJ
[0107]取K幅多航过SAR图像中第m行第n列像素点的Q个相互独立的像素对Κ?根据似然比假设检验理论，令似然比:
r ^ p{K?,^L...ep(x:w；e0)
[0108]Z = ^I—-T-J:-f = 1 I^7^f
[0109]#i?(X；,JH0), |^1醫1)代入上式化简得，
[0110]expi^Trlir；, ^f1
Unii I IvI)
[0111]取对数并忽略常数项，得检测统计量:
if A -l ΛI? Λ )
[0112]Z = Tr|r0 ^ri J￡XL?X；J'J1 = TrjriiGj
[0113]其中TrfdG}表示矩阵G的迹，也就是矩阵对角元素的和，
ΛA -1 Λ —1(>
Γ』=Γη -Γ, G = Yxi Xi Η ；
d OImn mi 5
i=l
[0114]S42:进行假设检验，定义一个和用来检测的图像一样大小的矩阵，记为R，选择门限T，判断门限T与检测统计量Z的大小:当Z>T时，判定假设H1成立，变化检测结果就是该像素对应的区域发生变化，令R对应的像素值为255 ;否则判定假设Htl成立，变化检测结果就是该像素对应的区域未发生变化，令R对应的像素值为O。
[0115]进一步地，所述的步骤S5中得到变化检测结果的方法为:依次选取多航过SAR图像像素对，重复步骤S4，直到确定矩阵R中所有的像素值，R即为变化检测结果。
[0116]本发明主要采用计算机仿真的方法进行验证，所有步骤、结论都在MATLAB_R2013a上验证正确，下面结合具体实施例进一步说明本发明的技术方案:本实施方式中以三航过变化检测为例，首先选取不同时间对同一地区多次观测获得的3幅SAR图像，如图2、图3和图4所示，经过配准后，记为{f\，f2, f3}，对比可知f2在图中黑色圆圈所表示区域发生较小的变化，对比可知f3在图中黑色椭圆所表示区域发生较大范围的变化。
[0117]依次选取f\，f2,&复图像对应的复矩阵第m行第η列的值分别记为/二./:,，./!，令复向量=ULJhCJ，m= I,2,-,451,η= 1，2，…，451，假设向量Xmn服从3维圆对称复高斯分布，Xmn的概率密度函数可表示为:
[0118]P(Xtfm) 二Χ,?,"ΓΧ?/)
[0119]假设Htl为目标区域未发生变化，此时对应的向量Xmn?CN (O, Γ 0)，假设H1为目标区域发生变化，此时对应的向量Xmn?CN(0，Γ\)。
[0120]对进行最大似然估计:根据最大似然估计理论，对于假设Htl，根据{f\，f2, f3}中明显未变化区域我们得到参数为Γ ^的随机向量Y的1681个相互独立的观测，每个观测Yi的概率密度函数为P (Yi I Ftl), q的似然函数IiHrtl)可以表示为:
Ifii
_ t.? -!>/% ?,
[0121]級11.气AI I鑛
[0122]令=0，求得协方差矩阵Γ ^的最大似然估计:
(.'Γ"Iuol ? ?
[0123]对Γ i进行最大似然估计:对于假设H1，根据{f\，f2, f3}中明显变化区域我们得到的参数为F1的随机向量X的121个相互独立的观测，每个观测Yi的概率密度函数为P(YiI T1), q的似然函数IiHr1)可以表示为:
|:丨.、 I次
[0124]lik(T^ = Up(YiITi) = -^7—H-exp.:
,=l^ IrlI
[0125]令^rtt(Ti)= Oi求出协方差矩阵r:的最大似然估计:€ =^TY1Y/1 O
4121
[0126]取3幅多航过SAR图像中第m行第η列像素点的9个相互独立的像素对xL?xL,..-Xm,根据似然比假设检验理论，令似然比:
p(xL,Xi?-Xl；H0) vp(x；w,(；H0)
[0128]将P(XLJH0)'P(XLJH1)代人上式得,
Z= Φ Cxpi^Trff0 1^fl 1IxxLxLn}]取对数并忽略常数项，得检测统计量:
U rO U V [V 知 JJ
=其中G = Xhh0
[0129]定义一个451 X451的矩阵，记为R。选择门限T = 0，当Z>0时，判定假设H1成立，令R对应的像素值为255，否则判定假设Htl成立，令R对应的像素值为O。
[0130]依次选取像素对，重复上述操作，直到确定矩阵R中所有的像素值，得到如图5所示的变化检测结果图像。
[0131]下面对本发明得到的变化检测方法进行性能分析:
[0132](I)计算一组检测概率、虚警概率:分别生成25个二维圆对称复高斯随机变量蒽、育’和三维圆对称复高斯随机变量X?、育'，其中Xf ?GV(^r20) , ~G￥(0，r’2i) , i = 1,2,…，25，
v., ,2.2686 0.45-2.2686]
刚 X'NCW(0,r30), X-CW(OXm), r2(i=10^45-12686 2.涵—，
|nf2.2686 O —
[0134]Γ? =10?
21 L O 0.95070
[0135]
—2.26860.45-2.2686 0.45-2.2686][2.2686 O O "
Γχ = 10 0.45-2.2686 2.2686 0.45-2.2686 ? Γ；0=10 O 1,0847 O 0
^0.45-2.26860.45-2.2686 2.2686 J[O O 0.95070—
[0136]可以看出Γ ’ 2(|，Γ ’ 3(|满足步骤2中的Htl假设，Γ ’ 21，Γ ’ 31满足步骤2中的H1假设。令门限T = O,由计算检测统计量Zf ,如果大于门限，则记为虚警。由1If计算检测统计量Zf.如果大于门限，则记为正确检测。同样地，由Xf计算检测统计量如果Zl:大于门限，则记为虚警。由 <'计算检测统计量Z#,如果大于门限，则记为正确检测。重复10000次上述步骤，统计虚警次数和正确检测的次数并分别除以10000，即得门限为O时两航过和三航过变化检测的虚警概率和检测概率。
[0137](2)用蒙特卡洛仿真方法求ROC曲线:在[-200，200]以I为间隔均匀改变步骤(I)中门限的值，计算401组检测概率和虚警概率。在以虚警概率为横轴，检测概率为纵轴的坐标系中，每个门限求得的虚警概率和检测概率对应于坐标轴中的点，用光滑曲线连接这些点得到R0C，如图6所示，可以看出和两航过相比，三航过相干变化检测在相同的虚警概率的条件下，检测概率明显提高。
【权利要求】
1.基于广义似然比的多航过SAR相干变化检测方法，其特征在于:包括以下步骤: 51:选取多航过SAR图像对:选取不同时间对同一地区多次观测获得的K幅SAR图像并进行配准后记为If1, f2,…，fK}, K彡3 ;52:选取多航过SAR图像像素对，依次选取f\，f2,…，fK对应的复矩阵第m行第η列的值，并记为向量Γ)，其中Γ =E[XmnXmnH]为协方差矩阵，设Htl表示目标区域未发生变化，得到协方差矩阵Γ 0,设H1表示目标区域发生变化，得到协方差矩阵Γ1; 53:对协方差矩阵进行最大似然估计:分别对协方差矩阵和^进行最大似然估计； 54:进行似然比假设检验； 55:依次选取多航过SAR图像像素对，重复步骤S4，得到变化检测结果。
2.根据权利要求1所述的基于广义似然比的多航过SAR相干变化检测方法，其特征在于:所述的步骤S2的具体实现方法为:由于图像存储为矩阵形式，因此f\，f2,…，fK均为复图像，依次选取f\，f2,…，&对应的复矩阵第m行第η列的值分别记为/K,,...，./=，令复向量 Κ/1 -..，/:]' m = 1，2，…，Μ，η = 1,2,…，N，其中[]Τ 表示转置运算，Μ、N为图像对应的复矩阵的大小； 设向量Xmn服从K维圆对称复高斯分布，即Χ_?0Ν(0，Γ)，其中协方差矩阵Γ =E [XnmXnmaI，因此Xmn的概率密度函数表示为:
= -TfH exp(-C:)
M 11 I 其中，E [XmnXmn11]表示求XmnXmn11的均值，Xmn11表示求复共轭转置，I Γ I表示Γ的行列式，exp表示指数运算； 设Htl表示目标区域未发生变化，此时对应的向量Xnm?CN(O，Γ J，其中，
σ,2P12^1Ct2 exp(-j#12)…ρικσ,σκ exp(-j#1A.)—Α2σ,σ2βχρ(]Φ)2)σ,2…/?2κσ2σκ exp(-jd>2K)I O 画'■..■tjA-σ.σΑ-exp(jflv) ρ2κστσκ exp(jfU2K)…σ/釋 其中，j 表示复数单元，P ab ?I，Oab ?O。，a= 1,2,..., K, b = 1,2,...,K, a<b; 设H1表示目标区域发生变化，此时对应的向量Xnm?CN(0，Γ J，其中，
cxpC-j#'I2)…exp(-沈)Γ = Pη)…/?'2Α.σ2σΑ.'6χρΗΦ'2Α.)y?rexp(j€)?)p?rexp(j.?)…遍 P' ab ^ O, Φ / ab# 0°，由于Γ。、Γι未知，因此为复合假设。
3.根据权利要求2所述的基于广义似然比的多航过SAR相干变化检测方法，其特征在于:所述的步骤S3对协方差矩阵Γ ^和Γ i进行最大似然估计的具体方法为: S31:对Γ ^进行最大似然估计，根据最大似然估计理论，对于假设Htl，根据多航过SAR图像对{f中明显未变化区域得到参数为h的复高斯随机向量Y的L个相互独立的观测Yi,每个观测Yi的概率密度函数为P (Yi Γ 0), Γ 0的似然函数Iik ( Γ 0)表示为:
級(Γ0)二.J.令' gr缺(T0)= O，求得协方差矩阵Γ。的最大似然估计:
T0=Itw ；
L i I S32:对^进行最大似然估计，对于假设H1，根据多航过SAR图像对{f\，f2，一，fK}中明显变化区域得到的参数为Γ i的复高斯随机向量Y的S个相互独立的观测Yi,每个观测Yi的概率密度函数为P (Yi I Γ\)，的似然函数IiHr1)可以表示为:
λ
s1-TviwI' X 令.M(F1) = Oi求出协方差矩阵Γ\的最大似然估计:
ΛI
Γ — 1 \ VV^
I —飞 Zmd ? ? Q
^ i I
4.根据权利要求3所述的基于广义似然比的多航过SAR相干变化检测方法，其特征在于:所述的步骤S4包括以下子步骤: S41:确定检验统计量，根据步骤S2的假设和步骤S3的参数估计得到:
I(Λ、
I Hq) =Λ~黑隱[O !顯
、 }
MXhw IH1) = exp?-XJ Γ? 1 XwwI
Z IΓ, 1、J 取K幅多航过SAR图像中第m行第η列像素点的Q个相互独立的像素对K…，疋，根据似然比假设检验理论，令似然比:
—p(xL>xL,-x￡,;H{I) ^ O ρ(χ;βι；Η0)
z— P(XL5XL3-XiiH1)"-.P(X^h1) 将P(X_ IH0)、p(Xm ] H,)代入上式化简得，exp -Tr|p? -Γ, jXxLxLH| 取对数并忽略常数项，得检测统计量:Z = TrIir0 -Γ, 1I^XLnXInnllI = TrirdG)
|vJ 1-]JiJ 其中Tr^dG丨表示矩阵Γ? G的迹，也就是矩阵对角元素的和，ΛΛ ^ Λ ""I,
Γ 署，f-ι /.、 X mr/ %ri ifa o1 山=L 臟 _ *
;=1 S42:进行假设检验，定义一个和用来检测的图像一样大小的矩阵，记为R，选择门限T，判断门限T与检测统计量Z的大小:当Z>T时，判定假设H1成立，变化检测结果就是该像素对应的区域发生变化，令R对应的像素值为255 ;否则判定假设Htl成立，变化检测结果就是该像素对应的区域未发生变化，令R对应的像素值为O。
5.根据权利要求4所述的基于广义似然比的多航过SAR相干变化检测方法，其特征在于:所述的步骤S5中得到变化检测结果的方法为:依次选取多航过SAR图像像素对，重复步骤S4，直到确定矩阵R中所有的像素值，R即为变化检测结果。
【文档编号】G01S13/90GK104166128SQ201410384298
【公开日】2014年11月26日 申请日期:2014年8月6日 优先权日:2014年8月6日
【发明者】黄钰林, 王园园, 武俊杰, 刘晓佳, 杨建宇 申请人:电子科技大学