一种MODIS影像云下海表面绿潮覆盖面积的估测方法技术

本发明专利技术公开了一种MODIS影像云下海表面绿潮覆盖面积的估测方法，包括以下步骤：对遥感影像数据进行预处理，得到研究海域的水体遥感反射率图像和DVI影像图；提取绿潮信息和周围有绿潮覆盖的云信息；利用矩形框确定周围有绿潮分布的云像元，假设其为中心像元，外围含绿潮的相同矩形框为其相邻像元，计算中心像元内的云覆盖面积和各相邻像元内的绿潮覆盖面积；计算相邻像元内的绿潮覆盖比例的均值，将其作为中心像元值，中心像元内的云覆盖面积与中心像元值的乘积即为该云下海表面绿潮的覆盖面积。本发明专利技术所公开的方法能够较好地估测MODIS影像云下海表面绿潮覆盖面积，降低了光学影像上云对绿潮面积估测的影响，可为绿潮的定量分析研究提供技术支持。析研究提供技术支持。

【技术实现步骤摘要】
一种MODIS影像云下海表面绿潮覆盖面积的估测方法


[0001]本专利技术涉及海洋监测
，特别涉及一种MODIS影像云下海表面绿潮覆盖面积的估测方法。

技术介绍

[0002]近年来，由漂浮大型藻类快速繁殖和聚集而引起的绿潮、金潮灾害，在在全球海洋中显著增加，成为一个全球性的海洋生态问题。自2007年以来，每年夏季黄海均会世界上最大的跨区域的绿潮灾害，造成了严重的生态环境和社会经济影响，引起了学者和大众的广泛关注。在诸多的研究中，卫星数据因其多视角、多尺度、大范围的优势发挥了重要作用。通过遥感数据，绿潮的起源及早期发展过程逐步清晰：南黄海绿潮是典型的跨区域生态灾害，起源于苏北浅滩紫菜养殖区，紫菜筏架回收过程中，附着在筏架上的绿藻被刮落；刮落绿藻入海后，在合适的温度、光照、营养盐环境下，迅速繁殖；并在风场和流场的驱动下，由南向北漂移至山东半岛；漂移过程中大规模暴发，至7、8月份开始衰退。
[0003]GOGI、MODIS、HJ-1A/1B、GF-1等光学影像是监测绿潮的主要卫星遥感图像。然而，受东亚夏季风的影响，绿潮发生期(即6月初至8月初)正处于黄海多雨季节，无云影像较少。光学影像容易受到云的影响，对于HJ-1A/1B、GF-1、Sentinel-2等高分辨率的少云影像，云覆盖区域的绿潮信息通过小窗口线性拉伸后，利用ENVI软件能够检测到；但对于MODIS影像，由于分辨率较低，少云影像上，云覆盖区域的绿潮信息容易漏检，导致绿潮面积估测的不准确，从而影响绿潮的定量分析结果。微波比可见光和红外辐射具有更大的波长，因此微波遥感图像受天气的影响较小。但微波遥感数据中的噪声信号对浒苔解译影响明显，利用微波遥感监测绿潮的方法尚处于起步阶段。
[0004]因此，目前关于绿潮面积和生物量的研究大多是先将云腌膜，再进行分析，没有考虑云下是否有漂浮大型藻类覆盖，这会增加估算的绿潮生物量与实际生物量之间的差异。

技术实现思路

