一种基于耀光反射差异的海面溢油光学遥感检测方法专利检索-溢油环境工程专利检索查询-专利查询网

积极推动地理标志专门立法

2022-03-10 地理标志，立法，知识产权
保护知识产权是对创新最大的激励

2022-03-10 保护知识产权，创新，激励
谢商华：加快制定知识产权基本法

2022-03-10 知识产权基本法
擦亮“双奥之城”品牌

2022-03-10 双奥，知识产权
让冰雪运动“热”力全开

2022-03-10 冰雪运动，知识产权
携手共奋进　走好强国路

2022-03-10 强国，知识产权
坚持创新引领　方能稳中求进

2022-03-10 创新，稳中求进，知识产权
答好“两张卷” 奋进新征程

2022-03-10 知识产权
专家解读政府工作报告中的创新和知识产权相关部署

2022-03-10 政府工作报告，创新，知识产权
今年政府工作报告指出：加强知识产权保护和运用

2022-03-10 政府工作报告，知识产权保护

一种基于耀光反射差异的海面溢油光学遥感检测方法

阅读：835发布：2021-02-23

IPRDB可以提供一种基于耀光反射差异的海面溢油光学遥感检测方法专利检索，专利查询，专利分析的服务。并且本发明涉及一种基于耀光反射差异的海面溢油光学遥感检测方法，利用海面溢油在不同太阳耀光反射下具有不同的光学信号特征，基于传感器观测方向与太阳入射光镜面反射方向的夹角来表示光学遥感影像上太阳耀光反射强弱，并结合临界角理论进行光学遥感影像分带，逐带扫描检测海面疑似溢油异常光学信号，提取出光学遥感影像中疑似海面溢油的像元。本发明所使用光学遥感影像仅需要进行辐射定标与瑞利校正，获得瑞利校正反射率影像；本发明适用于光学遥感影像因幅宽较大导致的耀光反射差异较大的情况，能够快速高效检测影像中的海面疑似溢油像元；该方法能满足海面溢油检测的需求，提高溢油污染治理的工作效率。，下面是一种基于耀光反射差异的海面溢油光学遥感检测方法专利的具体信息内容。

权利要求

1.一种基于耀光反射差异的海面溢油光学遥感检测方法，包括以下步骤：步骤1、数据预处理

对遥感影像数据进行预处理，包括辐射定标和瑞利校正，获取瑞利校正反射率图像Rrc；

步骤2、感兴趣区范围圈定

对于已知或疑似海面溢油事故发生区域，在经过步骤1预处理后的瑞利校正反射率图像Rrc上圈定出可能包含海面溢油的感兴趣区；

步骤3、去除云和陆地

检测云像元和陆地像元并去除，保留海面像元；

步骤4、基于临界角的影像分带

根据传感器观测方向与太阳入射光镜面反射方向的夹角θm对影像进行分带，θm≥130的

0 0 0

像元为弱耀光反射带，12＜θm＜13的像元为临界耀光反射带，θm≤12的像元为强耀光反射带；

步骤5、太阳耀光反射率计算

计算各海面像元的太阳耀光反射率，得到太阳耀光反射率图像LGN；

步骤6、针对影像各分带区域进行海面溢油异常检测

首先，利用下述公式计算归一化图像Rrc(NIR)normalRrc(NIR)normal＝Rrc(NIR)-LGN其中，Rrc(NIR)为步骤2中感兴趣区的近红外波段图像；

然后，以夹角θm为依据对归一化图像Rrc(NIR)normal进行分块扫描，检测溢油异常信号，具体步骤如下：

1)、分别针对弱耀光反射带和强耀光反射带，以夹角θm＝1.0°为间隔进行扫描窗口的划分，并逐一对归一化图像Rrc(NIR)normal进行扫描；

2)、对归一化图像Rrc(NIR)normal中处于同一扫描窗口的像元，统计置信区间为[1％，

99％]内的像元值的累计统计概率；

