不规则测网交叉点遍历搜索方法、存储介质和设备技术

提供一种不规则测网交叉点遍历搜索方法、存储介质和设备，包括：旋转该不规则测网使得主测线或副测线方向基本水平，所述不规则测网为海空重力不规则测网；确定平行于纵横坐标轴的、两个分别包围测网主副测线的最小外接矩形；取出两个外接矩形的重合矩形作为测网交叉点搜索范围，利用中心点将该重合矩形均分为四个子矩形；判定每个子矩形内主副测线包含的测点数乘积是否大于预定值，若是，则将该子矩形继续分为四个子矩形进行递归处理，若否，则输出该子矩形内主副测线包络矩形重合区域的测点索引。根据输出的测点索引，确定存在的交叉点和对应的测点信息。该方法可以在提高准确率的同时，提高搜索效率。提高搜索效率。提高搜索效率。

【技术实现步骤摘要】
不规则测网交叉点遍历搜索方法、存储介质和设备
[0001]本专利技术是中国专利申请号为ZL202010771152.7、专利技术名称为“海空重力测网交叉点搜索方法和装置”的分案申请，其内容通过全文参考结合于此。


[0002]本专利技术属于海空重力勘探领域，涉及一种不规则测网交叉点遍历搜索方法、存储介质和设备。

技术介绍

[0003]海空重力(海空重力为海面和航空重力的统称)测量中，数据往往是离散的一系列点，每条测线测点数量众多，用传统的逐一对比的遍历搜索法，求取交叉点效率非常低。
[0004]中国专利ZL201810100320.2公开了一种航空重力测网交叉点非遍历搜索方法，首先计算获得近似交叉点，并在测线和切割线上寻找离交叉点最近的测点；其次根据寻找到的目标测点，采用跳跃逼近法再次计算近似交叉点，第二次逼近交叉点；然后判断目标测点是否真正相交，如不相交又则对其进行微调滑动至完全相交。该方法避开了传统的遍历搜索方式，使计算效率大幅提升，有效缩短了航空重力测网交叉点搜索的耗时。
[0005]上述专利方法尽管在一定程度上提高了搜索效率，但是，对于大型测网而言，其搜索效率还有进一步提升的需求，并且，对于不规则的测网而言，该方法并不普遍适应，无法搜索到所有的交叉点，准确率会降低。
[0006]因此，有必要研究一种不规则测网交叉点遍历搜索方法、存储介质和设备来解决上述的一个或多个技术问题。

技术实现思路

