CN101587587A - 顾及多尺度马尔科夫场的合成孔径雷达图像分割方法 - Google Patents
顾及多尺度马尔科夫场的合成孔径雷达图像分割方法 Download PDFInfo
- Publication number
- CN101587587A CN101587587A CNA2009100631722A CN200910063172A CN101587587A CN 101587587 A CN101587587 A CN 101587587A CN A2009100631722 A CNA2009100631722 A CN A2009100631722A CN 200910063172 A CN200910063172 A CN 200910063172A CN 101587587 A CN101587587 A CN 101587587A
- Authority
- CN
- China
- Prior art keywords
- image
- segmentation result
- scale
- utilize
- segmentation
- 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
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V10/00—Arrangements for image or video recognition or understanding
- G06V10/40—Extraction of image or video features
- G06V10/54—Extraction of image or video features relating to texture
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Multimedia (AREA)
- Theoretical Computer Science (AREA)
- Image Analysis (AREA)
Abstract
本发明公开了一种顾及多特征多尺度MRF场的SAR图像分割方法,用于遥感影像处理,该方法包括:(1)对原始图像利用灰度共生矩阵提取纹理信息,在得到的纹理特征图像上进行多尺度分割;(2)利用小波变换检测原始图像的边缘信息,得到边缘特征图像;(3)利用线性加权法融合纹理多尺度分割结果与边缘特征,即直接对两幅图像的对应像素点进行加权叠加,得到最终分割结果。本发明方法减少了SAR图像的错误分割率,在消除斑点噪声的同时,很好的保留了图像的边缘等细节信息。
Description
技术领域
本发明涉及一种合成孔径雷达图像分割方法,具体地说是一种顾及多特征多尺度马尔科夫场的SAR图像分割方法,属于遥感影像处理领域。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)的全天候、大面积、高分辨率、能穿透地表云雾的特点,使其可在短时间内获得大量的信息,这也使得SAR图像应用成为遥感领域研究的热点。随着国内外星载SAR***的发展,获取了大量的星载SAR图像数据,与SAR数据源的快速增长趋势相比,作为SAR数据应用基础的SAR图像处理与分析相关技术的研究相对滞后。SAR图像分割是SAR图像应用中一种重要的图像技术,目的是把图像分成区域并提取出感兴趣目标,其分割精度的高低是对SAR图像中不同目标进行识别和提取的关键指标。但是SAR的后向散射机制决定了SAR图像中存在相干斑(Speckle)噪声,Speckle噪声可以认为是一种乘性噪声,且是一种非高斯分布的随机噪声,这导致了SAR图像数据模型的非线性化和图像数据分布的另类性。这些相干斑噪声降低了图像的空间分辨率,掩盖了图像的细节结构,使得SAR图像的稳定信息多存在于统计量而不是单纯的灰度值中,这也使得常规的图像分割技术难于直接应用于SAR图像中取得好结果。
目前,国内外研究较多的是含斑点噪声的SAR图像直接分割,其分割方法大致分为以下几类:基于边界提取的分割方法、基于区域的分割方法、基于特定理论工具的分割方法和混合的分割方法。基于边界提取的分割算法是人们最早研究的方法,它通过区域边缘上的像素灰度值的奇异性来解决不同区域间的分割问题。基于区域的分割方法侧重于利用区域内特征的相似性,主要包括直方图阈值化方法、MUM方法、***合并方法、流域分割方法和分类聚类方法等。结合特定理论工具的方法主要包括基于统计学的分割方法、基于数学形态学的分割方法、基于神经网络的分割方法、基于小波变换的分割方法、基于信息论的分割方法、基于水平集理论的分割方法和基于模糊集理论的分割方法等。以上方法各有优缺点,基于边缘检测的方法存在边界点的不连续和同质区域中的伪边界情况,后续处理难度大。基于阈值化的方法没有考虑图像的空间信息,对噪声很敏感。基于分类的方法可以较容易的同时利用多种特征,然而对空间信息利用不足,分割结果的连续性差;基于区域合并的分割方法,主要用于非纹理区域的分割。基于流域变换分割方法能否对纹理图像进行分割,决定于能否得到理想的纹理梯度图。基于MRF场的方法,对类别不多的图像分割精度高,充分利用了统计特性和纹理特性,但是MAP估计计算量大,费时。基于多尺度MRF场的方法,用一个新的多尺度随机场模型来替代MRF模型,并且用从一个新的估计标准推出连续的MAP(SMAP)估计来代替MAP估计,解决了MRF模型的计算问题,但是对斑点抑制的不够充分,且分割结果的边缘细节信息容易丢失。
发明内容
针对上述问题,本发明提出了一种顾及多特征多尺度马尔科夫场(MarkovRandom Field,MRF)的SAR图像分割方法,该方法利用灰度共生矩阵提取原始图像的纹理特征,生成纹理图像,并利用小波变换提取边缘特征,结合纹理和边缘特征对图像进行多尺度分割。该方法提高了SAR图像分割结果的精度,在保持图像大的区域性的同时,也很好的保留了图像的细节信息。
本发明的技术方案是:顾及多特征多尺度MRF场的SAR图像分割方法包括以下步骤:
(1)首先利用灰度共生矩阵对原始影像Y进行纹理分析,得到纹理特征图像Y′;
(2)利用模糊C均值方法对原始影像Y进行聚类分割,获得SAR图像的初始分割结果X(0);
(3)结合初始分割结果X(0),在纹理特征图像Y′上进行多尺度MRF分割,获得进一步分割结果X′(0)。
(4)利用小波变换检测原始影像Y边缘信息,得到边缘特征图像Y″。
(5)将边缘特征图像Y″与分割结果X′(0)进行融合,得到最终分割结果。
本发明的特点是:
(1)在得到的纹理特征图像上进行多尺度分割,保证了图像的整体性,使得分割结果区域性明显,受斑点噪声影响小;
(2)融合了小波提取的边缘特征,很好的保留了图像的边缘等细节信息,使得分割结果错误率小。
附图说明
图1为本发明方法的流程示意图。
图2为本发明方法中多尺度MRF模型层内与层间基团示意图,图中(a)为层内基团点s与其邻域,(b)为层内基团类型,(c)为点s的父节点及其邻域,(d)为层间基团类型。
图3为本发明方法的分割结果比较示意图。
具体实施方式
本发明顾及多特征多尺度MRF场的SAR图像分割方法的流程如图1所示,包括以下步骤:
(1)利用灰度共生矩阵对原始影像Y进行纹理分析,得到纹理特征图像Y′。设图像某一区域有N个灰度值,则对应该区域的灰度共生矩阵是一个N×N阶的矩阵,灰度共生矩阵中元素(i,j)的值表示沿某一方向,由距离d分开的分别具有灰度i和j的两相邻像素出现在图像上的相对频率。这样的矩阵是对称的,并且有一个两相邻像素之间角度关系的函数。既可以在整个图像上计算共生矩阵,也可以再扫描图像的小窗口中计算共生矩阵,本实施例采用的特征统计量有以下几种:
角二阶矩:
K阶逆差矩:
纹理相关性:
熵:
其中, Pij
是灰度共生矩阵中位置(i,j)处元素的值。
灰度共生矩阵提取图像纹理特征可以用以下过程描述:
①采用窗口扫描原图像,分别沿0°、45°、90°和135°四个方向,以两相邻像素距离为1生成共生矩阵,计算出每个像素每个方向相应的统计特征值;
②将四个方向的统计特征值求平均,以替代单个方向的统计特征值,并生成与原图像像素矩阵大小相同的特征值矩阵;
③找出特征值矩阵中的最大值与最小值,以特征值矩阵中各个元素与最小值之差同最大值与最小值之差的比值构成灰度值矩阵,进而转化生成特征图像;
④对这些特征图像像素进行直方图统计分析,发现这些特征量呈指数型分布,可在相应纹理特征图像上考虑选用Gauss分布来描述图像的纹理特征性质。
(2)利用模糊C均值方法对原始影像Y进行聚类分割,获得SAR图像的初始分割结果X(0);
(3)结合初始分割结果X(0),在纹理特征图像Y′上进行多尺度MRF分割,获得进一步分割结果X′(0)。
本实施例采用MLL模型描述标号场,如图2所示,只考虑2阶邻域***(8邻域)N={Ns,s∈S}和双像素集团(即只考虑β1,β2,β3,β4),并且令基团势能系数β1=β2=β3=β4=β,β取经验值。
在同一层内,点s和其邻域的局部相关性如图2(a、b)中的基团来表示,其MLL集团势能函数的计算如下:
层间的基团势能如图2中(c)和(d)所示,其MLL集团势能函数为:
本实施例选用Gauss模型来描述图像的分布,则可用均值μr和σr标准方差来表示其分布规律,似然函数可以表示为:
基于多尺度MRF模型的图像分割采用贝叶斯(Bayes)估计器实现,可由一个优化问题表示:
代价函数C给出了真实配置为X而实际分割结果为x时的代价。在已知图像y的情况下,最小化这一代价的期望,从而得到最佳的分割。本实施例采用序列最大后验概率对SAR图像进行多尺度分割,相应的贝叶斯估计式如下:
(4)小波变换提取边缘特征Y″,在用小波分析提取边缘信息时,采用的是小波模极大值的方法,其相应的模值和幅角公式如下:
小波变换的模值沿幅角方向的极大值点就是可能存在的边缘点,将图像的局部模极大值点连接成极值链,并选择阈值去除短的极值链,剩下的就是图像的边缘。
(5)利用线性加权法融合纹理多尺度分割结果X′(0)与边缘特征Y″,即直接对两幅图像的对应像素点进行加权叠加,得到最终分割结果,如图3所示。该图是对真实SAR图像农田进行的分割,并给出了几种分割算法进行比较,此农田图像大小为256×256像素。在图3中,(a)为原始图像,(b)为遗传算法KSW阈值分割结果,(c)为OSTU算法分割结果,(d)为单尺度MRF-ICM分割结果,(e)为多尺度MRF分割结果,(f)为本发明方法分割结果。(f)中圈A表示分割的整体性较好,圈B箭头所指,表示分割的细节性好。
Claims (2)
1.一种顾及多尺度马尔科夫场的合成孔径雷达图像分割方法,其特征在于包括以下步骤:
(1)对原始图像利用灰度共生矩阵提取纹理信息,在得到的纹理特征图像上进行多尺度分割;
(2)利用小波变换检测原始图像的边缘信息,得到边缘特征图像;
(3)利用线性加权法融合纹理多尺度分割结果与边缘特征,即直接对两幅图像的对应像素点进行加权叠加,得到最终分割结果。
2.根据权利要求1所述顾及多尺度马尔科夫场的合成孔径雷达图像分割方法,其特征在于:步骤(1)中首先利用灰度共生矩阵对原始影像进行纹理分析,得到纹理特征图像,并利用模糊C均值方法对原始影像进行聚类分割,获得SAR图像的初始分割结果,然后在纹理特征图像上结合所得到的初始分割结果,对图像进行多尺度MRF分割。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100631722A CN101587587A (zh) | 2009-07-14 | 2009-07-14 | 顾及多尺度马尔科夫场的合成孔径雷达图像分割方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNA2009100631722A CN101587587A (zh) | 2009-07-14 | 2009-07-14 | 顾及多尺度马尔科夫场的合成孔径雷达图像分割方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101587587A true CN101587587A (zh) | 2009-11-25 |
Family
ID=41371822
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNA2009100631722A Pending CN101587587A (zh) | 2009-07-14 | 2009-07-14 | 顾及多尺度马尔科夫场的合成孔径雷达图像分割方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101587587A (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101894365A (zh) * | 2010-07-13 | 2010-11-24 | 武汉大学 | 一种自适应变分遥感影像融合方法 |
CN102012511A (zh) * | 2010-10-13 | 2011-04-13 | 武汉大学 | 基于单时相单极化高分辨率sar图像的二次成像方法 |
CN102024258A (zh) * | 2010-12-15 | 2011-04-20 | 中国科学院遥感应用研究所 | 具有边界保持特性的遥感图像多尺度分割方法 |
CN102368331A (zh) * | 2011-10-31 | 2012-03-07 | 陈建裕 | 一种集成边缘信息的图像多尺度分割方法 |
CN102436642A (zh) * | 2011-10-24 | 2012-05-02 | 葛文英 | 结合mrf和神经网络的多尺度彩色纹理图像分割方法 |
CN102542570A (zh) * | 2011-12-30 | 2012-07-04 | 北京华航无线电测量研究所 | 一种微波图像中人体隐藏危险物体自动检测方法 |
CN102750700A (zh) * | 2012-06-05 | 2012-10-24 | 南京师范大学 | 一种结合邻域信息的快速鲁棒模糊c均值图像分割方法 |
CN104751192A (zh) * | 2015-04-24 | 2015-07-01 | 中国矿业大学(北京) | 基于图像块共生特征的煤岩识别方法 |
CN105139395A (zh) * | 2015-08-19 | 2015-12-09 | 西安电子科技大学 | 基于小波池化卷积神经网络的sar图像分割方法 |
CN106257498A (zh) * | 2016-07-27 | 2016-12-28 | 中南大学 | 基于异构纹理特征的锌浮选工况状态划分方法 |
CN109255334A (zh) * | 2018-09-27 | 2019-01-22 | 中国电子科技集团公司第五十四研究所 | 基于深度学习语义分割网络的遥感影像地物分类方法 |
CN109932370A (zh) * | 2017-12-15 | 2019-06-25 | 三星显示有限公司 | 利用改善的预处理的用于白色点斑检测的***和方法 |
CN111951282A (zh) * | 2020-08-12 | 2020-11-17 | 辽宁石油化工大学 | 一种基于马尔科夫随机场与区域合并的图像分割改进算法 |
CN112927257A (zh) * | 2021-03-22 | 2021-06-08 | 哈尔滨理工大学 | 一种基于改进mrf的冲压件缺陷图像分割方法 |
-
2009
- 2009-07-14 CN CNA2009100631722A patent/CN101587587A/zh active Pending
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101894365A (zh) * | 2010-07-13 | 2010-11-24 | 武汉大学 | 一种自适应变分遥感影像融合方法 |
CN102012511A (zh) * | 2010-10-13 | 2011-04-13 | 武汉大学 | 基于单时相单极化高分辨率sar图像的二次成像方法 |
CN102012511B (zh) * | 2010-10-13 | 2012-09-05 | 武汉大学 | 基于单时相单极化高分辨率sar图像的二次成像方法 |
CN102024258A (zh) * | 2010-12-15 | 2011-04-20 | 中国科学院遥感应用研究所 | 具有边界保持特性的遥感图像多尺度分割方法 |
CN102024258B (zh) * | 2010-12-15 | 2012-07-25 | 中国科学院遥感应用研究所 | 具有边界保持特性的遥感图像多尺度分割方法 |
CN102436642A (zh) * | 2011-10-24 | 2012-05-02 | 葛文英 | 结合mrf和神经网络的多尺度彩色纹理图像分割方法 |
CN102436642B (zh) * | 2011-10-24 | 2014-03-05 | 葛文英 | 结合mrf和神经网络的多尺度彩色纹理图像分割方法 |
CN102368331A (zh) * | 2011-10-31 | 2012-03-07 | 陈建裕 | 一种集成边缘信息的图像多尺度分割方法 |
CN102542570A (zh) * | 2011-12-30 | 2012-07-04 | 北京华航无线电测量研究所 | 一种微波图像中人体隐藏危险物体自动检测方法 |
CN102750700A (zh) * | 2012-06-05 | 2012-10-24 | 南京师范大学 | 一种结合邻域信息的快速鲁棒模糊c均值图像分割方法 |
CN104751192A (zh) * | 2015-04-24 | 2015-07-01 | 中国矿业大学(北京) | 基于图像块共生特征的煤岩识别方法 |
CN104751192B (zh) * | 2015-04-24 | 2017-12-22 | 中国矿业大学(北京) | 基于图像块共生特征的煤岩识别方法 |
CN105139395A (zh) * | 2015-08-19 | 2015-12-09 | 西安电子科技大学 | 基于小波池化卷积神经网络的sar图像分割方法 |
CN105139395B (zh) * | 2015-08-19 | 2018-03-06 | 西安电子科技大学 | 基于小波池化卷积神经网络的sar图像分割方法 |
CN106257498A (zh) * | 2016-07-27 | 2016-12-28 | 中南大学 | 基于异构纹理特征的锌浮选工况状态划分方法 |
CN106257498B (zh) * | 2016-07-27 | 2019-12-17 | 中南大学 | 基于异构纹理特征的锌浮选工况状态划分方法 |
CN109932370A (zh) * | 2017-12-15 | 2019-06-25 | 三星显示有限公司 | 利用改善的预处理的用于白色点斑检测的***和方法 |
CN109932370B (zh) * | 2017-12-15 | 2023-12-15 | 三星显示有限公司 | 利用改善的预处理的用于白色点斑检测的***和方法 |
CN109255334A (zh) * | 2018-09-27 | 2019-01-22 | 中国电子科技集团公司第五十四研究所 | 基于深度学习语义分割网络的遥感影像地物分类方法 |
CN109255334B (zh) * | 2018-09-27 | 2021-12-07 | 中国电子科技集团公司第五十四研究所 | 基于深度学习语义分割网络的遥感影像地物分类方法 |
CN111951282A (zh) * | 2020-08-12 | 2020-11-17 | 辽宁石油化工大学 | 一种基于马尔科夫随机场与区域合并的图像分割改进算法 |
CN112927257A (zh) * | 2021-03-22 | 2021-06-08 | 哈尔滨理工大学 | 一种基于改进mrf的冲压件缺陷图像分割方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101587587A (zh) | 顾及多尺度马尔科夫场的合成孔径雷达图像分割方法 | |
CN102722883B (zh) | 一种具有空间自适应性的极化sar图像分割方法 | |
CN101493935B (zh) | 基于剪切波隐马尔可夫模型的合成孔径雷达图像分割方法 | |
Berg et al. | Investigation of a hybrid algorithm for sea ice drift measurements using synthetic aperture radar images | |
CN102402685A (zh) | 基于Gabor特征的三马尔可夫场SAR图像分割方法 | |
CN102778680B (zh) | 基于参数化的匀加速运动刚体群目标成像方法 | |
CN103366371A (zh) | 基于k分布和纹理特征的sar图像分割方法 | |
CN107392887B (zh) | 一种基于同质像素点转化的异质遥感图像变化检测方法 | |
CN102867309A (zh) | 基于混合模型的sar图像变化检测方法 | |
CN104732552A (zh) | 基于非平稳条件场的sar图像分割方法 | |
CN104133983A (zh) | 基于加权负一阶rd和的多目标伯努利分布式融合方法 | |
Nazarinezhad et al. | A contextual-based segmentation of compact PolSAR images using Markov Random Field (MRF) model | |
Demir et al. | Extraction of coastlines with fuzzy approach using SENTINEL-1 SAR image | |
CN102496142B (zh) | 基于模糊的三马尔可夫场sar图像分割方法 | |
Wu et al. | Cost-sensitive latent space learning for imbalanced PolSAR image classification | |
CN101847259B (zh) | 基于加权信息熵和马尔可夫随机场的红外目标分割方法 | |
Srivastava et al. | Onboard detection of snow, ice, clouds and other geophysical processes using kernel methods | |
CN101877134B (zh) | 一种机场监视视频目标鲁棒跟踪方法 | |
CN105005987A (zh) | 基于广义gamma分布的SAR图像超像素生成方法 | |
CN104408731A (zh) | 基于区域图和统计相似性编码的sar图像分割方法 | |
CN104156955A (zh) | 一种高分辨率遥感影像的变化检测方法 | |
CN101872481B (zh) | 一种结合可见光图像信息的sar图像快速分割方法 | |
Li et al. | Automatic change detection of urban land-cover based on SVM classification | |
Sivapriya et al. | ViT-DexiNet: a vision transformer-based edge detection operator for small object detection in SAR images | |
CN104199006A (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 | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Open date: 20091125 |