3)、若扫描窗口位于强耀光反射带，则将累计统计概率大于70％所对应的像元作为溢油异常像元；若扫描窗口位于弱耀光反射带，则将累计统计概率小于30％所对应的像元作为溢油异常像元；

4)、针对每个扫描窗口内的溢油异常像元，对归一化后的反射率值进行自适应阈值分割，对于在强耀光反射带内的溢油异常像元，若归一化后的反射率值大于阈值，则该溢油异常像元判定为溢油像元；对于在弱耀光反射带内的溢油异常像元，若归一化后的反射率值小于阈值，则该溢油异常像元判定为溢油像元；

5)、将各扫描窗口的溢油像元进行合并，得到感兴趣区的溢油检测结果。

2.根据权利要求1所述的基于耀光反射差异的海面溢油光学遥感检测方法，其特征在于：步骤3中，对感兴趣区的瑞利校正反射率图像Rrc的每个像元计算云指数，以所有感兴趣区像元的云指数为样本，使用最大类间方差法求得用于分类的自适应阈值ε0，云指数≥ε0的像元判定为云像元。

3.根据权利要求1所述的基于耀光反射差异的海面溢油光学遥感检测方法，其特征在于：步骤4中，对感兴趣区的瑞利校正反射率图像Rrc的每个像元计算归一化水体指数，以所有感兴趣区像元的归一化水体指数为样本，使用最大类间方差法求得用于分类的自适应阈值ε1，归一化水体指数＜ε1的像元判定为陆地像元。

4.根据权利要求1所述的基于耀光反射差异的海面溢油光学遥感检测方法，其特征在于：步骤4中，根据夹角θm的大小生成分带图像Z：Z＝1的像元属于强耀光反射带，Z＝0的像元属于临界耀光反射带，Z＝-1的像元属于弱耀光反射带。

5.根据权利要求1所述的基于耀光反射差异的海面溢油光学遥感检测方法，其特征在于：所述步骤4中的夹角θm通过下式计算获得：其中，θ0为太阳天顶角，θ为卫星天顶角，为太阳与卫星的相对方位角。

6.根据权利要求1所述的基于耀光反射差异的海面溢油光学遥感检测方法，其特征在于：步骤5中，所述太阳耀光反射率通过下述公式计算，式中，θ0为太阳天顶角，θ为卫星天顶角，为太阳与卫星传感器的相对方位角，σ为海面粗糙度，ω为入射角，r为折射角，β为波浪小侧面倾斜角度；其中入射角ω和折射角r通过下述公式计算获得：式中，n1为光在海水中的折射率，n2为光在大气中的折射率。

7.根据权利要求6所述的基于耀光反射差异的海面溢油光学遥感检测方法，其特征在于：σ2＝0.003+0.00512W，其中W为每个像元的风速。

8.根据权利要求6所述的基于耀光反射差异的海面溢油光学遥感检测方法，其特征在于：n1＝1.34，n2＝1.0。

说明书全文

一种基于耀光反射差异的海面溢油光学遥感检测方法

技术领域

[0001] 本发明涉及遥感海洋溢油监测技术领域，特别是涉及一种基于耀光反射差异的海面溢油光学遥感检测方法。技术背景

[0002] 溢油污染是海洋环境监测的重要对象之一，遥感技术能为海洋环境监测与评估提供重要的技术支持。海洋溢油遥感面对复杂的来源、过程与对象：(1)海洋溢油污染的来源复杂。不同油种的污染，如不同产地的原油、各类油品的中间产品、各类型成品油等；不同平台类型，如海面的油轮、钻井平台，海底的钻井、输油管道，海岸带的码头港口等易发溢油事故对象(溢油的风化扩散过程会有差异)。(2)海洋溢油迁移与风化过程复杂。海洋溢油在风、浪、流等因素的作用下，其数量、组成、理化性质等方面不断地发生变化，主要过程包括扩散、漂移、乳化、蒸发、溶解、吸附沉淀、光氧化和生物降解等，会形成多种不同的溢油风化污染类型。(3)溢油污染遥感监测的对象复杂。这些不同的海洋溢油风化污染类型，对入射光具有不同的光学作用过程及光学响应特征(光谱特征、图型特征、纹理特征等)，是光学遥感鉴别与定量估算的理论基础。快速准确地识别海洋溢油污染的不同类型、估算溢油量，不仅有助于溢油污染的应急处理，也有助于海洋环境与社会经济的损害评估。

