航磁测量数据的检索提取方法与流程

1.本发明是关于航磁测量数据应用与服务领域，特别是关于航磁测量数据的检索提取方法。


背景技术：

2.航空磁力测量，又称航空磁测或航空磁力勘探，简称航磁测量，是将航空磁力仪及其配套的辅助设备装载在飞行器上，在测量地区(测区)上空按照预先设定的测线和高度对地磁场强度或梯度进行测量的地球物理方法。
3.近年来，随着对航空磁测勘查的需求不断加大，我国新增航磁测量数据持续增长。目前已基本实现我国陆域全覆盖，同一区域存在多次不同比例尺、高度、测量精度、收录方式等航磁数据重复覆盖情况。这些数据通常以测区为单元进行数据存储，在后期数据应用时，需要用户针对测量比例尺、高度、测量精度、收录方式等信息以及研究区域范围进行筛选，从而获得满足需要的最优测量数据。
4.此外，在数据提供时，由于用户无法直接识别原始测量数据，因此，往往将原始测量数据生成相应的网格化数据，当用户需求范围覆盖多个测区时，向用户提供数据时则根据用户需求范围信息分别对不同测区单元的网格化数据进行裁剪，不同测区间在边界处存在数据重叠区。
5.基于此，本技术的发明人发现，目前的航磁测量数据获取流程以人工筛选为主，过程复杂繁琐；数据提供以网格化数据为主，不是基于航磁测量获得的原始测量磁场值数据，受数据处理人为影响大；并且要同时提供多个测区时，重复工作量大，且不同测区边界处的重复覆盖数据，增加了用户数据处理难度。
6.公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解，而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。


技术实现要素：

