用于地形匹配的3D Zernike矩快速计算方法技术

本发明专利技术涉及一种用于地形匹配的3D Zernike矩快速计算方法，步骤：(1)利用本发明专利技术提出的几何矩动态存储技术，降低3D Zernike矩计算中的冗余计算量；(2)利用本发明专利技术提出的坐标分离计算技术，降低几何矩计算中的冗余计算量；(3)利用本发明专利技术提出的对称性技术，对几何矩分类计算，提高计算效率；(4)综合使用上述三种快速计算技术，实现对用于地形匹配的3D Zernike矩的快速计算。本发明专利技术与传统计算方法相比，在地形匹配中，可提高3D Zernike矩的计算速度约两个数量级。

A fast 3D Zernike moment computation method for terrain matching

The invention relates to a method for fast computation of terrain matching 3D Zernike moment steps: (1) using geometric moment dynamic storage technology provided by the invention reduces the computation complexity of the redundant 3D Zernike moment calculation; (2) separation using coordinate computing technology provided by the invention, to reduce the amount of redundant computing geometric moments in the calculation; (3) by the symmetry of the technology provided by the invention, the calculation of geometric moment classification, improve the computational efficiency; (4) the comprehensive use of the three fast computing technology, implementation for rapid calculation of terrain matching 3D Zernike moments. Compared with the traditional calculation method of the invention, the terrain matching, can improve the calculation speed of 3D Zernike moment is about two orders of magnitude.

【技术实现步骤摘要】
用于地形匹配的3DZernike矩快速计算方法
本专利技术涉及一种用于地形匹配的3DZernike矩快速计算方法，属于导航

