基于二分法的海洋重磁测线交点差计算方法

技术领域

本申请涉及海洋测量领域，尤其涉及基于二分法的海洋重磁测线交点差计算方法。

背景技术

当前海洋重磁测量中，往往需要根据具体项目的布设数量较多的测量线，由这些近似直线的直线段组成测线观测网，观测网中存在大量测量线与测量线之间的交叉点，交点不符值是重磁测量精度评估的重要依据。在数据处理过程中，需要根据测线交叉点的符合程度来评估项目或测区的测量精度，并根据交叉点的精度统计情况进行重力或磁力测量系统差的调整与补偿。

现有的计算交点不符值的主要方法有前后扩充法和测量点分段法，前后扩充法找到两直线的交点附近的测量点需耗费大量的计算机运算资源，不适用于大批量的测量线交点差求取。测量点分段法需要对单个测量线分段内的所有测量点进行搜索，当每一分段线较长时，则会比前后扩充法消耗更多的计算机运算资源，当每一分段线较短时，则可能会找不到相应的交叉点。

因此，亟需一种能快速获取海洋重磁测线的交点不符值，并且不依赖于分段线长度的方法。

发明内容

为克服相关技术中存在的问题，本申请提供基于二分法的海洋重磁测线交点差计算方法，该方法能够快速获取海洋重磁测线的交点不符值，并且不依赖于分段线长度。

本申请提供基于二分法的海洋重磁测线交点差计算方法，包括：

计算主测线和副测线的测线交点P