7.本发明的目的在于提供一种航磁测量数据的检索提取方法，其能够基于原始测量数据自动化的实现航磁测量数据的检索级提取。
8.为实现上述目的，本发明提供了一种航磁测量数据的检索提取方法，所述检索提取方法包括：获取用户检索条件，其中，所述用户检索条件包括空间检索范围以及提取精度，所述用户检索条件包括的空间检索范围为最初空间检索范围，所述提取精度包括测区级，测线级，测点级中的一个；获取与空间检索范围相交的全部测区的参数数据，并根据预设原则对与空间检索范围相交的全部测区进行打分；在所述与空间检索范围相交的全部测区中，获取得分最高的测区，作为第一测区；获取所述第一测区的空间范围与空间检索范围相交重叠的空间范围，作为与所述第一测区对应的重叠范围；从空间检索范围剔除与所述第一测区对应的重叠范围相交覆盖的部分，得到更新后的空间检索范围；在所述与空间检索范围相交的且没有作为或作为过第一测区的其他测区中，将打分最高的测区作为更新后
的第一测区，重复获取所述第一测区的空间范围与更新后的空间检索范围相交重叠的空间范围，作为所述更新后的第一测区对应的重叠范围的步骤，直至最初空间检索范围被全部覆盖或者与空间检索范围相交的全部测区对比完毕；根据所述提取精度以及所述重叠范围，提取数据，其中，所述重叠范围为所有的第一测区的空间范围与空间检索范围相交重叠的空间范围的和。
9.在本发明的一实施方式中，根据所述提取精度以及所述重叠范围，提取数据包括：当所述提取精度为测区级时，将所述重叠范围对应的所述第一测区的全部实体数据，作为所述第一测区的数据提取空间范围；根据所述第一测区的数据提取空间范围的和，提取数据。
10.在本发明的一实施方式中，根据所述提取精度以及所述重叠范围，提取数据包括：当所述提取精度为非测区级时，获取所述重叠范围对应的一个第一测区内的全部测线的信息；根据所述全部测线的信息，计算所述第一测区的数据提取空间范围；获取所述重叠范围对应的另一个第一测区内的全部测线的信息，并重复计算所述第一测区的数据提取空间范围的步骤，直至与所述重叠范围对应的全部测区计算完成；根据所述第一测区的数据提取空间范围的和，提取数据。
11.在本发明的一实施方式中，当所述提取精度为测线级时，根据所述全部测线的信息，计算所述第一测区的数据提取空间范围，包括：获取所述第一测区内全部测线信息中，与所述第一测区对应的重叠范围存在交点的测线对应的全部测点的信息；将所有与所述第一测区对应的重叠范围存在交点的测线对应的全部测点的信息作为所述第一测区的数据提取空间范围。
12.在本发明的一实施方式中，当所述提取精度为测点级时，根据所述全部测线信息，计算所述第一测区的数据提取空间范围，包括：在所述第一测区内全部测线信息中，若一条测线与所述第一测区对应的重叠范围存在交点，则计算所述测线与所述第一测区对应的重叠范围的交点；对所述测线上的测点进行判断，获取所述测线的测点在所述第一测区对应的重叠范围内的全部测点的信息，将所述全部测点的信息作为所述测线的提取空间范围；重复上述步骤，直至完成全部测线信息中与所述第一测区对应的重叠范围存在交点的测线的提取空间范围的计算；将所有与所述第一测区对应的重叠范围存在交点的测线的提取空间范围，作为所述第一测区的数据提取空间范围。
13.在本发明的一实施方式中，所述获取所述测线的测点在所述第一测区对应的重叠范围内的全部测点的信息，包括：根据以下公式分别计算所述测线两个交点附近的测点索引记index1、index2；
[0014][0015][0016]
其中，index1、index2取整数，(x
cross1
，y
cross1
)、(x
cross2
，y
cross2
)是测线与第一测区对
应的重叠范围存在两个交点坐标，(x
start，ystatr
)是测线起点坐标、(x
end
，y
end
)是终点坐标，index
start
是测线起点索引、index
end
是测线终点索引；获取所述测线上与所述index1、index2两个索引值对应的测点坐标(x
index1
，y
index1
)(x
index2
，y
index2
)，判断测点(x
index1
，y
index1
)(x
index2
，y
index2
)是否处于所述第一测区对应的重叠范围内；若处于，则逐一判断自(x
index1
，y
index1
)起，远离测点(x
index2
，y
index2
)，靠近测线起点或终点方向的连续相邻的测点是否处于第一测区对应的重叠范围内，直至找到最后一个处于第一测区对应的重叠范围内的测点，记录该测点的索引，作为测点数据提取的临界索引threshold1；同样逐一判断自测点(x
index2
，y
index2
)起，远离测点(x
index1
，y
index1
)，靠近测线起点或终点方向连续相邻的测点中最后一个处于第一测区对应的重叠范围内的测点索引，作为测点数据提取的另一个临界索引threshold2；将所述测线处于两个临界索引threshold1、threshold2之间的全部测点数据作为所述测线的测点在所述第一测区对应的重叠范围内的全部测点的信息。
[0017]
在本发明的一实施方式中，所述用户检索条件还包括属性约束条件；所述与空间检索范围相交的全部测区的参数数据为满足属性约束条件的且与空间检索范围相交的全部测区的参数数据；所述根据预设原则对与空间检索范围相交的全部测区进行打分包括：根据预设原则对所述满足属性约束条件并与空间检索范围相交的全部测区进行权值计算，对总权值进行排序，将总权值最大的测区作为得分最高的测区。
[0018]
在本发明的一实施方式中，在获取与空间检索范围相交的全部测区的参数数据之前，所述方法还包括：获取航磁原始编码数据；将一个测区的经航空地球物理数据处理后形成标准交换编码数据作为一个实体文件采用二进制形式进行存储，所述测区的实体文件中存储有所述测区的测点场值信息以及每个测点的索引信息。
[0019]
在本发明的一实施方式中，所述提取数据之后，还包括：通过图示向用户反馈提取到的数据。
[0020]
与现有技术相比，根据本发明的航磁测量数据的检索提取方法，能够实现全流程自动筛选并提取满足用户需求的指定区域单次覆盖最优航磁测量数据，提取效率高及提取精准。
[0021]
在数据提取时，可以按照测区、测线、测点三个数据提取精度级别进行数据获取，满足不同数据场景需求，特别是在测点级数据精度提取时，能够准确获取用户指定区域范围内单次覆盖的最优测量数据，并且解决来自不同测区的测量数据的边界处存在数据重复覆盖的问题，降低用户后期数据处理难度。
[0022]
此外本发明针对原始航磁测量场值和相关测量参数值定义统一的数据内容和数据组织模式，形成统一标准交换编码数据格式进行实体数据的组织存储，后期不同测区数据可基于该格式合并存储到同一个数据文件，减少用户多个测区文件进行处理的难度，同时基于原始航磁测量场值向用户提供数据服务，解决了当前基于网格化数据进行数据服务会存在虚假磁异常的问题。
附图说明
[0023]
图1是根据本发明一实施方式的航磁测量数据的检索提取方法的流程图；
[0024]
图2是根据本发明一实施方式的使用测区检索提取的实体数据提取空间范围示意图；
[0025]
图3是根据本发明一实施方式的使用测线检索提取的实体数据提取空间范围示意图；
[0026]
图4是根据本发明一实施方式的使用测点检索提取的实体数据提取空间范围示意图。
具体实施方式
[0027]
下面结合附图，对本发明的具体实施方式进行详细描述，但应当理解本发明的保护范围并不受具体实施方式的限制。
[0028]
除非另有其它明确表示，否则在整个说明书和权利要求书中，术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分，而并未排除其它元件或其它组成部分。
[0029]
如图1所示，其为根据本发明优选实施方式的一种航磁测量数据的检索提取方法的流程图。本实施例提供的一种航磁测量数据的检索提取方法，包括以下步骤：
[0030]
步骤1，获取用户检索条件，其中，所述用户检索条件包括空间检索范围以及提取精度，所述用户检索条件包括的空间检索范围为最初空间检索范围，所述提取精度包括测区级，测线级，测点级中的一个。
[0031]
步骤2，获取与空间检索范围相交的全部测区的参数数据，并根据预设原则对与空间检索范围相交的全部测区进行打分。其中，预设原则可根据用户需要自定义，可扩展。
[0032]
在一种实现方式中，所述用户检索条件还包括属性约束条件。例如，关于测量比例尺、测量时间、数据记录方式、平均飞行高度、测量总精度等属性信息的约束条件。
[0033]
所述与空间检索范围相交的全部测区的参数数据为满足属性约束条件的且与空间检索范围相交的全部测区的参数数据。
[0034]
相应的，步骤2中的根据预设原则对与空间检索范围相交的全部测区进行打分可以包括：根据预设原则对满足所述属性约束条件并与空间检索范围相交的全部测区进行权值计算，将权值最大的测区作为得分最高的测区。
[0035]
本实施例中，属性约束可包括但不限于测量比例尺、测量时间、数据记录方式、平均飞行高度、测量总精度等属性信息，属性约束可自定义，可扩展；检索条件中的空间范围信息基于2000国家大地坐标系。一个测区为一个数据存储单元，对应一个实体存储文件，文件内数据按照测线组织，测线则是由测点按空间位置、测量时间、测点采集顺序逐一存储，并建立测点索引。
[0036]
根据检索条件(包括空间检索范围和属性信息检索条件)在测区级元数据信息中进行空间-属性一体化检索，获取满足条件的测区编码信息及预设评分原则中涉及的属性元数据信息；
[0037]
最优检索原则按照用户从测区元数据中选定的测量比例尺、平均飞行高度、磁测总精度、数据记录方式等属性信息作为评分项并分别赋予权值，进行打分排序。
[0038]
步骤3，在所述与空间检索范围相交的全部测区中，获取得分最高的测区，作为第一测区；
[0039]
步骤4，获取所述第一测区的空间范围与空间检索范围相交重叠的空间范围，作为与所述第一测区对应的重叠范围。并记录所述测区的测区编码、测区优先级别和第一测区
对应的重叠范围。
[0040]
步骤5，从空间检索范围剔除与所述第一测区对应的重叠范围相交覆盖的部分，得到更新后的空间检索范围。
[0041]
步骤6，在所述与空间检索范围相交的且没有作为或作为过第一测区的其他测区中，将打分最高的测区作为更新后的第一测区，重复获取所述第一测区的空间范围与更新后的空间检索范围相交重叠的空间范围，作为所述更新后的第一测区对应的重叠范围的步骤，直至最初空间检索范围被全部覆盖或者与空间检索范围相交的全部测区对比完毕。
[0042]
具体的，可以获取满足检索条件且单次覆盖检索区域最优测区的测区编码和第一测区对应的重叠范围、实体数据存储路径信息列表。
[0043]
步骤7，根据所述提取精度以及所述重叠范围，提取数据，其中，所述重叠范围为所有的第一测区的空间范围与空间检索范围相交重叠的空间范围的和。
[0044]
其中，空间检索范围包括最初空间检索范围以及每次更新后的空间建设范围。
[0045]
由此，通过本实施例提供的航磁测量数据的检索提取方法，能够实现全流程自动筛选并提取满足用户需求的指定区域单次覆盖最优航磁测量数据，在数据提取时，可以按照测区、测线、测点三个数据提取精度级别进行数据获取，满足不同数据场景需求。
[0046]
在一种实现方式中，步骤7可以包括：
[0047]
步骤7011，当所述提取精度为测区级时，将所述重叠范围对应的所述第一测区的全部实体数据，作为所述第一测区的数据提取空间范围。
[0048]
步骤7012，根据所述第一测区的数据提取空间范围的和，提取数据。
[0049]
需要说明的是，如果用户希望获取全部最优测区编码的实体数据，即提取精度为测区级，则根据测区属性信息中的存储位置信息可以直接提取以测区为单元的实体数据以及相应的测区元数据信息，提取的数据信息参阅图2。
[0050]
在又一种实现方式中，步骤7可以包括：
[0051]
步骤7021，当所述提取精度为非测区级时，获取所述重叠范围对应的一个第一测区内的全部测线的信息；其中，测线的信息可以包括相应的测线编码、测线起点坐标、测线终点坐标、测线起点索引、测线终点索引和对应的测区编码列表。
[0052]
步骤7022，根据所述全部测线的信息，计算所述第一测区的数据提取空间范围。
[0053]
步骤7023，获取所述重叠范围对应的另一个第一测区内的全部测线的信息，并重复计算所述第一测区的数据提取空间范围的步骤，直至与所述重叠范围对应的全部测区计算完成。
[0054]
步骤7024，根据所述第一测区的数据提取空间范围的和，提取数据。
[0055]
在一种实现方式中，当所述提取精度为测线级时，步骤7022，可以包括：
[0056]
获取所述全部测线的信息中，与所述第一测区对应的重叠范围存在交点的测线信息，具体的，可以包括：得到符合条件的测线编码、所属测区编码、测线起点坐标、测线终点坐标、测线起点索引、测线终点索引；将所有与所述第一测区对应的重叠范围存在交点的测线对应的全部测点的信息作为所述第一测区的数据提取空间范围。
[0057]
需要说明的是，如果用户希望获取全部最优测区下与第一测区重叠范围相交的测线实体数据，即提取精度为测线级，则根据测区属性元数据信息中的实体数据存储位置信息以及测线元数据信息中的测线编码信息、测线起点索引、测线终点索引信息，直接提取以
测线为单元的实体数据，提取的数据信息参阅图3。
[0058]
在一种实现方式中，当所述提取精度为测点级时，步骤7022，可以包括：
[0059]
在所述全部测线信息中，若一条测线与所述第一测区对应的重叠范围存在交点，则计算所述测线与所述第一测区对应的重叠范围的交点(x
crosst
，y
crosst
)；
[0060]
当测线与第一测区对应的重叠范围存在两个交点坐标(x
cross1
，y
cross1
)、(x
cross2
，y
cross2
)，获取测线起点坐标(x
start，ystart
)、终点坐标(x
end
，y
end
)，测线起点索引index
start
、终点索引index
end
，则根据以下公式分别估算两个交点附近的测点索引记为index1、index2。
[0061][0062][0063]
其中，index1、index2取整数。
[0064]
根据index1、index2两个索引值获取该测线上相应索引的测点坐标(x
index1
，y
index1
)(x
index2
，y
index2
)，判断测点(x
index1
，y
index1
)(x
index2
，y
index2
)是否处于相应第一测区对应的重叠范围内；
[0065]
若处于，则逐一判断自(x
index1
，y
index1
)起，远离测点(x
index2
，y
index2
)，靠近测线起点或终点方向的连续相邻的测点是否处于第一测区对应的重叠范围内，直至找到最后一个处于第一测区对应的重叠范围内的测点，记录该测点的索引，作为测点数据提取的临界索引threshold1；同样逐一判断自测点(x
index2
，y
index2
)起，远离测点(x
index1
，y
index1
)，靠近测线起点或终点方向连续相邻的测点中最后一个处于第一测区对应的重叠范围内的测点索引，作为测点数据提取的另一个临界索引threshold2；
[0066]
若不处于，则逐一判断自(x
index1
，y
index1
)起，靠近测点(x
index2
，y
index2
)方向连续相邻的测点是否处于第一测区对应的重叠范围内，直至找到第一个处于第一测区对应的重叠范围内的测点，记录该测点的索引，作为测点数据提取的临界索引threshold1，同样对自测点(x
index2
，y
index2
)起靠近测点(x
index1
，y
index1
)方向的连续相邻的测点进行逐一判断，找到处于第一测区对应的重叠范围内的第一个测点的索引，作为测点数据提取的另一个临界索引threshold2。
[0067]
需要说明的是由于实际得到的测线只是近似直线，所以只能计算出交点附近的测点索引。计算测线与第一测区对应的重叠范围的交点时，用的都是测线元数据中的测线空间信息。测线空间信息由抽稀后的测点组成。
[0068]
当测线与第一测区对应的重叠范围仅有一个交点(x
cross1
，y
cross1
)时，按上述方法计算得到一个临界索引threshold1，同时进一步判断测线起点或终点是否处于第一测区对应的重叠范围内，若处于，则将该起点或终点的索引作为另一个临界索引threshold2。
[0069]
根据两个临界索引threshold1、threshold2提取所述测线处于两个索引范围的全部测点数据，将所述全部测点的信息作为所述测线的提取空间范围；具体可以提取，包括经纬度、椭球高、离地高度(真高)、原始测量的磁场值、补偿后的磁场值、正常地磁场值、磁日
变值、磁补偿值、δt磁异常值等信息中的多个组合值或全部值。由此，避免对测线上的全部测点进行判断是否处于第一测区对应的重叠范围内，速度快。
[0070]
重复上述步骤，直至完成全部测线信息中与所述第一测区对应的重叠范围存在交点的测线的提取空间范围的计算；将所有与所述第一测区对应的重叠范围存在交点的测线的提取空间范围，作为所述第一测区的数据提取空间范围。提取的数据信息参阅图4。
[0071]
由此，针对测点级数据精度提取时，能够准确获取用户指定区域范围内单次覆盖的最优测量数据，并且解决来自不同测区的测量数据的边界处存在数据重复覆盖的问题，降低用户后期数据处理难度。
[0072]
进一步的，在针对全国性等大范围数据应用时进行数据筛选和提取时，基于本方法可实现批量化、业务化的处理，优势更加明显。
[0073]
在一种实现方式中，在步骤2之前，所述方法还包括：
[0074]
获取航磁原始编码数据；将一个测区的航磁原始编码数据，经处理后形成标准交换编码数据作为一个实体文件进行存储，所述测区的实体文件中存储有所述测区的测点的磁场值信息以及测点索引。测线索引由测线的起点索引和终点索引组成。根据起点索引和终点索引可以提取测线的数据。同一测线上的测点按顺序存储，测区内的全部测点建立统一的测点索引。
[0075]
具体的，上述步骤可以通过以下方式实现：
[0076]
获取一个测区(调查区)的航空磁测数据；以测区为基本单元采集相关元数据信息，元数据信息包括空间信息和属性信息两部分，一个测区的航磁数据对应一条元数据信息。测区空间信息主要是航空磁测勘查所覆盖范围边界坐标信息，或称作测区范围信息，统一采用2000国家大地坐标系，以面要素进行表达。
[0077]
这里的航磁数据是指航磁勘查系统获取的航磁原始编码数据，经航空地球物理数据处理后形成标准交换编码数据作为最终数据存储文件。按照测区-》测线-》测点模式进行逐一存储。一个测区的实体数据作为一个存储文件，里面按照测线、测点进行组织存储。每个测点包括但不限于经纬度、椭球高、离地高度(真高)、原始测量的磁场值、补偿后的磁场值、正常地磁场值、磁日变值、磁补偿值、δt磁异常值等信息。
[0078]
测区属性信息包括但不限于测区编码、测区名称、测量时间、测量比例尺、平均飞行高度、磁测总精度、测线工作量、测量仪器名称、磁力仪灵敏度、数据记录方式、数据存储路径信息等，属性信息可根据实际需要进行扩展；数据存储路径信息包括原始编码数据存储路径和标准交换编码数据存储路径，后续实体数据提取则以标准交换编码数据为基础。采集每条测线的元数据信息，包括空间信息和属性信息两部分，一条测线对应一条元数据信息；测线空间信息以线要素进行表达，统一采用2000国家大地坐标系，测线空间信息由测线起始测点、终止测点的坐标以及抽稀后的测点坐标进行连线组成；属性信息包括测线编码、所属测区编码、架次号、测线方向、起始测点坐标、终止测点坐标、起始测点索引、终止测点索引，其中起始点索引和终止测点索引组成测线在实体文件中的索引信息。基于空间数据库存储测区元数据信息、测线元数据信息。保存经标准交换编码处理的航磁数据，并更新测区属性中的数据存储路径信息。
[0079]
进一步的，本实施例针对原始航磁测量场值和相关测量参数值定义统一的数据内容和数据组织模式，形成标准交换编码数据格式进行实体数据的组织存储，后期不同测区
数据可统一存储到同一个数据文件，减少用户多个测区文件进行处理的难度，同时基于原始航磁测量场值向用户提供数据服务，解决了当前基于网格化数据进行数据服务会存在虚假磁异常的问题。
[0080]
在一种实现方式中，标准交换编码数据在所述数据存储时为二进制编码形式存储，所述数据提取得到的数据为标准交换编码数据采用明码保存的xyz或csv格式数据。
[0081]
由此，相比于原始数据是飞机测量系统直接采集的数据，需要专业的软件才能识别，用户无法直接读取，本技术采用xyz或csv等格式进行分发，便于用户使用普通的数据处理软件就能识别，便于建立数据索引。
[0082]
在一种实现方式中，所述提取数据之后，还包括：通过图示向用户反馈提取到的数据，具体可以包括：
[0083]
数据图：将查询到的测点(反应磁场数据)按比例尺落到图面，以按颜色渲染的圆点表示，点的中心点为测点的经纬度位置，磁场(或其他值)按输入的色棒颜色(默认)进行分配，可按线性进行分配，不同颜色对应不同磁场值，同时在圆点上标注磁场数值(标注磁场异常值及离地高度(真高)，小数点位数最小0.01nt和0.1m)。
[0084]
工作程度图：以面要素表达，按前面最终优选的测区，按不同颜色区分形成工作程度图，作为附图，同时以测区的元数据信息标注图例；同时在其上绘制查询检索范围，以无填充的面表示。
[0085]
由此，输出仅以实际的磁场值表示的磁场数据图，不对数据进行任何改变。
[0086]
通过以下示例对检索部分进行详细说明。
[0087]
用户检索条件包括指定的空间检索范围和属性约束条件(如测量比例尺、数据收录方式、平均飞行高度、测量时间、测量仪器型号等测区属性元数据中的一个或多个)；对各个约束由用户确定相应的打分权值，总权值函数定义为m＝∑f(xi)。
[0088]
测区检索：
[0089]
步骤1、根据用户指定的空间检索范围和属性约束条件，基于测区元数据进行空间-属性一体化检索，检索出全部满足用户检索约束条件的测区，测区总数记为n；
[0090]
步骤2、根据预设的评分条件信息，对测区数据质量进行打分，评分条件可根据用户需求自定义，如可约定比例尺为1∶5万权重加40分、1∶10万权重加35分、1∶20万权重加30分、1∶50万权重加25分、1∶100万权重加20分；数据收录方式为数字记录，权重加30分，模拟记录加15分；测量飞行高度最高权重30分；最后统计各个条件的权重累加值，得到总权重值。
[0091]
例如，可以设置m＝f(x1) f(x2) f(x3)；
[0092]
其中，x1指比例尺的权重值，x2指数据收录方式的权重值，x3指测量飞行高度的权重值，当由该三项获得的权重值相同时，则进一步对比测区空间范围与检索范围的重叠率，重叠率大附加权重高。
[0093]
获取全部符合用户检索约束条件要求的测区及相应的总权值，按照总权值由大到小排序，确定优先级排序，总权值越大优先级越高，得到最终的测区优先级列表list1，包括测区编号、测区空间范围、总权值、优先级；
[0094]
步骤3、从测区优先级列表list1中获取优先级最高的测区，作为第一测区，得到该测区的测区编码、测区空间范围信息，并将该测区信息从list1中删除；
[0095]
步骤4、计算得到第一测区的测区空间范围与空间检索范围相交重叠的空间范围，记为该第一测区对应的重叠范围，添加至测区数据提取列表list2，包括测区编码、优先级、第一测区对应的重叠范围、实体数据存储路径；
[0096]
步骤5、计算从空间检索范围剔除与测区范围相交覆盖的部分，得到新的空间检索范围；
[0097]
步骤6、重复步骤3再得到list1中优先级最高的测区与步骤5中得到的新空间检索范围进行比较，如果两者存在重叠区域，则重复步骤4-5；如果两者不存在重叠区域，则重复步骤3继续从list1中选取下一个优先级最高的测区，作为更新后的第一测区，进行比较；
[0098]
步骤7、以此类推，直至用户约束条件中最初设置的空间检索范围被全部覆盖或者list1中全部测区对比完毕，最终得到满足用户约束条件的最优的测区数据提取列表list2，测区数量记为m，m≤n。
[0099]
测线检索：
[0100]
步骤1，从测区数据提取列表list2中获取一个测区的测区编码、该测区对应的重叠范围、实体数据存放路径；
[0101]
步骤2，根据测区编码、测区对应的重叠范围信息在测线元数据信息中进行空间-属性一体化检索，得到该测区内满足条件的全部测线信息，包括测线编码、所属测区编码、测线空间信息、测线起点坐标、测线终点坐标、测线起点索引、测线终点索引等信息；
[0102]
步骤3，获取测线一条测线空间信息并计算其与第一测区对应的重叠范围的两个交点(x
cross1
，y
cross1
)、(x
cross2
，y
cross2
)；
[0103]
步骤4，根据实体数据存放路径信息，以及步骤2获取的测线在实体文件中的起止索引信息，提取与第一测区对应的重叠范围存在交点的所有测线的全部测点信息。
[0104]
测点检索：重复测线检索的步骤1-3。
[0105]
步骤4，根据步骤2获取的测线起点坐标(x
start
，y
start
)、测线终点坐标(x
end
，y
end
)、测线起点索引index
start
、测线终点索引index
end
，步骤3获取的一条测线与第一测区的重叠范围的交点(x1，y1)(x2，y2)，计算两个索引值index1、index2；
[0106][0107][0108]
步骤5，获取测区中的实体数据存储路径信息，在实体数据文件中根据index1、index2两个索引获取对应的测点坐标(x
index1
，y
index1
)(x
index2
，y
index2
)，判读测点(x
index1
，y
index1
)(x
index2
，y
index2
)是否处于相应第一测区对应的重叠范围内；若处于，则进一步判断自测点(x
index1
，y
index1
)起，远离测点(x
index2
，y
index2
)且靠近测线起点或终点方向的连续相邻测点是否处于第一测区对应的重叠范围内，直至找到最后一个处于第一测区对应的重叠范围内测点，记录该测点索引，作为测点数据提取的临界索引threshold1；同样逐一判断自测点(x
index2
，y
index2
)起，远离测点(x
index1
，y
index1
)，靠近测线起点或终点方向连续相邻的测点中
最后一个处于第一测区对应的重叠范围内的测点索引，作为测点数据提取的另一个临界索引threshold2。
[0109]
若不处于，则逐一判断自(x
index1
，y
index1
)起，靠近测点(x
index2
，y
index2
)方向连续相邻的测点是否处于第一测区对应的重叠范围内，直至找到第一个处于第一测区对应的重叠范围内的测点，记录该测点的索引，作为测点数据提取的临界索引threshold1，同样对自测点(x
index2
，y
index2
)起靠近测点(x
index1
，y
index1
)方向的连续相邻的测点进行逐一判断，找到处于第一测区对应的重叠范围内的第一个测点的索引，作为测点数据提取的另一个临界索引threshold2。
[0110]
当测线与第一测区对应的重叠范围仅有一个交点时，按上述方法计算得到一个临界索引，同时进一步判断测线起点或终点是否处于第一测区对应的重叠范围内，若处于，则将该起点或终点的索引作为另一个临界索引。
[0111]
步骤6，根据步骤5获取的两个临界索引直接提取该条测线中两个临界索引范围的全部测点的磁场值记录信息；
[0112]
步骤7，重复步骤3-6，直至该测区内满足条件的全部测线的实体数据均被提取完毕；
[0113]
步骤8，重复步骤1，继续从测区数据提取列表list2中获取下一个测区的测区编码、第一测区对应的重叠范围、实体数据存放路径，并执行步骤2-7；
[0114]
以此类推，直至list2中全部测区数据被提取完毕，数据采用标准交换编码数据(xyz格式或csv格式)提供给用户，可同时提供数据图示。
[0115]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0116]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
[0117]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0118]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0119]
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述
并非想将本发明限定为所公开的精确形式，并且很显然，根据上述教导，可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用，从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。

