CN105005987A - 基于广义gamma分布的SAR图像超像素生成方法 - Google Patents
基于广义gamma分布的SAR图像超像素生成方法 Download PDFInfo
- Publication number
- CN105005987A CN105005987A CN201510348973.9A CN201510348973A CN105005987A CN 105005987 A CN105005987 A CN 105005987A CN 201510348973 A CN201510348973 A CN 201510348973A CN 105005987 A CN105005987 A CN 105005987A
- Authority
- CN
- China
- Prior art keywords
- sar image
- pixel
- region
- image
- sar
- 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
Classifications
-
- 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/20021—Dividing image into blocks, subimages or windows
-
- 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/20112—Image segmentation details
- G06T2207/20128—Atlas-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/30—Subject of image; Context of image processing
- G06T2207/30181—Earth observation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种基于广义gamma分布的SAR图像超像素生成方法。技术方案是:将SAR图像划分成若干区域,利用广义gamma分布描述SAR图像每个区域的区域PDF,PDF是指概率密度函数;结合SAR图像的区域PDF和空间距离信息构造局部聚类准则,对SAR图像的每一个像素进行聚类,更新SAR图像的区域划分;利用基于局部贝叶斯准则的边界演化方法消除SAR图像若干区域中的孤立小区域,则剩余的每个区域对应一个SAR图像超像素。本发明能够更好地利用SAR图像的上下文信息,不仅能消除孤立小区域,还能使得生成的超像素的边界与实际地物的边界更加吻合。
Description
技术领域
本发明涉及成像雷达图像分割技术领域,更具体地说,涉及一种SAR(Synthetic ApertureRadar,合成孔径雷达)图像超像素生成方法。
背景技术
SAR是一种主动微波成像传感器,具有全天时、全天候成像能力,这使得SAR图像在很多军事和民用领域中得到了广泛应用。SAR图像分割是SAR图像自动解译的重要内容,是大量其他SAR图像解译算法包括目标检测与识别、地物分类等的重要环节。近年来,随着SAR成像技术的发展,SAR图像的分辨率不断提高。与传统的中低分辨率SAR图像相比,高分辨率SAR图像能够获取地物和目标更为丰富的信息,如何有效地挖掘这些有效信息是当前SAR图像解译的重要研究内容。SAR图像超像素生成是一种通过对图像像素进行聚类,得到一系列大小相近的具有特定意义的图像区域的方法。SAR图像超像素生成本质上是一种图像分割方法,伴随着高分辨率SAR图像的出现而兴起,在基于区域的SAR图像解译方法研究和应用系统构建中具有广泛的应用前景。
当前,典型的图像超像素生成方法有基于图论的Normalized-Cut方法、凝聚式聚类(agglomerative clustering)方法、Turbopixel方法和SLIC(Simple Linear Iterative Clustering,简单线性迭代聚类)方法等,其研究主要针对光学图像开展。与其他方法相比,SLIC方法在边界拟合以及效率上都具有更优的性能。
SLIC方法主要包括两个步骤:(1)利用局部k-均值算法对图像像素聚类;(2)利用连通分量算法消除孤立的图像小区域。SLIC方法在聚类过程中,利用了各像素到聚类中心(即图像区域的均值)的欧氏距离进行聚类。然而,受SAR图像固有相干斑噪声以及SAR图像幅度值潜在的大动态变化特征的影响,SLIC方法难以生成有效的SAR图像超像素。另外,连通分量算法仅通过最近距离原则消除孤立小区域,没有充分利用图像本身的特征,难以处理步骤(1)生成的初始超像素的边界与实际地物边界偏离的问题。由于PDF(Probability DensityFunction,概率密度函数)比图像区域的均值能够更全面地描述SAR图像区域的统计特性,而广义gamma分布是一种用于描述SAR图像区域统计特性的重要PDF,因此,研究基于广义gamma分布的SAR图像超像素生成方法具有重要的理论意义和应用价值。
发明内容
本发明为了有效解决SAR图像超像素生成问题,提供了一种基于广义gamma分布的SAR图像超像素生成方法。本方法能够生成与实际地物的边界吻合较好的超像素,可以为基于区域的SAR图像解译方法研究和应用系统构建提供重要支撑。
本发明的技术方案是:一种基于广义gamma分布的SAR图像超像素生成方法,将SAR图像划分成若干区域,其特征在于,包括下述步骤:
利用广义gamma分布描述SAR图像的区域PDF;结合SAR图像的区域PDF和空间距离信息构造局部聚类准则,对SAR图像的每一个像素进行聚类,得到SAR图像的若干区域;
利用基于局部贝叶斯准则的边界演化方法消除SAR图像若干区域中的孤立小区域,则剩余的每个区域对应一个SAR图像超像素。
特别地,局部聚类准则是基于待分类像素和待聚类中心的幅度相似度和距离相似度的加权和进行判定。
本发明的有益效果是:结合SAR图像的区域PDF和空间距离信息构造局部聚类准则能够更全面地利用区域的统计信息,生成均匀性更好的初始图像区域。利用基于局部贝叶斯准则的边界演化方法能够更好地利用SAR图像的上下文信息,不仅能消除孤立小区域,还能使得生成的超像素的边界与实际地物的边界更加吻合。
附图说明
图1为本发明的SAR图像超像素生成方法流程图;
图2是利用实测图像进行实验的结果,其中(a)原始SAR图像;(b)采用标准SLIC的结果图;(c)是利用本发明的结果图;
图3是对应实测图像中区域1的实验结果,其中(a)原始SAR图像的区域1;(b)采用标准SLIC的结果图;(c)是利用本发明的结果图;
图4是对应实测图像中区域2的实验结果,其中(a)原始SAR图像的区域2;(b)采用标准SLIC的结果图;(c)是利用本发明的结果图;
具体实施方式
下面结合附图对本发明提供的SAR图像超像素生成方法进行详细说明。
图1是本发明的SAR图像超像素生成方法流程图。该流程的第一步为算法初始化,包括SAR图像标准化、SAR图像区域初始划分、SAR图像像素标签的初始化、SAR图像聚类中心初始化以及相似度矩阵初始化。第二步为基于SAR图像的区域PDF和空间距离信息的SAR图像局部聚类。对应每一个聚类中心,计算其搜索范围内的各像素与其之间的综合相似度,利用它对相似度矩阵进行更新,进而更新相应像素的标签。第三步为基于局部贝叶斯准则的SAR图像区域边界演化。由于受到相干斑噪声等因素的影响,第二步得到的SAR图像区域包含很多孤立的小区域,并且SAR图像区域的边界可能与实际地物的边界存在一定的偏差,通过采用基于局部贝叶斯准则的边界演化处理,不仅可以消除孤立的小区域,还使得生成的超像素的边界与地物的边界更加吻合。
下面具体说明上述各个步骤:
第一步:算法初始化。
步骤①,SAR图像标准化。设给定一幅原始SAR图像I0的幅度值为I0(x,y),1≤x≤M,1≤y≤N,其中(x,y)表示该SAR图像的像素坐标值。对原始SAR图像I0进行标准化可得标准化SAR图像I1,其幅度值为I1(x,y)=I0(x,y)/μ0,1≤x≤M,1≤y≤N,其中为原始SAR图像I0的均值。标准化的目的是为了压缩图像的动态范围,便于构造合理的像素聚类准则。
步骤②,SAR图像区域初始划分。将标准化SAR图像I1划分为一系列相邻的大小均为LS×LS像素的图像区域R1,R2,...,Rk,...,RNs,即其中为区域总数,为计算方便,假设M和N均为LS的整数倍,并且LS的取值根据需要确定。
步骤③,SAR图像像素标签初始化。对由步骤②划分得到的每个图像区域中的像素赋予一个与图像区域序号相同的标签,得到标签图像L,即:
步骤④,SAR图像聚类中心初始化。计算标准化SAR图像I1的梯度模图像G,其每一个像素G(x,y)由下式进行计算:
然后将图像区域Rk中具有最小梯度模的像素作为图像区域Rk的初始聚类中心Ck,即:若 则令Ck=I1(xk,yk)。
步骤⑤,相似度矩阵初始化。初始化相似度矩阵S为M行N列的全零矩阵,即令S={S(x,y)|S(x,y)=0,1≤x≤M,1≤y≤N}。相似度矩阵记录标准化SAR图像I1中的每一个像素点与某个聚类中心的相似度。
第二步:基于图像区域的PDF和空间距离信息的SAR图像局部聚类。
步骤①,估计标准化SAR图像I1中各图像区域的PDF。设图像各区域的PDF为广义gamma分布,第k个图像区域Rk的PDF的形状参数、能量参数和尺度参数为按照下述方式计算:
上式中:
其中: 为符号函数,Ψ(1,·)和Ψ(·)分别为第1阶polygamma函数和digamma函数,NRk表示图像区域Rk中的像素数目。
步骤②,基于图像区域的PDF和空间距离信息的局部聚类。
(a)聚类中心的搜索范围确定。令以Ck为中心的大小为2LS+1×2LS+1的窗口范围为聚类中心Ck的搜索范围,k=1,2,...,Ns。
(b)综合相似度计算。分别对应每个聚类中心Ck,依次计算Ck与其搜索范围内各像素之间的综合相似度为:
其中:表示在聚类中心Ck的搜索范围内的第m个像素在标准化SAR图像I1中的幅度值,表示Ck搜索范围内的像素数目;
(c)相似度矩阵S、SAR图像像素标签以及图像区域更新。
若
则令
将图像I1中标签L(x,y)=k的像素点划归到图像区域Rk中,同时更新
步骤③,聚类中心更新。
步骤④,重复迭代执行步骤①至③,直到聚类中心Ck的坐标值不再变化或者重复迭代次数达到设定的最大值
第三步:基于局部贝叶斯准则的SAR图像区域边界演化。
步骤①,边界像素最大似然演化。令标准化SAR图像I1的图像区域R1,R2,...,Rk,...,的边界像素集合为其中Nb为边界像素的总数,边界像素的坐标值为则按照如下最大似然准则对每一个边界像素的标签进行更新:
其中为边界像素的二阶邻域中的标签矩阵,即:
步骤②,重复迭代执行步骤①,直到迭代次数达到设定的最大次数
步骤③,边界像素最大后验概率演化。按照如下的最大后验概率准则对每一个边界像素的标签再次进行更新:
其中: 为权值矩阵,平滑参数α>0,典型地可取α=1.5;为的指示矩阵: 为的第i′行第j′列元素。
步骤④,重复执行步骤③,直到执行次数达到设定的最大次数然后将原始图像I0中标签L(x,y)=k的像素点划归到图像区域Rk,k=1,2,...,Ns,则图像区域Rk即为原始图像I0的第k个超像素。
图2(a)给出了一幅大小为300×450像素的原始SAR图像,图2(b)和图2(c)分别给出了采用标准SLIC和本发明方法生成的SAR图像超像素结果,其中本发明方法中参数设置为LS=15、w=0.6、和α=1.5。从图中可以看出,标准SLIC方法生成的超像素的边界能够较好吻合对比度较大的地物的边界,但是对于对比度较小的地物的边界吻合较差。相比之下,本发明方法生成的超像素的边界不仅与对比度大的地物的边界吻合很好,还与对比度相对较小的地物的边界吻合较好。
为了更好地观察标准SLIC方法和本发明方法生成的SAR图像超像素结果,图3和图4分别给出了利用两种方法对包含不同地物的典型区域(即分别对应图2(a)中左下角矩形框所示的区域1和右下角矩形框所示的区域2)的处理结果。图3(a)给出了与区域1对应的原始SAR图像,图3(b)和图3(c)分别给出了采用标准SLIC和本发明方法生成的与区域1对应的超像素结果,其中曲线表示超像素的边界;图4(a)给出了与区域2对应的原始SAR图像,图4(b)和图4(c)分别给出了采用标准SLIC和本发明方法生成的与区域2对应的原始SAR图像超像素结果,其中曲线表示超像素的边界。
从图3和图4可以清楚地看出,与标准SLIC方法相比,本发明方法能够生成与实际地物边界更加吻合的超像素,典型的如对应图3(a)中左侧中部椭圆所示的区域A和右下角椭圆所示区域B以及图4(a)中底部中间椭圆所示的区域C的超像素生成结果,从而验证了本发明方法的有效性。
Claims (2)
1.一种基于广义gamma分布的SAR图像超像素生成方法,SAR是指合成孔径雷达,将SAR图像划分成若干区域,其特征在于,还包括下述步骤:
利用广义gamma分布描述SAR图像每个区域的区域PDF,PDF是指概率密度函数;结合SAR图像的区域PDF和空间距离信息构造局部聚类准则,对SAR图像的每一个像素进行聚类,更新SAR图像的区域划分;利用基于局部贝叶斯准则的边界演化方法消除SAR图像若干区域中的孤立小区域,则剩余的每个区域对应一个SAR图像超像素。
2.根据权利要求1所述的基于广义gamma分布的SAR图像超像素生成方法,其特征在于,局部聚类准则是基于待分类像素和待聚类中心的幅度相似度和距离相似度的加权和进行判定。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510348973.9A CN105005987A (zh) | 2015-06-23 | 2015-06-23 | 基于广义gamma分布的SAR图像超像素生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510348973.9A CN105005987A (zh) | 2015-06-23 | 2015-06-23 | 基于广义gamma分布的SAR图像超像素生成方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN105005987A true CN105005987A (zh) | 2015-10-28 |
Family
ID=54378644
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510348973.9A Pending CN105005987A (zh) | 2015-06-23 | 2015-06-23 | 基于广义gamma分布的SAR图像超像素生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105005987A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105426914A (zh) * | 2015-11-19 | 2016-03-23 | 中国人民解放军信息工程大学 | 一种面向位置识别的图像相似度检测方法 |
CN105549009A (zh) * | 2015-12-09 | 2016-05-04 | 电子科技大学 | 一种基于超像素的sar图像cfar目标检测方法 |
CN106778821A (zh) * | 2016-11-25 | 2017-05-31 | 西安电子科技大学 | 基于slic和改进的cnn的极化sar图像分类方法 |
CN107016684A (zh) * | 2017-04-13 | 2017-08-04 | 中国人民解放军国防科学技术大学 | 一种极化sar图像的超像素快速生成方法 |
CN110717956A (zh) * | 2019-09-30 | 2020-01-21 | 重庆大学 | 一种有限角投影超像素引导的l0范数最优化重建方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130188869A1 (en) * | 2012-01-20 | 2013-07-25 | Korea Advanced Institute Of Science And Technology | Image segmentation method using higher-order clustering, system for processing the same and recording medium for storing the same |
WO2014166709A1 (en) * | 2013-04-12 | 2014-10-16 | Thomson Licensing | Superpixel generation with improved spatial coherency |
CN104361351A (zh) * | 2014-11-12 | 2015-02-18 | 中国人民解放军国防科学技术大学 | 一种基于区域统计相似度的合成孔径雷达图像分类方法 |
CN104463210A (zh) * | 2014-12-08 | 2015-03-25 | 西安电子科技大学 | 基于面向对象和谱聚类的极化sar图像分类方法 |
WO2015049118A1 (en) * | 2013-10-02 | 2015-04-09 | Thomson Licensing | Method and apparatus for generating superpixel clusters |
CN104680538A (zh) * | 2015-03-09 | 2015-06-03 | 西安电子科技大学 | 基于超像素的sar图像cfar目标检测方法 |
-
2015
- 2015-06-23 CN CN201510348973.9A patent/CN105005987A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130188869A1 (en) * | 2012-01-20 | 2013-07-25 | Korea Advanced Institute Of Science And Technology | Image segmentation method using higher-order clustering, system for processing the same and recording medium for storing the same |
WO2014166709A1 (en) * | 2013-04-12 | 2014-10-16 | Thomson Licensing | Superpixel generation with improved spatial coherency |
WO2015049118A1 (en) * | 2013-10-02 | 2015-04-09 | Thomson Licensing | Method and apparatus for generating superpixel clusters |
CN104361351A (zh) * | 2014-11-12 | 2015-02-18 | 中国人民解放军国防科学技术大学 | 一种基于区域统计相似度的合成孔径雷达图像分类方法 |
CN104463210A (zh) * | 2014-12-08 | 2015-03-25 | 西安电子科技大学 | 基于面向对象和谱聚类的极化sar图像分类方法 |
CN104680538A (zh) * | 2015-03-09 | 2015-06-03 | 西安电子科技大学 | 基于超像素的sar图像cfar目标检测方法 |
Non-Patent Citations (3)
Title |
---|
FACHAO QIN 等: "Superpixel Segmentation for Polarimetric SAR Imagery Using Local Iterative Clustering", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
XIANXIANG QIN 等: "SAR Image Segmentation via Hierarchical Region Merging and Edge Evolving With Generalized Gamma Distribution", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
张亚亚 等: "基于改进SLIC方法的彩色图像分割", 《计算机工程》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105426914A (zh) * | 2015-11-19 | 2016-03-23 | 中国人民解放军信息工程大学 | 一种面向位置识别的图像相似度检测方法 |
CN105426914B (zh) * | 2015-11-19 | 2019-03-15 | 中国人民解放军信息工程大学 | 一种面向位置识别的图像相似度检测方法 |
CN105549009A (zh) * | 2015-12-09 | 2016-05-04 | 电子科技大学 | 一种基于超像素的sar图像cfar目标检测方法 |
CN105549009B (zh) * | 2015-12-09 | 2017-12-15 | 电子科技大学 | 一种基于超像素的sar图像cfar目标检测方法 |
CN106778821A (zh) * | 2016-11-25 | 2017-05-31 | 西安电子科技大学 | 基于slic和改进的cnn的极化sar图像分类方法 |
CN107016684A (zh) * | 2017-04-13 | 2017-08-04 | 中国人民解放军国防科学技术大学 | 一种极化sar图像的超像素快速生成方法 |
CN110717956A (zh) * | 2019-09-30 | 2020-01-21 | 重庆大学 | 一种有限角投影超像素引导的l0范数最优化重建方法 |
CN110717956B (zh) * | 2019-09-30 | 2023-06-20 | 重庆大学 | 一种有限角投影超像素引导的l0范数最优化重建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Han et al. | A deep learning method for bias correction of ECMWF 24–240 h forecasts | |
CN108038445B (zh) | 一种基于多视角深度学习框架的sar自动目标识别方法 | |
Solomon et al. | Wasserstein propagation for semi-supervised learning | |
CN105005987A (zh) | 基于广义gamma分布的SAR图像超像素生成方法 | |
Yu et al. | A lightweight complex-valued DeepLabv3+ for semantic segmentation of PolSAR image | |
CN103955701B (zh) | 多层次结合的多视合成孔径雷达图像目标识别方法 | |
CN101551856A (zh) | 基于稀疏最小二乘支撑向量机的sar目标识别方法 | |
CN103258324B (zh) | 基于可控核回归和超像素分割的遥感图像变化检测方法 | |
CN110096994A (zh) | 一种基于模糊标签语义先验的小样本PolSAR图像分类方法 | |
US9449395B2 (en) | Methods and systems for image matting and foreground estimation based on hierarchical graphs | |
Wang et al. | Dynamic graph-level neural network for SAR image change detection | |
Park et al. | Bayesian manifold learning: the locally linear latent variable model (LL-LVM) | |
Liu et al. | Target recognition in synthetic aperture radar images via joint multifeature decision fusion | |
CN103366371A (zh) | 基于k分布和纹理特征的sar图像分割方法 | |
Wu et al. | Cost-sensitive latent space learning for imbalanced PolSAR image classification | |
Yang et al. | CEMA-LSTM: Enhancing contextual feature correlation for radar extrapolation using fine-grained echo datasets | |
CN114004998B (zh) | 基于多视张量积扩散的非监督极化sar图像地物分类方法 | |
Leinonen et al. | Unsupervised classification of snowflake images using a generative adversarial network and K-medoids classification | |
CN104123563B (zh) | 一种基于Cloude分解的极化SAR图像分类方法 | |
Wang et al. | RLPGB-Net: Reinforcement learning of feature fusion and global context boundary attention for infrared dim small target detection | |
CN103871060A (zh) | 基于平稳方向波域概率图模型的图像分割方法 | |
Du et al. | Aggregated-attention deformable convolutional network for few-shot SAR jamming recognition | |
Liu et al. | Pseudo-label growth dictionary pair learning for crowd counting | |
Chen et al. | Feature fusion based on convolutional neural network for SAR ATR | |
CN103295236A (zh) | 马尔可夫多特征随机场模型构建方法及其脑部mr图像分割技术 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20151028 |