技术介绍
3DZernike矩具有旋转平移尺度变换不变特性，以及适应性强、鲁棒性好等优点，在三维图形检索、分子生物学等诸多领域已得到广泛研究。地形在空间中表现为三维曲面，利用3DZernike矩可提取地形曲面的起伏特征，并构造特征向量，应用于地形匹配中。由于3DZernike矩的优良特性，使得应用3DZernike矩的地形匹配精度很高。制约3DZernike矩在应用的主要问题是大计算量问题，3DZernike矩的计算用时随地形区域的增大呈平方增长，随3DZernike矩阶次的提高呈指数增长。为提高3DZernike矩的计算效率，利用归一化单位球的体对称性降低计算量的许多方法被提出，并在计算三维实体和三维表面中取得较好效果。上述方法的针对对象，即三维实体和三维表面，与地形曲面之间存在差异，单位球的体对称性并不能对地形曲面的计算带来改善，相比之下，计算模板的面对称性更有应用价值，为3DZernike矩的计算简化提供了可能。按照传统算法，用于地形匹配的3DZernike矩的计算方法如下：(1)数据进行归一化处理。将地形高程数据归一化到单位球内。(2)计算几何矩。几何矩，定义为单位球内的积分Mrst＝∫|X|≤1f(X)xrysztdX,X＝(x,y,z)T(1)地形高程图在三维空间的表现为一曲面，在曲面上f(X)＝1，其它位置f(X)＝0。故几何矩Mrst可简化为其中N表示曲面包含的点数，(xi,yi,zi)表示三维空间中所有位于曲面上的点的坐标。(3)计算3DZernike矩。其中n为3DZernike矩的阶，l为重复度，且满足n≥1，n-l为偶数，m∈[-l,l]，为组合数，其与地形数据无关，仅取决于r,s,t,n,l,m六个参量，可以预先计算，使用时查表。当m<0时，有
技术实现思路
本专利技术技术解决问题：大幅提高了用于地形匹配的3DZernike矩的计算速度，降低了冗余计算。本专利技术的技术要点：1.几何矩动态存储技术由3DZernike矩的计算式(3)可知，在计算不同阶的时，几何矩Mrst存在大量重复使用的现象，即高阶次的3DZernike矩计算时会反复使用低阶的几何矩。因此，在第一次计算某组(r,s,t)的几何矩时，为几何矩分配动态存储空间，将(r,s,t)与Mrst保存下来。每当需要计算Mrst时，根据(r,s,t)值搜索存储空间中是否已存在，若存在则直接调用，若无则计算该值并存储。2.坐标分离计算技术由几何矩计算式(2)可知，在计算不同(r,s,t)组合时，xr、ys和zr存在大量重复使用的现象。由于r,s,t三个参数的数目比(r,s,t)组合的数目少的多，特别是当r,s,t值较大的时候，例如n取10时，独立的r,s,t组合有286种，而r,s,t三个参数的数目仅为33，因此在计算3DZernike矩时，可以根据r,s,t的取值计算一次xr、ys和zr的值并存储，从而在计算几何矩Mrst时直接调用，大幅降低冗余计算。3.对称性技术在计算3DZernike矩前，需进行归一化处理，从而产生了坐标的对称性。在用于地形匹配时，计算模板通常为方形或圆形，也具备对称性。利用归一化和计算模板对称的特性，可以提高几何矩的计算效率。具体的对称性表现如下：由归一化过程，有由模板对称性，当t＝0且r和s至少一个为奇数时，有本专利技术技术解决方案：用于地形匹配的3DZernike矩快速计算方法，步骤如下：(1)利用本专利技术提出的几何矩动态存储技术，降低3DZernike矩计算中的冗余计算量；(2)利用本专利技术提出的坐标分离计算技术，降低几何矩计算中的冗余计算量；(3)利用本专利技术提出的对称性技术，对几何矩分类计算，提高计算效率。所述步骤(1)中，所述几何矩动态存储技术为：通过动态分配存储空间，存储已计算的几何矩Mrst及其索引值(r,s,t)，从而在最终计算3DZernike矩时直接搜索调用，避免冗余计算。所述步骤(2)中，所述坐标分离计算技术为：几何矩中，xr、ys和zr拆开计算，且只计算一次，从而在计算不同r,s,t组合的几何矩Mrst时直接调用，避免冗余计算。所述步骤(3)中，所述对称性技术包括坐标归一化和计算模板对称性，由坐标归一化产生其中N为计算地形区域内的坐标点数；由模板对称性，当t＝0且r和s至少一个为奇数时，几何矩满足本专利技术与现有技术相比的优点在于：(1)本专利技术通过牺牲一定的存储空间，降低了3DZernike矩计算过程的冗余计算量。(2)结合地形曲面的特性，提高了几何矩的计算效率(3)大幅提升了用于地形匹配的3DZernike矩的计算速度。附图说明附图为本专利技术的3DZernike矩快速算法流程图；具体实施方式如图1所示，本专利技术具体实现如下：(1)参考传统算法，预计算组合数制表。(2)将计算区域的数字高程图坐标(x,y,z)归一化至单位球内。(3)根据设定的最高计算阶次，确定r,s,t的取值范围，计算xr、ys和zr并存储。(4)开始计算3DZernike矩(a)根据n确定r,s,t取值(b)搜索是否存在Mrst，若不存在，则转(c)；若存在，则转(d)(c)依据式(7)计算几何矩Mrst，并存储[r,s,t,Mrst]。(d)依据式(3)计算(5)判断n,l,m值是否达到最高阶，若未达到，则转(4)继续计算；若达到，则结束计算。提供以上实施例仅仅是为了描述本专利技术的目的，而并非要限制本专利技术的范围。本专利技术的范围由所附权利要求限定。不脱离本专利技术的精神和原理而做出的各种等同替换和修改，均应涵盖在本专利技术的范围之内。本文档来自技高网...

【技术保护点】
用于地形匹配的3D Zernike矩快速计算方法，其特征在于步骤如下：(1)利用本专利技术提出的几何矩动态存储技术，降低3D Zernike矩计算中的冗余计算量；(2)利用本专利技术提出的坐标分离计算技术，降低几何矩计算中的冗余计算量；(3)利用本专利技术提出的对称性技术，对几何矩分类计算，提高计算效率。

【技术特征摘要】
1.用于地形匹配的3DZernike矩快速计算方法，其特征在于步骤如下：(1)利用本发明提出的几何矩动态存储技术，降低3DZernike矩计算中的冗余计算量；(2)利用本发明提出的坐标分离计算技术，降低几何矩计算中的冗余计算量；(3)利用本发明提出的对称性技术，对几何矩分类计算，提高计算效率。2.根据权利要求1所述的用于地形匹配的3DZernike矩快速计算方法，其特征在于：所述步骤(1)中的几何矩动态存储技术为：通过动态分配存储空间，存储已计算的几何矩Mrst及其阶次(r,s,t)，从而在最终计算3DZernike矩时直接搜索调用，避免冗余计...

【专利技术属性】
技术研发人员：王可东，朱桐乾，王海涌，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京,11