[0003] 海洋溢油光学遥感理论研究，近年来取得了较好的发展，光学遥感在该领域的理论与应用研究积累表明：(1)光学遥感探测的分类目标已经明确，即溢油污染形成的海面油膜(不同厚度)、黑色浮油与油水混合物(不同类型与不同浓度)，海底烃渗漏形成的海面油膜与近海表大气碳氢化合物气体异常。(2)这些目标对入射光具有不同的光学作用过程(如反射、吸收、散射、干涉等)，会产生不同的光学响应特征，是光学遥感识别、分类与定量估算的理论基础。(3)在实际应用中，目标介质面(不同类型、折射率与粗糙度的油面与海面)的菲涅尔反射差异，有利于目标探测的同时，也给目标识别分类与定量估算带来诸多不确定性影响。光学遥感技术在本领域的应用特点与优势逐渐清晰，相关基础理论架构逐渐完善，充分展现了其对海洋溢油不同污染类型鉴别及定量估算的能力与潜力。

发明内容

[0004] 本发明要解决的技术问题是：针对溢油污染区域难以人为进入监测，水色背景噪声大，宽幅遥感影像溢油特征复杂的特点，基于太阳耀光反射差异对图像进行分带，检测不同分带内海面溢油的不同特征，对于海面溢油污染监测具有重要意义。

[0005] 为了解决上述技术问题，本发明提出的技术方案是：

[0006] 基于耀光反射差异的海面溢油光学遥感检测方法，包括以下步骤：

[0007] 步骤1、数据预处理

[0008] 对遥感影像数据进行预处理，包括辐射定标和瑞利校正，获取瑞利校正反射率图像Rrc；

[0009] 步骤2、感兴趣区范围圈定

[0010] 对于已知或疑似海面溢油事故发生区域，在经过步骤1预处理后的瑞利校正反射率图像Rrc上圈定出可能包含海面溢油的感兴趣区；

[0011] 步骤3、去除云和陆地

[0012] 检测云像元和陆地像元并去除，保留海面像元；

[0013] 步骤4、基于临界角的影像分带

[0014] 根据传感器观测方向与太阳入射光镜面反射方向的夹角θm对影像进行分带，θm≥13°的像元为弱耀光反射带，12°<θm<13°的像元为临界耀光反射带，θm≤12°的像元为强耀光反射带；

[0015] 步骤5、太阳耀光反射率计算

[0016] 计算各海面像元的太阳耀光反射率，得到太阳耀光反射率图像LGN；

[0017] 步骤6、针对影像各分带区域进行海面溢油异常检测

[0018] 首先，利用下述公式计算归一化图像Rrc(NIR)normal

[0019] Rrc(NIR)normal＝Rrc(NIR)-LGN

[0020] 其中，Rrc(NIR)为步骤2中感兴趣区的近红外波段图像；

[0021] 然后，以夹角θm为依据对归一化图像Rrc(NIR)normal进行分块扫描，检测溢油异常信号，具体步骤如下：

[0022] 1)、分别针对弱耀光反射带和强耀光反射带，以夹角θm＝1.0°为间隔进行扫描窗口的划分，并逐一对归一化图像Rrc(NIR)normal进行扫描；

[0023] 2)、对归一化图像Rrc(NIR)normal中处于同一扫描窗口的像元，统计置信区间为[1％,99％]内的像元值的累计统计概率；

[0024] 3)、若扫描窗口位于强耀光反射带，则将累计统计概率大于70％所对应的像元作为溢油异常像元；若扫描窗口位于弱耀光反射带，则将累计统计概率小于30％所对应的像元作为溢油异常像元；

