一种基于ODP的海底沙波特征自动识别方法技术

本发明专利技术公开了一种基于ODP的海底沙波特征自动识别方法。该方法根据水深曲面的最优剖面方向求导并判定极值，以提取海底沙波脊线和谷线特征点，用于海底沙波特征线的自动识别与提取，并通过1）构建数字水深矩阵，2）构建最优方向矩阵，3）提取海底沙波特征等三大步骤实现了海底沙波的自动识别。经实测多波束水深数据测试验证，采用该方法自动识别的海底沙波特征线与人工识别特征线的相关系数可达80%以上，与传统方法相比自动识别准确率平均提升约30%，且勿需设置阈值，大幅提升了工作效率。该技术方法可推广应用到其他类型海底地形以及陆地地形的自动识别，在海洋测绘、海洋地理信息系统、计算机图形学和海底科学研究等方面具有重要的实际应用价值。

【技术实现步骤摘要】
-种基于ODP的海底沙波特征自动识别方法
本专利技术涉及海底地形地貌制图、海洋测绘、海洋地理信息系统、计算机图形学和海 底科学等
。
技术介绍
海底沙波是发育在近海陆架上的一种极常见的海洋地貌，沙波的迁移和运动可表 征海洋沉积环境的特征，也可能影响甚至妨碍人类活动，如可能掏空或掩埋铺设的海底管 线，淤浅航道港口，导致风电场的构筑物基础失稳等，可造成严重的海底灾害。因此，海底沙 波的研究对了解浅海沉积物输运，合理设计规划海洋工程提供有力的信息和支持。沙波具 有似波的形态，波峰和波谷以及由其在平面方向上连结而成的特征线脊线和谷线是描述其 形态特征的最基本参量，也是研究沙波迁移运动规律的基础。 陆地山谷地形自动识别方面有多重方法，如：基于地表流水物理模拟的水文分析 法，基于地表几何形态的坡向分析法，以及综合运用分析法等。陆地山谷地形发育多受控于 构造运动，而海底沙波发育的主因是潮流和波浪，二者在几何尺度不可类比。因此难以简单 照搬这些方法进行海底沙波的识别。试验表明，应用上述方法提取海底沙波形态特征时对 于地形变化平缓和局部起伏差异大的区域，特征信息提取效果不佳。此外，分类阈值是这些 方法识别陆地地形特征的基础，但是目前阈值的选取还存在一些局限性，需要反复尝试且 易受主观因素影响，导致地形特征识别的自动化程度和工作效率降低。 在海底沙波研究方面，已有一些研究成果，如：一种基于MBES的海底沙波地貌运 动探测方法（ZL. 201310317429. 9) 公开了一种基于多波束测深数据来识别海底沙波运动 的方法，一种海底大型复杂沙波地貌的精确探测方法（ZL. 201310317430. 1) 公开了一种 基于剖面的海底沙波FFT分解与拟合方法。但这些方法均没提及如何自动识别海底沙波的 特征线，也就是海底沙波的脊线和谷线。 因此，目前的陆地地形特征识别方法难以简单移植到海底沙波特征自动识别，对 于海底沙波的形态自动识别尚缺乏有针对性的方法。
技术实现思路
针对现有地形识别方法强依赖阈值选取的局限性和识别率不高的问题，本专利技术公 开了一种基于0DP的海底沙波特征自动识别方法，具体而言是一种基于格网的海底沙波特 征线的最优方向剖面自动识别方法（〇DP，Optimally_Directional Profiling method):首先 基于水深曲面的最优剖面方向求导并判定极值，以提取海底沙波脊线和谷线特征点。 -种基于0DP的海底沙波特征自动识别方法，包括下列步骤： 1)测绘数字水深矩阵 1. 1)构建水深矩阵 1. 1. 1)若无多波束离散水深数据，采用全覆盖测量方式获取多波束水深数据，经 数据编辑和改正处理后获取多波束离散数据集合MBES，进入步骤1. 1. 2); 1. 1. 2)若有多波束离散水深数据集MBES，采用距离反比加权方法构建水深矩阵本文档来自技高网...

