1.一种基于星载激光雷达点云的水深反演方法,其特征在于,首先,获取ICESat‑2点云数据并进行预处理;
然后,基于卷积神经网络模型联合点云的邻域空间特征和光谱特征对预处理后的ICESat‑2点云数据提取水下地形光子;
结合点云的邻域空间特征利用空间密度差异及邻域位置关系对预处理后的点云数据进一步去噪,将预处理后的点云数据的任意光子均生成邻域空间特征张量,作为卷积神经网络模型的第一支输入,维度大小为128*128*3;
邻域空间特征张量的生成方法包括:
S1、以任意光子为中心建立矩形窗口,中心光子位于其周围四个网格的十字交织节点上,中心光子到矩形窗口的四个边界的距离相等,对每个矩形窗口在沿轨距离方向和高程方向设置网格分辨率划分网格;
S2、计算网格内的光子数,输出表达为第一通道特征张量,若网格内没有光子,则对应网格内的光子数为0;
S3、计算网格内光子高程相对于中心光子的网格协方差矩阵 ,求得 的特征值和 ,然后计算用于表征紧凑程度的高程全方差 ,高程全方差 填充对应网格作为第二通道特征张量,若网格内没有光子,则对应网格输出为0;
, ,式中, 为网格
内的总光子数, 为网格内第 个光子高程值, 为中心光子的高程值;
S4、计算网格内光子沿轨距离相对于中心光子的网格协方差矩阵 ,求得 的特征值 和 ,然后计算用于表征紧凑程度的沿轨距离全方差 ,沿轨距离全方差 填充对应网格作为第三通道特征张量,若网格内没有光子,则对应网格输出为0;
, ,式中, 为网格
内第 个光子对应的沿轨距离值, 为中心光子的沿轨距离值;
以任一光子匹配经预处理的WorldView‑2 2m高分辨率卫星影像的像素格为中央单元,获取卫星影像上邻域内9*9的像素单元遥感反射率;
将八个单波段的像元值作为八个通道光谱特征张量输出,每个通道光谱特征张量的尺寸为9*9,中央像元值为任一中心光子匹配卫星影像获得的单波段遥感反射率;
八通道的光谱特征张量作为卷积神经网络模型的第二支输入,输入大小维度为9*9*8;
卷积神经网络模型的训练及构建方法包括以下步骤:
A1、对预处理后的点云数据的同一光子分别构建邻域空间特征张量和光谱特征张量,进而对预处理后的所有点云光子生成数据集,将数据集按照8:2的比例划分为训练集和验证集;
A2、设置卷积神经网络模型结构:卷积神经网络模型为双分支网络,接受两个分支输入,卷积神经网络模型的结构分为输入特征提取层、特征融合层和预测分类层;
输入特征提取层的一个分支以稠密连接网络为主体架构,接受维度大小为128*128*3的邻域空间特征张量,另一个分支嵌入ECA注意力机制模块,接受维度大小为9*9*8光谱特征张量;
A3、导入训练集对卷积神经网络模型进行训练,将训练好的卷积神经网络模型应用到预处理后的点云数据,联合点云的邻域空间特征及光谱特征对预处理后的点云数据提取噪声光子、陆地光子、水面光子和水下地形光子;
最后,将水下地形光子代替原位测深点训练和构建星载激光雷达主被动融合水深反演模型,利用星载激光雷达主被动融合水深反演模型进行水深反演。
2.根据权利要求1所述的基于星载激光雷达点云的水深反演方法,其特征在于,点云数据的获取及预处理方法为:首先,下载h5文件格式的ICESat‑2原始点云数据集,将同一研究区域下的ATL03数据与对应的ATL08数据通过RGT号进行关联,获取有效原始ATL03数据;
然后,将h5文件格式的有效原始ATL03数据转化为mat文件格式的多轨点云原始数据;
随后在沿轨距离方向上每200米分段划分,每段区间根据高程直方图确定最大峰 ,即最大光子数,以及上下临近两个次级峰 和 ,次级峰 和 分别对应高程为和 ,且 > ,根据 向上延伸 ,根据 向下延伸 ,确定标高范围为大于且小于 ,标高范围以外的光子为噪声光子,标高范围内的光子置信度
小于1的光子为噪声光子,获得预处理点云数据。
3.根据权利要求2所述的基于星载激光雷达点云的水深反演方法,其特征在于,获取研究区域的多张WorldView‑2 2m高分辨率单景卫星影像,分别对多张WorldView‑2 2m高分辨率单景卫星影像进行地理校正、辐射校正、大气校正和耀斑校正,校正完成后对多张WorldView‑2 2m高分辨率单景卫星影像进行无缝拼接、去云和中值滤波操作,将预处理后的WorldView‑2 2m高分辨率卫星影像和ICESat‑2点云数据进行配准,获取每个波段像素单元遥感反射率。
4.根据权利要求3所述的基于星载激光雷达点云的水深反演方法,其特征在于,对提取的水下地形光子进行折射校正:,式中,
为折射校正后的水下地形光子高程, 为水面高度, 为
提取的水下地形光子的高程, 为水的折射率;
将折射校正后的水下地形光子高程与预处理后的WorldView‑2 2m高分辨率卫星影像的像素单元进行匹配,对匹配后的多个折射校正后的水下地形光子高程值取平均获得平均高程;
将匹配后的像素单元遥感反射率与对应的平均高程代入波段比值模型进行训练,构建星载激光雷达主被动融合水深反演模型,所述波段比值 模型表 示为:,式中, 表示水深值, 表示WorldView‑2 2m高分
辨率卫星影像中绿光波段的遥感反射率, 表示WorldView‑2 2m高分辨率卫星影像中蓝光波段的遥感反射率, 和 表示波段比值模型的系数。