技术特征：
1.一种航磁测量数据的检索提取方法，其特征在于，所述检索提取方法包括：获取用户检索条件，其中，所述用户检索条件包括空间检索范围以及提取精度，所述用户检索条件包括的空间检索范围为最初空间检索范围，所述提取精度包括测区级，测线级，测点级中的一个；获取与空间检索范围相交的全部测区的参数数据，并根据预设原则对与空间检索范围相交的全部测区进行打分；在所述与空间检索范围相交的全部测区中，获取得分最高的测区，作为第一测区；获取所述第一测区的空间范围与空间检索范围相交重叠的空间范围，作为与所述第一测区对应的重叠范围；从空间检索范围剔除与所述第一测区对应的重叠范围相交覆盖的部分，得到更新后的空间检索范围；在所述与空间检索范围相交的且没有作为或作为过第一测区的其他测区中，将打分最高的测区作为更新后的第一测区，重复获取所述第一测区的空间范围与更新后的空间检索范围相交重叠的空间范围，作为所述更新后的第一测区对应的重叠范围的步骤，直至最初空间检索范围被全部覆盖或者与空间检索范围相交的全部测区对比完毕；根据所述提取精度以及所述重叠范围，提取数据，其中，所述重叠范围为所有的第一测区的空间范围与空间检索范围相交重叠的空间范围的和。2.如权利要求1所述的检索提取方法，其特征在于，根据所述提取精度以及所述重叠范围，提取数据包括：当所述提取精度为测区级时，将所述重叠范围对应的所述第一测区的全部实体数据，作为所述第一测区的数据提取空间范围；根据所述第一测区的数据提取空间范围的和，提取数据。3.如权利要求1所述的检索提取方法，其特征在于，根据所述提取精度以及所述重叠范围，提取数据包括：当所述提取精度为非测区级时，获取所述重叠范围对应的一个第一测区内的全部测线的信息；根据所述第一测区内全部测线的信息，计算所述第一测区的数据提取空间范围；获取所述重叠范围对应的另一个第一测区内的全部测线的信息，并重复计算所述第一测区的数据提取空间范围的步骤，直至与所述重叠范围对应的全部测区计算完成；根据所述第一测区的数据提取空间范围的和，提取数据。4.如权利要求3所述的检索提取方法，其特征在于，当所述提取精度为测线级时，根据所述第一测区内全部测线的信息，计算所述第一测区的数据提取空间范围，包括：获取所述第一测区内全部测线的信息中，与所述第一测区对应的重叠范围存在交点的测线对应的全部测点的信息；将所有与所述第一测区对应的重叠范围存在交点的测线对应的全部测点的信息作为所述第一测区的数据提取空间范围。5.如权利要求3所述的检索提取方法，其特征在于，当所述提取精度为测点级时，根据所述全部测线信息，计算所述第一测区的数据提取空间范围，包括：在所述第一测区内的全部测线信息中，若一条测线与所述第一测区对应的重叠范围存
在交点，则计算所述测线与所述第一测区对应的重叠范围的交点；对所述测线上的测点进行判断，获取所述测线的测点在所述第一测区对应的重叠范围内的全部测点的信息，将所述全部测点的信息作为所述测线的提取空间范围；重复上述步骤，直至完成全部测线信息中与所述第一测区对应的重叠范围存在交点的测线的提取空间范围的计算；将所有与所述第一测区对应的重叠范围存在交点的测线的提取空间范围，作为所述第一测区的数据提取空间范围。6.如权利要求5所述的检索提取方法，其特征在于，所述获取所述测线的测点在所述第一测区对应的重叠范围内的全部测点的信息，包括：根据以下公式分别计算所述测线两个交点附近的测点索引记index1、index2；；其中，index1、index2取整数，(x
cross1
，y
cross1
)、(x
cross2
，y
cross2
)是测线与第一测区对应的重叠范围存在两个交点坐标，(x
start
，y
start
)是测线起点坐标、(x
end
，y
end
)是终点坐标，index
start
是测线起点索引、index
end
是测线终点索引；获取所述测线上与所述index1、index2两个索引值对应的测点坐标(x
index1
，y
index1
)(x
index2
，y
index2
)，判断测点(x
index1
，y
index1
)(x
index2
，y
index2
)是否处于所述第一测区对应的重叠范围内；若处于，则逐一判断自(x
index1
，y
index1
)起，远离测点(x
index2
，y
index2
)，靠近测线起点或终点方向的连续相邻的测点是否处于第一测区对应的重叠范围内，直至找到最后一个处于第一测区对应的重叠范围内的测点，记录该测点的索引，作为测点数据提取的临界索引threshold1；同样逐一判断自测点(x
index2
，y
index2
)起，远离测点(x
index1
，y
index1
)，靠近测线起点或终点方向连续相邻的测点中最后一个处于第一测区对应的重叠范围内的测点索引，作为测点数据提取的另一个临界索引threshold2；将所述测线处于两个临界索引threshold1、threshold2之间的全部测点数据作为所述测线的测点在所述第一测区对应的重叠范围内的全部测点的信息。7.如权利要求1所述的检索提取方法，其特征在于，所述用户检索条件还包括属性约束条件；所述与空间检索范围相交的全部测区的参数数据为满足属性约束条件的且与空间检索范围相交的全部测区的参数数据；所述根据预设原则对与空间检索范围相交的全部测区进行打分包括：根据预设原则对所述满足属性约束条件并与空间检索范围相交的全部测区进行权值计算，对总权值进行排序，将总权值最大的测区作为得分最高的测区。
8.如权利要求1所述的检索提取方法，其特征在于，在获取与空间检索范围相交的全部测区的参数数据之前，所述方法还包括：获取航磁原始编码数据；将一个测区的经航空地球物理数据处理后形成统一标准交换编码数据格式作为一个实体文件采用二进制形式进行存储，所述测区的实体文件中存储有所述测区的测点场值信息以及每个测点的索引信息。9.如权利要求1所述的检索提取方法，其特征在于，所述提取数据之后，还包括：通过图示向用户反馈提取到的数据。

技术总结
本发明公开了一种航磁测量数据的检索提取方法，所述检索提取方法包括：获取用户检索条件、与空间检索范围相交的全部测区的参数数据，并根据预设原则对与空间检索范围相交的全部测区进行打分；在所述与空间检索范围相交的全部测区中，获取得分最高的测区，作为第一测区；获取所述第一测区的空间范围与空间检索范围相交重叠的空间范围，作为与所述第一测区对应的重叠范围；从空间检索范围剔除与所述第一测区对应的重叠范围相交覆盖的部分，得到更新后的空间检索范围；根据所述提取精度以及所述重叠范围，提取数据。本发明的航磁测量数据的检索提取方法，可实现全自动化的数据检索、提取，提取效率高及提取精准。提取效率高及提取精准。提取效率高及提取精准。


技术研发人员：李敬敏 王林飞 骆遥 路文芬 冯磊
受保护的技术使用者：中国自然资源航空物探遥感中心
技术研发日：2022.02.18
技术公布日：2022/5/25