[0007]为了解决现有技术中的一个或多个技术问题，提供了一种海空重力测网交叉点搜索方法，其特征在于包括：
[0008]确定平行于纵横坐标轴的、两个分别包围测网主、副测线的最小外接矩形；
[0009]确定两个外接矩形的重合矩形作为测网交叉点搜索范围；
[0010]利用中心点将该重合矩形均分为四个子矩形；
[0011]判定每个子矩形内主副测线包含的测点数乘积是否大于预定值，若是，则将该子矩形继续分为四个子矩形进行递归处理，若否，则输出该子矩形内主副测线包络矩形重合区域内的测点索引。
[0012]可以理解的是，所述的测点索引是指用来存储测点位置等信息的数组下标。
[0013]根据所述测点索引，确定存在的交叉点和对应的测点。
[0014]根据本专利技术另一方面，所述预定值可以根据测网规模自适应选取，优选的，所述预定值为5000。
[0015]根据本专利技术另一方面，所述测网为不规则测网，旋转该不规则测网使得主测线或副测线方向为基本水平。
[0016]根据本专利技术另一方面，采用行列式法、投影法或矢量面积法确定存在的交叉点和对应的测点。
[0017]根据本专利技术另一方面，采用行列式法来求解，令两线段的四个端点坐标为(x1,y1),(x2,y2),(x3,y3),(x4,y4)，则两线段所在直线的参数方程分别为
[0018]x＝x1+t1*(x2‑
x1)
[0019]y＝y1+t1*(y2‑
y1)
[0020]x＝x3+t2*(x4‑
x3)
[0021]y＝y3+t2*(y4‑
y3)
[0022]若两直线相交，则x，y为交点坐标，t1和t2为待求参数，正好四个方程，四个未知数，写成矩阵形式如下：
[0023][0024]当0≤t1≤1且0≤t2≤1时，交叉点在线段上。
[0025]根据本专利技术另一方面，将子矩形内主副测线包络矩形重合区域内的测点索引予以保留，所述的包络矩形是指测线包含的所有线段的最小外接矩形的集合。
[0026]根据本专利技术另一方面，矩形递归剖分是自适应剖分，通过利用当前矩形内主、副测线测点数的乘积作为判定条件使得越是测点密集的区域剖分越细，加上利用包络矩形重合区域来选取索引，使得最终取出的索引和交叉点的量级相当，通过搜索查询和比较大小，避免了大量冗余的计算判定过程，极大地提高了遍历搜索交叉点的效率。
[0027]根据本专利技术又一方面，还提供了一种计算机可读存储介质，其特征在于，所述计算机可读存储介质中存储有计算机程序，所述计算机程序适于由处理器加载并执行，以使得具有所述处理器的计算机设备执行前述的方法。
[0028]根据本专利技术又一方面，还提供了一种计算机设备，其特征在于，包括：处理器和存储器；
[0029]所述处理器与存储器相连，其中，所述存储器用于存储计算机程序，所述处理器用于调用所述计算机程序，以使得所述计算机设备执行前述的方法。
[0030]根据本专利技术另一方面，还提供了一种海空重力测网交叉点搜索装置，其特征在于包括：
[0031]第一模块，用于确定平行于纵横坐标轴的、两个分别包围测网主、副测线的最小外接矩形；
[0032]第二模块，用于确定两个外接矩形的重合矩形作为测网交叉点搜索范围；
[0033]第三模块，用于利用中心点将该重合矩形均分为四个子矩形；
[0034]第四模块，用于判定每个子矩形内主副测线包含的测点数乘积是否大于预定值，若是，则将该子矩形继续分为四个子矩形进行递归处理，若否，则输出该子矩形内主副测线包络矩形重合区域内的测点索引；
[0035]第五模块，用于根据所述测点索引，确定存在的交叉点和对应的测点。
[0036]根据本专利技术另一方面，所述预定值可以根据测网规模自适应选取，优选的，所述预
定值为5000。
[0037]根据本专利技术另一方面，所述测网为不规则测网，旋转该不规则测网使得主测线或副测线方向为基本水平。
[0038]根据本专利技术另一方面，第五模块用于根据所述测点索引，采用行列式法、投影法或矢量面积法确定存在的交叉点和对应的测点。
[0039]与现有技术相比，本专利技术具有以下技术效果：提高了搜索效率；提高了搜索准确率；能够普遍适应于任何类型的海空重力测网；尤其是对于不规则测网而言，本专利技术相比于其它方法在准确率和/或效率方面具有明显优势。
附图说明
[0040]图1为根据本专利技术一种优选实施例的海空重力测网交叉点搜索方法的流程图；
[0041]图2为根据本专利技术一种优选实施例的海空重力测网交叉点搜索方法的应用示例；
[0042]图3为根据本专利技术一种优选实施例的重合判断(即确定每个子矩形中主测线的最小外包矩形和副测线的最小外包矩形是否重合)的示意图；
[0043]图4为根据本专利技术一种优选实施例的交叉点求解方法的示意图；
[0044]图5为图2中应用示例中的数据块的设计参数；
[0045]图6为根据本专利技术一种优选实施例的将子矩形内的重合区域的测点索引予以保留的方法示意图；
[0046]图7为根据本专利技术一种优选实施例的海空重力测网交叉点搜索方法的又一应用示例；
[0047]图8为图7中应用示例中未旋转测网的情形下的搜索过程(和结果)示意图；
[0048]图9为图7中应用示例中旋转测网后的搜索过程(和结果)示意图；
[0049]图10为图9中的搜索方法本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种不规则测网交叉点遍历搜索方法，其特征在于包括：旋转该不规则测网使得主测线或副测线方向基本水平，所述不规则测网为海空重力不规则测网；确定平行于纵横坐标轴的、两个分别包围测网主、副测线的最小外接矩形；确定两个外接矩形的重合矩形作为测网交叉点搜索范围；利用中心点将该重合矩形均分为四个子矩形；判定每个子矩形内主副测线包含的测点数乘积是否大于预定值，若是，则将该子矩形继续分为四个子矩形进行递归处理，若否，则取出该子矩形内主副测线包络矩形重合区域内的测点索引；根据所述测点索引，确定存在的交叉点和其对应的测点；其中，采用行列式法来求解确定存在的交叉点和对应的测点，令两线段的四个端点坐标为(x1，y1)，(x2，y2)，(x3，y3)，(x4，y4)，则两线段所在直线的参数方程分别为x＝x1+t1*(x2‑
x1)y＝y1+t1*(y2‑
y1)x＝x3+t2*(x4‑
x3)y＝y3+t2*(y4‑
y3)若两直线相交，则x，y为交点坐标，t1和t2为待求参数，正好四个方程，四个未知数，写...

【专利技术属性】
技术研发人员：徐光晶，舒晴，
申请(专利权)人：中国自然资源航空物探遥感中心，
类型：发明
国别省市：