[0005]为解决上述技术问题，本专利技术提供了一种MODIS影像云下海表面绿潮覆盖面积的估测方法，以达到降低光学影像上云对绿潮面积估测的影响，为绿潮的定量分析研究提供技术支持的目的。
[0006]为达到上述目的，本专利技术的技术方案如下：
[0007]一种MODIS影像云下海表面绿潮覆盖面积的估测方法，包括以下步骤：
[0008](1)对遥感影像数据进行预处理，得到研究海域的水体遥感反射率图像；
[0009](2)采用DVI算法对水体遥感反射率图像进行处理，得到研究海域的DVI影像图；
[0010](3)基于DVI影像图，利用动态阈值与假彩色图像的目视判读相结合的方法进行绿潮信息提取，检测绿潮像元数量；
[0011](4)基于MODIS影像红光波段的亮度动态阈值法，同时结合绿潮分布情况和假彩色影像的目视判读，提取周围有绿潮分布的云信息，检测周围有绿潮分布的云像元数量；
[0012](5)利用矩形框确定周围有绿潮分布的云像元，假设其为中心像元，外围含绿潮的相同矩形框为其相邻像元，计算中心像元内的云覆盖面积和各相邻像元内的绿潮覆盖面积；
[0013](6)计算相邻像元内的绿潮覆盖比例的均值，将其作为中心像元值，中心像元内的云覆盖面积与中心像元值的乘积即为该云下海表面绿潮的覆盖面积。
[0014]上述方案中，步骤(1)中，对遥感影像数据的预处理包括辐射定标和大气校正。
[0015]上述方案中，步骤(2)中，DVI算法公式如下：
[0016]DVI＝R
NIR-R
RED
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0017]其中，DVI代表影像上各个像元的差值植被指数值，R
NIR
和R
RED
分别代表影像上各个像元在近红外和红光波段的反射率。
[0018]上述方案中，步骤(5)中，云覆盖面积Area
_Cloud
计算公式如下：
[0019]Area
_Cloud
＝PS
×
n
_Cloud
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)
[0020]其中，PS代表卫星影像一个像元对应的地面面积的大小，单位为km2，n
_Cloud
代表检测到的云像元个数。
[0021]上述方案中，步骤(5)中，绿潮覆盖面积Area
_GT
计算公式如下：
[0022]Area
_GT
＝PS
×
N
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0023]其中，PS代表卫星影像一个像元对应的地面面积的大小，单位为km2，N代表检测到的绿潮像元个数。
[0024]上述方案中，步骤(6)中，相邻像元内的绿潮覆盖比例，即绿潮覆盖面积与对应的矩形框面积的比值。
[0025]通过上述技术方案，本专利技术提供的MODIS影像云下海表面绿潮覆盖面积的估测方法是基于低通滤波原理提出的，低通滤波是遥感上常用于图像平滑的滤波器，原理是使用3
×
3的变换核，每个变换核中的元素包含相同的权重，使用外围的均值来代替中心像元值。本方法能够较好地估测MODIS影像云下海表面绿潮覆盖面积，降低了光学影像上云对绿潮面积估测的影响，可为绿潮的定量分析研究提供技术支持。
附图说明
[0026]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。
[0027]图1为2015年7月1日观测到的南黄海海域的绿潮及云的覆盖区域；
[0028]图2为云四周均被绿藻完全包围的情况下，中心像元与相邻像元之间的关系示意图；
[0029]图3a为云的两侧有绿藻分布的情况下，中心像元与相邻像元之间的关系示意图；
[0030]图3b为云的一侧有绿藻分布的情况下，中心像元与相邻像元之间的关系示意图；
[0031]图3c为云的周围绿藻呈零散分布的情况下，中心像元与相邻像元之间的关系示意图；
[0032]图3d为云右侧无绿藻，但与之相接的A云、B云、C云周围有绿藻分布的情况下，中心像元与相邻像元之间的关系示意图。
具体实施方式
[0033]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述。
[0034]本专利技术提供了一种MODIS影像云下海表面绿潮覆盖面积的估测方法，具体实施例如下：
[0035](1)对遥感影像数据进行辐射定标和大气校正，得到研究海域的水体遥感反射率图像；
[0036](2)采用DVI算法对水体遥感反射率图像进行处理，得到研究海域的DVI影像图；DVI算法公式如下：
[0037]DVI＝R
NIR-R
RED
ꢀꢀꢀꢀ本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种MODIS影像云下海表面绿潮覆盖面积的估测方法，其特征在于，包括以下步骤：(1)对遥感影像数据进行预处理，得到研究海域的水体遥感反射率图像；(2)采用DVI算法对水体遥感反射率图像进行处理，得到研究海域的DVI影像图；(3)基于DVI影像图，利用动态阈值与假彩色图像的目视判读相结合的方法进行绿潮信息提取，检测绿潮像元数量；(4)基于MODIS影像红光波段的亮度动态阈值法，同时结合绿潮分布情况和假彩色影像的目视判读，提取周围有绿潮分布的云信息，检测周围有绿潮分布的云像元数量；(5)利用矩形框确定周围有绿潮分布的云像元，假设其为中心像元，外围含绿潮的相同矩形框为其相邻像元，计算中心像元内的云覆盖面积和各相邻像元内的绿潮覆盖面积；(6)计算相邻像元内的绿潮覆盖比例的均值，将其作为中心像元值，中心像元内的云覆盖面积与中心像元值的乘积即为该云下海表面绿潮的覆盖面积。2.根据权利要求1所述的一种MODIS影像云下海表面绿潮覆盖面积的估测方法，其特征在于，步骤(1)中，对遥感影像数据的预处理包括辐射定标和大气校正。3.根据权利要求1所述的一种MODIS影像云下海表面绿潮覆盖面积的估测方法，其特征在于，步骤(2)中，DVI算法公式如下：DVI＝R
NIR-R
RED
ꢀꢀꢀꢀꢀ
(1)其中，DVI代表影像上各个像元的...

【专利技术属性】
技术研发人员：安德玉，邢前国，禹定峰，
申请(专利权)人：山东省科学院海洋仪器仪表研究所，
类型：发明
国别省市：