CN113160264A - 基于融入空间信息的加权伽马混合模型的sar影像分割方法 - Google Patents
基于融入空间信息的加权伽马混合模型的sar影像分割方法 Download PDFInfo
- Publication number
- CN113160264A CN113160264A CN202110398274.0A CN202110398274A CN113160264A CN 113160264 A CN113160264 A CN 113160264A CN 202110398274 A CN202110398274 A CN 202110398274A CN 113160264 A CN113160264 A CN 113160264A
- Authority
- CN
- China
- Prior art keywords
- sar image
- model
- spatial information
- pixel
- weighted
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/143—Segmentation; Edge detection involving probabilistic approaches, e.g. Markov random field [MRF] modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10044—Radar image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20076—Probabilistic image processing
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Software Systems (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了基于融入空间信息的加权伽马混合模型的SAR影像分割方法,属于影像处理技术领域。读取待分割SAR影像;利用融入空间信息的加权伽马混合模型建立SAR影像统计模型,根据贝叶斯定理,得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型;设定迭代次数,并在每次迭代过程中,利用马尔可夫链蒙特卡洛方法求解形状参数,以及利用期望最大化方法求解分量权重和尺度参数,进而得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型的最优解;输出SAR影像分割结果。有效利用了像素光谱信息和空间信息,克服了由于形状参数结构复杂导致的参数求解困难、计算量大、效率低等问题,提高了SAR影像分割精度。
Description
技术领域
本发明涉及影像处理技术领域,尤其涉及一种基于融入空间信息的加权伽马混合模型的SAR影像分割方法。
背景技术
由于SAR(Synthetic Aperture Radar,合成孔径雷达)影像内像素光谱测度统计分布近似服从某概率分布,因此基于统计模型的SAR影像分割方法受到了学者们的广泛关注和研究,该方法假设像素光谱测度服从同一已知概率分布,通过最大似然估计将影像分割问题转化为模型参数求解问题,进而实现影像分割。
有限混合模型是一种可有效建模像素光谱测度统计分布的统计模型,由多个组份加权定义,各组份用于建模各区域内像素光谱测度统计分布,由同一概率分布定义,如伽马分布,进而形成伽马混合模型。伽马混合模型是由多个伽马分布加权定义,常用于SAR影像分割。基于伽马混合模型的SAR影像分割方法的具体实现过程是,采用伽马混合模型建模SAR影像像素光谱测度统计分布,并假设像素之间相互独立得到似然函数作为影像统计模型;由于SAR影像自身包含大量斑点噪声,采用Gibbs分布利用局部像素类属性建立组份权重的先验分布,以考虑到像素空间位置信息提高分割精度;根据贝叶斯定理,结合影像统计模型和组份权重的先验分布建立后验分布作为影像分割模型,通过最大化后验分布求解模型参数,进而实现影像分割。
目前,已提出的基于伽马混合模型的SAR影像分割方法存在以下几点问题:(1)伽马混合模型难以准确建立SAR影像统计模型。由于SAR影像自身固有的斑点噪声,导致其各区域内像素光谱测度统计分布呈现复杂的统计特性,如重尾、尖峰、双峰等,而伽马分布呈现单峰且右侧重尾特性,不能准确拟合复杂统计分布,进而伽马混合模型以伽马分布作为组份难以满足建立SAR影像统计模型的要求;(2)对像素空间信息统计建模增加了分割模型复杂性。为了避免SAR影像内斑点噪声对分割结果的影响,基于伽马混合模型的SAR影像分割方法通常利用邻域像素类属性建立组份权重的先验分布,根据贝叶斯定理,将影像空间信息引入分割模型,但同时增加了分割模型中参数结构的复杂性,给后续参数求解方法的设计带来了难题和挑战;(3)难以准确求解分割模型的参数。准确求解分割模型的参数是实现高精度影像分割的必要条件,但在基于伽马混合模型的SAR影像分割方法中,伽马分布的形状参数以伽马函数的形式存在其中,因此形状参数结构比较复杂,需采用参数优化方法求解形状参数,但导致该类方法计算量大,效率低。另外,由于组份权重先验分布中考虑到局部像素空间信息,通过最大似然估计不能推导出组份权重的解析式,同样需要采用参数优化方法求解组份权重,造成该类方法计算量大,效率低。
综上,由于伽马混合模型存在统计建模不准确、分割模型复杂、参数求解困难等问题,基于伽马混合模型的SAR影像分割方法的精度有待提高。
发明内容
针对上述现有技术存在的不足,本发明提供一种基于融入空间信息的加权伽马混合模型的SAR影像分割方法,旨在实现对SAR影像的精确分割。
本发明所采取的技术方案是:
一种基于融入空间信息的加权伽马混合模型的SAR影像分割方法,包括以下步骤:
步骤1:读取待分割的SAR影像,并在影像域上对其进行定义;
步骤2:在影像域上,利用融入空间信息的加权伽马混合模型建立SAR影像统计模型,根据贝叶斯定理,得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型;
步骤3:针对已建立的基于融入空间信息的加权伽马混合模型的SAR影像分割模型,设定迭代次数,并在每次迭代过程中,利用马尔可夫链蒙特卡洛方法求解形状参数,以及利用期望最大化方法求解分量权重和尺度参数,进而得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型的最优解;
步骤4:输出SAR影像分割结果。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,步骤1所述的在影像域上对所读取的待分割的SAR影像进行定义的具体内容为:将读取的SAR影像视为像素光谱测度集合,表示为x={xi;i=1,2,...,n},其中i为像素索引,xi为像素i的光谱测度,n为总像素数,将x定义为影像域D上随机场X={Xi;i=1,2,...,n}的一个实现,其中Xi为像素i光谱测度的随机变量。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,所述步骤2包括如下步骤:
步骤2.1:在影像域上,采用加权伽马混合模型建立像素光谱测度xi的概率分布,并通过对SAR影像像素光谱测度xi的概率分布连乘得到像素光谱测度联合概率分布,称为基于加权伽马混合模型的SAR影像统计模型;
步骤2.2:在基于加权伽马混合模型的SAR影像统计模型基础上,利用邻域像素类属性定义组份权重,将像素空间信息融入加权伽马混合模型,建立基于融入空间信息的加权伽马混合模型的SAR影像统计模型;
步骤2.3:根据贝叶斯理论,结合基于加权伽马混合模型的SAR影像统计模型和形状参数的先验分布建立基于融入空间信息的加权伽马混合模型的SAR影像分割模型。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,所述加权伽马混合模型的结构如下:采用伽马分布定义分量概率分布,用于建模各子类别像素光谱测度概率分布;通过分量概率分布加权定义组份概率分布,用于建模各类别像素光谱测度概率分布;通过组份概率分布加权定义加权伽马混合模型的概率分布,用于建模SAR影像内像素光谱测度概率分布。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,所述步骤2.2进一步包括如下步骤:
步骤2.2.1:令第t次迭代中加权伽马混合模型参数集为Ψ(t)={π(t),w(t),α(t),β(t)},根据贝叶斯定理,利用步骤2.1得到的组份概率分布和组份权重建立像素类属性后验概率,即像素隶属于各类别的概率;
步骤2.2.2:将像素空间位置关系视为马尔可夫随机场,利用邻域像素类属性后验概率定义中心像素组份权重;
步骤2.2.3:将步骤2.2.2得到的组份权重代入步骤2.1得到的基于加权伽马混合模型的SAR影像统计模型,得到基于融入空间信息的加权伽马混合模型的SAR影像统计模型。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,所述步骤2.3进一步包括如下步骤:
步骤2.3.1:建立形状参数的先验分布;
利用给定参数的高斯分布建立形状参数αlj的先验分布p(αlj),并假设各组份和分量的形状参数相互独立,通过该先验分布p(αlj)连乘建立形状参数的联合先验概率分布;
步骤2.3.2:根据贝叶斯理论,结合步骤2.2.3得到的基于融入空间信息的加权伽马混合模型的SAR影像统计模型和步骤2.3.1得到的形状参数先验分布建立后验分布,称为基于融入空间信息的加权伽马混合模型的SAR影像分割模型。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,步骤2.2.2中所述的利用邻域像素类属性后验概率定义中心像素组份权重,具体过程如下:对邻域像素类属性后验概率取均值,以平滑噪声或异常值的影响;为了满足组份权重的非负条件,对该均值取指数函数;为了满足组份权重对类别求和为1的条件,对该指数函数作归一化处理,进而建立满足约束条件的组份权重。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,所述步骤3进一步包括如下步骤:
步骤3.1:利用MCMC方法设计更新形状参数操作,选取候选形状参数集α*,并根据最大化后验分布准则计算候选形状参数集的接受率,若接受率等于1,则接受候选形状参数集α*;否则,保持形状参数集α不变;
步骤3.2:在步骤3.1更新形状参数操作的基础上,利用EM方法求解尺度参数和分量权重;
步骤3.2.1:对步骤2.3.2得到的基于融入空间信息的加权伽马混合模型的SAR影像分割模型取对数,并忽略与尺度参数和分量权重无关项,得到对数似然函数;
步骤3.2.2:在EM方法E步中,对步骤3.2.1得到的对数似然函数计算其条件期望函数;
步骤3.2.3:在EM方法M步中,利用步骤3.2.2得到的条件期望函数对尺度参数βlj求偏导,并令导数为0,得到尺度参数的解析式;
步骤3.2.4:在步骤3.2.2得到的条件期望函数的基础上,采用拉格朗日乘数法定义带有约束条件的目标函数,以满足分量权重的约束条件;
步骤3.2.5:利用步骤3.2.4得到的带有约束条件的目标函数对分量权重wlij和拉格朗日系数λw求偏导,并令导数为0,得到分量权重的解析式;
步骤3.3:按设定的迭代次数重复执行步骤3.1至步骤3.2,直到步骤3.2.1中的对数似然函数收敛,此时所对应的参数集{π,w,α,β}为基于融入空间信息的加权伽马混合模型的SAR影像分割模型的最优解。
进一步地,根据所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,步骤4所述的SAR影像分割结果为:利用步骤3.3得到的最优解计算类属性后验概率zli,通过最大化类属性后验概率得到像素标号集c={ci;i=1,2,...,n},即SAR影像最优分割结果,其中ci∈{1,2,...,k}为像素i类属性标号。
与现有技术相比,本发明具有如下有益效果:(1)本发明提出加权伽马混合模型以加权伽马分布作为组份,突破了传统伽马混合模型以单一伽马分布作为组份的传统思想,可更加准确建模SAR影像内像素光谱测度的复杂统计分布,进而有效利用SAR影像的光谱信息提高分割结果质量;(2)本发明利用邻域像素类属性建立组份权重的先验分布,以将像素空间信息融入加权伽马混合模型,在克服SAR影像噪声对分割结果影响的同时,避免了由于引入空间信息导致的分割模型结构复杂的问题;(3)本发明结合期望最大化方法(Exception Maximization,EM)和马尔可夫链蒙特卡洛方法(Markov Chain Monte Carlo,MCMC)求解模型参数,可求解形状参数的最优解,避免了由于形状参数结构复杂导致的参数求解困难、计算量大、效率低等问题,进而提高SAR影像分割精度。
附图说明
图1为本发明的基于融入空间信息的加权伽马混合模型的SAR影像分割方法的流程图;
图2为本发明实施例的待分割SAR影像,其中(a)为一SAR海冰影像,(b)为SAR城市影像;(c)为另一SAR海冰影像;
图3为本发明实施例的待分割SAR影像的分割结果影像,其中(a)为图2(a)示出的SAR海冰影像的分割结果;(b)为图2(b)示出的SAR城市影像的分割结果;(c)为图2(c)示出的SAR海冰影像的分割结果;
图4为本发明实施例的待分割SAR影像灰度直方图拟合结果,其中(a)为对图2(a)示出的SAR海冰影像灰度直方图的拟合结果;(b)为对图2(b)示出的SAR城市影像灰度直方图的拟合结果;(c)为对图2(c)示出的另一SAR海冰影像灰度直方图的拟合结果。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
图1是本发明的基于融入空间信息的加权伽马混合模型的SAR影像分割方法的流程图,所述基于融入空间信息的加权伽马混合模型的SAR影像分割方法,包括以下步骤:
步骤1:读取待分割的SAR影像,并在影像域上对其进行定义。
将读取的SAR影像视为像素光谱测度集合,表示为x={xi;i=1,2,...,n},其中i为像素索引,xi为像素i的光谱测度,n为总像素数,将x定义为影像域D上随机场X={Xi;i=1,2,...,n}的一个实现,其中Xi为像素i光谱测度的随机变量。
步骤2:在影像域上,利用融入空间信息的加权伽马混合模型建立SAR影像统计模型,根据贝叶斯理论,得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型。
步骤2.1:在影像域上,采用加权伽马混合模型对像素光谱测度xi的统计分布进行建模,得到基于加权伽马混合模型的SAR影像统计模型。
所述加权伽马混合模型的结构如下:采用伽马分布定义分量概率分布,用于建模各子类别像素光谱测度概率分布;通过分量概率分布加权定义组份概率分布,用于建模各类别像素光谱测度概率分布;通过组份概率分布加权定义加权伽马混合模型的概率分布,用于建模SAR影像内像素光谱测度概率分布。
在影像域上,采用加权伽马混合模型建立SAR影像像素光谱测度xi的概率分布,如下公式所示:
其中,Ψi={πi,wi,α,β}为像素i的加权伽马混合模型参数集,πi={πli;l=1,2,...,k}为像素i组份权重集,l为组份(类别)索引,k为总组份数,πli表示像素i隶属于类别l的先验概率,其满足条件0<πli<1和p(xi|wli,αl,βl)为组份概率分布,wi={wli;l=1,2,...,k}为像素i分量权重集,wli={wlij;j=1,2,...,m}为类别l内像素i分量权重集,j为分量(子类别)索引,m为总分量数,wlij表示像素i隶属于类别l中子类别j的概率,其满足条件0<wlij<1和G(xi|αlj,βlj)为分量概率分布,定义为伽马分布,α={αl;l=1,2,...,k}为形状参数集,αl={αlj;j=1,2,...,m}为组份l形状参数集,αlj为组份l内分量j的形状参数,β={βl;l=1,2,...,k}为尺度参数集,βl={βlj;j=1,2,...,m}为组份l尺度参数集,βlj为组份l内分量j的尺度参数,Γ(αlj)为形状参数αlj的伽马函数,如下公式所示:
其中,s为随机变量。
假设各像素之间相互独立,通过对SAR影像像素光谱测度xi的的概率分布连乘得到像素光谱测度联合概率分布,称为基于加权伽马混合模型的SAR影像统计模型,如下公式所示:
其中,Ψ={Ψi;i=1,2,...,n}为加权伽马混合模型参数集。
步骤2.2:在基于加权伽马混合模型的SAR影像统计模型基础上,利用邻域像素类属性定义组份权重,将像素空间信息融入加权伽马混合模型,建立基于融入空间信息的加权伽马混合模型的SAR影像统计模型。
步骤2.2.1:令第t次迭代中加权伽马混合模型参数集为Ψ(t)={π(t),w(t),α(t),β(t)},根据贝叶斯定理,利用步骤2.1得到的组份概率分布和组份权重建立像素类属性后验概率,即像素隶属于各类别的概率,如下公式所示:
同理,根据贝叶斯定理,利用步骤2.1得到的分量概率分布和分量权重建立像素子类属性后验概率,即像素隶属于各子类别的概率,如下公式所示:
步骤2.2.2:将像素空间位置关系视为马尔可夫随机场,利用邻域像素类属性后验概率定义中心像素组份权重,具体过程如下:对邻域像素类属性后验概率取均值,以平滑噪声或异常值的影响;为了满足组份权重的非负条件,对该均值取指数函数;为了满足组份权重对类别求和为1的条件,对该指数函数作归一化处理,进而建立满足约束条件的组份权重,如下公式所示:
其中,i'为像素i的邻域像素索引,Ni为像素i邻域像素索引集合,η为邻域像素对中心像素的平滑作用系数。
步骤2.2.3:将步骤2.2.2得到的组份权重代入步骤2.1得到的基于加权伽马混合模型的SAR影像统计模型,得到基于融入空间信息的加权伽马混合模型的SAR影像统计模型,如下公式所示:
步骤2.3:根据贝叶斯理论,结合基于融入空间信息的加权伽马混合模型的SAR影像统计模型和形状参数的先验分布建立基于融入空间信息的加权伽马混合模型的SAR影像分割模型。
步骤2.3.1:建立形状参数的先验分布;
利用给定参数的高斯分布建立形状参数αlj的先验分布p(αlj),并假设各组份和分量的形状参数相互独立,通过该先验分布p(αlj)连乘建立形状参数的联合先验概率分布,如下公式所示:
其中,μα=4为形状参数先验分布的均值,σα=2为形状参数先验分布的标准差。
步骤2.3.2:根据贝叶斯理论,结合步骤2.2.3得到的基于融入空间信息的加权伽马混合模型的SAR影像统计模型和步骤2.3.1得到的形状参数先验分布建立后验分布,称为基于融入空间信息的加权伽马混合模型的SAR影像分割模型,如下公式所示:
步骤3:针对已建立的基于融入空间信息的加权伽马混合模型的SAR影像分割模型,设定迭代次数,并在每次迭代过程中,利用马尔可夫链蒙特卡洛方法(Markov ChainMonte Carlo,MCMC)求解形状参数,利用期望最大化方法(Exception Maximization,EM)求解分量权重和尺度参数,进而得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型的最优解。
步骤3.1:利用MCMC方法设计更新形状参数操作,选取候选形状参数集,并根据最大化后验分布准则计算候选形状参数集的接受率。
令当前形状参数集为α={αlj;l=1,2,...,k,j=1,2,...,m},随机选择组份l∈{1,2,...,k}和分量j∈{1,2,...,m},所对应的αlj为待更新的形状参数,以αlj为均值、ε为标准差的高斯分布生成候选形状参数αlj *,进而候选形状参数集为α*={α1,...,αl-1,αl *,αl+1,...,αk},其中αl *={αl1,...,αlj-1,αlj *,αlj+1,...,αlm},根据最大化后验分布计算候选形状参数的接受率,如下公式所示:
若接受率等于1,则接受候选形状参数集α*;否则,保持形状参数集α不变。
步骤3.2:在步骤3.1更新形状参数操作的基础上,利用EM方法求解尺度参数和分量权重。EM方法包括E步和M步,其中在E步可计算出关于对数似然函数的条件期望函数,在M步通过最大化条件期望函数求解尺度参数和分量权重。
步骤3.2.1:对步骤2.3.2得到的基于融入空间信息的加权伽马混合模型的SAR影像分割模型取对数,并忽略与尺度参数和分量权重无关项,得到对数似然函数,如下公式所示:
步骤3.2.2:在EM方法E步中,对步骤3.2.1得到的对数似然函数计算其条件期望函数,如下公式所示:
步骤3.2.3:在EM方法M步中,利用步骤3.2.2得到的条件期望函数对尺度参数βlj求偏导,并令导数为0,得到尺度参数的解析式,如下公式所示:
步骤3.2.4:在步骤3.2.2得到的条件期望函数的基础上,采用拉格朗日乘数法定义带有约束条件的目标函数,以满足分量权重的约束条件,如下公式所示:
其中,λw为拉格朗日系数。
步骤3.2.5:利用步骤3.2.4得到的带有约束条件的目标函数对分量权重wlij和拉格朗日系数λw求偏导,并令导数为0,得到分量权重的解析式,如下公式所示:
步骤3.3:按设定的迭代次数重复执行步骤3.1至步骤3.2,直到步骤3.2.1中的对数似然函数收敛,此时所对应的参数集{π,w,α,β}为基于融入空间信息的加权伽马混合模型的SAR影像分割模型的最优解。
步骤4:输出SAR影像分割结果。
利用步骤3.3得到的最优解计算类属性后验概率zli,通过最大化类属性后验概率得到像素标号集c={ci;i=1,2,...,n},即SAR影像最优分割结果,其中ci∈{1,2,...,k}为像素i类属性标号。
图2中各幅图像示出的是待分割的SAR影像,其中(a)为一SAR海冰影像,(b)为SAR城市影像,(c)为另一SAR海冰影像。采用本发明的分割方法分别对(a)、(b)和(c)图进行分割,分割结果对应图3中的(a)、(b)和(c)图,即(a)为图2(a)示出的SAR海冰影像的分割结果;(b)为图2(b)示出的SAR城市影像的分割结果;(c)为图2(c)示出的SAR海冰影像的分割结果。从图中可以看出本发明方法能够将SAR影像内各区域准确分割开,获得高质量的SAR影像分割结果。采用本发明的分割方法分别对图2(a)、(b)和(c)所示SAR影像的灰度直方图进行拟合,拟合结果对应图4中的(a)、(b)和(c)图,图中横纵坐标分别为像素光谱测度值及其对应频数,实线代表加权伽马混合模型的拟合曲线,虚线代表加权伽马混合模型组份的拟合曲线,从组份概率分布拟合曲线可看出,加权伽马混合模型组份通过分量加权结构可准确拟合灰度直方图内各峰值,进而加权伽马混合模型通过各组份概率分布加权结构可准确拟合SAR影像灰度直方图,因此本发明方法可对SAR影像内复杂的光谱测度统计分布进行准确建模。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (9)
1.一种基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,包括:
步骤1:读取待分割的SAR影像,并在影像域上对其进行定义;
步骤2:在影像域上,利用融入空间信息的加权伽马混合模型建立SAR影像统计模型,根据贝叶斯理论,得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型;
步骤3:针对已建立的基于融入空间信息的加权伽马混合模型的SAR影像分割模型,设定迭代次数,并在每次迭代过程中,利用马尔可夫链蒙特卡洛方法求解形状参数,以及利用期望最大化方法求解分量权重和尺度参数,进而得到基于融入空间信息的加权伽马混合模型的SAR影像分割模型的最优解;
步骤4:输出SAR影像分割结果。
2.根据权利要求1所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,步骤1所述的在影像域上对所读取的待分割的SAR影像进行定义的具体内容为:将读取的SAR影像视为像素光谱测度集合,表示为x={xi;i=1,2,...,n},其中i为像素索引,xi为像素i的光谱测度,n为总像素数,将x定义为影像域D上随机场X={Xi;i=1,2,...,n}的一个实现,其中Xi为像素i光谱测度的随机变量。
3.根据权利要求1或2所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,所述步骤2包括如下步骤:
步骤2.1:在影像域上,采用加权伽马混合模型对像素光谱测度xi的统计分布进行建模,并通过对SAR影像像素光谱测度xi的概率分布连乘得到像素光谱测度联合概率分布,称为基于加权伽马混合模型的SAR影像统计模型;
步骤2.2:在基于加权伽马混合模型的SAR影像统计模型基础上,利用邻域像素类属性定义组份权重,将像素空间信息融入加权伽马混合模型,建立基于融入空间信息的加权伽马混合模型的SAR影像统计模型;
步骤2.3:根据贝叶斯理论,结合基于融入空间信息的加权伽马混合模型的SAR影像统计模型和形状参数的先验分布建立基于融入空间信息的加权伽马混合模型的SAR影像分割模型。
4.根据权利要求3所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,所述加权伽马混合模型的结构如下:采用伽马分布定义分量概率分布,用于建模各子类别像素光谱测度概率分布;通过分量概率分布加权定义组份概率分布,用于建模各类别像素光谱测度概率分布;通过组份概率分布加权定义加权伽马混合模型的概率分布,用于建模SAR影像内像素光谱测度概率分布。
5.根据权利要求3所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,所述步骤2.2进一步包括如下步骤:
步骤2.2.1:令第t次迭代中加权伽马混合模型参数集为Ψ(t)={π(t),w(t),α(t),β(t)},根据贝叶斯定理,利用步骤2.1得到的组份概率分布和组份权重建立像素类属性后验概率,即像素隶属于各类别的概率;
步骤2.2.2:将像素空间位置关系视为马尔可夫随机场,利用邻域像素类属性后验概率定义中心像素组份权重;
步骤2.2.3:将步骤2.2.2得到的组份权重代入步骤2.1得到的基于加权伽马混合模型的SAR影像统计模型,得到基于融入空间信息的加权伽马混合模型的SAR影像统计模型。
6.根据权利要求5所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,所述步骤2.3进一步包括如下步骤:
步骤2.3.1:建立形状参数的先验分布;
利用给定参数的高斯分布建立形状参数αlj的先验分布p(αlj),并假设各组份和分量的形状参数相互独立,通过该先验分布p(αlj)连乘建立形状参数的联合先验概率分布;
步骤2.3.2:根据贝叶斯理论,结合步骤2.2.3得到的基于融入空间信息的加权伽马混合模型的SAR影像统计模型和步骤2.3.1得到的形状参数先验分布建立后验分布,称为基于融入空间信息的加权伽马混合模型的SAR影像分割模型。
7.根据权利要求5所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,步骤2.2.2中所述的利用邻域像素类属性后验概率定义中心像素组份权重,具体过程如下:对邻域像素类属性后验概率取均值,以平滑噪声或异常值的影响;为了满足组份权重的非负条件,对该均值取指数函数;为了满足组份权重对类别求和为1的条件,对该指数函数作归一化处理,进而构建出满足约束条件的组份权重。
8.根据权利要求6所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,所述步骤3进一步包括如下步骤:
步骤3.1:利用MCMC方法设计更新形状参数操作,选取候选形状参数集α*,并根据最大化后验分布准则计算候选形状参数集的接受率,若接受率等于1,则接受候选形状参数集α*;否则,保持形状参数集α不变;
步骤3.2:在步骤3.1更新形状参数操作的基础上,利用EM方法求解尺度参数和分量权重;
步骤3.2.1:对步骤2.3.2得到的基于融入空间信息的加权伽马混合模型的SAR影像分割模型取对数,并忽略与尺度参数和分量权重无关项,得到对数似然函数;
步骤3.2.2:在EM方法E步中,对步骤3.2.1得到的对数似然函数计算其条件期望函数;
步骤3.2.3:在EM方法M步中,利用步骤3.2.2得到的条件期望函数对尺度参数βlj求偏导,并令导数为0,得到尺度参数的解析式;
步骤3.2.4:在步骤3.2.2得到的条件期望函数的基础上,采用拉格朗日乘数法定义带有约束条件的目标函数,以满足分量权重的约束条件;
步骤3.2.5:利用步骤3.2.4得到的带有约束条件的目标函数对分量权重wlij和拉格朗日系数λw求偏导,并令导数为0,得到分量权重的解析式;
步骤3.3:按设定的迭代次数重复执行步骤3.1至步骤3.2,直到步骤3.2.1中的对数似然函数收敛,此时所对应的参数集{π,w,α,β}为基于融入空间信息的加权伽马混合模型的SAR影像分割模型的最优解。
9.根据权利要求8所述的基于融入空间信息的加权伽马混合模型的SAR影像分割方法,其特征在于,步骤4所述的SAR影像分割结果为:利用步骤3.3得到的最优解计算类属性后验概率zli,通过最大化类属性后验概率得到像素标号集c={ci;i=1,2,...,n},即SAR影像最优分割结果,其中ci∈{1,2,...,k}为像素i类属性标号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110398274.0A CN113160264A (zh) | 2021-04-14 | 2021-04-14 | 基于融入空间信息的加权伽马混合模型的sar影像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110398274.0A CN113160264A (zh) | 2021-04-14 | 2021-04-14 | 基于融入空间信息的加权伽马混合模型的sar影像分割方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113160264A true CN113160264A (zh) | 2021-07-23 |
Family
ID=76890446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110398274.0A Pending CN113160264A (zh) | 2021-04-14 | 2021-04-14 | 基于融入空间信息的加权伽马混合模型的sar影像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113160264A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113723775A (zh) * | 2021-08-16 | 2021-11-30 | 国网上海市电力公司 | 一种基于电力大数据的企业及行业运营风险的评估方法 |
-
2021
- 2021-04-14 CN CN202110398274.0A patent/CN113160264A/zh active Pending
Non-Patent Citations (3)
Title |
---|
王玉等: "多尺度曲波分解下的可变类SAR图像分割", 《信号处理》 * |
石雪: "基于层次化混合模型的高分辨率遥感影像分割方法", 《中国优秀博硕士学位论文全文数据库(博士)基础科学辑》 * |
石雪等: "层次混合模型快速遥感影像分割算法", 《控制与决策》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113723775A (zh) * | 2021-08-16 | 2021-11-30 | 国网上海市电力公司 | 一种基于电力大数据的企业及行业运营风险的评估方法 |
CN113723775B (zh) * | 2021-08-16 | 2023-09-12 | 国网上海市电力公司 | 一种基于电力大数据的企业及行业运营风险的评估方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Liu et al. | Incorporating adaptive local information into fuzzy clustering for image segmentation | |
Wagner et al. | On expected-improvement criteria for model-based multi-objective optimization | |
Hoogerheide et al. | Bayesian forecasting of value at risk and expected shortfall using adaptive importance sampling | |
CN107993237A (zh) | 一种基于窄带约束的几何活动轮廓模型图像局部分割方法 | |
WO2022041598A1 (zh) | 一种遥感影像分割方法、系统、终端以及存储介质 | |
CN111275724B (zh) | 一种基于八叉树和边界优化的机载点云屋顶平面分割方法 | |
CN110211126A (zh) | 基于直觉模糊c均值聚类的图像分割方法 | |
CN104346814B (zh) | 基于层次视觉语义的sar图像分割方法 | |
CN113160264A (zh) | 基于融入空间信息的加权伽马混合模型的sar影像分割方法 | |
US10846559B2 (en) | Image detection method, system and non-volatile computer readable medium | |
CN106600611B (zh) | 基于稀疏三重马尔可夫场的sar图像分割方法 | |
CN110136143A (zh) | 基于admm算法的马氏场下多分辨率遥感图像分割方法 | |
CN117437421A (zh) | 基于空间约束加权伽马混合模型的模糊聚类图像分割方法 | |
CN111553926A (zh) | 一种基于二维Tsallis灰度熵快速迭代的阈值分割方法 | |
Sfikas et al. | Majorization-minimization mixture model determination in image segmentation | |
De Bortoli et al. | Efficient stochastic optimisation by unadjusted langevin monte carlo. application to maximum marginal likelihood and empirical bayesian estimation | |
CN115937465A (zh) | 一种顾及像素区域预测的深度图融合算法 | |
Hong et al. | Optimal visual tracking using Wasserstein transport proposals | |
CN113486904B (zh) | 一种特征自适应激光扫描投影图形控制点提取及优化方法 | |
CN108932520A (zh) | 结合先验概率估计的sar影像水体概率制图方法 | |
CN114782765A (zh) | 一种基于IoU改进损失函数的目标检测方法 | |
CN114359151A (zh) | 组合SegNet网络和局部Nakagami分布拟合模型的超声图像分割方法 | |
CN111008977B (zh) | 基于多节点拓扑重叠测度的高阶mrf模型的图像分割方法 | |
CN109934825B (zh) | 一种结合层次化高斯混合模型和m-h的遥感影像分割方法 | |
CN112419310A (zh) | 一种基于交并融合边框优选的目标检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210723 |
|
WD01 | Invention patent application deemed withdrawn after publication |