[发明专利]一种航空重力测网交叉点非遍历搜索方法

权利要求书

1.一种航空重力测网交叉点非遍历搜索方法，其特征在于：包括以下步骤：

(1)加载测网数据，包括以下步骤：

读取测网数据；

提取测线和切割线；

(2)逐一计算交叉点，直至完成所有交叉点的计算；

所述计算交叉点的方法包括如下步骤：

初步计算交叉点位置；

寻找测线、切割线上距离交叉点位置最近的4个目标测点；

通过获得的4个目标测点，采用跳跃逼近法第二次计算交叉点位置；

更新测线、切割线上距离交叉点位置最近的4个目标测点；

判断目标测点之间是否存在交叉点；若否，则微调4个目标测点，并再次判断，直至4个目标测点之间存在交叉点；

计算并缓存交叉点位置信息和重力场值信息；

所述采用跳跃逼近法第二次计算交叉点位置的方法为：

通过选定的测线上距离交叉点位置最近的2个目标测点的投影坐标，构建第三直线；

通过选定的切割线上距离交叉点位置最近的2个目标测点的投影坐标，构建第四直线；

计算第三直线和第四直线的相交点，即为采用跳跃逼近法第二次计算得到的交叉点位置；

(3)输出交叉点的信息。

2.根据权利要求1所述的航空重力测网交叉点非遍历搜索方法，其特征在于：所述初步计算交叉点位置的方法为：

对测线和切割线两端进行外扩；

通过测线两端点的投影坐标构建第一直线；通过切割线两端点的投影坐标构建第二直线；

计算第一直线和第二直线的相交点，即为交叉点位置。

3.根据权利要求2所述的航空重力测网交叉点非遍历搜索方法，其特征在于：所述寻找测线、切割线上距离交叉点位置最近的4个目标测点的方法为：

通过测线两端点的投影坐标，计算测线总长度及测点之间的平均点距；通过测线一个端点至交叉点位置之间的距离，选定测线上距离交叉点位置最近的2个测点；

通过切割线两端点的投影坐标，计算切割线总长度及测点之间的平均点距；通过切割线一个端点至交叉点位置之间的距离，选定切割线上距离交叉点位置最近的2个测点；

所述测线上距离交叉点位置最近的2个测点和切割线上距离交叉点位置最近的2个测点，即为4个目标测点。

4.根据权利要求3所述的航空重力测网交叉点非遍历搜索方法，其特征在于：所述判断目标测点之间是否存在交叉点的方法为：

判断选定的测线上距离交叉点位置最近的2个目标测点，是否真正位于切割线两侧；若否，则判定4个目标测点之间不存在交叉点；

若是，则进一步判断选定的切割线上距离交叉点位置最近的2个目标测点，是否真正位于测线两侧；若否，则判定4个目标测点之间不存在交叉点；

若是，则判定4个目标测点之间存在交叉点。

5.根据权利要求4所述的航空重力测网交叉点非遍历搜索方法，其特征在于：所述微调4个目标测点的方法为：

根据测线或切割线上的目标测点与交叉点位置之间的距离及方向，采取+1或-1个测点位置的方式，滑动微调测线或切割线上2个相邻的目标测点。

6.根据权利要求5所述的航空重力测网交叉点非遍历搜索方法，其特征在于：所述判断选定的测线上距离交叉点位置最近的2个目标测点是否真正位于切割线两侧的方法为：

记选定的测线上距离交叉点位置最近的2个测点分别为L1、L2，记选定的切割线上距离交叉点位置最近的2个测点分别为T1、T2；计算ΔT1T2L1和ΔT1T2L2的向量面积的正、负符号；

若ΔT1T2L1和ΔT1T2L2的向量面积符号相反，则判定选定的测线上距离交叉点位置最近的2个测点真正位于切割线两侧。

7.根据权利要求6所述的航空重力测网交叉点非遍历搜索方法，其特征在于：所述判断选定的切割线上距离交叉点位置最近的2个目标测点是否真正位于测线两侧的方法为：

记选定的测线上距离交叉点位置最近的2个测点分别为L1、L2，记选定的切割线上距离交叉点位置最近的2个测点分别为T1、T2；计算ΔL1L2T1和ΔL1L2T2的向量面积的正、负符号；

若ΔL1L2T1和ΔL1L2T2的向量面积符号相反，则判定选定的切割线上距离交叉点位置最近的2个测点真正位于测线两侧。