CN110097101A - 一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 - Google Patents
一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 Download PDFInfo
- Publication number
- CN110097101A CN110097101A CN201910319782.8A CN201910319782A CN110097101A CN 110097101 A CN110097101 A CN 110097101A CN 201910319782 A CN201910319782 A CN 201910319782A CN 110097101 A CN110097101 A CN 110097101A
- Authority
- CN
- China
- Prior art keywords
- image
- sar
- optical
- indicate
- pixel
- 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
- 238000000034 method Methods 0.000 title claims abstract description 25
- 230000004927 fusion Effects 0.000 title claims abstract description 20
- 230000003287 optical effect Effects 0.000 claims abstract description 54
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 25
- 238000005381 potential energy Methods 0.000 claims abstract description 24
- 239000000284 extract Substances 0.000 claims abstract description 17
- 238000012549 training Methods 0.000 claims abstract description 17
- 238000001514 detection method Methods 0.000 claims abstract description 6
- 238000000605 extraction Methods 0.000 claims description 13
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 241001269238 Data Species 0.000 abstract description 4
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 11
- 238000010276 construction Methods 0.000 description 9
- 230000003595 spectral effect Effects 0.000 description 7
- 230000000694 effects Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003325 tomography Methods 0.000 description 2
- 239000003643 water by type Substances 0.000 description 2
- 240000006541 Dactyloctenium aegyptium Species 0.000 description 1
- 238000007635 classification algorithm Methods 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000003550 marker Substances 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000004304 visual acuity Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/253—Fusion techniques of extracted features
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/25—Fusion techniques
- G06F18/254—Fusion techniques of classification results, e.g. of results related to same input data
Landscapes
- Engineering & Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于改进可靠性因子的遥感图像融合与海岸带分类方法，主要包括：读取SAR图像与光学图像并配准；选用FLF线检测算法提取海岸线；将分类区域分割为均匀区域与不均匀区域；进行图像融合并提取灰度共生矩阵纹理并选定训练数据；计算当前像素所属类别的条件概率及每个像素属于每一类的势能，取势能最小的类别标签赋给当前像素作为分类结果；由所述分类结果的每一类别中提取20％的数据作为训练集直到每个像素属于每一类的势能不再变化，则输出最终分类结果。本发明针对不同的地物区域，对不同的传感器数据给出不同的可靠性因子，进而实现海岸带的精确分类。
Description
技术领域
本发明涉及海岸带图像提取技术领域，具体而言，尤其涉及一种基于改 进可靠性因子的遥感图像融合与海岸带分类方法。
背景技术
随着地理、海洋、地球物理以及气象监测等领域研究的深入开展，依靠 卫星传感技术为研究不断提供丰富数据支持成为必然。然而，在许多应用 中，单个传感器所提供的数据存在信息缺失、一致性差，导致提取的数据并 不精确。选自多种传感器的图像源间存在一定的互补关系，同时合并多源数 据不仅能够提供某一个场景的一致性解译，而且可以相应地减轻数据类别不 确定性的影响。因此，图像融合对于遥感图像解译具有重要意义。海岸带区 域作为社会与经济发展的核心地带，其土地利用以及海域使用情况变化剧 烈，因此实现海岸带区域的环境资源的有效检测，必然有助于海岸带地区的 可持续发展。
目前，对于海岸带区域的遥感图像的解译分类大致分为基于像元的分类 方法与基于对象的分类方法。由于受到海陆共同作用，海岸带地物分布情况 复杂，提高了分类算法的解译难度，难以达到理想的应用效果。因此，有许 多学者提出了不同的改进的海岸带分类算法，例如结合地学知识、DEM数 据、图斑空间信息、水文与气象数据等进行分类。但由于这种分类方法对海 岸带一般是分层对各个地物进行分类，流程繁杂，缺乏自动性。另外基于面 向对象的分类方法在很大程度上取决于对象分割的准确度，若是某个对象包 含不同的地物类别，而非单一的地物类别，则很容易错判，从而导致整个对 象中的像素被错判。
现有基于SAR图像与光学图像融合的遥感图像分类算法往往仅考虑不同 传感器数据的不确定性，没有充分利用不同传感器数据对不同地物具有不同 分辨力的特性，从而造成分类效果不佳。
发明内容
根据上述提出的使用不确定性因子作为传感器数据的可靠性度量的不准 确性问题，而提供基于改进可靠性因子的遥感图像融合与海岸带分类方法， 针对不同的地物区域，对不同的传感器数据给出不同的可靠性因子，进而实 现海岸带的精确分类。
本发明采用的技术手段如下：
一种基于改进可靠性因子的遥感图像融合与海岸带分类方法，其特征在 于包括以下步骤：
步骤S1、读取SAR图像与光学图像，并对所述SAR图像与所述光学图 像进行配准，其中所述SAR图像为Sentinel-1图像，所述光学图像为Landsat-8 图像第4、3、2波段合成的真彩色图像；
步骤S2、对Landsat-8图像的第5波段图像提取海岸线，以海岸线为界， 向陆地方向扩300个像素宽的区域作为分类区域，其中线状物的提取选用 FLF线检测算法；
步骤S3、将所述分类区域分割为均匀区域不均匀区域，包括根据SAR 图像的熵纹理信息提取出初步的不均匀区域，综合光学图像的灰度值信息得 到的得到最终的不均匀区域标记场；
步骤S4、对SAR图像与光学图像进行融合分类，提取SAR图像的灰度 共生矩阵纹理并选定训练数据；
步骤S5、计算SAR图像与光学图像融合后的当前像素所属类别的条件 概率，并计算每个像素属于每一类的势能，取势能最小的类别标签赋给当前 像素作为分类结果；
步骤S6、由所述分类结果的每一类别中提取20％的数据作为训练集，包 括对像素所属类别的条件概率排序，提取最高的20％作为新的训练集；
步骤S7、反复执行步骤S5-步骤S6，直到每个像素属于每一类的势能不 再变化，则输出最终分类结果。
较现有技术相比，本发明具有以下优点：
本发明提供了一种均匀性度量可靠性因子遥感图像融合与海岸带分类算 法，在不确定性因子的基础上，引入一种均匀性度量算子，该算子可将图像 分为均匀区域与不均匀区域。然后利用不同传感器对不同地物辨别力的差别 (例如SAR图像对于不均匀区域的辨别力较强，而光学图像含有丰富的光谱 信息，对于没有细节信息的地物类别辨别力较强)，针对不同的地物区域， 对不同的传感器数据给出不同的可靠性因子，进而定义了一种改进的势能函 数，最终实现海岸带的精确分类。综上，应用本发明的技术方案解决了现有技术中仅使用不确定性因子作为传感器数据的可靠性度量而影响分类精度的 问题，可在遥感图像分类领域广泛推广。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实 施例或现有技术描述中所需要使用的附图做以简单地介绍，显而易见地，下 面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在 不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明方法流程图。
图2a为本发明实施例中厦门地区Sentinel-1图像预览。
图2b为本发明实施例中厦门地区谷歌地图图像。
图3a-3l为本发明实施例中本发明算法与对比算法结果对比图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发 明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述， 显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。 以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本 发明及其应用或使用的任何限制。基于本发明中的实施例，本领域普通技术 人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明 保护的范围。
一种基于改进可靠性因子的遥感图像融合与海岸带分类方法，其特征在 于包括以下步骤：
步骤S1、读取SAR图像与光学图像，并对所述SAR图像与所述光学图 像进行配准，其中所述SAR图像为Sentinel-1图像，所述光学图像为Landsat-8 图像第4、3、2波段合成的真彩色图像，其具有丰富的光谱信息，可用于分 类。
步骤S2、对Landsat-8图像的第5波段图像提取海岸线，以海岸线为界， 向陆地方向扩300个像素宽的区域作为分类区域，其中线状物的提取选用 FLF线检测算法。在海岸线提取时，不需要丰富的光谱信息，而且432波段 的干扰较多，第5波段对比度高，利于海岸线提取，因此选用第5波段提取 海岸线。
步骤S3、将所述分类区域分割为均匀区域与不均匀区域，包括根据SAR 图像的熵纹理信息提取出初步的不均匀区域，以及综合光学图像的灰度值信 息得到的得到最终的不均匀区域。
其中根据SAR图像的熵纹理信息提取出初步的不均匀区域包括：
根据公式
计算SAR图像的熵纹理信息，其中，pij表示灰度共生矩阵中的像素概 率，K表示灰度的级数；以及
设定第一提取阈值，将SAR图像的熵纹理信息大于所述第一提取阈值的 SAR图像记为初步的不均匀区域。
综合光学图像的灰度值信息得到的得到最终的不均匀区域记场包括：
计算光学图像的灰度值；
设定第二提取阈值，提取灰度值大于第二提取阈值的二值图像与初步的 不均匀区域相与运算得到最终的不均匀区域标记场。
步骤S4、对SAR图像与光学图像进行融合分类，提取SAR图像的灰度 共生矩阵纹理并选定训练数据；
步骤S5、计算SAR图像与光学图像融合后的当前像素所属类别的条件 概率，并计算每个像素属于每一类的势能，取势能最小的类别标签赋给当前 像素作为分类结果。
计算SAR图像与光学图像融合后的当前像素所属类别的条件概率包括根 据公式
计算SAR图像与光学图像融合后的条件概率，其中Maski＝1表示当前像 素i处于不均匀区域，而Maski≠1表示当前的像素i处于均匀区域，ωB表示人工 建筑区域标签，ωB′表示非人工建筑区域的标签，λ′SAR,i表示SAR图像中第i个 像素的不确定性因子，λSAR,i表示λ′SAR,i归一化后的不确定性因子，λ′Optical,i表示光 学图像中第i个像素的不确定性因子，λOptical,i表示λ′Optical,i归一化后的不确定性 因子，XSAR,i表示SAR图像中第i个像素，XOptical,i表示光学图像中第i个像素， XFused,i表示SAR图像中第i个像素与XOptical,i表示光学图像中第i个像素联合后 的多维度数据，λe表示常数1，λe′表示常数0，ωj表示类别标签，ep表示极小值0.00001。
提取均匀区域非建筑区类别的可靠性因子：
αSAR,i＝λSAR,i+λe，αOptical,i＝λOptical,,i+λe′
均匀区域建筑区类别的可靠性因子：
不均匀区域非建筑区类别的可靠性因子：
αs,i＝λs,i+λe
不均匀区域建筑区类别的可靠性因子为：
αs,i＝λs,i+λe′
并计算每个像素属于每一类的势能包括
根据公式
Udata(XFused)+Usp(C)
＝-{(λSAR,i+λ′e)log(P(XSAR,i|ωB))+(λOptical,i+λ′e)log(P(XOptical,i|ωB))}+Usp(C)
计算若当前像素处于不均匀区域，且判定其属于建筑类别的势能；
根据公式
计算若当前像素处于不均匀区域，且判定其属于非建筑类别的势能；
根据公式
计算若当前像素处于均匀区域，且判定其属于建筑类别的势能；
根据公式
Udata(XFused)+Usp(C)
＝-{(λSAR,i+λe)log(P(XSAR,i|ωB′))+(λOptical,i+λ′)log(P(XOptical,i|ωB′))}+Usp(C)
计算若当前像素处于均匀区域，且判定其属于非建筑类别的势能；
其中，Udata表示数据项Udata(XS)＝log(P(XS|C))，Usp表示空间项 Usp(C)＝log(P(C))，C表示类别标签的集合，C＝{C(i,j)；1≤i≤M,1≤j≤N}是 对应的所有像素的标签集，其中C(i,j)∈{ω1,ω2,...,ωk}。
步骤S6、由所述分类结果的每一类别中提取20％的数据作为训练集，包 括对像素所属类别的条件概率排序，提取最高的20％作为新的训练集；
步骤S7、反复执行步骤S5-步骤S6，直到每个像素属于每一类的势能不 再变化，则输出最终分类结果。
本发明旨在研究一种基于SAR图像与光学图像海岸带融合分类的方法， 以解决现有技术中使用不确定性因子作为传感器数据的可靠性度量的不准确 性。提出了一种均匀性度量可靠性因子遥感图像融合与海岸带分类算法，其 技术思路是：在不确定性因子的基础上，引入一种均匀性度量算子，该算子 可将图像分为均匀区域与不均匀区域。然后利用不同传感器对不同地物辨别 力的差别(例如SAR图像对于建筑区的辨别力较强，而光学图像含有丰富的 光谱信息，对于没有细节信息的地物类别辨别力较强)，针对不同的地物区域，对不同的传感器数据给出不同的可靠性因子，定义了一种改进的势能函 数，进而实现海岸带的精确分类。
下面通过具体的应用与对比实例，对发明进行验证：
本发明将对中国厦门地区进行实验。厦门地区图像大小为762×805(该尺 寸为Sentinel-1图像大小，Landsat-8图像需要配准并上采样后达到该尺寸)， 该区域的分类类别主要有:人工建筑区、山林以及水域。如图2所示，图2(a) 为sentinel-1的预览图像，图中深灰色矩形区域为实验中使用的厦门区域。图 2(b)为谷歌地球中分辨率为2.17m的图像。本发明将参考该光学图像对厦门 实验地区标注Ground Truth图。
本发明算法参数设置：海岸线检测时采用Landsat-8图像第5波段，厦门 种子点数分别为750，最大迭代次数为10。高斯混合模型采用g＝3混合高斯， 图像配准时，控制点对的误差不超过2个像素。海岸带提取选取海岸线陆地 方向上300米宽区域分类。分类时，厦门SAR图像的灰度共生矩阵窗大小为 9，马尔科夫随机场的空间光滑项权重β值为0.01(根据参数性能分析所得)。 SAR图像的阈值选取为0.4(根据参数性能分析所得)。训练样本选择和停止策 略采用与对比算法一致的方式，每次迭代时，选取分类精度前20％的数据作为训练数据集，为了减少算法运行时间并保证算法的准确度，我们选择标签 变化率小于5％时，即停止迭代。
仅使用光学图像分类算法参数设置：采用Landsat-8图像432波段的真彩 色图像分类，厦门的马尔科夫随机场的空间光滑项权重β值分别为0.01。每 次迭代时，选取分类精度前20％的数据作为训练数据集，当标签变化率小于 5％时，即停止迭代。
仅使用SAR图像分类算法参数设置：采用Sentinel-1图像灰度共生矩阵 纹理分类，厦门的SAR图像的灰度共生矩阵窗大小为9。马尔科夫随机场的 空间光滑项权重β值为分别为0.01。每次迭代时，选取分类精度前20％的数 据作为训练数据集，当标签变化率小于5％时，即停止迭代。
对比算法采用基于像素融合分类算法,其参数设置为：多孔小波变换(à trouswavelet transform,ATWT)分解的层数为3，经验模式分解(empirical modedecomposition,EMD)分解的层数为3。融合后得到含SAR图像信息的 高分辨率多光谱图像(基于Landsat-8图像4、3、2波段真彩色图像与全色波 段锐化(pansharpening)后图像)。厦门的马尔科夫随机场的空间光滑项权重β值 为0.01。每次迭代时，选取分类精度前20％的数据作为训练数据集，当标签 变化率小于5％时，即停止迭代。
本发明基于可靠性因子的决策层融合分类算法，其参数设置为：厦门的 马尔科夫随机场的空间光滑项权重β值分别为0.01。每次迭代时，选取分类 精度前20％的数据作为训练数据集，当标签变化率小于5％时，即停止迭代。
图3(a)-(l)分别为厦门地区Landsat-8图像的4、3、2真彩色图像与第8 波段pansharpening图像、Sentinel-1SAR图像、中国专利【201810546924.X】 中所述算法的海岸线检测结果、光学图像提取的海岸带区域、SAR图像提取 的海岸带区域、真实标签(Groundtruth)、本发明算法结果、用于对比的仅使 用光学图像分类结果、用于对比的仅使用SAR图像分类结果、用于对比的基 于不确定性因子的图像融合分类结果、用于对比的基于经验模式分解的图像 融合结果、用于对比的基于像素融合后的分类结果。如图3(g)所示,，本发明算法与真实数据标签的一致性是最高的，其中人工建筑区域与提取的人工建 筑区域基本吻合，因此，本发明发明算法分类结果对于人工建筑区域的提取 结果具有较好的效果。在图3(g)中水域中间比较窄的地方有断层，这是因为 马尔科夫随机场有空间光滑的效果，即数据项的能量本身分辨性能不高，而 空间光滑项的能量差距较大时，便能造成这种断层。如图3(h)，仅仅基于光 学图像的分类结果可以实现非人工建筑区域分类，但是植被区也会出现一些 点状的错误标签。而人工建筑区域则不间隔的错分为植被区，这是因为基于 光学图像的分类时，仅仅采用了光学图像的光谱特性，建筑区域中往往夹杂 着阴影、黑色的道路以及植被引起的浅灰色板块。而仅仅采用SAR图像分类 时，如图3(i)所示，对于建筑区域的分类可以达到很好的效果，而那些不在 建筑区域的类别，由于单纯的使用电磁波的反射信息，而缺乏细节纹理信息， 使在非人工建筑区域产生错分的现象。
对于仅仅使用不确定性作为传感器可靠性度量的图像融合分类结果，如 图3(j)所示，由于不能排除在建筑区时，光学图像的不确定性更小，即光学 图像对图像融合分类的目标函数贡献比SAR图像的贡献更大。而在非人工建 筑区域时，也有可能是SAR图像的不确定性更小，即SAR图像对图像融合 分类的目标函数贡献比光学图像的贡献更大。也就是说，该方法是一种没有 一种导向性的缺乏目的性的融合。本发明算法由于有人工区域标记场的存在， 使得融合具有很强的导向性。即在建筑区域时，发挥出SAR图像的优势(利 用SAR图像对于建筑物的识别率较高的特性)，而在非人工建筑区域，发挥 出光学图像的优势(利用光学图像的光谱特征)。对于基于像素融合的分类方 法，如图3(l)所示,可以看出虽然融合图像(图3(k))相比光学图像而言，有了 SAR图像的特性，即在建筑区域存在SAR图像的电磁反射信息，但是在非人 工区域也带来了一些相干斑噪声。因此，该方法相比仅仅使用光学图像分类 的结果，在植被区明显多了一些建筑区的分类标签，即SAR图像的噪声带来 的效果。也就是说基于经验模式分解的图像融合算法，也应该有一种导向性， 才能避免该问题。即在人工建筑区，SAR图像的特性更突出，而在非人工区 域时，融合的图像中光学图像的光谱特性更加突出。
表1厦门地区实验性能对比
本发明算法检测出结果的精度度量为：制图精度(Product’s Accuraccy, PA)、用户精度(User’s Accuraccy,UA)、总体精度(Overall Accuraccy,OA)以 及Kappa系数。这些参数的值越大，分类效果越好。表1为厦门地区实验性能 的对比，相比于仅使用光学图像分类结果、仅适用SAR图像分类结果、基于不 确定性因子的图像融合分类结果、基于经验模式分解的图像融合结果、基于像 素融合后的分类结果，本发明算法的总体精度与Kappa系数最大(表1第二行所 示)，分别达到93.61％与0.8717，说明本算法的分类效果最优。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对 其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通 技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修 改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替 换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (5)
1.一种基于改进可靠性因子的遥感图像融合与海岸带分类方法，其特征在于包括以下步骤：
步骤S1、读取SAR图像与光学图像，并对所述SAR图像与所述光学图像进行配准，其中所述SAR图像为Sentinel-1图像，所述光学图像为Landsat-8图像第4、3、2波段合成的真彩色图像；
步骤S2、对Landsat-8图像的第5波段图像提取海岸线，以海岸线为界，向陆地方向扩300个像素宽的区域作为分类区域，其中线状物的提取选用FLF线检测算法；
步骤S3、将所述分类区域分割为均匀区域与不均匀区域，包括根据SAR图像的熵纹理信息提取出初步的不均匀区域，综合光学图像的灰度值信息得到的得到最终的不均匀区域标记场；
步骤S4、对SAR图像与光学图像进行融合分类，提取SAR图像的灰度共生矩阵纹理并选定训练数据；
步骤S5、计算SAR图像与光学图像融合后的当前像素所属类别的条件概率，并计算每个像素属于每一类的势能，取势能最小的类别标签赋给当前像素作为分类结果；
步骤S6、由所述分类结果的每一类别中提取20％的数据作为训练集，包括对像素所属类别的条件概率排序，提取最高的20％作为新的训练集；
步骤S7、反复执行步骤S5-步骤S6，直到每个像素属于每一类的势能不再变化，则输出最终分类结果。
2.根据权利要求1所述的方法，其特征在于所述根据SAR图像的熵纹理信息提取出初步的不均匀区域包括：
根据公式
计算SAR图像的熵纹理信息，其中，pij表示灰度共生矩阵中的像素概率，K表示灰度的级数；以及
设定第一提取阈值，将SAR图像的熵纹理信息大于所述第一提取阈值的SAR图像记为初步的不均匀区域。
3.根据权利要求1所述的方法，其特征在于所述综合光学图像的灰度值信息得到的得到最终的不均匀区域标记场包括：
计算光学图像的灰度值；
设定第二提取阈值，提取灰度值大于第二提取阈值的二值图像与初步的不均匀区域相与运算得到最终的不均匀区域标记场。
4.根据权利要求1所述的方法，其特征在于，所述所述计算SAR图像与光学图像融合后的当前像素所属类别的条件概率包括根据公式
计算SAR图像与光学图像融合后的条件概率，其中Maski＝1表示当前像素i处于不均匀区域，而Maski≠1表示当前的像素i处于均匀区域，ωB表示人工建筑区域标签，ωB′表示非人工建筑区域的标签，λ′SAR,i表示SAR图像中第i个像素的不确定性因子，λSAR,i表示λ′SAR,i归一化后的不确定性因子，λ′Optical,i表示光学图像中第i个像素的不确定性因子，λOptical,i表示λ′Optical,i归一化后的不确定性因子，XSAR,i表示SAR图像中第i个像素，XOptical,i表示光学图像中第i个像素，XFused,i表示SAR图像中第i个像素与XOptical,i表示光学图像中第i个像素联合后的多维度数据，λe表示常数1，λe′表示常数0，ωj表示类别标签，ep表示极小值0.00001。
5.根据权利要求1所述的方法，其特征在于，所述计算每个像素属于每一类的势能包括
根据公式
Udata(XFused)+Usp(C)
＝-{(λSAR,i+λ′e)log(P(XSAR,i|ωB))+(λOptical,i+λ′e)log(P(XOptical,i|ωB))}+Usp(C)
计算若当前像素处于不均匀区域，且判定其属于建筑类别的势能；
根据公式
计算若当前像素处于不均匀区域，且判定其属于非建筑类别的势能；
根据公式
计算若当前像素处于均匀区域，且判定其属于建筑类别的势能；
根据公式
Udata(XFused)+Usp(C)
＝-{(λSAR,i+λe)log(P(XSAR,i|ωB′))+(λOptical,i+λ′)log(P(XOptical,i|ωB′))}+Usp(C)
计算若当前像素处于均匀区域，且判定其属于非建筑类别的势能；
其中，Udata表示数据项Udata(XS)＝log(P(XS|C))，Usp表示空间项Usp(C)＝log(P(C))，C表示类别标签的集合，C＝{C(i,j)；1≤i≤M,1≤j≤N}是对应的所有像素的标签集，其中C(i,j)∈{ω1,ω2,...,ωk}。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910319782.8A CN110097101B (zh) | 2019-04-19 | 2019-04-19 | 一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910319782.8A CN110097101B (zh) | 2019-04-19 | 2019-04-19 | 一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110097101A true CN110097101A (zh) | 2019-08-06 |
| CN110097101B CN110097101B (zh) | 2022-09-13 |
Family
ID=67445343
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910319782.8A Expired - Fee Related CN110097101B (zh) | 2019-04-19 | 2019-04-19 | 一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110097101B (zh) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110705578A (zh) * | 2019-10-15 | 2020-01-17 | 中船（浙江）海洋科技有限公司 | 一种海岸线年际演变分析系统及方法 |
| CN111079847A (zh) * | 2019-12-20 | 2020-04-28 | 郑州大学 | 一种基于深度学习的遥感影像自动标注方法 |
| CN111339959A (zh) * | 2020-02-28 | 2020-06-26 | 西南交通大学 | 基于sar和光学影像融合的近海浮筏养殖区提取方法 |
| CN112016441A (zh) * | 2020-08-26 | 2020-12-01 | 大连海事大学 | 基于Radon变换多特征融合的Sentinel-1图像海岸带养殖池提取方法 |
| CN112364289A (zh) * | 2020-11-02 | 2021-02-12 | 首都师范大学 | 一种通过数据融合提取水体信息的方法 |
| CN113112533A (zh) * | 2021-04-15 | 2021-07-13 | 宁波甬矩空间信息技术有限公司 | 基于多分辨率分析的sar-多光谱-高光谱一体化融合方法 |
| CN113538536A (zh) * | 2021-07-21 | 2021-10-22 | 中国人民解放军国防科技大学 | 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 |
| CN113538306A (zh) * | 2021-06-15 | 2021-10-22 | 西安电子科技大学 | 一种sar图像与低分辨率光学图像多图融合方法 |
| CN113838107A (zh) * | 2021-09-23 | 2021-12-24 | 哈尔滨工程大学 | 一种基于稠密连接的异源图像自动配准方法 |
| CN116129145A (zh) * | 2023-04-14 | 2023-05-16 | 广东海洋大学 | 一种高分辨率遥感影像的砂质海岸线提取方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150334315A1 (en) * | 2009-03-02 | 2015-11-19 | Flir Systems, Inc. | Infrared imaging enhancement with fusion |
| CN106372590A (zh) * | 2016-08-29 | 2017-02-01 | 江苏科技大学 | 一种基于机器视觉的海面船只智能跟踪系统及其方法 |
| WO2017071160A1 (zh) * | 2015-10-28 | 2017-05-04 | 深圳大学 | 一种大幅面遥感图像海陆分割的方法及系统 |
| CN107256399A (zh) * | 2017-06-14 | 2017-10-17 | 大连海事大学 | 一种基于Gamma分布超像素算法和基于超像素TMF 的SAR 图像海岸线检测算法 |
| CN109448016A (zh) * | 2018-11-02 | 2019-03-08 | 三亚中科遥感研究所 | 一种基于面向对象及其隶属规则的遥感影像海岸线提取方法 |
-
2019
- 2019-04-19 CN CN201910319782.8A patent/CN110097101B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150334315A1 (en) * | 2009-03-02 | 2015-11-19 | Flir Systems, Inc. | Infrared imaging enhancement with fusion |
| WO2017071160A1 (zh) * | 2015-10-28 | 2017-05-04 | 深圳大学 | 一种大幅面遥感图像海陆分割的方法及系统 |
| CN106372590A (zh) * | 2016-08-29 | 2017-02-01 | 江苏科技大学 | 一种基于机器视觉的海面船只智能跟踪系统及其方法 |
| CN107256399A (zh) * | 2017-06-14 | 2017-10-17 | 大连海事大学 | 一种基于Gamma分布超像素算法和基于超像素TMF 的SAR 图像海岸线检测算法 |
| CN109448016A (zh) * | 2018-11-02 | 2019-03-08 | 三亚中科遥感研究所 | 一种基于面向对象及其隶属规则的遥感影像海岸线提取方法 |
Non-Patent Citations (1)
| Title |
|---|
| 庄翠蓉: "一种基于面向对象的海岸带信息提取技术研究", 《三峡环境与生态》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110705578A (zh) * | 2019-10-15 | 2020-01-17 | 中船（浙江）海洋科技有限公司 | 一种海岸线年际演变分析系统及方法 |
| CN111079847A (zh) * | 2019-12-20 | 2020-04-28 | 郑州大学 | 一种基于深度学习的遥感影像自动标注方法 |
| CN111079847B (zh) * | 2019-12-20 | 2023-05-02 | 郑州大学 | 一种基于深度学习的遥感影像自动标注方法 |
| CN111339959A (zh) * | 2020-02-28 | 2020-06-26 | 西南交通大学 | 基于sar和光学影像融合的近海浮筏养殖区提取方法 |
| CN112016441A (zh) * | 2020-08-26 | 2020-12-01 | 大连海事大学 | 基于Radon变换多特征融合的Sentinel-1图像海岸带养殖池提取方法 |
| CN112016441B (zh) * | 2020-08-26 | 2023-10-13 | 大连海事大学 | 基于Radon变换多特征融合的Sentinel-1图像海岸带养殖池提取方法 |
| CN112364289A (zh) * | 2020-11-02 | 2021-02-12 | 首都师范大学 | 一种通过数据融合提取水体信息的方法 |
| CN113112533B (zh) * | 2021-04-15 | 2022-05-03 | 宁波甬矩空间信息技术有限公司 | 基于多分辨率分析的sar-多光谱-高光谱一体化融合方法 |
| CN113112533A (zh) * | 2021-04-15 | 2021-07-13 | 宁波甬矩空间信息技术有限公司 | 基于多分辨率分析的sar-多光谱-高光谱一体化融合方法 |
| CN113538306A (zh) * | 2021-06-15 | 2021-10-22 | 西安电子科技大学 | 一种sar图像与低分辨率光学图像多图融合方法 |
| CN113538306B (zh) * | 2021-06-15 | 2024-02-13 | 西安电子科技大学 | 一种sar图像与低分辨率光学图像多图融合方法 |
| CN113538536A (zh) * | 2021-07-21 | 2021-10-22 | 中国人民解放军国防科技大学 | 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 |
| CN113538536B (zh) * | 2021-07-21 | 2022-06-07 | 中国人民解放军国防科技大学 | 一种sar图像信息辅助的遥感光学图像浓云检测方法及系统 |
| CN113838107A (zh) * | 2021-09-23 | 2021-12-24 | 哈尔滨工程大学 | 一种基于稠密连接的异源图像自动配准方法 |
| CN113838107B (zh) * | 2021-09-23 | 2023-12-22 | 哈尔滨工程大学 | 一种基于稠密连接的异源图像自动配准方法 |
| CN116129145B (zh) * | 2023-04-14 | 2023-06-23 | 广东海洋大学 | 一种高分辨率遥感影像的砂质海岸线提取方法及系统 |
| CN116129145A (zh) * | 2023-04-14 | 2023-05-16 | 广东海洋大学 | 一种高分辨率遥感影像的砂质海岸线提取方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110097101B (zh) | 2022-09-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110097101A (zh) | 一种基于改进可靠性因子的遥感图像融合与海岸带分类方法 | |
| Liu et al. | Multiscale morphological compressed change vector analysis for unsupervised multiple change detection | |
| Huang et al. | A new building extraction postprocessing framework for high-spatial-resolution remote-sensing imagery | |
| Yin et al. | Hot region selection based on selective search and modified fuzzy C-means in remote sensing images | |
| Wen et al. | A novel automatic change detection method for urban high-resolution remotely sensed imagery based on multiindex scene representation | |
| Li et al. | A review of remote sensing image classification techniques: The role of spatio-contextual information | |
| Wang et al. | Optimal segmentation of high-resolution remote sensing image by combining superpixels with the minimum spanning tree | |
| Huang et al. | Morphological building/shadow index for building extraction from high-resolution imagery over urban areas | |
| Zhou et al. | Multiscale water body extraction in urban environments from satellite images | |
| CN103632363B (zh) | 基于多尺度融合的对象级高分辨率遥感影像变化检测方法 | |
| CN101840581B (zh) | 一种从卫星遥感影像中提取建筑物轮廓的方法 | |
| Chen et al. | Building detection in an urban area using lidar data and QuickBird imagery | |
| CN109063754B (zh) | 一种基于OpenStreetMap的遥感影像多特征联合分类方法 | |
| Liu et al. | Hierarchical semantic model and scattering mechanism based PolSAR image classification | |
| Zhao et al. | Waterbody information extraction from remote-sensing images after disasters based on spectral information and characteristic knowledge | |
| Peng et al. | Object-based change detection from satellite imagery by segmentation optimization and multi-features fusion | |
| CN107330875A (zh) | 基于遥感图像正反向异质性的水体周边环境变化检测方法 | |
| CN106611422B (zh) | 基于素描结构的随机梯度贝叶斯sar图像分割方法 | |
| CN106408030A (zh) | 基于中层语义属性和卷积神经网络的sar图像分类方法 | |
| CN103839267A (zh) | 一种基于形态学建筑物指数的建筑物提取方法 | |
| CN106991686A (zh) | 一种基于超像素光流场的水平集轮廓跟踪方法 | |
| CN109753916A (zh) | 一种归一化差分植被指数尺度转换模型构建方法及装置 | |
| Kurtz et al. | Multi-resolution region-based clustering for urban analysis | |
| CN107992856A (zh) | 城市场景下的高分遥感建筑物阴影检测方法 | |
| Zhu et al. | Loess terrain segmentation from digital elevation models based on the region growth method |
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: 20220913 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |