CN103049653B - 基于em算法的g0分布参数最大似然估计方法 - Google Patents
基于em算法的g0分布参数最大似然估计方法 Download PDFInfo
- Publication number
- CN103049653B CN103049653B CN201210546164.5A CN201210546164A CN103049653B CN 103049653 B CN103049653 B CN 103049653B CN 201210546164 A CN201210546164 A CN 201210546164A CN 103049653 B CN103049653 B CN 103049653B
- Authority
- CN
- China
- Prior art keywords
- beta
- parameter
- sigma
- distribution
- prime
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000007476 Maximum Likelihood Methods 0.000 title claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 48
- 238000012804 iterative process Methods 0.000 description 8
- 238000013179 statistical model Methods 0.000 description 3
- 238000009795 derivation Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000002940 Newton-Raphson method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000009472 formulation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 1
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于EM算法的G0分布参数最大似然估计方法,首先使用矩估计法对G0分布参数进行估计,然后是将矩估计法得到的参数估计值作为初始值,采用EM算法,以迭代求解的方式估计G0分布参数。本发明所设计的基于EM算法的G0分布参数最大似然估计方法具有较高的参数估计精度。
Description
技术领域
本发明属于合成孔径雷达图像解译领域,涉及一种基于EM算法的G0分布参数最大似然估计方法。
背景技术
由于合成孔径雷达(syntheticapertureradar,SAR)图像数据具有噪声严重、背景杂波复杂的特点,因此基于统计模型来展开SAR图像解译工作得到了广泛的关注。统计模型对实测SAR图像数据统计特性描述的准确性会在很大程度上影响SAR图像解译的性能,为此各国学者发展了很多用于描述SAR图像数据的统计模型。其中,G0分布,具有适用范围广、建模能力强优点,近年来被广泛地应用于SAR图像的解译中。
G0分布在SAR图像解译中的应用离不开对其进行参数估计。所谓参数估计,就是利用观测得到的SAR图像数据来估计G0分布的参数。设观测得到的SAR图像幅度数据为:y={yt,t=1,…,T},且独立同分布,则yt的G0分布表达式为:
其中,参数N为SAR图像数据视数的2倍,可以先验获得不需要估计。因此,参数估计即为利用观测得到的SAR图像幅度数据y={yt,t=1,…,T}来估计G0分布中的参数β和σ。
目前,G0分布最常用的参数估计方法是矩估计法(methodofmoments,MoM)。矩估计法,计算相对比较简单,但是精度不高,并且由于其方法自身固有的限制,不能实现G0分布全范围的参数估计。G0分布另外一种参数估计方法是基于Mellin变化的,该方法可以实现G0分布全范围的参数估计,但是参数估计的精度同样不是很高。
最大似然(MaximumLikelihood,ML)估计是统计意义上最优的参数估计方法,但是由于G0分布表达式的复杂性,最大似然估计方法给出的方程组具有高度的非线性,肯定没有解析解。而如果运用数值方法来求解最大似然估计方程组,则面临着计算量巨大和可能无法收敛到正确解的困境。这些都限制了最大似然估计方法在G0分布上的运用。但是,作为统计意义上最优的方法,最大似然估计法可能会消耗更多的计算量,理论上却是一定会给出更准确的参数估计。本专利给出了一种基于EM算法的G0分布参数最大似然估计方法。EM算法是一种迭代的方法来寻找统计模型的最大似然估计,常常用于最大似然估计的方程无法直接求解的情况。
发明内容
本发明所要解决的技术问题是提供一种具有较高的参数估计精度的基于EM算法的G0分布最大似然参数估计方法。
本发明为解决上述技术问题采用如下技术方案:本发明设计了一种基于EM算法的G0分布参数最大似然估计方法,包括如下具体步骤:
步骤(1):采用矩估计法,为观测得到的SAR图像生成参数估计初始值;
步骤(2):根据步骤(1)中的参数估计初始值,采用迭代算法估计G0分布参数。
作为本发明的一种优化方法:所述步骤(1)还包括如下具体步骤:
步骤(11):设观测得到的SAR图像幅度数据为:y={yt,t=1,…,T},且独立同分布,则yt的G0分布表达式为:
其中;N定义为等效视数,σ定义为形状参数,β定义为尺度参数,t定义为SAR图像幅度数据的序号,T定义为SAR图像幅度数据的个数;
步骤(12):利用矩估计法来估计参数初始值,采用如下公式:
其中,定义为k阶样本矩,即有
得到尺度参数估计初始值和形状参数估计初始值
作为本发明的一种优化方法:所述步骤(2)包括如下具体处理:
步骤(21):令建立如下函数变量A和G:
其中,Ψ(·)定义为digamma函数;
步骤(22):令利用步骤(21)中的函数变量A和G与公式进行迭代,当满足迭代终止条件后,结束迭代过程,令参数估计值等于最后一次迭代所得的则σ的估计值由公式求得,并将迭代终止时的和作为最终的参数估计值;若不满足迭代终止条件,则返回步骤(21)。
本发明与现有技术先比具有如下优点:
本发明所设计的基于EM算法的G0分布最大似然参数估计方法的优点为:通过EM算法可以实现对G0分布的参数最大似然估计,该方法比现有的矩估计法和基于Mellin变化的方法具有更高的参数估计精度。
附图说明
图1为本发明所设计的基于EM算法的G0分布最大似然参数估计方法的流程示意图。
具体实施方式
下面结合附图对本发明做进一步的详细说明:
本发明提出一种基于EM算法的G0分布参数最大似然估计方法。由于SAR图像数据具有噪声严重、背景杂波复杂的特点,因此基于统计模型来展开SAR图像解译工作得到了广泛的关注。Frery等人给出了一种新的统计分布模型G分布,G分布的特殊形式G0分布,具有适用范围广、参数估计容易的优点。其中,参数估计是G0分布研究的一个核心问题。矩估计法以及基于Mellin变换的参数估计方法是目前常用的参数估计方法。但是最大似然估计作为统计意义上最优的参数估计方法由于表达式的复杂性一直没有得到应用。为了解决这个问题,本专利提出了基于EM算法的G0分布参数估计方法,该方法是通过迭代的方式来寻找G0分布参数的最大似然估计。
如图1所示,本发明设计了一种基于EM算法的G0分布参数最大似然估计方法,包括以下步骤:
步骤1:采用矩估计法,生成参数估计的初值,具体方法如下:
其中表示k阶样本矩,即有
步骤2:根据步骤1的矩估计法给出的初始值,采用EM算法以迭代求解的方式估计G0分布参数,具体方法如下:
a)令 根据下面的(4)式和(5)式求取函数变量A和G,
其中,Ψ(·)表示digamma函数;
b)令利用(6)式进行迭代,当满足迭代终止条件后,结束迭代过程,令β的估计值等于最后一次迭代所得的
c)由(7)式求取σ的估计值
d)满足迭代终止条件则结束迭代过程,将迭代终止时的和作为最终的参数估计值;若不满足迭代终止条件,则回到a)。
为了进一步说明本发明所提出的方法,下面给出本发明方法的理论推导过程。
G0分布模型是由Frery等人根据经典的乘积模型发展而得。就SAR图像幅度数据而言,可以表示为符合单位均值方根Gamma分布Γ1/2(n,n)的斑点噪声分量与符合逆方根Gamma分布Γ-1/2(α,γ)的表征地物RCS起伏特性的后向散射幅度的乘积。设观测得到SAR图像幅度数据为y={yt,t=1,…,T},且独立同分布,则yt的G0分布模型表达式为:
其中,n表示等效视数,α为形状参数反映了被测区域的均匀度,γ为尺度参数与被测区域的平均能量有关,t定义为SAR图像幅度数据的序号,T定义为SAR图像幅度数据的个数,式(8)是当前关于G0分布的一种常用的表达式。
本专利的方法是基于一种新的G0分布表达式及其推导过程。由贝叶斯理论,我们可以将G0分布写成
假设变量yt为N个独立同分布的高斯随机变量(均值为零,方差为ωt)的平方和的平方根,则有p(yt|ωt)符合generalizedRayleigh分布,即
假设参数ωt符合inverseGamma分布,即有
式中λ=(β,σ)为inverseGamma分布的参数。将(10)和(11)带入(9)式,积分后可以得到
将(12)式与传统的G0分布表达式(8)式相比较,得
由此可见(12)式也是一种G0分布的表达式,本专利中的参数估计是基于(12)式的。上述的G0分布的推导也相当于提供了一种新的方式来解释符合G0分布的SAR数据。
由最大似然估计法可知,G0分布的最大似然参数估计使得
(14)
由于没有(14)式的解析解,我们采用EM算法来求解(14)式的最大似然估计值。EM算法是一种迭代的算法,根据当前的参数估计值产生一个新的参数估计值,并且新的参数估计值相比当前参数估计值具有更大的似然。EM算法的特性保证了,在达到稳定点之前,新的参数估计值总是比当前的参数估计值具有更大的似然。
EM算法由最大化下面的辅助函数得到:
λ′为当前的参数估计值。由(15)式可以得到:
将(11)式代入(16)式,将其等号右部分别对σ和β进行求导并令其为0,可得
为了求解(17),(18)两式中的积分,首先求解下式:
将(10)、(11)、(12)式带入可得
令(20)式中s=-1,则可得(17)式所需的
利用关系式 可以得到(18)式所需的
将(21),(22)的结果代入式(17)和(18)可得
其中,
式(23)~(26)构成了EM算法的迭代求解过程来求解G0分布参数,在给出初始的σ′和β′值时(本专利中,初始的σ′和β′值为矩估计法得到参数估计值),迭代求解过程就可以开始。迭代求解过程可以在满足某些设定的条件的情况下停止。方程(23)没有解析解,本专利采用了不动点迭代(fixed-pointiteration)的方法来求解方程(23),这种方法较常用的Newton-Raphson方法简单,而且可以避免求解digamma函数的导数。对于一个给定的ln(AG)值,不动点迭代法给出的迭代求解方程为:
该迭代过程的起始值可以为当前的估计值β′。由(27)的迭代过程求得后,可以直接由(24)式求得,然后将求得和值代替β′和σ′,完成一次EM迭代。不难看出,整个参数求解过程有两个迭代过程,其中内部的迭代为用式(27)来求解方程(23),而外部的迭代则为EM算法的迭代过程。迭代过程的终止条件可以设定为连续的参数估计值相差不超过某个设定的门限值或者达到设定的最大迭代次数。
Claims (2)
1.一种基于EM算法的G0分布参数最大似然估计方法,其特征在于,包括如下具体步骤:
步骤(1):采用矩估计法,为观测得到的SAR图像生成参数估计初始值;
步骤(2):根据步骤(1)中的参数估计初始值,采用迭代算法估计G0分布参数;
所述步骤(1)还包括如下具体步骤:
步骤(11):设观测得到的SAR图像幅度数据为:y={yt|t=1,…,T},且独立同分布,为实数集,则yt的G0分布表达式为:
其中;N定义为等效视数,σ定义为形状参数,β定义为尺度参数,t定义为SAR图像幅度数据的序号,T定义为SAR图像幅度数据的个数,Γ(·)为Gamma函数;
步骤(12):利用矩估计法来估计参数初始值,采用如下公式:
其中,定义为k阶样本矩,即有
得到尺度参数估计初始值和形状参数估计初始值
2.根据权利要求1所述的基于EM算法的G0分布参数最大似然估计方法,其特征在于,所述步骤(2)包括如下具体处理:
步骤(21):令建立如下函数变量A和G:
其中,Ψ(·)定义为digamma函数;
步骤(22):令利用步骤(21)中的函数变量A和G与公式进行迭代,当满足迭代终止条件后,结束迭代过程,令参数估计值等于最后一次迭代所得的则σ的估计值由公式求得,并将迭代终止时的和作为最终的参数估计值;若不满足迭代终止条件,则返回步骤(21)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210546164.5A CN103049653B (zh) | 2012-12-17 | 2012-12-17 | 基于em算法的g0分布参数最大似然估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210546164.5A CN103049653B (zh) | 2012-12-17 | 2012-12-17 | 基于em算法的g0分布参数最大似然估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103049653A CN103049653A (zh) | 2013-04-17 |
CN103049653B true CN103049653B (zh) | 2016-04-06 |
Family
ID=48062287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210546164.5A Expired - Fee Related CN103049653B (zh) | 2012-12-17 | 2012-12-17 | 基于em算法的g0分布参数最大似然估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103049653B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105184062A (zh) * | 2015-08-25 | 2015-12-23 | 中国人民解放军后勤工程学院 | 在群智感知网络中基于置信区间的用户感知质量评估方法 |
CN105743593B (zh) * | 2016-01-25 | 2018-06-05 | 重庆邮电大学 | 一种基于双对数累积量期望的Gamma-Gamma分布参数估计方法 |
CN107808380B (zh) * | 2016-12-28 | 2021-05-25 | 中国测绘科学研究院 | 一种基于G0和Gamma联合分布的多尺度SAR影像水体分割方法 |
CN109323876B (zh) * | 2018-09-17 | 2020-10-16 | 中国人民解放军海军工程大学 | 一种估计伽玛型单元可靠性参数的方法 |
CN109145502B (zh) * | 2018-09-17 | 2023-05-12 | 中国人民解放军海军工程大学 | 一种威布尔型单元寿命分布参数估计方法 |
CN110967184B (zh) * | 2019-12-03 | 2021-06-11 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542561A (zh) * | 2011-11-23 | 2012-07-04 | 浙江工商大学 | 基于Fisher分布的活动轮廓SAR图像分割方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7624006B2 (en) * | 2004-09-15 | 2009-11-24 | Microsoft Corporation | Conditional maximum likelihood estimation of naïve bayes probability models |
-
2012
- 2012-12-17 CN CN201210546164.5A patent/CN103049653B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102542561A (zh) * | 2011-11-23 | 2012-07-04 | 浙江工商大学 | 基于Fisher分布的活动轮廓SAR图像分割方法 |
Non-Patent Citations (2)
Title |
---|
Maximum likelihood from incomplete data via the EM algorithm;A. P. DEMPSTER等;《Journal of the Royal Statistical Society, Series B》;19770430;第39卷(第1期);第1-38页 * |
基于Mellin变换的G0分布参数估计方法;时公涛等;《自然科学进展 》;20090619;第19卷(第6期);第677-689页 * |
Also Published As
Publication number | Publication date |
---|---|
CN103049653A (zh) | 2013-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103049653B (zh) | 基于em算法的g0分布参数最大似然估计方法 | |
CN105719023A (zh) | 一种基于混合高斯分布的风电功率实时预测误差分析方法 | |
CN103942457A (zh) | 基于关联向量机回归的水质参数时间序列预测方法 | |
CN104899380A (zh) | 一种基于蒙特卡洛模拟的边坡稳定可靠度敏感性分析方法 | |
CN103426030A (zh) | 计及老化因素的电力设备故障率预测方法 | |
CN107038292A (zh) | 一种基于自适应多变量非参数核密度估计的多风电场出力相关性建模方法 | |
CN106354695A (zh) | 一种仅输出线性时变结构模态参数辨识方法 | |
WO2020134388A1 (zh) | 一种基于随机等几何分析的叶片高刚度设计方法 | |
CN107844627A (zh) | 一种仅输出时变结构模态参数贝叶斯估计方法 | |
CN106330197A (zh) | 一种建筑风洞测压试验数据压缩方法 | |
CN103927556A (zh) | 一种基于小波包和近似熵的心电信号分类方法 | |
CN104573876A (zh) | 基于时序长记忆模型的风电场短期风速预测方法 | |
CN106203723A (zh) | 基于rt重构eemd‑rvm组合模型的风功率短期区间预测方法 | |
CN103136239A (zh) | 一种基于张量重建的交通数据丢失恢复方法 | |
CN105320809A (zh) | 一种针对风电场空间相关性的风速预测方法 | |
CN104050604A (zh) | 基于概率潮流的电力系统静态安全评估方法 | |
CN106154243A (zh) | 海杂波Pareto分布模型的参数估计范围拓展方法 | |
Halilu et al. | Enhanced matrix-free method via double step length approach for solving systems of nonlinear equations | |
CN104795063A (zh) | 一种基于声学空间非线性流形结构的声学模型构建方法 | |
CN104166777A (zh) | 计及多重相关性的风速矢量数据模拟生成方法 | |
CN106228468A (zh) | 一种潮汐流能发电场输出功率的概率模拟方法 | |
CN105939026A (zh) | 基于混合Laplace分布的风电功率波动量概率分布模型建立方法 | |
CN105046057A (zh) | 基于Morlet小波核的LSSVM脉动风速预测方法 | |
CN104408525A (zh) | 作业车间调度风险的量化评估与控制方法 | |
CN105183997A (zh) | 一种基于双层嵌套不确定性传播的热传导模型校准方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160406 Termination date: 20171217 |