[发明专利]一种绘制高精度水底地形与淤泥厚度双界面分布图的方法

权利要求书

1.一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于，包括以下步骤：

步骤S1：将单频测深仪放入水中发出多波次声波信号；

步骤S2：将多频测深仪放入水中发出多波次声波信号；

步骤S3：采集单频测深仪所发出的声波信号，处理之后得点位i水深数据记为H_i1，采集多频测深仪所发出的声波信号，处理之后得点位i水深数据记为H_i2,处理之后得到点位i淤泥厚度数据记为N_i；

步骤S4：根据H_i1和H_i2数据，校正H_i1数据；

步骤S5：重复步骤S1～S4，根据各个空间点的H_i1计算得水深网格文件；

步骤S6：重复步骤上述S1～S4得到淤泥下界面深度数据H_i3＝H_i+N_i，根据各个空间点的H_i1计算得到淤泥下界面网格文件；

步骤S7：使用软件将水底网格文件和淤泥厚度网格文件拟合，得到水深泥厚双界面分布图。

2.根据权利要求1所述的一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于：

步骤S4包括以下步骤：

a、计算收敛因子K_d，根据公式计算所有点位H1和H2的差值平均值记为A1(d表示迭代次数)；

b、根据收敛因子k_d，若满足收敛因子k_d条件，无需校正水深数据H_i1，若不满足收敛因子k_d条件，校正得新的水深数据H_i1；

c、将新的水深数据H_i重新迭代计算收敛因子k_d，直到收敛因子k_d满足要求，输出水深数据H_i1。

3.根据权利要求2所述的一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于：在步骤S4中的a步骤中：设置收敛因子其中σ_d为H_i1与H_i2的标准差，η为±1cm+0.4％×h，h为该点位的水深值。

4.根据权利要求3所述的一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于：若k_d≤1,则无需修正；若k_d1，H_i1与H_i2差值大于平均值A1，将H_i1的值修正±0.5σi，得新的H_i1。

5.根据权利要求4所述的一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于：在步骤S4中的c步骤中：将新的H_i1，重新带入到a步骤，重复循环步骤，直到k_d≤1为止。

6.根据权利要求1所述的一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于：在步骤S6中，H1和H3进行数值插值计算，设置相同的数据范围和网格大小，从而获得网格尺寸。

7.根据权利要求6所述的一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于：在步骤S6中，制定网格尺寸之后，分别导出中心坐标点位一致的记录水深的grd文件与记录淤泥厚度的grd文件。

8.根据权利要求7所述的一种绘制高精度水底地形与淤泥厚度双界面分布图的方法，其特征在于：在步骤S7中，使用matlab分别提取记录水深的grd文件中具体的水深数据和记录淤泥厚度的grd文件淤泥厚度数据，再使用matlab将同一坐标点的水深数据和淤泥厚度数据合并，之后重复数据提取和合并过程，直到所有坐标点的数据处理完成，将所有坐标点的数据导出绘图软件中，绘制出水深泥厚双界面分布图。