基于固有光学参数的线性水深反演方法与流程

本发明涉及海洋探测领域，特别地涉及一种基于固有光学参数的线性水深反演方法。

背景技术：

水深是了解海洋环境的一个关键参数，也是水文测量和海图制作的重要组成部分。对海岸研究和海洋工程建设非常重要，也为海上运输和航行提供重要信息，同时有利于珊瑚礁的管理监测与生态系统的保护。传统的船载水深测量方法存在不足之处，对于船只无法到达的复杂区域不能适用，而且测量过程耗时耗力。

随着一系列地球观测卫星的升空，遥感平台的多样化和可用遥感数据的富足极大的推进了光学水深遥感的发展。机载激光雷达测深是一种主动式遥感测深技术，可以测得某些船舶无法到达的区域如岛礁、暗礁的有效水深，但是这种方法成本高，测量范围较小，很难广泛应用。被动光学遥感水深探测因其覆盖范围广、更新时间快、成本低的优势，逐渐成为获取水深信息的一种重要方法。

被动光学遥感技术包括多光谱遥感及高光谱遥感。基于水体辐射传输理论，水深与水体的固有光学特性(如水柱的扩散衰减系数、吸收系数及后向散射系数)等相关，据此发展了大量的水深被动光学反演模型。根据模型的参数化程度，基本可以分为理论解析模型与半理论半经验模型两类。理论解析模型在构建过程中所需的水体光学参数多，计算复杂。半理论半经验模型是理论模型与经验算法的结合，其中对数线性模型和对数转换比值模型应用相对广泛。

技术实现要素：

有鉴于此，本发明提供一种基于固有光学参数的线性水深反演方法，以解决现有技术中的水深探测方法适用范围小、成本高、计算复杂的问题。

为达上述目的，本发明提供的一种基于固有光学参数的线性水深反演方法的技术方案是：

一种基于固有光学参数的线性水深反演方法，所述反演方法包括：

对遥感影像进行预处理得到预处理影像；

选取所述预处理影像中预定数量的水深控制点；

分别提取各个水深控制点的蓝、绿波段反射率值；

根据各个水深控制点的蓝、绿波段反射率值计算得到各个所述水深控制点的蓝、绿波段的固有光学参数；

根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型；

利用所述线性水深反演模型对水深进行反演。

优选地，所述固有光学参数与吸收系数和后向散射系数相关。

优选地，所述固有光学参数u(λ)通过如下公式得到：

其中，p0和p1为随水体光学特性变化的模型常量；

rrs(λ)为将相应波段的反射率rrs(λ)转换为恰好在水面之下的遥感反射率。

优选地，p0＝0.0895；

p1＝0.1247。

优选地，rrs(λ)通过如下公式得到：

优选地，根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型的方法包括：

构建水深值z与蓝、绿波段的固有光学参数的关系式：

其中，a和b为回归系数；

u(λi)为蓝波段i的固有光学参数；

u(λj)为绿波段j的固有光学参数；

利用各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数的比值进行回归得到所述关系式中的回归系数a和b，从而构建出所述线性水深反演模型。

优选地，所述对遥感影像进行预处理得到预处理影像包括：

对所述遥感影像进行辐射定标；

对辐射定标后的遥感影像进行大气校正处理；

对大气校正处理后的遥感影像中的陆地区域进行掩膜，得到所述预处理影像。

优选地，所述遥感影像为worldview-2影像，对所述遥感影像进行辐射定标的转换公式为：

其中，l(λi)为第i波段的辐亮度值；

abscalfactori为第i波段的绝对定标系数；

dni为第i波段图像像元灰度值；

δλi为第i波段的等效波段宽度。

优选地，采用flaash模块进行大气校正。

优选地，所述遥感影像为多光谱遥感影像。

本发明提供的基于固有光学参数的线性水深反演方法的有益效果为：

在缺乏现场信息与大量水深实测深度的条件下能够获得准确的水深，适用范围广、运算量小、成本低。

附图说明

附图用于更好地理解本发明，不构成对本发明的不当限定。其中：

图1是本发明具体实施例提供的基于固有光学参数的线性水深反演方法的流程图；

图2是本发明具体实施例提供的ioplm散点拟合图；

图3是本发明具体实施例提供的水深反演结果图。

具体实施方式

下面结合实施例对本发明作进一步说明。

针对现有技术中的水深探测方法适用范围小、成本高、计算复杂的问题，本实施例提供了一种基于固有光学参数的线性水深反演方法，尤其是对于海面的水深反演方法，如图1所示，其包括如下步骤：

s100、对遥感影像进行预处理得到预处理影像；

s200、选取所述预处理影像中预定数量的水深控制点；

s300、分别提取各个水深控制点的蓝、绿波段反射率值；

s400、根据各个水深控制点的蓝、绿波段反射率值计算得到各个所述水深控制点的蓝、绿波段的固有光学参数；

s500、根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型；

s600、利用所述线性水深反演模型对水深进行反演。

本发明提供的基于固有光学参数的线性水深反演方法，能够获得如下有益效果：

在缺乏现场信息与大量水深实测深度的条件下能够获得准确的水深，适用范围广、运算量小、成本低，尤其适用于岛礁、暗礁等船舶难以到达区域的水深反演。