[0025] 4)、针对每个扫描窗口内的溢油异常像元，对归一化后的反射率值进行自适应阈值分割，对于在强耀光反射带内的溢油异常像元，若归一化后的反射率值大于阈值，则该溢油异常像元判定为溢油像元；对于在弱耀光反射带内的溢油异常像元，若归一化后的反射率值小于阈值，则该溢油异常像元判定为溢油像元；

[0026] 5)、将各扫描窗口的溢油像元进行合并，得到感兴趣区的溢油检测结果。

[0027] 本发明基于临界角理论将影像进行分带，临界角理论表述为：海面溢油在不同的太阳耀光反射下表现不同的光学特征，以θm来表示太阳耀光反射的强弱时，存在临界角：当θm小于该临界角时，海面油膜相对于海水表现为亮特征，当θm大于该临界角时海面油膜相对于海水则表现为暗特征，本发明中12.0°～13.0°为潜在临界角范围；基于此，将光学影像分为强耀光反射带、临界耀光反射带和弱耀光反射带，需要注意的是，本发明中的强弱是相对于临界角进行划分的相对概念。

[0028] 本发明具有如下有益效果：

[0029] 本发明针对海面溢油在不同太阳耀光反射下表现出的不同特征，利用传感器观测方向与镜面反射方向的夹角θm来表示太阳耀光反射差异，并结合临界角理论对光学遥感影像进行分带检测。实际效果表明，本发明能够降低溢油区检测的时间、人力和物力成本，保证检测精度，提高实时检测效率。

附图说明

[0030] 图1为本发明流程示意图。

[0031] 图2为传感器观测方向与太阳入射光镜面反射方向夹角θm示意图。

[0032] 图3为实施例MODIS遥感影像。

[0033] 图4-a为影像云检测图。

[0034] 图4-b为影像海陆标识图。

[0035] 图4-c为掩膜图像。

[0036] 图5为圈定的感兴趣区影像。

[0037] 图6为传感器观测方向与太阳入射光镜面反射方向夹角θm图像。

[0038] 图7为太阳耀光反射率LGN图像。

[0039] 图8-a为未进行LGN归一化影像。

[0040] 图8-b为LGN归一化后影像。

[0041] 图9为海面溢油异常检测结果图。

[0042] 下面结合附图对本发明的较佳实施例进行详细阐述，以使本发明的优点和特征能更易于被本领域技术人员理解，从而对本发明的保护范围做出更为清楚明确的界定。

[0043] 本实施例应用于1景MODIS多光谱影像数据。

[0044] 如图1所示，为本发明实施例基于耀光反射差异的海面溢油光学遥感检测方法流程示意图，包括以下步骤：

[0045] 步骤1、数据预处理

[0046] 对多光谱遥感影像数据进行预处理，包括辐射定标和瑞利校正，获取瑞利校正反射率图像Rrc。如图3所示，MODIS影像因幅宽较大，影像中存在明显的明暗差异，这种差异的存在为海面溢油检测造成了一定的挑战。

[0047] 步骤2、感兴趣区范围圈定

[0048] 对于已知或疑似海面溢油事故发生区域，通过目视(或其他辅助手段)在经过步骤1预处理后的瑞利校正反射率图像Rrc上圈定出可能包含海面溢油的感兴趣区。在本实施例中，所圈定感兴趣如图5所示。

[0049] 步骤3、云检测

[0050] 对感兴趣区的瑞利校正反射率图像Rrc的每个像元，计算云指数(Cloud Index)，获得云指数影像CI。云指数(Cloud Index)采用下述公式计算：

[0051]

[0052] 其中，Rrc(Blue)max、Rrc(Blue)min为Rrc影像的蓝光波段反射率的最大值与最小值；Rrc(IR)max、Rrc(IR)min为Rrc影像的红外波段反射率的最大值与最小值；Rrc(IR)为Rrc影像的红外波段反射率值，Rrc(Blue)为Rrc影像的蓝光波段反射率值。

[0053] 在此基础上基于OSTU(最大类间方差法)求得自适应阈值ε0，并对CI影像进行自动化分割，获得云检测结果影像Cloud_mask，公式如下：

[0054]

