基于固有光学参数的线性水深反演方法技术

本发明专利技术提供一种基于固有光学参数的线性水深反演方法，反演方法包括：对遥感影像进行预处理得到预处理影像；选取所述预处理影像中预定数量的水深控制点；分别提取各个水深控制点的蓝、绿波段反射率值；根据各个水深控制点的蓝、绿波段反射率值计算得到各个所述水深控制点的蓝、绿波段的固有光学参数；根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型；利用所述线性水深反演模型对水深进行反演。该方法在缺乏现场信息与大量水深实测深度的条件下能够获得准确的水深，适用范围广、运算量小、成本低。

【技术实现步骤摘要】
基于固有光学参数的线性水深反演方法
本专利技术涉及海洋探测领域，特别地涉及一种基于固有光学参数的线性水深反演方法。
技术介绍
水深是了解海洋环境的一个关键参数，也是水文测量和海图制作的重要组成部分。对海岸研究和海洋工程建设非常重要，也为海上运输和航行提供重要信息，同时有利于珊瑚礁的管理监测与生态系统的保护。传统的船载水深测量方法存在不足之处，对于船只无法到达的复杂区域不能适用，而且测量过程耗时耗力。随着一系列地球观测卫星的升空，遥感平台的多样化和可用遥感数据的富足极大的推进了光学水深遥感的发展。机载激光雷达测深是一种主动式遥感测深技术，可以测得某些船舶无法到达的区域如岛礁、暗礁的有效水深，但是这种方法成本高，测量范围较小，很难广泛应用。被动光学遥感水深探测因其覆盖范围广、更新时间快、成本低的优势，逐渐成为获取水深信息的一种重要方法。被动光学遥感技术包括多光谱遥感及高光谱遥感。基于水体辐射传输理论，水深与水体的固有光学特性(如水柱的扩散衰减系数、吸收系数及后向散射系数)等相关，据此发展了大量的水深被动光学反演模型。根据模型的参数化程度，基本可以分为理论解析模型与半理论半经验模型两类。理论解析模型在构建过程中所需的水体光学参数多，计算复杂。半理论半经验模型是理论模型与经验算法的结合，其中对数线性模型和对数转换比值模型应用相对广泛。
技术实现思路
有鉴于此，本专利技术提供一种基于固有光学参数的线性水深反演方法，以解决现有技术中的水深探测方法适用范围小、成本高、计算复杂的问题。为达上述目的，本专利技术提供的一种基于固有光学参数的线性水深反演方法的技术方案是：一种基于固有光学参数的线性水深反演方法，所述反演方法包括：对遥感影像进行预处理得到预处理影像；选取所述预处理影像中预定数量的水深控制点；分别提取各个水深控制点的蓝、绿波段反射率值；根据各个水深控制点的蓝、绿波段反射率值计算得到各个所述水深控制点的蓝、绿波段的固有光学参数；根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型；利用所述线性水深反演模型对水深进行反演。优选地，所述固有光学参数与吸收系数和后向散射系数相关。优选地，所述固有光学参数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个控制点，并保证水深控本文档来自技高网...

【技术保护点】
1.一种基于固有光学参数的线性水深反演方法，其特征在于，所述反演方法包括：/n对遥感影像进行预处理得到预处理影像；/n选取所述预处理影像中预定数量的水深控制点；/n分别提取各个水深控制点的蓝、绿波段反射率值；/n根据各个水深控制点的蓝、绿波段反射率值计算得到各个所述水深控制点的蓝、绿波段的固有光学参数；/n根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型；/n利用所述线性水深反演模型对水深进行反演。/n

【技术特征摘要】
1.一种基于固有光学参数的线性水深反演方法，其特征在于，所述反演方法包括：
对遥感影像进行预处理得到预处理影像；
选取所述预处理影像中预定数量的水深控制点；
分别提取各个水深控制点的蓝、绿波段反射率值；
根据各个水深控制点的蓝、绿波段反射率值计算得到各个所述水深控制点的蓝、绿波段的固有光学参数；
根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波段的固有光学参数构建线性水深反演模型；
利用所述线性水深反演模型对水深进行反演。


2.根据权利要求1所述的基于固有光学参数的线性水深反演方法，其特征在于，所述固有光学参数与吸收系数和后向散射系数相关。


3.根据权利要求2所述的基于固有光学参数的线性水深反演方法，其特征在于，所述固有光学参数u(λ)通过如下公式得到：



其中，p0和p1为随水体光学特性变化的模型常量；
rrs(λ)为将相应波段的反射率Rrs(λ)转换为恰好在水面之下的遥感反射率。


4.根据权利要求3所述的基于固有光学参数的线性水深反演方法，其特征在于，p0＝0.0895；
p1＝0.1247。


5.根据权利要求3所述的基于固有光学参数的线性水深反演方法，其特征在于，rrs(λ)通过如下公式得到：





6.根据权利要求3至5之一所述的基于固有光学参数的线性水深反演方法，其特征在于，根据各个水深控制点的实际深度值以及各个水深控制点的蓝、绿波...

【专利技术属性】
技术研发人员：张靖宇，张雪纯，马毅，
申请(专利权)人：自然资源部第一海洋研究所，
类型：发明
国别省市：山东;37