法律状态公告日
法律状态信息
法律状态
2020-07-24
授权
授权
2017-11-10
实质审查的生效 IPC(主分类):G06K9/00 申请日:20170527
实质审查的生效
2017-10-13
公开
公开
技术领域
本发明涉及一种深圳海域初级生产力遥感反演方法及系统。
背景技术
深圳市位于广东省南部沿海,陆域范围为22°26'59"——22°51'49"N,113°45'44"——114°37'21",东临大亚湾和大鹏湾,西濒珠江口和伶仃洋,南隔深圳河与香港相望,北接东莞和惠州;海域连接南海与太平洋,主要划分为珠江口、深圳湾、大鹏湾和大亚湾四个海区,海岸线全长257公里,海域面积1145平方公里。属南亚热带海洋性季风气候,沿岸海区年均表层水温18~28℃,近岸海域主要污染物为石油类、无机氮和活性磷酸盐,近12年共发生赤潮62次,累积发生面积621.85km2,主要发生在2~6月,已成为该海域日益突出的生态问题。通过现有技术方法对深圳市大亚湾和珠江口等海区进行现场调查后发现,叶绿素和初级生产力的周年变化明显呈现双周期型,且珠江口的初级生产力远高于大亚湾。大亚湾叶绿素a与初级生产力在春夏季变化趋势较一致,而大鹏澳海区却没有明显的季节规律,虽然其叶绿素a含量在近岸养殖区较高,但季节上的变化却十分不明显。通过以往相关技术方法的研究可知,深圳市海域的净初级生产力在自然区域上呈现出较大的差别,因而采用同一标准进行深圳市海域的整体研究往往不能得到更为可靠的结果,所以适当划分研究海区,将深圳市海域化整为零不失为一种有效的处理手段。
现有技术对海洋净初级生产力的估算主要集中在其估算模式、时空分布和影响因素等几个方面,相应的技术方法主要有较为传统的黑白瓶法、同位素固碳法、基于海洋卫星遥感的经验模型、半经验模型、BPM(Bedford Productivity Model)、LPCM(Laboratoire dePhysique et Chimie Marines)、VGPM(Vertically Generalized Production Model)等估算模型。
现有的方法大多采用传统的现场测定方法来研究海洋净初级生产力,并且其研究尺度多锁定在大亚湾等局部海湾,采样耗时长且难以做到大面积采样;另一方面,为了解决传统现场测定方法的效率问题,也有一些新的技术方法涌现,譬如利用遥感数据来反演海洋净初级生产力,虽然这些遥感反演技术的时间连续性较好且便于分析大尺度空间内海洋净初级生产力的分布特点,但是以VGPM模型为首的遥感反演手段的适用尺度往往差强人意,在面临小尺度的近海海域问题上往往不能得到很好的结果,再加上不但近海海域上情况复杂,连近海海域的不同海区也可能面临遥感反演上的差异性,除此之外,长时间序列的海洋净初级生产力估算难度较大,参与运算的参数繁多,模型十分复杂,因而对于以深圳海域为例的研究区域而言,现有的技术都存在着一定的局限性。基于以上原因,以往的相似技术方法也多选择东海、太湖、珠江口等区域为主要研究区,而关于深圳近岸海域的相关NPP遥感反演技术方法则较少。
发明内容
有鉴于此,有必要提供一种深圳海域初级生产力遥感反演方法及系统。
本发明提供一种深圳海域初级生产力遥感反演方法,该方法包括如下步骤:a.输入Landsat TM遥感影像数据;b.根据输入的Landsat TM遥感影像数据,计算深圳海域的归一化差分植被指数数据NDVI;c.将计算得到的NDVI数据进行线性压缩,得到NDVI′;d.计算上述得到的NDVI′数据的梯度值
其中,所述Landsat TM遥感影像数据包括7个谱段:B1蓝绿谱段、B2绿谱段、B3红谱段、B4近红外谱段、B5近短波红外谱段、B6热红外谱段、B7近短波红外谱段。
所述的步骤b具体包括:
采用如下公式计算深圳海域的归一化差分植被指数数据NDVI:
NDVI=(B4-B3)/(B4+B3)
其中,NDVI表示归一化差分植被指数;B3表示Landsat>4表示Landsat>
所述的步骤c具体包括:
将NDVI数据的数据范围由NDVI∈[-1,1]调整为NDVI∈[0,1],线性压缩公式为:
NDVIi′=(χi-χmin)/(χmax-χmin),i=1、2、3、…、n
其中,NDVI′表示线性压缩后的NDVI值;χi表示每个像元的NDVI值;χmin表示所有像元中的NDVI的最小值;χmsx表示所有像元中的NDVI的最大值;n表示该NDVI数据的像元数量。
所述的步骤d具体包括:
利用公式:
dx(i,j)=NDVI′(i+1,j)-NDVI′(i,j)
dy(i,j)=NDVI′(i,j+1)-NDVI′(i,j)
其中,
本发明提供一种深圳海域初级生产力遥感反演系统,该系统包括该系统包括输入模块、归一化模块、线性压缩模块、梯度值计算模块、控制点模块、插值运算模块以及反演结果模块,其中:所述输入模块用于输入Landsat TM遥感影像数据;所述归一化模块用于根据输入的Landsat TM遥感影像数据,计算深圳海域的归一化差分植被指数数据NDVI;所述线性压缩模块用于将计算得到的NDVI数据进行线性压缩,得到NDVI′;所述梯度值计算模块用于计算上述得到的NDVI′数据的梯度值
其中,所述的Landsat TM遥感影像数据包括7个谱段:B1蓝绿谱段、B2绿谱段、B3红谱段、B4近红外谱段、B5近短波红外谱段、B6热红外谱段、B7近短波红外谱段。
所述的归一化模块具体用于:
采用如下公式计算深圳海域的归一化差分植被指数数据NDVI:
NDVI=(B4-B3)/(B4+B3)
其中,NDVI表示归一化差分植被指数;B3表示Landsat>4表示Landsat>
所述的线性压缩模块具体用于:
将NDVI数据的数据范围由NDVI∈[-1,1]调整为NDVI∈[0,1],线性压缩公式为:
NDVIi′=(χi-χmin)/(χmax-χmin),i=1、2、3、…、n
其中,NDVI′表示线性压缩后的NDVI值;χi表示每个像元的NDVI值;χmin表示所有像元中的NDVI的最小值;χmax表示所有像元中的NDVI的最大值;n表示该NDVI数据的像元数量。
所述的梯度值计算模块具体用于:
利用公式:
dx(i,j)=NDVI′(i+1,j)-NDVI′(i,j)
dy(i,j)=NDVI′(i,j+1)-NDVI′(i,j)
其中,
本发明针对深圳海域净初级生产力的时空分布特点,以Landsat TM数据作为插值控制点的选择背景,生成含有小尺度纹理信息的控制点集,在此基础上对MODIS/AQUA卫星数据及其产品进行反距离加权插值操作,实现对粗空间分辨率遥感数据进行降尺度的处理,之后在VGPM模型的基础上对其必要参数进行了优化,构造出了一个在深圳海域尺度上更为精准的OVGPM模型,将只适用于一类水体遥感反演的MODIS/AQUA数据产品扩展至以深圳近海海域为例的二级水体中去,首次尝试将MODIS/AQUA数据产品和小尺度的海区联系在一起,不仅扩展了MODIS/AQUA数据产品的适用范围,更有效地利用了现有技术手段和资源。除此之外,本发明还将深圳海域区域进行分区反演处理,最大程度的适应深圳海域净初级生产力“西高东低,北强南弱”的分布特点,最终实现更为精准、简易、快速的深圳市海域净初级生产力的反演。本发明反演结果数据精度较高,反演流程简单易操作,虽然所用参数较多,但都是可以通过MODIS/AQUA数据产品直接获取的遥感数据。
附图说明
图1为本发明深圳海域初级生产力遥感反演方法的流程图;
图2为本发明深圳海域初级生产力遥感反演系统的硬件架构图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细的说明。
参阅图1所示,是本发明深圳海域初级生产力遥感反演方法较佳实施例的作业流程图。
步骤S1,输入Landsat TM遥感影像数据。其中:
所述Landsat TM遥感影像数据包括7个谱段,分别为:3个可见光谱段、1个近红外谱段、2个近短波红外谱段、1个热红外谱段,所述3个可见光谱段包括:蓝绿谱段、绿谱段、红谱段;并将所述7个谱段依次命名为:B1(蓝绿谱段)、B2(绿谱段)、B3(红谱段)、B4(近红外谱段)、B5(近短波红外谱段)、B6(热红外谱段)、B7(近短波红外谱段)。
步骤S2,根据输入的Landsat TM遥感影像数据,计算研究海域的归一化差分植被指数数据NDVI。具体而言:
通过Landsat TM遥感影像数据,计算研究海域的归一化差分植被指数数据(NDVI),计算公式如下:
NDVI=(B4-B3)/(B4+B3)
其中,NDVI表示归一化差分植被指数;B3表示Landsat>4表示Landsat>
步骤S3,将计算得到的NDVI数据进行线性压缩,得到NDVI′。具体而言:
将NDVI数据的数据范围由NDVI∈[-1,1]调整为NDVI∈[0,1],线性压缩公式为:
NDVI′′=(χi-χmin)/(χmax-χmin),i=1、2、3、…、n
其中,NDVI′表示线性压缩后的NDVI值;χi表示每个像元的NDVI值;χmin表示所有像元中的NDVI的最小值;χmax表示所有像元中的NDVI的最大值;n表示该NDVI数据的像元数量。
步骤S4,计算上述得到的NDVI′数据的梯度值
计算NDVI′数据的梯度grad(NDVI),计算公式如下:
dx(i,j)=NDVI′(i+1,j)-NDVI′(i,j)
dy(i,j)=NDVI′(i,j+1)-NDVI′(i,j)
其中,
步骤S5,在深圳的四大主体海域随机布设插值控制点作为一般控制点。其中:所述四大主体海域分别是珠江口海区、深圳湾海区、大鹏湾海区和大亚湾海区。具体而言:
本实施例中,依据深圳海域净初级生产力“西高东低,北强南弱”的分布特点,在梯度数据grad(NDVI)上划定处理区域,处理区域包括了深圳市海域的四大主体海域,分别是珠江口海区、深圳湾海区、大鹏湾海区和大亚湾海区。在划定好的四大主体海域,按照珠江口海区60、深圳湾海区45、大鹏湾海区100、大亚湾海区95的分配方式在各自区域上随机布设共计300个插值控制点来作为一般控制点。
步骤S6,按照
本实施例中,为了保留小尺度数据上的纹理信息,在划定好的四大主体海域上选择
步骤S7,对输入的MODIS/AQUA数据中的月平均海表温度数据、月平均海水漫衰减系数数据、月平均海水光合有效辐射数据、月平均海水叶绿素a浓度数据,分别进行反距离加权插值运算。具体包括:
步骤S71,对输入的MODIS/AQUA数据中的月平均海表温度数据,进行反距离加权插值运算,得到水体最大光合速率:
首先,输入MODIS/AQUA三级数据产品中的月平均海表温度数据(SST),空间分辨率为4km×4km,根据数据头文件将数据亮度值转化为摄氏度(℃);
然后,将输入的MODIS/AQUA月平均海表温度数据(SST)进行IDW插值运算,将该数据的空间分辨率由4km×4km提高至30m×30m;
最后,依据公式计算水体最大光合速率
f(SST)=-3.27×10-8SST7+3.41×10-6SST6-1.35×10-4SST5
+2.46×10-3SST4-0.0205SST3+0.0617SST2
+0.2749SST+C3
其中,
步骤S72,对输入的MODIS/AQUA数据中的月平均海水漫衰减系数数据,进行反距离加权插值运算,得到真光层深度。
首先,输入MODIS/AQUA三级数据产品中的月平均海水漫衰减系数数据(Kd490),空间分辨率为4km×4km,单位是m-1;
然后,将输入的MODIS/AQUA月平均海水漫衰减系数数据(Kd490)进行IDW插值运算,将该数据的空间分辨率由4km×4km提高至30m×30m;
最后,依据公式计算真光层深度Zeu:
Zeu=C1/Kd490
ρ=Rrs(488)/Rrs(547)
其中,Zeu表示真光层深度,单位是m;Kd490表示海水漫衰减系数,单位是m-1;Rrs(488)表示MODIS/AQUA数据488nm波段亮度值;Rrs(547)表示MODIS/AQUA数据547nm波段亮度值;C1表示常系数,通常可以取值2.7834,可以根据用户实测数据进行调整;C2表示常系数,通常可以取值0.0166,可以根据用户实测数据进行调整。
步骤S73,对输入的MODIS/AQUA数据中的月平均海水光合有效辐射数据,进行反距离加权插值运算,得到表面日光合有效辐射。
首先,输入MODIS/AQUA三级数据产品中的月平均海水光合有效辐射数据(PAR)来代替PPeu反演公式中的E0,其空间分辨率为4km×4km,单位是Ein/(m2·d);
然后,将输入的MODIS/AQUA月平均海水光合有效辐射数据(PAR)进行IDW插值运算,将该数据的空间分辨率由4km×4km提高至30m×30m。
步骤S74,对输入的MODIS/AQUA数据中的月平均海水叶绿素a浓度数据,进行反距离加权插值运算,得到垂直叶绿素a浓度。
首先,输入MODIS/AQUA三级数据产品中的月平均海水叶绿素a浓度数据(Chl-a)来代替PPeu反演公式中的Copt,其空间分辨率为4km×4km,单位是mg/m3;
然后,将输入的MODIS/AQUA月平均海水叶绿素a浓度数据(Chl-a)进行IDW插值运算,将该数据的空间分辨率由4km×4km提高至30m×30m。
步骤S8,对输入的研究区域反演时间当日的光照周期数据,进行反距离加权插值运算。具体而言:
输入研究区域反演时间当日的光照周期数据(Dirr),该数据为离散常数数值,单位是小时(h),该数据可由当地天文台查询获得。
步骤S9,根据步骤S7、S8中进行反距离加权插值运算的结果,利用OVGPM模型计算得到深圳海域初级生产力遥感反演结果PPeu。具体而言:
利用OVGPM模型计算PPeu,公式如下:
其中,PPeu是海水表层到真光层的初级生产力,以每平方米产生的碳的毫克数计,单位是mg/m2;
参阅图2所示,是本发明深圳海域初级生产力遥感反演系统10的硬件架构图。该系统包括:输入模块101、归一化模块102、线性压缩模块103、梯度值计算模块104、控制点模块105、插值运算模块106以及反演结果模块107。
所述输入模块101用于输入Landsat TM遥感影像数据。其中:
所述Landsat TM遥感影像数据包括7个谱段,分别为:3个可见光谱段、1个近红外谱段、2个近短波红外谱段、1个热红外谱段,所述3个可见光谱段包括:蓝绿谱段、绿谱段、红谱段;并将所述7个谱段依次命名为:B1(蓝绿谱段)、B2(绿谱段)、B3(红谱段)、B4(近红外谱段)、B5(近短波红外谱段)、B6(热红外谱段)、B7(近短波红外谱段)。
所述归一化模块102用于根据输入的Landsat TM遥感影像数据,计算研究海域的归一化差分植被指数数据NDVI。具体而言:
所述归一化模块102通过Landsat TM遥感影像数据,计算研究海域的归一化差分植被指数数据(NDVI),计算公式如下:
NDVI=(B4-B3)/(B4+B3)
其中,NDVI表示归一化差分植被指数;B3表示Landsat>4表示Landsat>
所述线性压缩模块103用于将计算得到的NDVI数据进行线性压缩,得到NDVI′。具体而言:
所述线性压缩模块103将NDVI数据的数据范围由NDVI∈[-1,1]调整为NDVI∈[0,1],线性压缩公式为:
NDVIi′=(χi-χmin)/(χmax-χmin),i=1、2、3、…、n
其中,NDVI′表示线性压缩后的NDVI值;χi表示每个像元的NDVI值;χmin表示所有像元中的NDVI的最小值;χmax表示所有像元中的NDVI的最大值;n表示该NDVI数据的像元数量。
所述梯度值计算模块104用于计算上述得到的NDVI′数据的梯度值
所述梯度值计算模块104计算得到NDVI′数据的梯度grad(NDVI),计算公式如下:
dx(i,j)=NDVI′(i+1,j)-NDVI′(i,j)
dy(i,j)=NDVI′(i,j+1)-NDVI′(i,j)
其中,
所述控制点模块105用于在深圳的四大主体海域随机布设插值控制点作为一般控制点。其中:所述四大主体海域分别是珠江口海区、深圳湾海区、大鹏湾海区和大亚湾海区。具体而言:
本实施例中,依据深圳海域净初级生产力“西高东低,北强南弱”的分布特点,在梯度数据grad(NDVI)上划定处理区域,处理区域包括了深圳市海域的四大主体海域,分别是珠江口海区、深圳湾海区、大鹏湾海区和大亚湾海区。在划定好的四大主体海域,按照珠江口海区60、深圳湾海区45、大鹏湾海区100、大亚湾海区95的分配方式在各自区域上随机布设共计300个插值控制点来作为一般控制点。
所述控制点模块105还用于按照
本实施例中,为了保留小尺度数据上的纹理信息,在划定好的四大主体海域上选择
所述插值运算模块106用于对输入的MODIS/AQUA数据中的月平均海表温度数据、月平均海水漫衰减系数数据、月平均海水光合有效辐射数据、月平均海水叶绿素a浓度数据,分别进行反距离加权插值运算。具体包括:
对输入的MODIS/AQUA数据中的月平均海表温度数据,进行反距离加权插值运算,得到水体最大光合速率;
对输入的MODIS/AQUA数据中的月平均海水漫衰减系数数据,进行反距离加权插值运算,得到真光层深度;
对输入的MODIS/AQUA数据中的月平均海水光合有效辐射数据,进行反距离加权插值运算,得到表面日光合有效辐射;
对输入的MODIS/AQUA数据中的月平均海水叶绿素a浓度数据,进行反距离加权插值运算,得到垂直叶绿素a浓度。
所述插值运算模块106还用于对输入的研究区域反演时间当日的光照周期数据,进行反距离加权插值运算。具体而言:
输入研究区域反演时间当日的光照周期数据(Dirr),该数据为离散常数数值,单位是小时(h),该数据可由当地天文台查询获得。
所述反演结果模块107用于根据所述插值运算模块106进行反距离加权插值运算的结果,利用OVGPM模型计算得到深圳海域初级生产力遥感反演结果PPeu。具体而言:
利用OVGPM模型计算PPeu,公式如下:
其中,PPeu是海水表层到真光层的初级生产力,以每平方米产生的碳的毫克数计,单位是mg/m2;
本发明对VGPM模型的主要输入参数进行优化,将只适用于一类水体遥感反演的MODIS/AQUA数据产品扩展至以深圳近海海域为例的二级水体中去,首次尝试将MODIS/AQUA数据产品和小尺度的海区联系在一起,不仅扩展了MODIS/AQUA数据产品的适用范围,更有效地利用了现有技术手段和资源,为深圳市近海海域环境监测提供了廉价、高效、简易的技术服务支持。
虽然本发明参照当前的较佳实施方式进行了描述,但本领域的技术人员应能理解,上述较佳实施方式仅用来说明本发明,并非用来限定本发明的保护范围,任何在本发明的精神和原则范围之内,所做的任何修饰、等效替换、改进等,均应包含在本发明的权利保护范围之内。
机译: 船舶的实际海域性能评估方法,实际海域性能评估程序以及船舶的实际海域性能评估系统
机译: 用于将泥浆从邻近海域运输到邻近海域的位置的上升系统和配置Said系统的方法
机译: 遥感器,遥感系统和遥感方法