CN107609577B - 一种利用随机森林的极化sar海面油膜提取方法 - Google Patents
一种利用随机森林的极化sar海面油膜提取方法 Download PDFInfo
- Publication number
- CN107609577B CN107609577B CN201710729485.1A CN201710729485A CN107609577B CN 107609577 B CN107609577 B CN 107609577B CN 201710729485 A CN201710729485 A CN 201710729485A CN 107609577 B CN107609577 B CN 107609577B
- Authority
- CN
- China
- Prior art keywords
- oil film
- node
- random forest
- follows
- coefficient
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 46
- 238000007637 random forest analysis Methods 0.000 title claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims abstract description 47
- 238000000605 extraction Methods 0.000 claims abstract description 39
- 230000001427 coherent effect Effects 0.000 claims abstract description 30
- 238000011160 research Methods 0.000 claims abstract description 19
- 238000012549 training Methods 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 12
- 239000013535 sea water Substances 0.000 claims abstract description 5
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 238000003066 decision tree Methods 0.000 claims description 41
- 238000004364 calculation method Methods 0.000 claims description 34
- 238000001514 detection method Methods 0.000 claims description 32
- 238000010586 diagram Methods 0.000 claims description 27
- 230000008569 process Effects 0.000 claims description 17
- 238000011156 evaluation Methods 0.000 claims description 16
- 238000005259 measurement Methods 0.000 claims description 9
- 230000010287 polarization Effects 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000010845 search algorithm Methods 0.000 claims description 6
- 230000007246 mechanism Effects 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 3
- 229940125944 selective estrogen receptor degrader Drugs 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 239000000284 extract Substances 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003305 oil spill Substances 0.000 description 1
- 230000003746 surface roughness Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Investigating Or Analysing Materials By Optical Means (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种利用随机森林的极化SAR海面油膜提取方法,步骤如下:获得原始的极化SAR影像,对极化SAR影像预处理,选取极化SAR影像的研究区域,对研究区域进行滤波处理,然后进行相干矩阵的提取,得到相干矩阵和相干矩阵的特征值;利用相干矩阵的特征值对相干矩阵进行计算,得到若干油膜特征;选取油膜样本、类油膜样本和海水样本构成原始样本集,在原始样本集中选取若干子样本集,利用油膜特征对子样本集进行训练,得到随机森林模型;利用随机森林模型中的基尼系数提取油膜特征的贡献度;油膜特征构成油膜特征向量集,基于所述随机森林模型对油膜特征向量集进行分类,获得油膜提取结果。本发明能有效提高油膜提取结果精确度、减小噪声影响。
Description
技术领域
本发明涉及图像处理技术领域,尤其涉及一种利用随机森林的极化SAR海面油膜提取方法。
背景技术
石油在开采运输等过程中发生的泄漏会对海洋生态环境造成严重的破坏,同时造成巨大的经济损失。因此准确的油膜检测和提取在溢油事件发生后能为污染清理和灾后评价提供决策支持。同时,海面自然油膜的位置对海底石油勘探产生指导作用,从而大面积缩小海底石油勘探的作业范围。极化SAR(Polarimetric Synthetic Aperture Radar,PolSAR)影像由于其具有覆盖范围广、穿透性强、全天时全天候的工作优点,使得极化SAR在油膜检测中有独特的优势,成为海面油膜检测的主要途径。
针对极化SAR海面油膜提取方法,现有的提取方法包括:
(1)基于单特征的极化SAR海面油膜提取方法:该方法主要利用某一个极化特征对海面油膜进行提取,在特定的研究区内有较好的效果,但是在不同入射角、不同海况下,得到的检测结果存在较大的差异,存在不适用的情况,单特征的极化SAR海面油膜提取方法无法全面的利用到极化信息;
(2)基于多特征的极化SAR海面油膜提取方法:该方法利用了多个特征,但是存在不知哪个特征在提取中有最大的贡献的问题,对后续的验证和油膜检测指导作用微弱;
(3)利用K均值的极化SAR海面油膜提取方法:该方法利用在单特征和多特征下的海面油膜提取中,是目前较为常用的提取方法,但是利用该方法时不同迭代次数不同阈值的设定会对提取结果产生非常大的影响,在不同的影像上需要不同的阈值设定,适用性较差,同时其噪声和异质的容忍度较差,在提取结果中常出现误分、错分的现象,容易出现较大的错误提取,难以控制其准确度。
发明内容
有鉴于此,本发明提供了一种能有效提高油膜提取结果精确度、减小噪声影响的利用随机森林的极化SAR海面油膜提取方法。
本发明提供一种利用随机森林的极化SAR海面油膜提取方法,包括以下步骤:
S101,获得原始的极化SAR影像,对所述极化SAR影像进行预处理,选取所述极化SAR影像的研究区域,对所述研究区域进行滤波处理,然后利用PolSARpro对研究区域进行相干矩阵的提取,得到相干矩阵和相干矩阵的特征值;
S102,利用所述相干矩阵的特征值对相干矩阵进行计算,得到若干油膜特征;
S103,选取油膜样本、类油膜样本和海水样本构成原始样本集,在所述原始样本集中选取若干子样本集,利用所述油膜特征对子样本集进行训练,得到随机森林模型;
S104,利用所述随机森林模型中的基尼系数提取所述油膜特征的贡献度,所述油膜特征的贡献度用于为提取油膜时的特征选择提供支持;
S105,所述油膜特征构成油膜特征向量集,基于所述随机森林模型对油膜特征向量集进行分类,获得油膜提取结果。
进一步地,步骤S101中,所述研究区域为油膜与油膜类似物同时存在的海面,所述滤波处理的方式选择精细Lee滤波和Sigma Lee滤波中的任一种。
进一步地,步骤S102中,所述油膜特征包括熵、反熵、基准高度、同极化相位差的标准差、一致性系数、几何强度、布拉格散射占比、单次反射特征值差异度、自相似系数和镜面相似系数;
所述熵用以整体描述研究区域内不同散射类型在统计意义上的无序性,所述熵的计算公式为:
式中,H为熵,pi(i=1,2,3)为由相干矩阵的特征值获取的伪概率;
所述反熵用以反映研究区域内不同散射类型在统计意义上的无序性,所述反熵的计算公式为:
式中,A12为反熵,λ1、λ2分别为相干矩阵的第一特征值和第二特征值;
所述基准高度的计算公式为:
式中,PH为基准高度,λ3为相干矩阵的第三特征值,λ1>λ2>λ3;
所述同极化相位差的标准差能反映海面后向散射机制的差异,所述同极化相位差的标准差的计算公式为:
式中,σijCO为同极化相位差的标准差,φhh为水平通道下的相位,φvv为垂直通道下的相位;
所述一致性系数的计算公式为:
式中,μ为一致性系数,T13为相干矩阵中第一行第三列的值,T22为相干矩阵中第二行第二列的值,Span为极化SAR影像的总能量;
所述几何强度的计算公式为:
V=(det(T))1/d
式中,V为几何强度,T为相干矩阵,d为相干矩阵的秩;
所述布拉格散射占比的计算公式为:
其中,η为布拉格散射占比,T11为相干矩阵中第一行第一列的值,T12为相干矩阵中第一行第二列的值;
所述单次反射特征值差异度的计算公式为:
式中,SERD为单次反射特征值差异度,i=1,2,3;
所述自相似系数的计算公式为:
式中,rrrs为自相似系数;
所述镜面相似系数的计算公式为:
式中,rrrm为镜面相似系数。
进一步地,步骤103中,所述得到随机森林模型的具体过程为:
3.1通过Bagging抽样方法对所述原始样本集进行有放回的随机抽样n次,n的取值大于2,得到n个大小为N的子样本集;
3.2将所述每个子样本集作为一棵决策树的训练样本,所述训练样本用作对应决策树的根节点,然后由所述根节点开始,执行一个自顶向下的贪婪搜索算法,得到n棵训练完成的决策树;
3.3将所述n棵训练完成的决策树组合,得到随机森林模型。
进一步地,所述贪婪搜索算法的过程为:
3.2.1计算当前决策树的节点的基尼系数,所述基尼系数的计算公式为:
式中,I(A)为基尼系数,A为节点,P(w)为节点处类别w的样本数占总样本数的比例,w=1、2…Nw-1、Nw,所述节点处的每个样本对应一个类别,若节点处所有样本的类别一致,到步骤3.2.4;若节点处所有样本的类别不一致,到步骤3.2.2;
3.2.2利用所述油膜特征对节点进行划分;
3.2.3对所述节点进行划分的过程中,通过比较所述决策树的基尼系数,选择出使基尼系数下降最快的油膜特征,将该油膜特征作为该节点处的最佳划分,同时记录该油膜特征;判断基尼系数下降的参数为度量参数,所述度量参数的计算公式为:
ΔI(A)=I(A)-PLI(AL)-(1-PL)I(AR)
式中,ΔI(A)为度量参数,PL为划分到左分支节点AL中的样本占节点A中样本数的比例,AR为右分支节点;ΔI(A)的值越大说明节点A处的基尼系数下降越快;
3.2.4确定节点处的最佳划分后,首先寻找节点处是否存在右兄弟节点,若存在,则该右兄弟节点作为新的当前节点;若不存在,则寻找决策树中节点处的下一层最左边节点作为当前节点;确定新的当前节点后,重复执行步骤3.2.1-3.2.4;如果节点处不存在右兄弟节点,且决策树中节点处的下一层也没有其他节点,说明每个节点处的样本的类别一致,则决策树的训练过程结束,得到训练完成的决策树。
进一步地,步骤S104中,所述提取油膜特征的贡献度的过程为:通过遍历所有决策树的节点,计算每个油膜特征对应的基尼系数下降总和,所述基尼系数下降总和为该油膜特征的贡献度。
进一步地,步骤S105中,所述获得油膜提取结果的过程为:
5.1所述油膜特征向量集为:
Q=(q1,q2,q3,…,qN)
式中,Q为油膜特征向量集,qj为油膜特征向量集中第j个像元的特征向量,j=1,2,…,N;N为像元总数;
5.2将所述油膜特征向量集投入到所述随机森林模型中;
5.3对于油膜特征向量集中的每一个像元,随机森林模型中的n棵决策树各有一个投票结果,统计n棵决策树的投票结果,其中得票数最多的类别为该像元的类别,则实现了对每一个像元的分类,从而得到油膜提取结果。
进一步地,利用所述步骤S105中的油膜提取结果提取油膜检测范围图,将所述油膜检测范围图与真实油膜范围图叠加,定义所述油膜检测范围图中与真实油膜范围图的重叠部分的油膜的像元数与真实油膜范围图的油膜的像元数的比值为油膜检测率,定义油膜检测范围图中未与真实油膜范围图重叠的部分的油膜的像元数与真实油膜范围图的油膜的像元数的比值为油膜虚警率。
进一步地,根据所述油膜检测率和油膜虚警率定义评价参数,利用所述评价参数评价油膜检测结果,定义所述评价参数为:
式中,F1为评价参数,DR为油膜检测率,FAR为油膜虚警率,F1∈[0,1],评价参数越接近1表明油膜的检测精度越高。
本发明提供的技术方案带来的有益效果是:本发明利用随机森林算法提取海面油膜,训练速度快且能简单实现,不仅有效提高了油膜提取的效率、增加了油膜和油膜类似物的分辨能力,而且减少了油膜提取中的噪声,同时利用本发明的方法能得到不同油膜特征的贡献度,该贡献度可以为油膜提取时的特征选择提供指导,进一步克服了现有技术中信息利用不全、噪声大、误分严重的情况。
附图说明
图1是本发明一种利用随机森林的极化SAR海面油膜提取方法的流程示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地描述。
请参考图1,本发明的实施例提供了一种利用随机森林的极化SAR海面油膜提取方法,包括以下步骤:
步骤S101,获得原始的极化SAR(Polarimetric Synthetic Aperture Radar,PolSAR)影像,对极化SAR影像进行预处理,选取极化SAR影像的研究区域,对研究区域进行滤波处理,然后利用PolSARpro对研究区域进行相干矩阵的提取,得到相干矩阵和相干矩阵的特征值。
通常选取的研究区域为油膜与油膜类似物同时存在的海面,由于极化SAR影像存在大量的相干斑,为了降低噪声的影响,一般需要对原始的极化SAR影像进行滤波处理,滤波处理的方式包括精细Lee滤波和Sigma Lee滤波,一实施例中,采用精细Lee滤波对极化SAR影像进行滤波处理;当影像分辨率较低时,为了保留影像信息也可以不做滤波处理。
步骤S102,利用相干矩阵的特征值对相干矩阵进行计算,得到若干油膜特征。
具体地,步骤S102中,油膜特征包括熵、反熵、基准高度、同极化相位差的标准差、一致性系数、几何强度、布拉格散射占比、单次反射特征值差异度、自相似系数和镜面相似系数。
熵用以整体描述研究区域内不同散射类型在统计意义上的无序性,其计算公式为:
式中,H为熵,Pi(i=1,2,3)为由相干矩阵的特征值获取的伪概率。
反熵用以反映研究区域内不同散射类型在统计意义上的无序性,其计算公式为:
式中,A12为反熵,λ1、λ2分别为相干矩阵的第一特征值和第二特征值。
基准高度的计算公式为:
式中,PH为基准高度,λ3为相干矩阵的第三特征值,λ1>λ2>λ3。
同极化相位差的标准差能反映海面后向散射机制的差异,其计算公式为:
式中,σijCO为同极化相位差的标准差,φhh为水平通道下的相位,φvv为垂直通道下的相位。
一致性系数的计算公式为:
式中,μ为一致性系数,T13为相干矩阵中第一行第三列的值,T22为相干矩阵中第二行第二列的值,Span为原始的极化SAR影像的总能量。
几何强度的计算公式为:
V=(det(T))1/d
式中,V为几何强度,T为相干矩阵,d为相干矩阵的秩。
布拉格散射占比的计算公式为:
其中,η为布拉格散射占比,T11为相干矩阵中第一行第一列的值,T12为相干矩阵中第一行第二列的值;海面是以布拉格散射为主的散射机制,油膜是非布拉格散射为主的散射机制,因此计算原始的极化SAR影像中的布拉格散射占比可以辨别油膜和海水。
单次反射特征值差异度的计算公式为:
式中,SERD为单次反射特征值差异度,i=1,2,3;单次反射特征值差异度对海表面粗糙度敏感,利用其可以辨别油膜和非油膜。
自相似系数的计算公式为:
式中,rrrs为自相似系数。
镜面相似系数的计算公式为:
式中,rrrm为镜面相似系数。
步骤S103,选取油膜样本、类油膜样本和海水样本构成原始样本集,在原始样本集中选取若干子样本集,利用油膜特征对子样本集进行训练,得到随机森林模型。
具体地,步骤S103中,得到随机森林模型的具体过程为:
3.1通过Bagging抽样方法对原始样本集进行有放回的随机抽样n次,n的取值大于2,得到n个大小为N的子样本集;
3.2将每个子样本集作为一棵决策树的训练样本,训练样本用作对应决策树的根节点,然后由根节点开始,执行一个自顶向下的贪婪搜索算法,得到n棵训练完成的决策树;
贪婪搜索算法的过程为:
3.2.1计算当前决策树的节点的基尼系数,基尼系数的计算公式为:
式中,I(A)为基尼系数,A为节点,P(w)为节点A处类别w的样本数占总样本数的比例,w=1、2…Nw-1、Nw,节点A处的每个样本对应一个类别,由于每个子样本集的生成具有随机性,以及节点A在决策树中层次位置不同,因此节点A处有些类别的样本数可能为0;若节点A处所有样本的类别一致,到步骤3.2.4;若节点A处所有样本的类别不一致,到步骤3.2.2;
3.2.2利用步骤S102中得到的油膜特征对节点A进行划分;
3.2.3对节点A进行划分的过程中,通过比较决策树的基尼系数,选择出使基尼系数下降最快的油膜特征,将该油膜特征作为该节点处的最佳划分,同时记录该油膜特征;判断基尼系数下降的参数为度量参数,度量参数的计算公式为:
ΔI(A)=I(A)-PLI(AL)-(1-PL)I(AR)
式中,ΔI(A)为度量参数,PL为划分到左分支节点AL中的样本占节点A中样本数的比例,AR为右分支节点;ΔI(A)的值越大说明节点A处的基尼系数下降越快;
3.2.4确定节点A处的最佳划分后,首先寻找节点A处是否存在右兄弟节点,若存在,则该右兄弟节点作为新的当前节点;若不存在,则寻找决策树中节点A处的下一层最左边节点作为当前节点;确定新的当前节点后,重复执行步骤3.2.1-3.2.4;如果节点A处不存在右兄弟节点,且决策树中节点A处的下一层也没有其他节点,说明每个节点处的样本的类别一致,则决策树的训练过程结束,得到训练完成的决策树;
3.3将n棵训练完成的决策树组合,得到随机森林模型。
一实施例中,决策树的算法采用GART算法。
步骤S104,利用随机森林模型中的基尼系数提取油膜特征的贡献度,并将油膜特征的贡献度可视化表现,油膜特征的贡献度用于为提取油膜时的特征选择提供支持。
具体地,步骤S104中,提取油膜特征的贡献度的过程为:通过遍历所有决策树的节点,计算每个油膜特征对应的基尼系数下降总和,基尼系数下降总和为该油膜特征的贡献度。
步骤S105,油膜特征构成油膜特征向量集,基于随机森林模型对油膜特征向量集进行分类,获得油膜提取结果。
具体地,步骤S105中,获得油膜提取结果的过程为:
5.1油膜特征向量集为:
Q=(q1,q2,q3,…,qN)
式中,qj为油膜特征向量集Q中第j个像元的特征向量,j=1,2,…,N,N为像元总数;
5.2将油膜特征向量集Q投入到随机森林模型中;
5.3对于油膜特征向量集Q中的每一个像元,随机森林模型中的n棵决策树各有一个投票结果,统计n棵决策树的投票结果,其中得票数最多的类别为该像元的类别,则实现了对每一个像元的分类,从而得到油膜提取结果。
本发明利用步骤S105中的油膜提取结果提取油膜检测范围图,将油膜检测范围图与真实油膜范围图叠加,定义油膜检测范围图中与真实油膜范围图的重叠部分的油膜的像元数与真实油膜范围图的油膜的像元数的比值为油膜检测率,定义油膜检测范围图中未与真实油膜范围图重叠的部分的油膜的像元数与真实油膜范围图的油膜的像元数的比值为油膜虚警率;根据油膜检测率和油膜虚警率定义评价参数,利用评价参数评价油膜检测结果,定义评价参数为:
式中,F1为评价参数,DR为油膜检测率,FAR为油膜虚警率,F1∈[0,1],评价参数越接近1表明油膜的检测精度越高。
本发明利用随机森林算法提取海面油膜,训练速度快且能简单实现,不仅有效提高了油膜提取的效率、增加了油膜和油膜类似物的分辨能力,而且减少了油膜提取中的噪声,同时利用本发明的方法能得到不同油膜特征的贡献度,该贡献度可以为油膜提取时的特征选择提供指导,进一步克服了现有技术中信息利用不全、噪声大、误分严重的情况。
在不冲突的情况下,本文中上述实施例及实施例中的特征可以相互结合。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种利用随机森林的极化SAR海面油膜提取方法,其特征在于,包括以下步骤:
S101,获得原始的极化SAR影像,对所述极化SAR影像进行预处理,选取所述极化SAR影像的研究区域,对所述研究区域进行滤波处理,然后利用PolSARpro对研究区域进行相干矩阵的提取,得到相干矩阵和相干矩阵的特征值;
S102,利用所述相干矩阵的特征值对相干矩阵进行计算,得到若干油膜特征;
S103,选取油膜样本、类油膜样本和海水样本构成原始样本集,在所述原始样本集中选取若干子样本集,利用所述油膜特征对子样本集进行训练,得到随机森林模型;
S104,利用所述随机森林模型中的基尼系数提取所述油膜特征的贡献度,所述油膜特征的贡献度用于为提取油膜时的特征选择提供支持;
S105,所述油膜特征构成油膜特征向量集,基于所述随机森林模型对油膜特征向量集进行分类,获得油膜提取结果。
2.如权利要求1所述的利用随机森林的极化SAR海面油膜提取方法,其特征在于,步骤S101中,所述研究区域为油膜与油膜类似物同时存在的海面,所述滤波处理的方式选择精细Lee滤波和Sigma Lee滤波中的任一种。
3.如权利要求1所述的利用随机森林的极化SAR海面油膜提取方法,其特征在于,步骤S102中,所述油膜特征包括熵、反熵、基准高度、同极化相位差的标准差、一致性系数、几何强度、布拉格散射占比、单次反射特征值差异度、自相似系数和镜面相似系数;
所述熵用以整体描述研究区域内不同散射类型在统计意义上的无序性,所述熵的计算公式为:
式中,H为熵,pi为由相干矩阵的特征值获取的伪概率,i=1,2,3;
所述反熵用以反映研究区域内不同散射类型在统计意义上的无序性,所述反熵的计算公式为:
式中,A12为反熵,λ1、λ2分别为相干矩阵的第一特征值和第二特征值;
所述基准高度的计算公式为:
式中,PH为基准高度,λ3为相干矩阵的第三特征值,λ1>λ2>λ3;
所述同极化相位差的标准差能反映海面后向散射机制的差异,所述同极化相位差的标准差的计算公式为:
式中,σijCO为同极化相位差的标准差,φhh为水平通道下的相位,φvv为垂直通道下的相位;
所述一致性系数的计算公式为:
式中,μ为一致性系数,T13为相干矩阵中第一行第三列的值,T22为相干矩阵中第二行第二列的值,Span为极化SAR影像的总能量;
所述几何强度的计算公式为:
V=(det(T))1/d
式中,V为几何强度,T为相干矩阵,d为相干矩阵的秩;
所述布拉格散射占比的计算公式为:
其中,η为布拉格散射占比,T11为相干矩阵中第一行第一列的值,T12为相干矩阵中第一行第二列的值;
所述单次反射特征值差异度的计算公式为:
式中,SERD为单次反射特征值差异度,i=1,2,3;
所述自相似系数的计算公式为:
式中,rrrs为自相似系数;
所述镜面相似系数的计算公式为:
式中,rrrm为镜面相似系数。
4.如权利要求1所述的利用随机森林的极化SAR海面油膜提取方法,其特征在于,步骤103中,所述得到随机森林模型的具体过程为:
3.1通过Bagging抽样方法对所述原始样本集进行有放回的随机抽样n次,n的取值大于2,得到n个大小为N的子样本集;
3.2将每个所述子样本集作为一棵决策树的训练样本,所述训练样本用作对应决策树的根节点,然后由所述根节点开始,执行一个自顶向下的贪婪搜索算法,得到n棵训练完成的决策树;
3.3将所述n棵训练完成的决策树组合,得到随机森林模型。
5.如权利要求4所述的利用随机森林的极化SAR海面油膜提取方法,其特征在于,所述贪婪搜索算法的过程为:
3.2.1计算当前决策树的节点的基尼系数,所述基尼系数的计算公式为:
式中,I(A)为基尼系数,A为节点,P(w)为节点处类别w的样本数占总样本数的比例,w=1、2…Nw-1、Nw,所述节点处的每个样本对应一个类别,若节点处所有样本的类别一致,到步骤3.2.4;若节点处所有样本的类别不一致,到步骤3.2.2;
3.2.2利用所述油膜特征对节点进行划分;
3.2.3对所述节点进行划分的过程中,通过比较所述决策树的基尼系数,选择出使基尼系数下降最快的油膜特征,将该油膜特征作为该节点处的最佳划分,同时记录该油膜特征;判断基尼系数下降的参数为度量参数,所述度量参数的计算公式为:
ΔI(A)=I(A)-PLI(AL)-(1-PL)I(AR)
式中,ΔI(A)为度量参数,PL为划分到左分支节点AL中的样本占节点A中样本数的比例,AR为右分支节点;ΔI(A)的值越大说明节点A处的基尼系数下降越快;
3.2.4确定节点处的最佳划分后,首先寻找节点处是否存在右兄弟节点,若存在,则该右兄弟节点作为新的当前节点;若不存在,则寻找决策树中节点处的下一层最左边节点作为当前节点;确定新的当前节点后,重复执行步骤3.2.1-3.2.4;如果节点处不存在右兄弟节点,且决策树中节点处的下一层也没有其他节点,说明每个节点处的样本的类别一致,则决策树的训练过程结束,得到训练完成的决策树。
6.如权利要求5所述的利用随机森林的极化SAR海面油膜提取方法,其特征在于,步骤S104中,所述提取油膜特征的贡献度的过程为:通过遍历所有决策树的节点,计算每个油膜特征对应的基尼系数下降总和,所述基尼系数下降总和为该油膜特征的贡献度。
7.如权利要求6所述的利用随机森林的极化SAR海面油膜提取方法,其特征在于,步骤S105中,所述获得油膜提取结果的过程为:
5.1所述油膜特征向量集为:
Q=(q1,q2,q3,…,qN)
式中,Q为油膜特征向量集,qj为油膜特征向量集中第j个像元的特征向量,j=1,2,…,N;N为像元总数;
5.2将所述油膜特征向量集投入到所述随机森林模型中;
5.3对于油膜特征向量集中的每一个像元,随机森林模型中的n棵决策树各有一个投票结果,统计n棵决策树的投票结果,其中得票数最多的类别为该像元的类别,则实现了对每一个像元的分类,从而得到油膜提取结果。
8.如权利要求1所述的利用随机森林的极化SAR海面油膜提取方法,其特征在于,利用所述步骤S105中的油膜提取结果提取油膜检测范围图,将所述油膜检测范围图与真实油膜范围图叠加,定义所述油膜检测范围图中与真实油膜范围图的重叠部分的油膜的像元数与真实油膜范围图的油膜的像元数的比值为油膜检测率,定义油膜检测范围图中未与真实油膜范围图重叠的部分的油膜的像元数与真实油膜范围图的油膜的像元数的比值为油膜虚警率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710729485.1A CN107609577B (zh) | 2017-08-23 | 2017-08-23 | 一种利用随机森林的极化sar海面油膜提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710729485.1A CN107609577B (zh) | 2017-08-23 | 2017-08-23 | 一种利用随机森林的极化sar海面油膜提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107609577A CN107609577A (zh) | 2018-01-19 |
CN107609577B true CN107609577B (zh) | 2020-05-01 |
Family
ID=61065864
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710729485.1A Expired - Fee Related CN107609577B (zh) | 2017-08-23 | 2017-08-23 | 一种利用随机森林的极化sar海面油膜提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107609577B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108932520B (zh) * | 2018-04-26 | 2021-12-07 | 武汉大学 | 结合先验概率估计的sar影像水体概率制图方法 |
CN108805186B (zh) * | 2018-05-29 | 2020-11-17 | 北京师范大学 | 一种基于多维显著特征聚类的sar图像圆形油库检测方法 |
CN110646795B (zh) * | 2019-09-16 | 2021-09-24 | 武汉大学 | 一种简缩极化sar的海洋溢油检测方法和*** |
CN111025291B (zh) * | 2019-11-06 | 2021-07-16 | 中国石油大学(华东) | 一种基于全极化sar新特征的海洋溢油检测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609709A (zh) * | 2012-02-03 | 2012-07-25 | 清华大学 | 一种基于极化sar数据融合的海面溢油分割方法 |
CN103092645A (zh) * | 2012-12-31 | 2013-05-08 | 武汉中地数码科技有限公司 | 一种基于微内核技术的地理空间信息应用***及其实现方法 |
CN105866775A (zh) * | 2016-03-28 | 2016-08-17 | 南京信息工程大学 | 一种基于全极化合成孔径雷达图像的海面溢油检测方法 |
CN105931256A (zh) * | 2016-06-03 | 2016-09-07 | 中国地质大学(武汉) | 基于cuda的大幅面遥感影像快速分割方法 |
CN106845347A (zh) * | 2016-12-19 | 2017-06-13 | 航天恒星科技有限公司 | 海上油膜识别提取方法及*** |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
BR112015003444A2 (pt) * | 2012-08-20 | 2017-07-04 | Behavioral Recognition Sys Inc | método e sistema para detectar óleo em superfície de mar |
US10037477B2 (en) * | 2015-08-31 | 2018-07-31 | Massachusetts Institute Of Technology | Combined intensity and coherent change detection in images |
-
2017
- 2017-08-23 CN CN201710729485.1A patent/CN107609577B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102609709A (zh) * | 2012-02-03 | 2012-07-25 | 清华大学 | 一种基于极化sar数据融合的海面溢油分割方法 |
CN103092645A (zh) * | 2012-12-31 | 2013-05-08 | 武汉中地数码科技有限公司 | 一种基于微内核技术的地理空间信息应用***及其实现方法 |
CN105866775A (zh) * | 2016-03-28 | 2016-08-17 | 南京信息工程大学 | 一种基于全极化合成孔径雷达图像的海面溢油检测方法 |
CN105931256A (zh) * | 2016-06-03 | 2016-09-07 | 中国地质大学(武汉) | 基于cuda的大幅面遥感影像快速分割方法 |
CN106845347A (zh) * | 2016-12-19 | 2017-06-13 | 航天恒星科技有限公司 | 海上油膜识别提取方法及*** |
Non-Patent Citations (2)
Title |
---|
《Oil spill detection analyses based on small patch mergence algorithm of SAR image 》;Liu Xiuguo et al;;《The 1st International Conference on Information Science and Engineering》;20091231;第1352-1355页; * |
《面向高分辨率极化SAR四分量分解的相干矩阵自适应估计》;杨帅 等;;《***工程与电子技术》;20160131;第38卷(第1期);第53-59页; * |
Also Published As
Publication number | Publication date |
---|---|
CN107609577A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107609577B (zh) | 一种利用随机森林的极化sar海面油膜提取方法 | |
CN111291639B (zh) | 基于分层变分自编码的跨源舰船特征融合学习与识别方法 | |
CN110533631A (zh) | 基于金字塔池化孪生网络的sar图像变化检测方法 | |
CN103440508B (zh) | 基于视觉词袋模型的遥感图像目标识别方法 | |
CN111429496B (zh) | 一种顾及统计特性的时序PolSAR影像无监督变化检测方法 | |
CN109977968A (zh) | 一种深度学习分类后比较的sar变化检测方法 | |
CN110991257A (zh) | 基于特征融合与svm的极化sar溢油检测方法 | |
Shanjiang et al. | Classification of sea and land waveforms based on deep learning for airborne laser bathymetry | |
Zhang et al. | Automatic recognition of oil industry facilities based on deep learning | |
Xiu et al. | 3D semantic segmentation for high-resolution aerial survey derived point clouds using deep learning | |
CN113111706B (zh) | 一种面向方位角连续缺失的sar目标特征解缠与识别方法 | |
CN116206203B (zh) | 基于SAR与Dual-EndNet的溢油检测方法 | |
Alohali et al. | Automated fault detection in the Arabian Basin | |
Chen et al. | Class-aware domain adaptation for coastal land cover mapping using optical remote sensing imagery | |
Zhao et al. | An anchor-free method for arbitrary-oriented ship detection in SAR images | |
CN113191996A (zh) | 一种遥感影像变化检测方法、装置及其电子设备 | |
CN115456957B (zh) | 一种全尺度特征聚合的遥感影像变化检测的方法 | |
Suidong et al. | Improving SSD for detecting small target in Remote Sensing Image | |
Li et al. | A new automatic oceanic mesoscale eddy detection method using satellite altimeter data based on density clustering | |
Lei et al. | Multi-feature fusion sonar image target detection evaluation based on particle swarm optimization algorithm | |
Yu et al. | A lightweight ship detection method in optical remote sensing image under cloud interference | |
Wang et al. | Point cloud classification and accuracy analysis based on feature fusion | |
Tan et al. | Divide to attend: A multiple receptive field attention module for object detection in remote sensing images | |
Liu et al. | Peaks fusion assisted early-stopping strategy for overhead imagery segmentation with noisy labels | |
Concea et al. | Texture Analysis for Images with Forested Areas |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200501 |