根据所述交点坐标(X

计算主测线测量点数N

判断是否满足主副测线测量点配对条件，若是，对所述主测线和所述副测线的测量点进行配对，得到测量点组，寻找最小距离测量点组；

根据所述最小距离测量点组寻找相交临近点，根据所述相交临近点计算交点不符值；

所述主副测线测量点配对条件为：

所述主测线测量点数N

或，所述主测线首末测量点和所述副测线首末测量点是否不变。

所述计算主测线和副测线的测线交点P

所述主测线首末测量点包括：主测线起始点P

所述副测线首末测量点包括：副测线起始点P

根据所述主测线首末测量点和所述副测线首末测量点计算所述交点坐标(X

Y

其中，X

所述根据所述交点坐标(X

比较所述交点坐标(X

若所述交点坐标在所述主测线起始点的坐标与所述主测线中点的坐标的范围内，则所述测线交点P

若所述交点坐标在所述主测线中点的坐标与所述主测线终止点的坐标的范围内，则所述测线交点P

若所述交点坐标不在所述主测线起始点的坐标与所述主测线终止点的坐标范围内，则所述主测线和所述副测线无交点，保持所述主测线起始点P

所述主测线中点为所述主测线上序号为i2/2的点。

所述根据所述交点坐标(X

比较所述交点坐标(X

若所述交点坐标在所述副测线起始点的坐标与所述副测线中点的坐标的范围内，则所述测线交点P

若所述交点坐标在所述副测线中点的坐标与所述副测线终止点的坐标的范围内，则所述测线交点P

若所述交点坐标不在所述副测线起始点的坐标与所述副测线终止点的坐标范围内，则所述主测线和所述副测线无交点，保持所述副测线起始点P

所述副测线中点为所述副测线上序号为j2/2的点。

所述计算主测线测量点数N

通过以下公式计算所述主测线测量点数N

N

N

其中，i2为所述主测线终止点P

进一步地，所述判断是否满足主副测线测量点配对条件，还包括：

若否，则根据所述交点坐标(X

所述对所述主测线和所述副测线的测量点进行配对，包括：

将所述主测线起始点P

将所述主测线的测量点P

所述根据所述最小距离测量点组寻找相交临近点，包括：

所述最小距离测量点组包括：主测线最小距离测量点P

取所述主测线最小距离测量点P

所述主测线最小距离测量点P

进一步地，所述根据所述最小距离测量点组寻找相交临近点之后，所述根据所述相交临近点计算交点不符值之前，还包括：

连接所述相交临近点P

连接所述相交临近点P

计算临近线交点P

所述根据所述相交临近点计算交点不符值，包括：

根据所述相交临近点P

Z

其中，Z

本申请提供的技术方案可以包括以下有益效果：

本申请技术方案提供了基于二分法的海洋重磁测线交点差计算方法，本申请技术方案首先计算主测线和副测线的交点坐标，然后根据交点坐标的位置使用测量线中点坐标更新主测线和副测线的起始点坐标或终止点坐标，将测量点的搜索范围缩小，在小范围的主测线和副测线的测量点集合里寻找最小距离测量点组，能够快速根据最小距离测量点组寻找相交临近点，从而能够快速计算出交点不符值。另外，由于不需要人工对测量线进行分段，避免了人工分段带来的丢失测线交点的问题。

应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。

附图说明

通过结合附图对本申请示例性实施方式进行更详细的描述，本申请的上述以及其它目的、特征和优势将变得更加明显，其中，在本申请示例性实施方式中，相同的参考标号通常代表相同部件。

图1是本申请实施例示出的基于二分法的海洋重磁测线交点差计算方法的流程示意图；

图2是本申请实施例示出的海洋重磁测量线交点的示意图；

图3是本申请实施例示出的海洋重磁测量线网的示意图。

具体实施方式

下面将参照附图更详细地描述本申请的优选实施方式。虽然附图中显示了本申请的优选实施方式，然而应该理解，可以以各种形式实现本申请而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本申请更加透彻和完整，并且能够将本申请的范围完整地传达给本领域的技术人员。

在本申请使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本申请。在本申请和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。

应当理解，尽管在本申请可能采用术语“第一”、“第二”、“第三”等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本申请范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本申请的描述中，“多个”的含义是两个或两个以上，除非另有明确具体的限定。

实施例一

在海洋重磁测量中，需要通过计算测量线之间的交点的交点不符值来评估重磁测量的精度，现有的前后扩充法计算复杂度高，会消耗大量的计算机运算资源，现有的测量点分段法依赖于每一分段线的长度。

针对上述问题，本申请实施例提供基于二分法的海洋重磁测线交点差计算方法，以下结合附图详细描述本申请实施例的技术方案。

图1是本申请实施例示出的基于二分法的海洋重磁测线交点差计算方法的流程示意图，图2是本申请实施例示出的海洋重磁测量线交点的示意图，图3是本申请实施例示出的海洋重磁测量线网的示意图，请参阅图1、图2和图3。基于二分法的海洋重磁测线交点差计算方法的具体内容如下：

A1、计算主测线和副测线首末测量点连线的交点坐标。

根据事先选择好的主测线和副测线确定主测线和副测线的首末测量点：主测线起始点P

根据所述主测线首末测量点和所述副测线首末测量点计算所述测线交点P

Y

其中，X

当R

A2、根据交点坐标在主测线上的位置更新主测线首末测量点。

比较所述交点坐标的横坐标X

当X

当X

当X

当Y

更新主测线首末测量点具体为更新主测线计算段的首末测量点。

A3、根据交点坐标在副测线上的位置更新副测线首末测量点。

比较所述交点坐标的横坐标X

当X

当X

当X

当Y

更新副测线首末测量点具体为更新副测线计算段的首末测量点。

A4、判断主测线和副测线的首末测量点是否不变。

判断主测线和副测线的首末测量点i，j，i2，j2是否不变，若是，则执行A7，若否，则执行A5。

A5、计算主测线测量点数N

N

N

其中，i2为所述主测线终止点P

A6、判断N

判断主测线测量点数N

示例性，测量点数阈值设置为300。测量点数阈值根据实际情况而定，此处不作限定。

A7、寻找最小距离测量点组。

对所述主测线和所述副测线的测量点进行配对，得到测量点组，从测量点组中寻找最小距离测量点组。

将主测线上的每一个测量点和副测线上的每一个测量点进行配对，得到(i2-i)×(j2-j)个测量点组，计算每一个测量点组的距离，得到最小距离测量点组。

示例性，主测线上序号为i3的测量点和副测线上的测量点进行配对，得到的测量点组为(i3，j)，(i3，j+1)，...，(i3，j2-1)，(i3，j2)。其中，i＜i3＜i2。

示例性，测量点组(i3，j)的距离为

A8、计算交点不符值。

根据所述最小距离测量点组寻找相交临近点，再求相交临近线的测线交点的交点不符值。

本申请实施例提供了基于二分法的海洋重磁测线交点差计算方法，本申请实施例首先计算主测线和副测线的交点坐标，然后根据交点坐标的位置使用测量线中点坐标更新主测线和副测线的起始点坐标或终止点坐标，将测量点的搜索范围缩小，在小范围的主测线和副测线的测量点集合里寻找最小距离测量点组，能够快速根据最小距离测量点组寻找相交临近点，从而能够快速计算出交点不符值。另外，由于不需要人工对测量线进行分段，避免了人工分段带来的丢失测线交点的问题。

实施例二

计算交点不符值需要先根据所述最小距离测量点组寻找相交临近点，然后求相交临近线交点处的主副测线重磁测量值的差值。计算交点不符值的具体过程如下：

最小距离测量点组为(P

取主测线最小距离测量点P

根据所述相交临近点计算交点不符值，连接相交临近点P

计算临近线交点P

根据所述相交临近点P

Z

其中，Z

本申请实施例提供了计算交点不符值的具体过程，根据主测线和副测线上相交临近点的距离和相交临近线交点和相交临近点的距离计算出在测线交点处的主测线重磁测量值和副测线重磁测量值，再计算出主测线重磁测量值和副测线重磁测量值的差值，得到交点不符值。本申请实施例根据已知的各点处的重磁测量值，仅需计算不同点的距离即可求出交点不符值，计算过程较为简单，能够快速得到交点不符值。

上文中已经参考附图详细描述了本申请的方案。在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详细描述的部分，可以参见其他实施例的相关描述。本领域技术人员也应该知悉，说明书中所涉及的动作和模块并不一定是本申请所必须的。另外，可以理解，本申请实施例方法中的步骤可以根据实际需要进行顺序调整、合并和删减，本申请实施例装置中的模块可以根据实际需要进行合并、划分和删减。

附图中的流程图和框图显示了根据本申请的多个实施例的系统和方法的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现中，方框中所标记的功能也可以以不同于附图中所标记的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或操作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

以上已经描述了本申请的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术的改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。