CN107169533B - 一种超像素的概率因子tmf的sar图像海岸线检测算法 - Google Patents

一种超像素的概率因子tmf的sar图像海岸线检测算法 Download PDF

Info

Publication number
CN107169533B
CN107169533B CN201710448649.3A CN201710448649A CN107169533B CN 107169533 B CN107169533 B CN 107169533B CN 201710448649 A CN201710448649 A CN 201710448649A CN 107169533 B CN107169533 B CN 107169533B
Authority
CN
China
Prior art keywords
superpixel
representing
pixel
super
texture
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
Application number
CN201710448649.3A
Other languages
English (en)
Other versions
CN107169533A (zh
Inventor
史晓非
王智罡
冯建德
丁星
刘玲
马海洋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dalian Maritime University
Original Assignee
Dalian Maritime University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dalian Maritime University filed Critical Dalian Maritime University
Priority to CN201710448649.3A priority Critical patent/CN107169533B/zh
Publication of CN107169533A publication Critical patent/CN107169533A/zh
Application granted granted Critical
Publication of CN107169533B publication Critical patent/CN107169533B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

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/13—Edge detection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F18/00—Pattern recognition
    • G06F18/20—Analysing
    • G06F18/24—Classification techniques
    • G06F18/241—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches
    • G06F18/2415—Classification techniques relating to the classification model, e.g. parametric or non-parametric approaches based on parametric or probabilistic models, e.g. based on likelihood ratio or false acceptance rate versus a false rejection rate
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00—Image analysis
    • G06T7/40—Analysis of texture
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Probability & Statistics with Applications (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Image Analysis (AREA)

Abstract

本发明公开了一种超像素的概率因子TMF的SAR图像海岸线检测算法读取合成孔径雷达采集到的SAR图像，生成对应的纹理图像Itexture均匀选取种子点，遍历所述的纹理图像Itexture，将无边缘的局部窗作为邻域点的集合C中的特征点，计算种子点的特征均值和方差，根据均值法更新每一类标签的特征和每一个超像素的特征，并计算出ps,i，进而更新二维向量的标记场X，重复循环纸质超像素保持不变，进行条件迭代，得到海岸的检测结果。本发明不再采用传统的增加特征维数的方法，而是由液体扩散系数的启发提出了用统计量表示图像的方法，提出这种纹理，使用该统计量可以在一定程度上降低SAR图像中海洋和陆地的纹理复杂性，同时在一定程度上增强海洋和陆地的对比度。

Description

一种超像素的概率因子TMF的SAR图像海岸线检测算法
技术领域
本发明涉及一种超像素的概率因子TMF的SAR图像海岸线检测算法，属于海岸线检测领域。
背景技术
海岸线检测是海域使用和管理中的一个重要的环节，也是国家管理海洋和陆地资源的一个重要的手段，因此研究海岸线的变化有着重要的意义，导致海岸线变化的原因主要由以下几方面：(1)地壳的运动，这是海岸线变化的最主要原因。(2)冰川的融化和扩展。合成孔径雷达(Synthetic Aperture Radar,SAR)是一种主动式微波探测器，利用合成孔径原理，信号处理方法和脉冲压缩技术，通过较小尺寸的真实天线孔径合成较大的等效天线孔径来成像。相对于光学成像，合成孔径雷达具有全天候成像等优势，因而使得SAR图像在战略目标识别与探测，灾害控制，国土资源监测，海域使用管理，地图测绘，船舰目标识别等领域起到了重要的作用，虽然SAR图像虽然SAR图像有全天候的特性，但是SAR图像的解析是较为困难的，从海岸线检测角度来看，其原因是由于SAR图像有很强的相干斑，使得较精确提取海岸线成为较困难的问题，因此研究SAR图像的海岸线检测具有重要的意义。
目前为止，常见的用于海岸线检测的技术有：基于邻域相似性局部窗的超像素算法和三重马尔科夫随机场(TMF)算法。但是传统的超像素算法在计算像素特征时采用的是恒定大小的矩形窗作为邻域窗，将邻域窗内所有像素特征的均值作为中心点的特征，当中心点位于边缘处时，此时矩形窗内含有边缘，就会导致计算出的特征很模糊，进而使边缘处超像素的边缘贴合度较低。由于在已有基于超像素算法文献中，未发现基于超像素的区域合并海岸线检测算法。而已有马尔科夫随机场算法具有贴合度不佳的特点，在海岸线检测时特别是低对比度边界容易脱离海岸边缘。但是现有的算法尚不能有效的解决这一问题。
发明内容
本发明针对以上问题的提出，提出一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征在于包括如下步骤：
S1：读取合成孔径雷达采集到的SAR图像，生成对应的纹理图像Itexture；
S2：根据图像的大小输入种子点数k，其取值范围为所述纹理图像Itexture大小的10％-20％，并在所述纹理图像中均匀选取种子点；
S3：遍历所述的纹理图像Itexture，每一个点的7像素×7像素的局部窗内计算局部纹理图像Ipart，判断当前局部窗内是否含有边缘并根据所述局部窗内是否含有边缘确定用于计算特征的邻域点的集合C，将无边缘的局部窗作为邻域点的集合C中的特征点；
S4：根据所述步骤S3所得的邻域点的集合C计算所述步骤S2中所选取的种子点所有点的特征均值和方差；
S5：计算每个邻域点和选取的种子点之间的特征距离dci,j，邻域点和种子点之间位置的欧氏距离dsi,j及种子点与邻域点空间距离和特征距离的加权距离di,j，通过比较di,j的大小将当前所述di,j，将di,j较小的邻域点j归入较小的种子点中形成超像素；根据均值法更新每个选取的种子点特征；
S6：重复步骤S3-S5直至前后两次超像素的边缘不再发生变化，超像素恒定后进行步骤S7；
S7：初始化超像素标签，同时初始化辅助场U；遍历所述合成孔径雷达采集到的SAR图像计算出超像素S属于第i类的概率ps,i初始化二维向量的标记场X并计算第i类超像素的均值μi和第i类超像素方差σi；
S8：根据步骤S7计算出的第i类超像素的均值μi计算Gamma分布
Figure BDA0001321890030000021
S9：根据计算的二维向量的标记场X，辅助场U计算势能函数，并更新超像素和辅助场U；
S10：根据均值法更新每一类标签的特征和每一个超像素的特征，并计算出ps,i，进而更新二维向量的标记场X；
S11：重复步骤S7-S10，直至标记场不变时停止更新；更新后使用ICM算法和SG算法对参数集θ＝{αx,αua,αub}更新，最后通过ICM算法进行条件迭代，得到海岸的检测结果。
进一步的，步骤S3中：当局部窗不含有边缘时，则使用全部的像素值计算特征；当局部窗含有边缘时，则使用和中心点属于同一类的像素集合计算特征。
进一步的，每一个像素点的统计量st形成的图像即所述纹理图像Itexture的计算过程包括：
Figure BDA0001321890030000031
其中，μi表示中心点i局部窗内的均值，Ij表示中心点i局部窗内的第j像素的像素值，Ni表示局部窗内的像素的集合，M表示Ni中像素的个数。
进一步的，所述邻域点和种子点之间的距离dc1,2：
Figure BDA0001321890030000032
其中，μ1表示种子点的统计量st的均值，μ2表示邻域点的统计量st的均值，var1表示种子点的统计量st的方差，var2表示邻域点的统计量st的方差，Wvar表示均值项的权值，W2表示方差项的权值；
所述均值项的权值Wvar：
Figure BDA0001321890030000033
Figure BDA0001321890030000034
Figure BDA0001321890030000035
其中，邻域点和中心点统计量方差的相似程度varmin，μmin表示中心点局部窗内统计量均值和邻域点局部窗内统计量均值的比值的较小值；
所述方差项的权值W2：
Figure BDA0001321890030000036
进一步的，所述邻域点和种子点之间位置的欧氏距离dsi,j：
Figure BDA0001321890030000037
所述邻域点和种子点空间距离和特征距离的加权距离di,j：
Figure BDA0001321890030000038
其中，m表示dci,j和dsi,j之间的权重系数，S表示第一次选取种子点时，两个种子点之间的距离。
进一步的，二维向量的标记场X每一维表示当前超像素属于两类之一的概率为：
Figure BDA0001321890030000041
其中，i表示超像素的标签，s表示超像素的位置；
所述第s位置的超像素和第i类的相似程度第s位置的超像素的特征和第i类特征之间的距离ds,i为：
Figure BDA0001321890030000042
其中，i表示超像素的标签，Isp表示当前超像素的均值，σsp表示当前超像素的方差，μi表示第i类超像素的均值，σi表示第i类超像素方差。
进一步的，将ps,ii∈{1,2}形成一个向量作为标记场X的实现，计算μi和σi：
Figure BDA0001321890030000043
Figure BDA0001321890030000044
标记场的后验概率服从均值为μi，视数为L的Gamma分布为：
Figure BDA0001321890030000045
其中，i表示标签，i的取值范围为i∈{1,2}，L表示SAR图像的视数，Isp表示超像素内像素均值，S表示的是整个图像中所有超像素位置的集合。
进一步的，基于标记场X的势能函数模型为：
Figure BDA0001321890030000046
其中，ps,1表示第s位置超像素属于第一类的概率，ps,2表示第s位置超像素属于第二类的概率，t是s的邻域超像素，pt,1表示超像素t属于第一类的概率，pt,2表示超像素t属于第二类的概率，{αx,αua,αub}表示参数集；
Wedge(xs,xt)计算如下：
Figure BDA0001321890030000051
将初始的超像素图像通过初始聚类将图像中的超像素分成两类即标签值是1和2两类当超像素s和t都属于第一类即1时，
Figure BDA0001321890030000052
进一步的，当超像素s和超像素t的特征均介于第一类和第二类中间且偏向于第一类略大于第二类时，根据辅助场U的值判断其种类；当超像素s和超像素t中有一个与第一/二类相似，所述与第一类相似就是指ps,1/pt,1远大于ps,2/pt,2，与第二类非常相似就是指ps,2/pt,2远大于ps,1/pt,1，另一个与该类比较相似时，根据标签值和辅助场U的值联合判断其种类。
本发明的优点在于：本专利不再采用传统的增加特征维数的方法，而是由液体扩散系数的启发提出了用统计量表示图像的方法，提出这种纹理，使用该统计量可以在一定程度上降低SAR图像中海洋和陆地的纹理复杂性，同时在一定程度上增强海洋和陆地的对比度。同时，在本专利模型中，标记场不再仅仅是一个标签值而是一个概率因子本专利模型中将标记场变成一个二维向量，能够有效的解决贴合度不佳的特点，在海岸线检测时特别是低对比度边界容易脱离海岸边缘的问题。
附图说明
为了更清楚的说明本发明的实施例或现有技术的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的整体流程图。
图2(a)-(d)为本发明得到的海岸线检测图。
具体实施方式
为使本发明的实施例的目的、技术方案和优点更加清楚，下面结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚完整的描述：
如图1所示，一种超像素的概率因子TMF的SAR图像海岸线检测算法，作为优选的实施方式，包括如下步骤：
S1：读取合成孔径雷达采集到的SAR图像，生成对应的纹理图像Itexture；
S2：根据图像的大小输入种子点数k，其取值范围为纹理图像Itexture大小的10％-20％，并在纹理图像Itexture中均匀选取种子点；
S3：遍历的纹理图像Itexture，在选取的每一个点的7像素×7像素的局部窗内计算局部纹理图像Ipart，判断当前局部窗内是否含有边缘并根据局部窗内是否含有边缘确定用于计算特征的邻域点的集合C；
S4：根据步骤S3所得的邻域点的集合C计算步骤S2中所选取的种子点所有点的特征均值和方差；
S5：计算每个邻域点和选取的种子点之间的距离dci,j，邻域点和种子点之间位置的欧氏距离dsi,j及di,j，通过比较di,j的大小将当前di,j，将di,j较小的邻域点j归入较小的种子点中形成超像素；根据均值法更新每个选取的种子点特征；
S6：重复步骤S3-S5直至前后两次超像素的边缘不再发生变化，超像素恒定后进行步骤S7；
S7：初始化超像素标签，同时初始化辅助场U；遍历原图像计算出ps,i，初始化二维向量的标记场X并计算第i类超像素的均值μi和第i类超像素方差σi；
S8：根据步骤S7计算出的第i类超像素的均值μi计算Gamma分布
Figure BDA0001321890030000061
S9：根据计算的二维向量的标记场X，辅助场U计算势能函数，并对超像素和辅助场U进行更新；
S10：根据均值法更新每一类标签的特征和每一个超像素的特征，并计算出ps,i，进而更新二维向量的标记场X；
S11：重复步骤S7-S10，直至标记场不变时停止更新；更新后使用ICM算法和SG算法对参数集θ＝{αx,αua,αub}更新，最后通过ICM算法进行条件迭代，得到海岸的检测结果。
在本实施方式中，步骤S3中：当局部窗不含有边缘则使用全部的像素值计算特征；当局部窗含有边缘时，则使用和中心点属于同一类的像素集合计算特征。
在本实施方式中，每一个像素点的统计量st形成的图像即纹理图像Itexture的计算过程包括：
Figure BDA0001321890030000071
其中，μi表示中心点i局部窗内的均值，Ij表示中心点i局部窗内的第j像素的像素值，Ni表示局部窗内的像素的集合，M表示Ni中像素的个数。可以理解为在其它实施方式中，纹理图像Itexture的可以用其它的物理量代替或表示，只要能够满足能够清楚的区别开所需的纹理图像并且能够有效的表示图像即可。
在本实施方式中，邻域点和种子点之间的距离dc1,2：
Figure BDA0001321890030000072
其中，μ1表示种子点的统计量st的均值，μ2表示邻域点的统计量st的均值，var1表示种子点的统计量st的方差，var2表示邻域点的统计量st的方差，Wvar表示均值项的权值，W2表示方差项的权值；
均值项的权值Wvar：
Figure BDA0001321890030000073
Figure BDA0001321890030000074
Figure BDA0001321890030000075
其中，邻域点和中心点统计量方差的相似程度varmin，μmin表示种子点的统计量均值和邻域点的统计量均值的比值的较小值；方差项的权值W2：
Figure BDA0001321890030000076
在本实施方式中对于dsi,j的计算公式如下：
Figure BDA0001321890030000077
di,j的计算方式为：
Figure BDA0001321890030000078
二维向量的标记场X每一维表示当前超像素属于两类之一的概率为：
Figure BDA0001321890030000079
其中，i表示超像素的标签，s表示超像素的位置；第s位置的超像素和第i类的相似程度。
ds,i的计算方式为：
Figure BDA0001321890030000081
其中，i表示超像素的标签，Isp表示当前超像素的均值，σsp表示当前超像素的方差，μi表示第i类超像素的均值，σi表示第i类超像素方差。
作为优选的实施方式，将ps,ii∈{1,2}形成一个向量作为标记场X的实现，计算μi和σi：
Figure BDA0001321890030000082
Figure BDA0001321890030000083
假设标记场的后验概率服从均值为μi，视数为L的Gamma分布为：
Figure BDA0001321890030000084
其中，i表示标签，i的取值范围为i∈{1,2}，L表示SAR图像的视数，Isp表示超像素内像素均值。
在本实施方式中，基于标记场X的势能函数模型为：
Figure BDA0001321890030000085
其中，ps,1表示第s位置超像素属于第一类的概率，ps,2表示第s位置超像素属于第二类的概率，t是s的邻域超像素，pt,1表示超像素t属于第一类的概率，pt,2表示超像素t属于第二类的概率，{αx,αua,αub}表示参数集；
在本实施方式中，Wedge(xs,xt)计算如下：
Figure BDA0001321890030000086
将初始的超像素图像通过初始聚类将图像中的超像素分成两类即标签值是1和2两类当超像素s和t都属于第一类即1时，
Figure BDA0001321890030000091
在本实施方式中，当超像素s和超像素t的特征均介于第一类和第二类中间且偏向于第一类略大于第二类时，根据辅助场U的值判断其种类；当超像素s和超像素t中有一个与第一/二类相似，与第一类相似就是指ps,1/pt,1远大于ps,2/pt,2，与第二类非常相似就是指ps,2/pt,2远大于ps,1/pt,1，另一个与该类比较相似时，根据标签值和辅助场U的值联合判断其种类。
实施例一：Envisat图像检测结果对比
取种子数k分别为：150，150，250，250。每个邻域点和选取的种子点之间的距离dc与邻域点和种子点之间位置的欧氏距离ds之间的权重系数m为0.5，局部窗大小是5×5，最大迭代次数为5次。对于TMF部分αx的初始值为1，αua的值为0.6，αub的值为0.6，外迭代的迭代次数是10，内迭代的迭代次数是10，用于判断不再发生变化的门限值是0.0001。
对比试验GMRF的参数设置如下：条件迭代模式的迭代次数为5，滤波窗的大小是3×3。
对比试验基于超像素的TMF的参数设置如下：种子数k取值为150，150，250，250，dc与ds之间的权重系数m为0.5，局部窗大小是3×3，最大迭代次数为5次，TMF部分αx的初始值为1，αua的值为0.6，αub的值为0.6，外迭代的迭代次数是10，内迭代的迭代次数是10，用于判断不再发生变化的门限值是0.0001。得到图2(a)-(d)所示为本具体实施方式得到的海岸线检测图。
表1针对Envisat图像的三种算法RMSE对比
Figure BDA0001321890030000092
实施例二：Terra图像检测结果对比与分析
种子数k分别是150，150，150，150。dc与ds之间的权重系数m为0.5，局部窗大小是5×5，最大迭代次数为5次。对于TMF部分αx的初始值为1，αua的值为0.6，αub的值为0.6，外迭代的迭代次数是10，内迭代的迭代次数是10，用于判断不再发生变化的门限值是0.0001。对比算法GMRF的参数设置如下：条件迭代模式的迭代次数为5，滤波窗的大小是3×3。对比算法基于超像素的TMF的参数设置如下：种子数k取值为150，150，150，150，dc与ds之间权重系数m为0.5，局部窗大小是3×3，最大迭代次数为5次。对于TMF部分αx的初始值为1，αua的值为0.6，αub的值为0.6，外迭代的迭代次数是10，内迭代的迭代次数是10，用于判断不再发生变化的门限值是0.0001。
表2针对Terra图像的三种算法RMSE对比
Figure BDA0001321890030000101
实施例三：Radasat图像检测结果对比与分析
种子数k分别是250，250，250，250。dc与ds之间的权重系数m为0.5，局部窗大小是5×5，最大迭代次数为5次。对于TMF部分αx的初始值为1，αua的值为0.6，αub的值为0.6，外迭代的迭代次数是10，内迭代的迭代次数是10，用于判断不再发生变化的门限值是0.0001。对比算法GMRF的参数设置如下：条件迭代模式的迭代次数为5，滤波窗的大小是3×3。对比算法基于超像素的TMF的参数设置如下：种子数k取值为250，250，250，250，dc与ds之间的权重系数m为0.5，局部窗大小是3×3，最大迭代次数为5次。对于TMF部分αx的初始值为1，αua的值为0.6，αub的值为0.6，外迭代的迭代次数是10，内迭代的迭代次数是10，用于判断不再发生变化的门限值是0.0001。
以上，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。

Claims (9)

1.一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征在于包括如下步骤：
S1：读取合成孔径雷达采集到的SAR图像，生成对应的纹理图像Itexture；
S2：根据所述纹理图像Itexture的大小输入种子点数k，其取值范围为所述纹理图像Itexture大小的10％-20％，并在所述纹理图像Itexture中均匀选取种子点；
S3：遍历所述的纹理图像Itexture，在选取的种子点的7像素×7像素的局部窗内计算局部纹理图像Ipart，判断当前局部窗内是否含有边缘并根据所述局部窗内是否含有边缘确定用于计算特征的邻域点的集合C，将无边缘的局部窗作为邻域点的集合C中的特征点；
S4：根据所述步骤S3所得的邻域点的集合C计算所有点的特征均值和方差；
S5：计算每个邻域点和选取的种子点之间的距离dci,j，邻域点和种子点之间位置的欧氏距离dsi,j及种子点与邻域点空间距离和特征距离的加权距离di,j，通过比较di,j的大小将di,j较小的邻域点j归入种子点i形成超像素；根据均值法更新每个选取的种子点特征；
S6：重复步骤S3-S5直至前后两次超像素的边缘不再发生变化，超像素恒定后进行步骤S7；
S7：初始化超像素标签，同时初始化辅助场U；遍历所述合成孔径雷达采集到的SAR图像计算出超像素S属于第i类的概率ps,i，初始化二维向量的标记场X并计算第i类的超像素的均值μi和第i类的超像素方差σi；
S8：根据步骤S7计算出的第i类的超像素的均值μi计算Gamma分布pμi,L(Isp|xi)；其中，L表示SAR图像的视数，Isp表示超像素内像素均值，xi表示超像素i的标记场；
S9：根据计算的二维向量的标记场X，辅助场U计算势能函数，并更新超像素和辅助场U；
S10：根据均值法更新每一类标签的特征和每一个超像素的特征，并计算出ps,i，进而更新二维向量的标记场X；
S11：重复步骤S7-S10，直至超像素标签的大小不变停止更新；更新后使用ICM算法和SG算法对参数集θ＝{αx,αua,αub}更新，最后通过ICM算法进行条件迭代，得到海岸的检测结果。
2.根据权利要求1所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
步骤S3中：当图像不含有边缘时，则使用全部的像素值计算特征；当图像含有边缘时，则使用和中心点属于同一类的像素集合计算特征。
3.根据权利要求1所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
每一个像素点的统计量st形成的图像即所述纹理图像Itexture的计算过程包括：
Figure FDA0002709930420000021
其中，μi表示中心点i局部窗内的均值，Ij表示中心点i局部窗内的第j像素的像素值，Ni表示局部窗内的像素的集合，M表示Ni中像素的个数。
4.根据权利要求1所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
所述邻域点和种子点之间的距离dc1,2：
Figure FDA0002709930420000022
其中，μ1表示种子点的统计量st的均值，μ2表示邻域点的统计量st的均值，var1表示种子点的统计量st的方差，var2表示邻域点的统计量st的方差，Wvar表示均值项的权值，W2表示方差项的权值；
所述均值项的权值Wvar：
Figure FDA0002709930420000023
Figure FDA0002709930420000024
Figure FDA0002709930420000031
其中，邻域点和中心点统计量方差的相似程度varmin，μmin表示表示中心点局部窗内统计量均值和邻域点局部窗内统计量均值的比值的较小值；
所述方差项的权值W2：
Figure FDA0002709930420000032
5.根据权利要求1所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
所述邻域点和种子点之间位置的欧氏距离dsi,j
Figure FDA0002709930420000033
所述邻域点和种子点空间距离和特征距离的加权距离di,j：
Figure FDA0002709930420000034
其中，m表示dci,j和dsi,j之间的权重系数，S表示第一次选取种子点时，两个种子点之间的距离。
6.根据权利要求1所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
二维向量的标记场X每一维表示当前超像素属于两类之一的概率为：
Figure FDA0002709930420000035
其中，i表示超像素的标签，s表示超像素的位置；
所述第s位置的超像素和第i类的相似程度第s位置的超像素的特征和第i类特征之间的距离ds,i为：
Figure FDA0002709930420000036
其中，i表示超像素的标签，Isp表示当前超像素的均值，σsp表示当前超像素的方差，μi表示第i类的超像素的均值，σi表示第i类的超像素方差。
7.根据权利要求1所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
将ps,ii∈{1,2}形成一个向量作为标记场X的实现，计算μi和σi：
Figure FDA0002709930420000041
Figure FDA0002709930420000042
标记场的后验概率服从均值为μi，视数为L的Gamma分布为：
Figure FDA0002709930420000043
其中，i表示标签，i的取值范围为i∈{1,2}，L表示SAR图像的视数，Isp表示超像素内像素均值，S表示的是整个图像中所有超像素位置的集合。
8.根据权利要求1所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
基于标记场X的势能函数模型为：
Figure FDA0002709930420000044
其中，ps,1表示第s位置超像素属于第一类的概率，ps,2表示第s位置超像素属于第二类的概率，t是s的邻域超像素，pt,1表示超像素t属于第一类的概率，pt,2表示超像素t属于第二类的概率，{αx,αua,αub}表示参数集；
Wedge(xs,xt)计算如下：
Figure FDA0002709930420000045
将初始的超像素图像通过初始聚类将图像中的超像素分成两类即标签值是1和2两类当超像素s和t都属于第一类即1时，
Figure FDA0002709930420000051
9.根据权利要求8所述的一种超像素的概率因子TMF的SAR图像海岸线检测算法，其特征还在于：
当超像素s和超像素t的特征均介于第一类和第二类中间且偏向于第一类略大于第二类时，根据辅助场U的值判断其种类；当超像素s和超像素t中有一个与第一/二类相似，所述与第一类相似就是指ps,1/pt,1远大于ps,2/pt,2，与第二类非常相似就是指ps,2/pt,2远大于ps,1/pt,1，另一个与该类比较相似时，根据标签值和辅助场U的值联合判断其种类。
CN201710448649.3A 2017-06-14 2017-06-14 一种超像素的概率因子tmf的sar图像海岸线检测算法 Expired - Fee Related CN107169533B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710448649.3A CN107169533B (zh) 2017-06-14 2017-06-14 一种超像素的概率因子tmf的sar图像海岸线检测算法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710448649.3A CN107169533B (zh) 2017-06-14 2017-06-14 一种超像素的概率因子tmf的sar图像海岸线检测算法

Publications (2)

Publication Number Publication Date
CN107169533A CN107169533A (zh) 2017-09-15
CN107169533B true CN107169533B (zh) 2021-01-15

Family

ID=59818542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710448649.3A Expired - Fee Related CN107169533B (zh) 2017-06-14 2017-06-14 一种超像素的概率因子tmf的sar图像海岸线检测算法

Country Status (1)

Country Link
CN (1) CN107169533B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108537116B (zh) * 2018-03-05 2021-09-24 中国地质大学（武汉） 一种基于多尺度特征的海岸线二级类型提取方法及系统
CN108789431A (zh) * 2018-06-11 2018-11-13 深圳万发创新进出口贸易有限公司 一种智能引导机器人

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102903102A (zh) * 2012-09-11 2013-01-30 西安电子科技大学 基于非局部的三马尔可夫随机场sar图像分割方法
CN104036515A (zh) * 2014-06-30 2014-09-10 西安电子科技大学 基于单幅sar图像的车辆目标高度估计方法
CN104680549A (zh) * 2015-03-24 2015-06-03 西安电子科技大学 基于高阶邻域tmf模型的sar图像变化检测方法
CN106156758A (zh) * 2016-07-29 2016-11-23 辽宁工程技术大学 一种sar海岸图像中海岸线提取方法
CN106446914A (zh) * 2016-09-28 2017-02-22 天津工业大学 基于超像素和卷积神经网络的道路检测
CN106709948A (zh) * 2016-12-21 2017-05-24 浙江大学 一种基于超像素分割的快速双目立体匹配方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102903102A (zh) * 2012-09-11 2013-01-30 西安电子科技大学 基于非局部的三马尔可夫随机场sar图像分割方法
CN104036515A (zh) * 2014-06-30 2014-09-10 西安电子科技大学 基于单幅sar图像的车辆目标高度估计方法
CN104680549A (zh) * 2015-03-24 2015-06-03 西安电子科技大学 基于高阶邻域tmf模型的sar图像变化检测方法
CN104680549B (zh) * 2015-03-24 2017-09-26 西安电子科技大学 基于高阶邻域tmf模型的sar图像变化检测方法
CN106156758A (zh) * 2016-07-29 2016-11-23 辽宁工程技术大学 一种sar海岸图像中海岸线提取方法
CN106446914A (zh) * 2016-09-28 2017-02-22 天津工业大学 基于超像素和卷积神经网络的道路检测
CN106709948A (zh) * 2016-12-21 2017-05-24 浙江大学 一种基于超像素分割的快速双目立体匹配方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
SAR Image Multiclass Segmentation Using a;Peng Zhang等;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20121231;第8卷(第6期);全文 *
SAR图像海岸线检测与地物分类;魏明月;《中国优秀硕士学位论文全文数据库信息科技辑》;20170315;全文 *

Also Published As

Publication number Publication date
CN107169533A (zh) 2017-09-15

Similar Documents

Publication Publication Date Title
CN111723654B (zh) 基于背景建模、YOLOv3与自优化的高空抛物检测方法及装置
CN108573276A (zh) 一种基于高分辨率遥感影像的变化检测方法
CN111027446B (zh) 一种高分辨率影像的海岸线自动提取方法
CN110309781B (zh) 基于多尺度光谱纹理自适应融合的房屋损毁遥感识别方法
Aytekın et al. Unsupervised building detection in complex urban environments from multispectral satellite imagery
CN112598713A (zh) 一种基于深度学习的近岸海底鱼类检测、跟踪统计方法
CN110569782A (zh) 一种基于深度学习目标检测方法
CN106778635A (zh) 一种基于视觉显著性的人体区域检测方法
CN107123150A (zh) 全局颜色对比度检测和分割显著图的方法
CN109035196B (zh) 基于显著性的图像局部模糊检测方法
WO2016165064A1 (zh) 基于多视角学习的鲁棒性前景检测方法
CN106991686A (zh) 一种基于超像素光流场的水平集轮廓跟踪方法
Javed et al. SAR image segmentation based on active contours with fuzzy logic
CN107301649B (zh) 一种基于超像素的区域合并sar图像海岸线检测算法
CN102903102A (zh) 基于非局部的三马尔可夫随机场sar图像分割方法
CN107452013A (zh) 基于Harris角点检测和Sugeno模糊积分的显著性检测方法
CN109345536B (zh) 一种图像超像素分割方法及其装置
CN107169533B (zh) 一种超像素的概率因子tmf的sar图像海岸线检测算法
Gao et al. Agricultural image target segmentation based on fuzzy set
CN115393635A (zh) 一种基于超像素分割以及数据增强的红外小目标检测方法
CN107704864B (zh) 基于图像对象性语义检测的显著目标检测方法
CN110211106B (zh) 基于分段Sigmoid带宽的均值漂移SAR图像海岸线检测方法
Vukadinov et al. An algorithm for coastline extraction from satellite imagery
CN106971402A (zh) 一种基于光学辅助的sar图像变化检测方法
CN107146230A (zh) 基于k‑s距离合并代价的sar图像分割方法

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
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20210115

Termination date: 20210614

CF01 Termination of patent right due to non-payment of annual fee