【技术保护点】
一种基于ODP的海底沙波特征自动识别方法，其特征在于，包括下列步骤：1)测绘数字水深矩阵1.1)构建水深矩阵1.1.1)若无多波束离散水深数据，采用全覆盖测量方式获取多波束水深数据，经数据编辑和改正处理后获取多波束离散数据集合MBES，进入步骤1.1.2)；1.1.2)若有多波束离散水深数据集合MBES，采用距离反比加权方法构建水深矩阵其中G(i，j)为第i行第j列水深矩阵点，nx为矩阵总行数，ny为矩阵总列数，i、j、nx、ny均为自然数；1.2)构建水深矩阵的数据索引以水深矩阵的行i为顺序，依次建立每行的数据索引Indexi＝{id[m]}m＝1，n，id[m]中存储着水深矩阵Depth第i行的第m个有水深数值的格网所在列，n表示水深矩阵Depth第i行有水深数值的格网的个数，其中m和n均为自然数；2)构建最优剖面方向矩阵2.1)建立最优剖面方向矩阵令a(i，j)＝0；2.2)建立步骤1.1.2)所述的G(i，j)的四邻域窗口矩阵A=G(i,j-1)G(i,j+1)G(i-1,j)G(i+1,j);]]>2.3)计算G(i，j)的最优剖面方向值其中dx＝[G(i，j‑1)‑G(i，j+1)]，dy＝[G(i+1，j)‑G(i‑1，j)]，根据步骤2.3.1)至2.3.3)三种情况计算θ值：2.3.1)当dy＝0时：a)若dx＞0，则θ＝90°；b)若dx＜0，则θ＝270°；2.3.2)当dy＞0时：a)若dx≥0，则θ＝θ；b)若dx＜0，则θ＝360°‑abs(θ)；2.3.3)当dy＜0时：a)若dx≥0，则θ＝180°‑θ；b)若dx＜0，则θ＝180°+abs(θ)；2.3.4)将步骤2.3.1)至2.3.3)中计算的最优剖面方向值θ赋予a(i，j)＝θ；2.4)根据步骤1.2)所述的数据索引Indexi，提取对应的G(i，j)，循环步骤2.2)至2.3)，建立完整的最优剖面方向矩阵Aspect；3)提取海底沙波特征线3.1)确定求导方向提取G(i，j)对应的最优剖面方向a(i，j)，a(i，j)如步骤2.3)所述，根据步骤3.1.1)至3.1.4)确定对应G(i，j)的求导方向：3.1.1)若a(i，j)在区间I或区间V，则沿y＝0方向求导；3.1.2)若a(i，j)在区间II或区间VI，则沿y＝x方向求导；3.1.3)若a(i，j)在区间III或区间VII，则沿x＝0方向求导；3.1.4)若a(i，j)在区间IV或区间VIII，则沿y＝‑x方向求导；其中，区间I＝(0°，22.5°]∪(337.5°，360°]，区间II＝(22.5°，67.5°]，区间III＝(67.5°，112.5°]，区间IV＝(112.5°，157.5°]，区间V＝(157.5°，202.5°]，区间VI＝(202.5°，247.5°]，区间VII＝(247.5°，292.5°]，区间VIII＝(292.5°，337.5°]；3.2)判定极值3.2.1)构建G(i，j)的八邻域窗口矩阵B，：B=G(i-1,j-1)G(i-1,j)G(i-1,j+1)G(i,j-1)G(i,j)G(i,j+1)G(i+1,j-1)G(i+1,j)G(i+1,j+1);]]>3.2.2)建立极小值矩阵令c(i，j)＝0，Crest也称为沙波脊线矩阵；3.2.3)建立极大值矩阵令t(i，j)＝0，Trough也称为沙波谷线矩阵；3.2.4)根据下述步骤a)至d)来判断沙波特征：a)当G(i，j)的沿y＝0方向求导时：i)若[G(i，j)‑G(i‑1，j)]＜0且[G(i，j)‑G(i+1，j)]＜0，则G(i，j)为极小值点，c(i，j)＝1；ii)b)若[G(i，j)‑G(i‑1，j)]＞0且[G(i，j)‑G(i+1，j)]＞0，则G(i，j)为极大值点，t(i，j)＝1；b)当G(i，j)的沿y＝x方向求导时：i)若[G(i，j)‑G(i‑1，j+1)]＜0且[G(i，j)‑G(i+1，j‑1)]＜0，则G(i，j)为极小值点，c(i，j)＝1；ii)b)若[G(i，j)‑G(i‑1，j+1)]＞0且[G(i，j)‑G(i+1，j‑1)]＞0，则G(i，j)为极大值点，t(i，j)＝1；c)当G(i，j)的沿x＝0方向求导时：i)若[G(i，j)‑G(i，j‑1)]＜0且[G(i，j)‑G(i，j+1)]＜0，则G(i，j)为极小值点，c(i，j)＝1；ii)b)若[G(i，j)‑G(i，j‑1)]＞0且[G(i，j)‑G(i，j+1)]＞0，则G(i，j)为极大值点，t(i，j)＝1；d)当G(i，j)的沿y＝‑x方向求导时：i)若[G(i，j)‑G(i‑1，j‑1)]＜0且[G(i，j)‑G(i+1，j+1)]＜0，则G(i，j)为极小值点，c(...

【技术特征摘要】
1. 一种基于ODP的海底沙波特征自动识别方法，其特征在于，包括下列步骤： 1) 测绘数字水深矩阵 1. 1)构建水深矩阵 1. 1. 1)若无多波束离散水深数据，采用全覆盖测量方式获取多波束水深数据，经数据 编辑和改正处理后获取多波束离散数据集合MBES，进入步骤1. 1. 2); 1. 1. 2)若有多波束离散水深数据集合MBES，采用距离反比加权方法构建水深矩阵'其中G(i, j)为第i行第j列水深矩阵点，nx为矩阵总行数，ny为 矩阵总列数，i、j、nx、ny均为自然数； 1. 2)构建水深矩阵的数据索引 以水深矩阵的行i为顺序，依次建立每行的数据索引IndeXi = id[m]中存储着水深矩阵D印th第i行的第m个有水深数值的格网所在列， η表示水深矩阵D印th第i行有水深数值的格网的个数，其中m和η均为自然数；2) 构建最优剖面方向矩阵 2. 1)建立最优剖面方向矩阵 2. 2)建立步骤1. 1.2)所述的G(i，j)的四邻域窗口矩阵2. 3)计算G(i，j)的最优剖面方向值其中dx = [G(i，j-l)-G(i，j+1)]， dy = [G(i+1，j)-G(i-l，j)]，根据步骤 2. 3. 1)至 2. 3. 3)三种情况计算 Θ 值：2. 3. 4)将步骤2. 3. 1)至2. 3. 3)中计算的最优剖面方向值Θ赋予a(i，j) = Θ ; 2.4)根据步骤1.2)所述的数据索引111(16^，提取对应的6(1，」），循环步骤2.2)至 2. 3)，建立完整的最优剖面方向矩阵Aspect ; 3) 提取海底沙波特征线 3. 1)确定求导方向 提取G(i，j)对应的最优剖面方向a(i，j)，a(i，j)如步骤2. 3)所述，根据步骤3. 1. 1) 至3. 1.4)确定对应G(i，j)的求导方向： 3. 1. 1)若a(i，j)在区间I或区间V，则沿y = 0方向求导； 3. 1. 2)若a(i，j)在区间II或区间VI，则沿y = X方向求导； 3. 1.3)若a(i，j)在区间III或区间VII，则沿x = 0方向求导； 3. 1.4)若a(i，j)在区间IV或区间VIII，则沿y = -x方向求导； 其中，区间1=(0°,22.5。]1^337.5。，360。]，区间11=(22.5。，67.5。]，区间 111=(67.5。，112.5。]，区间1￥=(112.5。，157.5。]，区间￥=(157.5。，202.5。]， 区间 VI =(202.5°，247.5° ...

【专利技术属性】
技术研发人员：吴自银，周洁琼，赵荻能，李守军，尚继宏，梁裕扬，周勐佳，
申请(专利权)人：国家海洋局第二海洋研究所，
类型：发明
国别省市：浙江;33