一种利用星载单光子激光雷达海面点云反演水深的方法.pdf

(19)国家知识产权局(12)发明专利申请(10)申请公布号CN115508854A(43)申请公布日2022.12.23(21)申请号202211094879.1(22)申请日2022.09.05(71)申请人山东科技大学地址266590山东省青岛市黄岛区前湾港路579号(72)发明人马跃王霖张婷朱凯郭忠磊杨坚郑慧莹宿殿鹏阳凡林(74)专利代理机构青岛智地领创专利代理有限公司37252专利代理师王鸣鹤(51)Int.Cl.G01S17/88(2006.01)G01S7/48(2006.01)G01C13/00(2006.01)权利要求书2页说明书5页附图3页(54)发明名称一种利用星载单光子激光雷达海面点云反演水深的方法(57)摘要本发明涉及利用星载单光子激光雷达海面点云反演水深的方法，属于激光遥感技术领域。对星载单光子激光雷达获取的原始光子点云进行海面光子信号提取，通过数据平滑算法拟合海面高程轮廓，提取出海面沿卫星飞行方向不同位置的波浪波长参数。基于线性波理论建立海面波浪波长、周期与水深的关系，利用深水区域波浪波长计算出波浪周期信息；利用波浪周期在深水、浅水区域不变的原理，利用单光子激光雷达计算的浅水波浪波长、波浪周期信息，计算浅水区域的水深结果。本发明通过星载单光子激光雷达获取的光子点云数据，快速准确地计算当地水深；利用所观测的海面波浪信息计算水深结果，对于填补使用激光雷达测量浑浊水体区域的水深具有重要的作用。CN115508854ACN115508854A权利要求书1/2页1.一种利用星载单光子激光雷达海面点云反演水深的方法，其特征在于，包括：S1.处理和提取星载单光子激光雷达飞越水面时所获取的水面光子点云信号，并拟合水面的波浪轮廓；S1.1使用点云密度提取算法，提取所述水面光子点云信号；每个空间网格大小在沿卫星飞行方向长度的间隔为Δl，在高程方向的间隔为Δh，水面光子点云信号沿卫星飞行方向的总长度为l，高程方向总长度为h，沿卫星飞行方向的空间网格数量nl和在高程方向的空间网格数量nh分别表示为：nl＝ceil(l/Δl)，nh＝ceil(h/Δh)，ceil表示向上取整函数；将研究区域内的水面光子点云信号进行空间网格划分，在沿卫星飞行方向的第i个且沿高度方向第j个网格的空间坐标范围F(i,j)满足公式(1)；其中，lstart和lend分别是研究区域内沿卫星飞行方向的起始坐标和终止坐标，单位为米，hmin和hmax分别是研究区域内光子点云高程方向的记录范围，单位为米；统计每个空间网格坐标范围F(i,j)内的光子点云数量，即单位空间内的光子点云密度，记为N(i,j)，计算所有空间网格点云密度的均值μ和标准差σ，设定有效信号阈值TH＝μ+ea×σ，ea是比例系数，在研究区域内，所有空间网格点云密度N(i,j)大于鉴别阈值TH的空间范围内点云都分类为水面光子点云信号，其余为噪声点云；S1.2使用水面光子点云信号拟合得到水面高程轮廓H(x)；选取沿卫星飞行方向总长度为l的水面光子点云信号，沿飞行方向每隔距离Δx，计算每个距离间隔Δx内所有水面信号光子的平均高程Hphotonavg(xn)，即：其中，xn表示第n个距离间隔Δx，M表示当前第n个距离间隔Δx内包含的水面信号光子个数，h(m)表示该距离间隔内第m个水面光子点云的高程；对计算得到的平均高程序列Hphotonavg(x)再进行Savitzky‑Golay中值滤波，得到水面高程轮廓序列H(x)，其中，中值滤波窗宽同样使用Δx，P为距离间隔倍数，取10～15之间；S2.提取远海深水区的海面波浪波长，利用深水区的水面波浪波长计算波浪周期；对于海面重力波，依据线性波理论，波浪波长L和周期T在水深D时，频散关系满足公式(2)：其中，g是重力加速度常数，π是圆周率常数，tanh是双曲正切函数；在深水区域，即当D>L/2时，有tanh(2πD/L)≈1，(2)式可以简化为：给出深水区波浪波长Ldeep，即可以求取波浪周期T；对深水处的水面高程轮廓序列H(x1)进行快速傅里叶变换，得到自变量为波数k的幅度谱分布Y(k)＝FFT(H(x1))，选取幅度谱Y(k)峰值位置处对应波数k作为深水区域波数kdeep的大小；2CN115508854A权利要求书2/2页根据波数与波长之间的关系，2π＝k×L计算出深水处的波长Ldeep，则波浪周期T可以用公式(3)计算得到：S3.提取近岸浅水区的海面波浪波长，利用浅水区的水面波长反演对应位置的水深；对浅水处的水面高程轮廓序列H(x2)进行短时傅里叶变换STFT，STFT窗宽选择为2×Ldeep，沿飞行X方向移动步长选择为0.2×Ldeep，每个步长对应沿卫星飞行方向的位置为xp，通过提取出每个窗宽范围内幅度谱峰值位置的波数作为当前位置xp的波数ks