在步骤s100中，遥感影像优选为多光谱遥感影像，对遥感影像进行预处理使得后续的运算数据更加精确，对遥感影像进行预处理的过程优选包括如下步骤：

s110、对所述遥感影像进行辐射定标；

s120、对辐射定标后的遥感影像进行大气校正处理；

s130、对大气校正处理后的遥感影像中的陆地区域进行掩膜，得到所述预处理影像。

其中，步骤s110中，对遥感影像进行辐射定标即将影像上无量纲的灰度值转换为具有实际意义的辐亮度值的过程，在一个具体的实施例中，所述遥感影像为worldview-2影像，对所述遥感影像进行辐射定标的转换公式为：

其中，l(λi)为第i波段的辐亮度值；

abscalfactori为第i波段的绝对定标系数；

dni为第i波段图像像元灰度值；

δλi为第i波段的等效波段宽度。

当然，可以理解的是，遥感影像也可以是其他多光谱影像，转换公式也相应有所变化。

大气的衰减作用对不同波长的光是有选择性的，因而大气对不同波段影像的影响是不同的，针对这一问题，采用步骤s120对辐射定标后的遥感影像进行大气校正处理，从而去除大气散射、气溶胶等的影响，尽可能还原地表真实反射率。在一个优选的实施例中，采用flaash模块对影像进行大气校正。最后，采用步骤s130对陆地区域进行掩膜，得到处理后的最终影像。

步骤s200中，每米间隔约选取4至6个控制点，并保证水深控制点在空间上大致分布均匀。

步骤s400中，固有光学参数与吸收系数和后向散射系数相关。例如，在一个优选的实施例中，所述固有光学参数u(λ)通过如下公式得到：

其中，p0和p1为随水体光学特性变化的模型常量，例如，在一个具体的实施例中，p0＝0.0895，p1＝0.1247；

rrs(λ)为将相应波段的反射率rrs(λ)转换为恰好在水面之下的遥感反射率。

进一步优选地，rrs(λ)通过如下公式得到：

进一步地，根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型的方法包括：

构建水深值z与蓝、绿波段的固有光学参数的关系式：

其中，a和b为回归系数；

u(λi)为蓝波段i的固有光学参数；

u(λj)为绿波段j的固有光学参数；

利用各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数的比值进行回归得到所述关系式中的回归系数a和b，从而构建出所述线性水深反演模型(ioplm，inherentopticalparameterslinearmodel)，回归方法例如是采用最小二程法进行回归。

下面选取中国西沙群岛某珊瑚岛礁进行水深反演研究，以worldview-2多光谱遥感影像为例，具体说明本申请提供的线性水深反演模型(ioplm，inherentopticalparameterslinearmodel)的构建方法。

1.遥感影像预处理

首先对遥感影像进行辐射定标，即将影像上无量纲的灰度值转换为具有实际意义的辐亮度值的过程，对于worldview-2影像，转换公式如下：

其中，l(λi)为第i波段的辐亮度值，单位为w/(m²·sr·μm)；

abscalfactori为第i波段的绝对定标系数；

dni为第i波段图像像元灰度值；

δλi为第i波段的等效波段宽度。

接着采用flaash模块对影像进行大气校正处理，大气的衰减作用对不同波长的光是有选择性的，因而大气对不同波段影像的影响是不同的，大气校正可以去除大气散射、气溶胶等的影响，尽可能还原地表真实反射率。影像经过上述处理步骤后，再对陆地区域进行掩膜，得到处理后的最终影像。

2.基于固有光学参数构建线性水深反演模型

选取一定数量的水深控制点，提取经过预处理步骤影像的蓝、绿波段反射率值，并计算出各自的固有光学参数u(λ)，利用实测水深值与蓝、绿波段固有光学参数的比值进行回归，求得回归参数a和b。本例中的参数为a＝31.734；b＝-30.729，水深值与固有光学参数比值的拟合散点图如图2所示。

3.浅海水深反演结果

利用上节构建的线性水深反演模型(ioplm)对整个研究区进行水深反演，得到浅海(≤20m)水深反演结果如图3所示。

4.水深反演精度评价

选取一定数量的检查点对模型进行精度评价，分别计算模型整体以及分水深段的平均绝对误差(meanabsoluteerror)、平均相对误差(meanrelativeerror)及均方根误差(rootmeansquarederror)。其中，平均绝对误差(meanabsoluteerror)的计算公式如下：

平均相对误差(meanrelativeerror)的计算公式如下：

平均绝对误差(meanabsoluteerror)的计算公式如下：

具体结果见表1。

表1水深反演精度评价

由表1可知，本申请构建的线性水深反演模型(ioplm)整体反演精度较高，平均绝对误差为1.17m，平均相对误差为17.5％，均方根误差为1.49m。

本领域内的技术人员应明白，本公开的实施例可提供为方法、系统、或计算机程序产品。因此，本公开可采用硬件实施例、软件实施例、或结合软件和硬件方面的实施例的形式。而且，本公开可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。

本公开是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，所述的存储介质可为磁碟、光盘、只读存储记忆体(read-onlymemory，rom)或随机存储记忆体(randomaccessmemory，ram)等。

上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本公开保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。