[0055] Cloud_mask影像中为1的像元表示云，为0的像元表示非云。如图4-a所示，白色像元为云。

[0056] 步骤4、海陆标识

[0057] 对感兴趣区的瑞利校正反射率图像Rrc的每个像元，计算归一化水体指数(NDWI)，获得水体指数影像WI。归一化水体指数(NDWI)采用下述公式计算：

[0058]

[0059] 其中，Rrc(Green)为Rrc影像的绿光波段反射率；Rrc(NIR)为Rrc影像的近红外波段反射率。

[0060] 在此基础上基于OSTU(最大类间方差法)求得自适应阈值ε1，并对WI影像进行自动化分割，获得海陆标识结果影像Water_mask，公式如下：

[0061]

[0062] 影像中为1的像元表示水体，为0的像元表示陆地。如图4-b所示，白色像元为陆地。

[0063] 步骤5、数据掩膜

[0064] 利用步骤3、4生成的Cloud_mask影像与Water_mask影像将步骤2中圈定出的瑞利校正反射率图像Rrc中的云像元与陆地像元剔除。本例中，生成的掩膜是在步骤3和步骤4生成的云检测和海陆标识的基础上扩展1-2像元的缓冲区合成得来，减少云和陆地对海面疑似溢油异常检测的影响。

[0065] 如图4-c所示，利用掩膜将图像中的云和陆地像元剔除。

[0066] 步骤6、θm计算与基于临界角的影像分带

[0067] θm为传感器观测方向与太阳入射光镜面反射方向的夹角，θm的含义如图2所示。该角度能够有效地表达出影像的太阳耀光反射的强弱，并基于临界角理论进行影像分带。该步骤输入数据为：步骤2中圈定的瑞利校正反射率图像Rrc对应的各像元角度数据(包括太阳天顶角、太阳方位角、传感器天顶角与传感器方位角)。计算获得θm图像。夹角θm通过下式计算获得：

[0068]

[0069] 其中，θ0为太阳天顶角，θ为卫星天顶角，为太阳与卫星的相对方位角。

[0070] 然后根据角度θm的大小，将图像分为强耀光反射带，临界耀光反射带和弱耀光反射带，分别赋值为1，0，-1，生成的θm分带图像Z为：

[0071]

[0072] Z＝1的像元属于强耀光反射带，Z＝0的像元属于临界耀光反射带，Z＝-1的像元属于弱耀光反射带。如图6为卫星观测方向与镜面反射方向的夹角图像θm，其中每个像元的值单位为度。

[0073] 步骤7、太阳耀光反射率计算

[0074] 输入数据为：步骤2中圈定影像对应的各像元风速数据(每个像元风速W，单位m/s)和角度数据(包括太阳天顶角、太阳方位角、传感器天顶角与传感器方位角)，计算太阳耀光反射率图像LGN，与Rrc图像像元一一对应，为下一步的溢油检测提供信息。如图7所示为所计算的太阳耀光反射率图像LGN。

[0075] 本实施例中，太阳耀光反射率通过下述公式计算，

[0076]

[0077]

[0078]

[0079]

[0080] 式中，θ0为太阳天顶角，θ为卫星天顶角，为太阳与卫星传感器的相对方位角，σ为海面粗糙度，ω为入射角，r为折射角，β为波浪小侧面倾斜角度；其中入射角ω和折射角r通过下述公式计算获得：

[0081]

[0082]

[0083] 式中，n1为光在海水中的折射率，n2为光在大气中的折射率。

[0084] 本例中，n1＝1.34，n2＝1.0，σ2＝0.003+0.00512W，其中W为风速。

[0085] 步骤8、不同影像分带的海面溢油异常检测

[0086] 根据步骤6计算获得的θm图像、太阳耀光反射强弱分带图像Z以及步骤7计算的LGN图像，针对不同分带(图像Z中不同的区域)采用不同的海面疑似溢油异常检测方法：

[0087] 首先，利用LGN图像与步骤2中圈定的Rrc图像的近红外波段(本实施例MODIS影像中选择865nm波段)进行归一化计算：

