CN103049653A - 基于em算法的g0分布参数最大似然估计方法 - Google Patents
基于em算法的g0分布参数最大似然估计方法 Download PDFInfo
- Publication number
- CN103049653A CN103049653A CN2012105461645A CN201210546164A CN103049653A CN 103049653 A CN103049653 A CN 103049653A CN 2012105461645 A CN2012105461645 A CN 2012105461645A CN 201210546164 A CN201210546164 A CN 201210546164A CN 103049653 A CN103049653 A CN 103049653A
- Authority
- CN
- China
- Prior art keywords
- beta
- parameter
- sigma
- algorithm
- distribution
- 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.)
- Granted
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于EM算法的G0分布参数最大似然估计方法,首先使用矩估计法对G0分布参数进行估计,然后是将矩估计法得到的参数估计值作为初始值,采用EM算法,以迭代求解的方式估计G0分布参数。本发明所设计的基于EM算法的G0分布参数最大似然估计方法具有较高的参数估计精度。
Description
技术领域
本发明属于合成孔径雷达图像解译领域,涉及一种基于EM算法的G0分布参数最大似然估计方法。
背景技术
由于合成孔径雷达(synthetic aperture radar,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分布最常用的参数估计方法是矩估计法(method of moments,MoM)。矩估计法,计算相对比较简单,但是精度不高,并且由于其方法自身固有的限制,不能实现G0分布全范围的参数估计。G0分布另外一种参数估计方法是基于Mellin变化的,该方法可以实现G0分布全范围的参数估计,但是参数估计的精度同样不是很高。
最大似然(Maximum Likelihood,ML)估计是统计意义上最优的参数估计方法,但是由于G0分布表达式的复杂性,最大似然估计方法给出的方程组具有高度的非线性,肯定没有解析解。而如果运用数值方法来求解最大似然估计方程组,则面临着计算量巨大和可能无法收敛到正确解的困境。这些都限制了最大似然估计方法在G0分布上的运用。但是,作为统计意义上最优的方法,最大似然估计法可能会消耗更多的计算量,理论上却是一定会给出更准确的参数估计。本专利给出了一种基于EM算法的G0分布参数最大似然估计方法。EM算法是一种迭代的方法来寻找统计模型的最大似然估计,常常用于最大似然估计的方程无法直接求解的情况。
发明内容
本发明所要解决的技术问题是提供一种具有较高的参数估计精度的基于EM算法的G0分布最大似然参数估计方法。
本发明为解决上述技术问题采用如下技术方案:本发明设计了一种基于EM算法的G0分布参数最大似然估计方法,包括如下具体步骤:
步骤(1):采用矩估计法,为观测得到的SAR图像生成参数估计初始值;
步骤(2):根据步骤(1)中的参数估计初始值,采用迭代算法估计G0分布参数。
作为本发明的一种优化方法:所述步骤(1)还包括如下具体步骤:
其中;N定义为等效视数,σ定义为形状参数,β定义为尺度参数,t定义为SAR图像幅度数据的序号,T定义为SAR图像幅度数据的个数;
步骤(12):利用矩估计法来估计参数初始值,采用如下公式:
得到尺度参数估计初始值和形状参数估计初始值
作为本发明的一种优化方法:所述步骤(2)包括如下具体处理:
其中,Ψ(·)定义为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:采用矩估计法,生成参数估计的初值,具体方法如下:
步骤2:根据步骤1的矩估计法给出的初始值,采用EM算法以迭代求解的方式估计G0分布参数,具体方法如下:
其中,Ψ(·)表示digamma函数;
为了进一步说明本发明所提出的方法,下面给出本发明方法的理论推导过程。
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)符合generalized Rayleigh分布,即
假设参数ωt符合inverse Gamma分布,即有
式中λ=(β,σ)为inverse Gamma分布的参数。将(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-point iteration)的方法来求解方程(23),这种方法较常用的Newton-Raphson方法简单,而且可以避免求解digamma函数的导数。对于一个给定的ln(AG)值,不动点迭代法给出的迭代求解方程为:
Claims (3)
1.一种基于EM算法的G0分布参数最大似然估计方法,其特征在于,包括如下具体步骤:
步骤(1):采用矩估计法,为观测得到的SAR图像生成参数估计初始值;
步骤(2):根据步骤(1)中的参数估计初始值,采用迭代算法估计G0分布参数。
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 true CN103049653A (zh) | 2013-04-17 |
CN103049653B 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) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105184062A (zh) * | 2015-08-25 | 2015-12-23 | 中国人民解放军后勤工程学院 | 在群智感知网络中基于置信区间的用户感知质量评估方法 |
CN105743593A (zh) * | 2016-01-25 | 2016-07-06 | 重庆邮电大学 | 一种基于双对数累积量期望的Gamma-Gamma分布参数估计方法 |
CN107808380A (zh) * | 2016-12-28 | 2018-03-16 | 中国测绘科学研究院 | 一种基于G0和Gamma联合分布的多尺度高分辨率SAR影像水体分割方法 |
CN109145502A (zh) * | 2018-09-17 | 2019-01-04 | 中国人民解放军海军工程大学 | 一种威布尔型单元寿命分布参数估计方法 |
CN109323876A (zh) * | 2018-09-17 | 2019-02-12 | 中国人民解放军海军工程大学 | 一种估计伽玛型单元可靠性参数的方法 |
CN110967184A (zh) * | 2019-12-03 | 2020-04-07 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060074630A1 (en) * | 2004-09-15 | 2006-04-06 | Microsoft Corporation | Conditional maximum likelihood estimation of naive bayes probability models |
CN102542561A (zh) * | 2011-11-23 | 2012-07-04 | 浙江工商大学 | 基于Fisher分布的活动轮廓SAR图像分割方法 |
-
2012
- 2012-12-17 CN CN201210546164.5A patent/CN103049653B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060074630A1 (en) * | 2004-09-15 | 2006-04-06 | Microsoft Corporation | Conditional maximum likelihood estimation of naive bayes probability models |
CN102542561A (zh) * | 2011-11-23 | 2012-07-04 | 浙江工商大学 | 基于Fisher分布的活动轮廓SAR图像分割方法 |
Non-Patent Citations (2)
Title |
---|
A. P. DEMPSTER等: "Maximum likelihood from incomplete data via the EM algorithm", 《JOURNAL OF THE ROYAL STATISTICAL SOCIETY, SERIES B》, vol. 39, no. 1, 30 April 1977 (1977-04-30), pages 1 - 38 * |
时公涛等: "基于Mellin变换的G0分布参数估计方法", 《自然科学进展 》, vol. 19, no. 6, 19 June 2009 (2009-06-19), pages 677 - 689 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105184062A (zh) * | 2015-08-25 | 2015-12-23 | 中国人民解放军后勤工程学院 | 在群智感知网络中基于置信区间的用户感知质量评估方法 |
CN105743593A (zh) * | 2016-01-25 | 2016-07-06 | 重庆邮电大学 | 一种基于双对数累积量期望的Gamma-Gamma分布参数估计方法 |
CN105743593B (zh) * | 2016-01-25 | 2018-06-05 | 重庆邮电大学 | 一种基于双对数累积量期望的Gamma-Gamma分布参数估计方法 |
CN107808380A (zh) * | 2016-12-28 | 2018-03-16 | 中国测绘科学研究院 | 一种基于G0和Gamma联合分布的多尺度高分辨率SAR影像水体分割方法 |
CN109145502A (zh) * | 2018-09-17 | 2019-01-04 | 中国人民解放军海军工程大学 | 一种威布尔型单元寿命分布参数估计方法 |
CN109323876A (zh) * | 2018-09-17 | 2019-02-12 | 中国人民解放军海军工程大学 | 一种估计伽玛型单元可靠性参数的方法 |
CN109323876B (zh) * | 2018-09-17 | 2020-10-16 | 中国人民解放军海军工程大学 | 一种估计伽玛型单元可靠性参数的方法 |
CN110967184A (zh) * | 2019-12-03 | 2020-04-07 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和*** |
CN110967184B (zh) * | 2019-12-03 | 2021-06-11 | 合肥工业大学 | 基于振动信号分布特征识别的变速箱故障检测方法和*** |
Also Published As
Publication number | Publication date |
---|---|
CN103049653B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103049653A (zh) | 基于em算法的g0分布参数最大似然估计方法 | |
CN105975645B (zh) | 一种基于多步的含激波区域飞行器流场快速计算方法 | |
Furukawa et al. | Convolution quadrature time-domain boundary element method for 2-D and 3-D elastodynamic analyses in general anisotropic elastic solids | |
Mazzi et al. | Fractal-generated turbulence | |
CN107247687B (zh) | 一种基于特征正交分解的非平稳随机过程快速模拟方法 | |
Gao et al. | Probability density function for wave elevation based on Gaussian mixture models | |
CN106330197B (zh) | 一种建筑风洞测压试验数据压缩方法 | |
CN103136239A (zh) | 一种基于张量重建的交通数据丢失恢复方法 | |
CN106646417A (zh) | 广义帕累托分布参数的迭代最大似然估计方法 | |
CN104318593A (zh) | 一种雷达海杂波的仿真方法和*** | |
CN106842163A (zh) | 一种弹道目标回波信号时频特性估计方法 | |
CN103984844A (zh) | 一种不等长序列相似度度量算法 | |
Szmelter et al. | An unstructured-mesh atmospheric model for nonhydrostatic dynamics: Towards optimal mesh resolution | |
CN104795063A (zh) | 一种基于声学空间非线性流形结构的声学模型构建方法 | |
CN105939026A (zh) | 基于混合Laplace分布的风电功率波动量概率分布模型建立方法 | |
CN107064893A (zh) | 基于对数矩的广义帕累托分布参数估计方法 | |
Loseille et al. | Anisotropic adaptive simulations in aerodynamics | |
Ma et al. | Improved weighted compact nonlinear scheme for implicit large-eddy simulations | |
Wang et al. | Study on a new aerodynamic model of HAWT based on panel method and Reduced Order Model using Proper Orthogonal Decomposition | |
Lim et al. | Design wave method for the extreme horizontal slow-drift motion of moored floating platforms | |
Li et al. | Performance study of a simplified shape optimization strategy for blended-wing-body underwater gliders | |
CN102928827B (zh) | 一种基于past的快速降维空时自适应处理方法 | |
CN104731762A (zh) | 基于循环移位的立方相位信号参数估计方法 | |
EP2416266A2 (en) | Computer-assisted method for the cost-optimised calculation of variable distributions on an aerodynamic surface | |
Akbarzade et al. | Application of variational iteration method to partial differential equation systems |
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 |
Granted publication date: 20160406 Termination date: 20171217 |
|
CF01 | Termination of patent right due to non-payment of annual fee |