[0088] Rrc(865)normal＝Rrc(865)-LGN。其中，Rrc(865)为步骤2中865nm波段图像，Rrc(865)normal为归一化图像。在本实施例中，如图8-a所示为未进行归一化前图像，图8-b为进行归一化处理后的图像，可见影像中的太阳耀光反射率差异得到有效减弱，为海面溢油检测提供更为精确的信息。

[0089] Rrc(865)normal影像对应于不同太阳耀光反射分带检测不同溢油异常特征。θm<12°的强耀光反射分带内(图像Z中值为1的区域)海面溢油(油膜)与水呈正对比，检测亮特征作为海面溢油异常特征；θm>13°的弱耀光反射分带内(图像Z中值为-1的区域)油膜与水为负对比，检测暗特征；12°<θm<13°的临界耀光反射区内(图像Z中值为0的区域)，做掩膜忽略掉，因为该区域内溢油与海水信号差异较小。

[0090] 考虑到MODIS等遥感影像幅宽较大，影像中存在明显的太阳耀光反射差异，故以θm＝1.0°为间隔对影像进行分块扫描，有益于更精细的海面溢油检测。

[0091] 故本实施例中，以夹角θm为依据对归一化图像Rrc(865)normal进行分块扫描，检测溢油异常信号，具体步骤如下：

[0092] 1)、分别针对弱耀光反射带和强耀光反射带，以夹角θm＝1.0°为间隔进行扫描窗口的划分，并逐一对归一化图像Rrc(865)normal进行扫描；；

[0093] 2)、对归一化图像Rrc(865)normal中处于同一扫描窗口的像元，统计置信区间为[1％,99％]内的像元值的累计统计概率；

[0094] 3)、若扫描窗口位于强耀光反射带(分带图像Z的值为1的区域)，则将累计统计概率大于70％所对应的像元作为溢油异常像元；若扫描窗口位于弱耀光反射带(分带图像Z的值为-1的区域)，则将累计统计概率小于30％所对应的像元作为溢油异常像元；

[0095] 4)、针对每个扫描窗口内的溢油异常像元，对归一化后的反射率值进行最大类间方差法自适应阈值分割，对于在强耀光反射带内的溢油异常像元，若归一化后的反射率值大于阈值，则该溢油异常像元判定为溢油像元，对溢油像元赋值为1，其余像元赋值为0；对于在弱耀光反射带内的溢油异常像元，若归一化后的反射率值小于阈值，则该溢油异常像元判定为溢油像元，对溢油像元赋值为1，其余像元赋值为0；

[0096] 5)、将各扫描窗口的溢油像元进行合并，值为1的像元即为海面溢油异常检测结果。如图9即为本方法海面溢油异常检测结果。

[0097] 以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换，或直接或间接运用在其他相关的技术领域，均同理包括在本发明的专利保护范围内。

[0098] 除上述实施例外，本发明还可以有其他实施方式。凡采用等同替换或等效变换形成的技术方案，均落在本发明要求的保护范围。

标题	发布/更新时间	阅读量
一种溢油切断阀-专利编号CN103358893A	2020-05-13	477
喷油泵的溢油阀-专利编号CN107762698A	2020-05-12	706
溢油锅盖-专利编号CN104188546B	2020-05-11	864
智能溢油器-专利编号CN104352183A	2020-05-12	922
智能溢油器-专利编号CN104352183B	2020-05-12	440
汽车溢油循环器-专利编号CN103241119A	2020-05-13	830
溢油锅盖装置-专利编号CN104188551B	2020-05-13	95
溢油锅盖-专利编号CN104188546A	2020-05-11	385
一种溢油回收船-专利编号CN209836975U	2020-05-11	170
一种海湾溢油清理装置-专利编号CN210310803U	2020-05-11	63

一种基于耀光反射差异的海面溢油光学遥感检测方法

一种基于耀光反射差异的海面溢油光学遥感检测方法

技术领域

发明内容

附图说明

IPRDB

热门服务

关于我们

友情链接

联系方式