CN114528454A - Retrieval and Extraction Method of Aeromagnetic Survey Data - Google Patents

Abstract

The invention discloses a retrieval and extraction method for aeromagnetic survey data. The retrieval and extraction method includes: acquiring user retrieval conditions and parameter data of all survey areas intersecting with the spatial retrieval range, and aligning the spatial retrieval range with the parameter data according to a preset principle. All the intersecting survey areas are scored; among all the survey areas intersecting with the spatial retrieval range, the survey area with the highest score is obtained as the first survey area; the spatial scope of the obtained first survey area intersects with the spatial retrieval scope The overlapping spatial range is taken as the overlapped range corresponding to the first survey area; the part overlapping and covered by the overlapped range corresponding to the first survey area is removed from the spatial retrieval range to obtain the updated spatial retrieval range; according to the Extract the precision as well as the overlap range, extract the data. The retrieval and extraction method of aeromagnetic measurement data of the present invention can realize fully automatic data retrieval and extraction, and has high extraction efficiency and extraction accuracy.

Description

Retrieval and Extraction Method of Aeromagnetic Survey Data

technical field

The invention relates to the application and service field of aeromagnetic survey data, in particular to a retrieval and extraction method of aeromagnetic survey data.

Background technique

Aeromagnetic survey, also known as aeromagnetic survey or aeromagnetic survey, referred to as aeromagnetic survey, is to load the aeromagnetic instrument and its supporting auxiliary equipment on the aircraft, and follow the preset survey line over the survey area (measurement area). and altitude geophysical methods for measuring the strength or gradient of the geomagnetic field.

In recent years, with the increasing demand for aeromagnetic surveying, my country's new aeromagnetic survey data continues to grow. At present, the full coverage of my country's land area has been basically achieved, and the same area has repeated coverage of aeromagnetic data of different scales, altitudes, measurement accuracy, and recording methods. These data are usually stored in units of survey areas. In the later data application, users are required to filter information such as measurement scale, height, measurement accuracy, recording method, and the scope of the research area, so as to obtain the optimal measurement data that meets the needs.

In addition, when the data is provided, since the user cannot directly identify the original measurement data, the corresponding gridded data is often generated from the original measurement data. The demand range information cuts the gridded data of different survey area units respectively, and there are data overlapping areas at the boundaries of different survey areas.

Based on this, the inventor of the present application found that the current process of obtaining aeromagnetic measurement data is mainly based on manual screening, and the process is complicated and tedious; , which is greatly affected by data processing; and when multiple survey areas are to be provided at the same time, the repetitive workload is large, and the repeated coverage of data at the boundaries of different survey areas increases the difficulty of user data processing.

The information disclosed in this Background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person of ordinary skill in the art.

SUMMARY OF THE INVENTION

The purpose of the present invention is to provide a retrieval and extraction method of aeromagnetic measurement data, which can automatically realize retrieval-level extraction of aeromagnetic measurement data based on the original measurement data.

In order to achieve the above object, the present invention provides a retrieval and extraction method for aeromagnetic survey data, the retrieval and extraction method includes: obtaining user retrieval conditions, wherein the user retrieval conditions include a spatial retrieval range and an extraction accuracy, and the user retrieval conditions include: The spatial retrieval scope included in the retrieval conditions is the initial spatial retrieval scope, and the extraction precision includes one of the survey area level, the survey line level, and the survey point level; the parameter data of all survey areas intersecting with the spatial retrieval scope are obtained, and based on the prediction Suppose the principle is to score all the survey areas intersecting with the spatial retrieval range; in all the survey areas intersecting with the spatial retrieval range, obtain the survey area with the highest score as the first survey area; obtain the first survey area. The spatial range where the spatial range intersects and overlaps with the spatial retrieval range is taken as the overlapping range corresponding to the first survey area; the part that intersects and covers the overlapped range corresponding to the first survey area is removed from the spatial retrieval range to obtain the updated Spatial search range; among the other survey areas that intersect with the spatial search range and that have not acted as or have been used as the first survey area, take the survey area with the highest score as the updated first survey area, and repeat the acquisition of the first survey area. The spatial range where the spatial range of the survey area intersects with the updated spatial search range is taken as the step of the overlapped range corresponding to the updated first survey area, until the initial spatial search range is completely covered or intersects with the spatial search range. All survey areas are compared; data is extracted according to the extraction accuracy and the overlap range, wherein the overlap range is the sum of the overlapped spatial ranges of all the spatial ranges of the first survey area and the spatial retrieval range.

In an embodiment of the present invention, according to the extraction accuracy and the overlap range, extracting data includes: when the extraction accuracy is at the survey area level, extracting all the first survey areas corresponding to the overlap range The entity data is used as the data extraction space range of the first survey area; the data is extracted according to the sum of the data extraction space ranges of the first survey area.

In an embodiment of the present invention, according to the extraction accuracy and the overlap range, extracting data includes: when the extraction accuracy is a non-measurement area level, acquiring the data in a first survey area corresponding to the overlap range information of all survey lines; calculate the data extraction space range of the first survey area according to the information of all survey lines; obtain information of all survey lines in another first survey area corresponding to the overlapping range, and The step of calculating the data extraction spatial range of the first survey area is repeated until the calculation of all the survey areas corresponding to the overlapping range is completed; and data is extracted according to the sum of the data extraction spatial ranges of the first survey area.

In an embodiment of the present invention, when the extraction accuracy is at the survey line level, calculating the data extraction space range of the first survey area according to the information of all survey lines includes: acquiring the first survey line. In the information of all survey lines in the area, the information of all survey points corresponding to survey lines whose overlapping ranges corresponding to the first survey area have intersection points; all survey lines whose overlapping ranges corresponding to the first survey area have intersection points The information of all the corresponding measuring points is used as the data extraction spatial range of the first measuring area.

In an embodiment of the present invention, when the extraction accuracy is at the measuring point level, calculating the data extraction space range of the first survey area according to all the survey line information, including: in the first survey area In all survey line information, if there is an intersection point between a survey line and the overlapping range corresponding to the first survey area, the intersection point of the survey line and the overlap range corresponding to the first survey area is calculated; The measurement points on the measurement line are judged, and the information of all measurement points in the overlapping range corresponding to the first measurement area of the measurement points of the measurement line is obtained, and the information of all measurement points is used as the extraction space of the measurement line. range; repeat the above steps until the calculation of the extraction space range of the survey lines where the overlapping range corresponding to the first survey area in all survey line information has an intersection point is completed; all overlapping ranges corresponding to the first survey area exist The extraction space range of the survey line at the intersection is taken as the data extraction space range of the first survey area.

In an embodiment of the present invention, the acquiring the information of all the measuring points of the measuring point of the survey line in the overlapping range corresponding to the first survey area includes: calculating the two points of the survey line respectively according to the following formulas Indexes of measuring points near each intersection are marked INDEX ₁ , INDEX ₂ ;

Wherein, INDEX ₁ and INDEX ₂ are integers, (X _cross1 , Y _cross1 ), (X _cross2 , Y _cross2 ) are two intersection coordinates in the overlapping range corresponding to the survey line and the first survey area, (X _start, Y _statr ) are the coordinates of the starting point of the survey line, (X _end , Y _end ) are the coordinates of the end point, INDEX _start is the index of the starting point of the survey line, and INDEX _end is the index of the end point of the survey line; obtain two indexes on the survey line and the INDEX ₁ and INDEX ₂ The coordinates of the measuring point corresponding to the value (X _index1 , Y _index1 ) (X _index2 , Y _index2 ), determine whether the measuring point (X _index1 , Y _index1 ) (X _index2 , Y _index2 ) is in the overlapping range corresponding to the first measuring area If it is within, then judge one by one from (X _index1 , Y _index1 ), away from the measuring point (X _index2 , Y _index2 ), and close to the starting point or end point of the measuring line. Whether the adjacent measuring points are in the corresponding first measuring area Within the overlapping range of the first measurement area, until the last measurement point in the corresponding overlapping range of the first measurement area is found, record the index of this measurement point as the critical index THRESHOLD ₁ extracted from the measurement point data; also judge the self-measurement points one by one (X _index2 , Y _index2 ), away from the measuring point (X _index1 , Y _index1 ), close to the starting point or end point of the measuring line, the last measuring point index in the overlapping range corresponding to the first measuring area is used as the measuring point index. Another critical index THRESHOLD ₂ for point data extraction; take all the measurement point data of the survey line between the two critical indexes THRESHOLD ₁ and THRESHOLD ₂ as the measurement points of the survey line corresponding to the first survey area. Information of all measuring points within the overlapping range.

In an embodiment of the present invention, the user search conditions further include attribute constraints; the parameter data of all survey areas intersecting with the spatial search range are all survey areas that meet the attribute constraints and intersect the spatial search range the parameter data; the scoring according to the preset principle for all the survey areas intersecting with the spatial retrieval range includes: according to the preset principle, performing weight calculation on all the survey areas that meet the attribute constraints and intersect with the spatial retrieval scope, Sort the total weights, and take the test area with the largest total weight as the test area with the highest score.

In an embodiment of the present invention, before acquiring the parameter data of all survey areas intersecting with the spatial retrieval range, the method further includes: acquiring aeromagnetic original encoded data; processing the aeronautical geophysical data of one survey area after processing The standard exchange coded data is stored as an entity file in binary form, and the entity file of the survey area stores the field value information of the survey points of the survey area and the index information of each survey point.

In an embodiment of the present invention, after the data is extracted, the method further includes: feeding back the extracted data to the user through an illustration.

Compared with the prior art, according to the method for retrieving and extracting aeromagnetic measurement data of the present invention, the whole process can be automatically screened and the optimal single-coverage aeromagnetic measurement data of a designated area that meets the needs of users can be extracted, and the extraction efficiency is high and the extraction is accurate. .

During data extraction, data can be acquired according to the three data extraction accuracy levels of survey area, survey line, and survey point to meet the needs of different data scenarios, especially when data accuracy at the measurement point level is extracted, it can accurately obtain the user-specified area range. Optimal measurement data for single coverage, and solve the problem of repeated data coverage at the boundary of measurement data from different measurement areas, reducing the difficulty of user post-processing data.

In addition, the present invention defines a unified data content and data organization mode for the original aeromagnetic measurement field value and related measurement parameter values, and forms a unified standard exchange coding data format for the organization and storage of entity data. Later, data from different survey areas can be combined and stored based on this format. To the same data file, it reduces the difficulty of processing multiple survey area files for users, and provides data services to users based on the original aeromagnetic survey field values, which solves the problem of false magnetic anomalies in current data services based on gridded data. .

Description of drawings

FIG. 1 is a flowchart of a method for retrieving and extracting aeromagnetic survey data according to an embodiment of the present invention;

2 is a schematic diagram of a spatial range of entity data extraction using survey area retrieval and extraction according to an embodiment of the present invention;

3 is a schematic diagram of a spatial range of entity data extraction using survey line retrieval and extraction according to an embodiment of the present invention;

FIG. 4 is a schematic diagram of extracting a spatial range of entity data extracted using measuring point retrieval according to an embodiment of the present invention.

Detailed ways

The specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings, but it should be understood that the protection scope of the present invention is not limited by the specific embodiments.

Unless expressly stated otherwise, throughout the specification and claims, the term "comprising" or its conjugations such as "comprising" or "comprising" and the like will be understood to include the stated elements or components, and Other elements or other components are not excluded.

As shown in FIG. 1 , it is a flowchart of a method for retrieving and extracting aeromagnetic survey data according to a preferred embodiment of the present invention. A method for retrieving and extracting aeromagnetic survey data provided by this embodiment includes the following steps:

Step 1: Obtain user retrieval conditions, wherein the user retrieval conditions include spatial retrieval range and extraction precision, the spatial retrieval range included in the user retrieval conditions is the initial spatial retrieval range, and the extraction precision includes survey area level, survey line level, one of the measuring point levels.

Step 2: Acquire parameter data of all survey areas intersecting with the spatial retrieval range, and score all survey areas intersecting with the spatial retrieval scope according to a preset principle. Among them, the preset principle can be customized and extensible according to user needs.

In an implementation manner, the user retrieval conditions further include attribute constraints. For example, constraints on attribute information such as measurement scale, measurement time, data recording method, average flying height, and total measurement accuracy.

The parameter data of all survey areas intersecting with the spatial retrieval range are parameter data of all survey areas that meet the attribute constraints and intersect the spatial retrieval scope.

Correspondingly, scoring all survey areas intersecting with the spatial retrieval range according to the preset principle in step 2 may include: according to the preset principle, assigning weights to all survey areas that meet the attribute constraints and intersect with the spatial retrieval scope. Calculate, and take the test area with the largest weight as the test area with the highest score.

In this embodiment, the attribute constraints may include, but are not limited to, attribute information such as measurement scale, measurement time, data recording method, average flight height, and total measurement accuracy. The attribute constraints can be customized and expandable; the spatial range information in the retrieval conditions is based on 2000 National Geodetic Coordinate System. A survey area is a data storage unit, corresponding to an entity storage file. The data in the file is organized according to the survey line, and the survey line is stored one by one by the measuring points according to the spatial position, measurement time, and measuring point collection order, and the measuring point index is established.

According to the retrieval conditions (including the spatial retrieval scope and attribute information retrieval conditions), the spatial-attribute integrated retrieval is performed in the metadata information of the survey area, and the coding information of the survey area that meets the conditions and the attribute metadata information involved in the preset scoring principle are obtained. ;

The optimal retrieval principle is based on the measurement scale, average flight height, total magnetic measurement accuracy, data recording method and other attribute information selected by the user from the metadata of the survey area as the scoring items and assigns weights respectively to score and sort.

Step 3, in all the survey areas intersecting with the spatial retrieval range, obtain the survey area with the highest score as the first survey area;

Step 4: Acquire a spatial range in which the spatial range of the first survey area and the spatial retrieval range intersect and overlap, as the overlapped range corresponding to the first survey area. And record the measurement area code of the measurement area, the priority level of the measurement area, and the overlapping range corresponding to the first measurement area.

Step 5: Eliminate the overlapped part of the overlapping range corresponding to the first survey area from the spatial retrieval range to obtain an updated spatial retrieval range.

Step 6: In the other survey areas that intersect with the spatial retrieval range and have not acted as or have been used as the first survey area, take the survey area with the highest score as the updated first survey area, and repeat the acquisition of the first survey area. The spatial range in which the spatial range of the area intersects and overlaps with the updated spatial retrieval range is taken as the step of the overlapped range corresponding to the updated first survey area, until the initial spatial retrieval range is completely covered or all areas that intersect with the spatial retrieval range. The comparison of the measurement area is completed.

Specifically, a survey area code that satisfies the retrieval conditions and covers the optimal survey area in a single search area, and the overlapping range and entity data storage path information list corresponding to the first survey area can be obtained.

Step 7: Extract data according to the extraction precision and the overlapping range, wherein the overlapping range is the sum of the overlapping spatial ranges of all the spatial ranges of the first survey area and the spatial retrieval range.

The spatial retrieval scope includes the initial spatial retrieval scope and the spatial construction scope after each update.

Therefore, through the retrieval and extraction method of aeromagnetic measurement data provided in this embodiment, the whole process can be automatically screened and the optimal single coverage of aeromagnetic measurement data in a designated area that meets the needs of users can be extracted. The three data extraction accuracy levels of area, survey line and survey point are used for data acquisition to meet the needs of different data scenarios.

In one implementation, step 7 may include:

Step 7011, when the extraction accuracy is at the survey area level, use all the entity data of the first survey area corresponding to the overlapping range as the data extraction spatial range of the first survey area.

Step 7012: Extract data according to the sum of the spatial ranges of the data in the first survey area.

It should be noted that if the user wants to obtain all the entity data encoded by the optimal survey area, that is, the extraction accuracy is the survey area level, the entity data with the survey area as a unit can be directly extracted according to the storage location information in the survey area attribute information. Refer to Figure 2 for the corresponding metadata information of the survey area and the extracted data information.

In yet another implementation, step 7 may include:

Step 7021, when the extraction accuracy is non-measurement area level, obtain the information of all the survey lines in a first survey area corresponding to the overlapping range; Line start point coordinates, survey line end point coordinates, survey line start point index, survey line end point index and the corresponding survey area code list.

Step 7022: Calculate the data extraction spatial range of the first survey area according to the information of all the survey lines.

Step 7023: Acquire information of all survey lines in another first survey area corresponding to the overlap range, and repeat the step of calculating the data extraction space range of the first survey area until all the data corresponding to the overlap range are The calculation of the survey area is completed.

Step 7024: Extract data according to the sum of the spatial ranges of the data in the first survey area.

In an implementation manner, when the extraction accuracy is line level, step 7022 may include:

In acquiring the information of all the survey lines, the information of survey lines with intersections in the overlapping range corresponding to the first survey area may include: obtaining the code of the survey line that meets the conditions, the code of the survey area to which it belongs, and the starting point of the survey line. coordinates, coordinates of the end point of the survey line, index of the start point of the survey line, and index of the end point of the survey line; the information of all the survey points corresponding to the survey lines whose overlapping ranges corresponding to the first survey area have intersections are taken as the information of the first survey area. Data extraction spatial extent.

It should be noted that, if the user wishes to obtain the entity data of the survey lines that intersect with the overlapping range of the first survey region under all the optimal survey areas, that is, the extraction accuracy is the survey line level, the entity data in the metadata information of the survey area attributes will be stored. The location information and the survey line code information, the survey line start point index, and the survey line end index information in the survey line metadata information directly extract the entity data with the survey line as a unit. See Figure 3 for the extracted data information.

In an implementation manner, when the extraction accuracy is at the measuring point level, step 7022 may include:

In all the survey line information, if there is an intersection point between a survey line and the overlapping range corresponding to the first survey area, then calculate the intersection point of the survey line and the overlap range corresponding to the first survey area (X _crosst , Y _crosst );

When there are two intersection coordinates (X _cross1 , Y _cross1 ) and (X _cross2 , Y _cross2 ) in the overlapping range corresponding to the survey line and the first survey area, obtain the coordinates of the starting point (X _start, Y _start ) and the coordinates of the end point (X _end , Y _end ), the starting point index INDEX _start of the survey line, and the end point index INDEX _end , then according to the following formulas, the indexes of the survey points near the two intersection points are respectively estimated and recorded as INDEX ₁ and INDEX ₂ .

Wherein, INDEX ₁ and INDEX ₂ are integers.

According to the two index values of INDEX ₁ and INDEX ₂ , obtain the coordinates of the measuring point (X _index1 , Y _index1 ) (X _index2 , Y _index2 ) of the corresponding index on the survey line, and judge the measuring point (X _index1 , Y _index1 ) (X _index2 , Y _index2 ) is within the overlapping range corresponding to the corresponding first survey area;

If it is, then judge one by one from (X _index1 , Y _index1 ), away from the measuring point (X _index2 , Y _index2 ), and close to the starting point or end point of the measuring line. Whether the adjacent measuring points are in the corresponding overlap of the first measuring area Within the range, until the last measurement point in the overlapping range corresponding to the first measurement area is found, record the index of the measurement point as the critical index THRESHOLD ₁ for extraction of measurement point data; also judge the self-measurement points one by one (X _index2 , Y _index2 ), away from the measuring point (X _index1 , Y _index1 ), close to the starting point or end point of the measuring line, the last index of the measuring point that is in the overlapping range corresponding to the first measuring area is used as the measuring point data Extracted another critical index THRESHOLD ₂ ;

If not, then judge one by one from (X _index1 , Y _index1 ), whether the adjacent measuring points in the direction of the measuring point (X _index2 , Y _index2 ) are within the overlapping range corresponding to the first measuring area, until the first measuring point is found. For a measuring point in the overlapping range corresponding to the first measuring area, record the index of the measuring point as the critical _{index THRESHOLD 1} for extracting the measuring point data. X _index1 , Y _index1 ) direction continuous adjacent measuring points are judged one by one, find the index of the first measuring point in the overlapping range corresponding to the first measuring area, as another critical index THRESHOLD ₂ for measuring point data extraction .

It should be noted that since the actually obtained measuring line is only an approximate straight line, only the index of the measuring point near the intersection point can be calculated. When calculating the intersection point of the survey line and the overlapping range corresponding to the first survey area, the survey line space information in the survey line metadata is used. The spatial information of the survey line consists of the thinned survey points.

When the overlapping range corresponding to the survey line and the first survey area has only one intersection (X _cross1 , Y _cross1 ), a critical index THRESHOLD ₁ is obtained by calculating the above method, and at the same time, it is further judged whether the starting point or the end point of the survey line is in the first survey area. In the corresponding overlapping range, if it is within the range, the index of the start or end point is used as another critical index THRESHOLD ₂ .

According to the two critical indexes THRESHOLD ₁ and THRESHOLD ₂ , extract all the measuring point data of the survey line in the two index ranges, and use the information of all the survey points as the extraction space range of the survey line; specifically, it can be extracted, including longitude and latitude , ellipsoid height, height above ground (true height), original measured magnetic field value, compensated magnetic field value, normal geomagnetic field value, magnetic diurnal variation value, magnetic compensation value, ΔT magnetic anomaly value, etc. value or all values. Therefore, it is avoided to judge whether all the measuring points on the survey line are within the overlapping range corresponding to the first survey area, and the speed is fast.

Repeat the above steps until the calculation of the extraction space range of the survey lines in which the overlapping ranges corresponding to the first survey area have intersections in all survey line information is completed; The extraction space range of the survey line is taken as the data extraction space range of the first survey area. The extracted data information is shown in Figure 4.

In this way, when extracting the accuracy of measurement point-level data, the optimal measurement data for a single coverage within the user-specified area can be accurately obtained, and the problem of repeated data coverage at the boundary of measurement data from different measurement areas can be solved, reducing the need for users to Difficulty in post data processing.

Further, when performing data screening and extraction for national and other large-scale data applications, batch and business-oriented processing can be realized based on this method, and the advantages are more obvious.

In an implementation manner, before step 2, the method further includes:

Obtain the aeromagnetic original encoded data; the aeromagnetic original encoded data of a survey area is processed to form standard exchange encoded data for storage as an entity file, and the entity file of the survey area stores the survey points of the survey area The magnetic field value information and measuring point index. The survey line index consists of the start point index and the end point index of the survey line. The data of the survey line can be extracted according to the start point index and the end point index. The measuring points on the same survey line are stored in order, and all the measuring points in the survey area establish a unified measuring point index.

Specifically, the above steps can be implemented in the following ways:

Obtain the aeromagnetic survey data of a survey area (survey area); take the survey area as the basic unit to collect relevant metadata information. The metadata information includes two parts: spatial information and attribute information. The aeromagnetic data of one survey area corresponds to one piece of metadata information. . The spatial information of the survey area is mainly the boundary coordinate information of the coverage area of the aerial magnetic survey, or the information of the survey area, which is expressed by the 2000 National Geodetic Coordinate System and expressed by surface elements.

The aeromagnetic data here refers to the aeromagnetic original coded data obtained by the aeromagnetic survey system, which is processed by the aeronautical geophysical data to form the standard exchange coded data as the final data storage file. Store them one by one according to the surveying area -> surveying line -> surveying point mode. The entity data of a survey area is used as a storage file, which is organized and stored according to survey lines and survey points. Each measuring point includes but is not limited to latitude and longitude, ellipsoid height, height above ground (true height), original measured magnetic field value, compensated magnetic field value, normal geomagnetic field value, magnetic diurnal variation value, magnetic compensation value, ΔT magnetic field value outliers, etc.

The attribute information of the survey area includes but is not limited to the survey area code, the survey area name, the measurement time, the measurement scale, the average flying height, the total magnetic measurement accuracy, the survey line workload, the name of the measuring instrument, the sensitivity of the magnetometer, the data recording method, and the data storage. Path information, etc., attribute information can be expanded according to actual needs; data storage path information includes original encoded data storage path and standard exchange encoded data storage path, and subsequent entity data extraction is based on standard exchange encoded data. Collect the metadata information of each survey line, including spatial information and attribute information. One survey line corresponds to one metadata information; the spatial information of the survey line is expressed by line elements, and the 2000 national geodetic coordinate system is uniformly used, and the spatial information of the survey line is used. It is composed of the starting measuring point of the survey line, the coordinates of the ending measuring point and the thinned measuring point coordinates. The attribute information includes the survey line code, the code of the survey area, the flight number, the survey line direction, and the coordinates of the starting survey point. , end point coordinates, start point index, end point index, where the start point index and the end point index constitute the index information of the survey line in the entity file. Based on the spatial database, the metadata information of the survey area and the survey line metadata information are stored. Save the aeromagnetic data processed by standard exchange coding, and update the data storage path information in the property of the survey area.

Further, this embodiment defines a unified data content and data organization mode for the original aeromagnetic survey field value and related measurement parameter values, and forms a standard exchange coding data format for the organization and storage of entity data. Later, the data of different survey areas can be uniformly stored in the The same data file reduces the difficulty of processing multiple survey area files for users, and provides data services to users based on the original aeromagnetic survey field values, which solves the problem of false magnetic anomalies in current data services based on gridded data.

In an implementation manner, the standard exchange coded data is stored in binary coding form when the data is stored, and the data obtained by the data extraction is XYZ or CSV format data stored in plain code as the standard exchange coded data.

Therefore, compared with the original data, which is directly collected by the aircraft measurement system, professional software is required to identify it, and users cannot directly read it. This application is distributed in XYZ or CSV formats, which is convenient for users to use ordinary data processing software. It can be identified and easy to establish data index.

In an implementation manner, after the data is extracted, the method further includes: feeding back the extracted data to the user through an icon, which may specifically include:

Data map: drop the queried measuring points (reaction magnetic field data) on the graph according to the scale, and represent them as dots rendered by color. The color bar color (default) is assigned, which can be assigned linearly. Different colors correspond to different magnetic field values. At the same time, the magnetic field value is marked on the dots (marking the abnormal value of the magnetic field and the height above the ground (true height), and the minimum number of decimal places is 0.01nT). and 0.1m).

Work degree map: Expressed by surface elements, according to the final preferred survey area, the work degree map is formed by different colors, as a drawing, and the legend is marked with the metadata information of the survey area; at the same time, the query and retrieval range is drawn on it, Represented by unfilled faces.

As a result, a magnetic field data map represented only by the actual magnetic field value is output without any change in the data.

The retrieval section is explained in detail with the following example.

User retrieval conditions include the specified spatial retrieval range and attribute constraints (such as one or more of the measurement area attribute metadata, such as measurement scale, data recording method, average flight height, measurement time, measurement instrument model, etc.); The user determines the corresponding scoring weight, and the total weight function is defined as M=∑f(x _i ).

Survey area search:

Step 1. According to the spatial retrieval range and attribute constraints specified by the user, based on the metadata of the survey area, perform an integrated spatial-attribute search, and retrieve all the survey areas that satisfy the user's search constraints, and the total number of survey areas is recorded as N;

Step 2. According to the preset scoring condition information, score the data quality of the survey area. The scoring conditions can be customized according to user needs. For example, the scale can be agreed as 1:50,000 weight plus 40 points, 1:100,000 weight plus 35 points , 1:200,000 weight plus 30 points, 1:500,000 weight plus 25 points, 1:1,000,000 weight plus 20 points; data recording method is digital record, weight plus 30 points, analog record plus 15 points; the highest measured flight altitude The weight is 30 points; finally, the accumulated weight value of each condition is counted to obtain the total weight value.

For example, M=f(x ₁ )+f(x ₂ )+f(x ₃ ) can be set;

Among them, x ₁ refers to the weight value of the scale, x ₂ refers to the weight value of the data recording method, and x ₃ refers to the weight value of measuring the flight height. The overlap rate of the retrieval range, the larger the overlap rate, the higher the additional weight.

Obtain all survey areas that meet the requirements of the user's retrieval constraints and the corresponding total weights, sort the total weights from large to small, and determine the priority order. The larger the total weight, the higher the priority, and the final survey area priority is obtained. List List1, including survey area number, spatial range of survey area, total weight, priority;

Step 3. Obtain the survey area with the highest priority from the survey area priority list List1 as the first survey area, obtain the survey area code and spatial range information of the survey area, and delete the survey area information from List1 ;

Step 4. Calculate the spatial range where the spatial range of the survey area and the spatial retrieval range of the first survey area intersect and overlap, which is denoted as the overlapped range corresponding to the first survey area, and added to the survey area data extraction list List2, including the survey area code, Priority, overlapping range corresponding to the first survey area, entity data storage path;

Step 5. Calculate the part that intersects with the survey area from the spatial retrieval range to obtain a new spatial retrieval range;

Step 6. Repeat step 3 to obtain the survey area with the highest priority in List1 and compare it with the new spatial search range obtained in step 5. If there is an overlapping area between the two, repeat steps 4-5; if there is no overlapping area between the two , then repeat step 3 to continue to select the next highest priority survey area from List1 as the updated first survey area for comparison;

Step 7, and so on, until the spatial retrieval range originally set in the user constraints is completely covered or all the survey areas in List1 are compared, and finally the optimal survey area data extraction list List2 that satisfies the user constraints is obtained, and the number of survey areas is obtained. Denoted as M, M≤N.

Line search:

Step 1, obtain the survey area code of a survey area, the corresponding overlapping range of the survey area, and the entity data storage path from the survey area data extraction list List2;

Step 2, according to the survey area code and the overlapping range information corresponding to the survey area, perform a space-attribute integrated search in the survey line metadata information, and obtain all survey line information that meets the conditions in the survey area, including the survey line code, the survey line to which it belongs. Area code, survey line space information, survey line start point coordinates, survey line end point coordinates, survey line start point index, survey line end point index and other information;

Step 3, obtains a survey line space information of survey line and calculates two intersections (X _cross1 , Y _cross1 ), (X _cross2 , Y _cross2 ) of its overlapping range corresponding to the first survey area;

Step 4, according to the entity data storage path information and the start and end index information of the survey line in the entity file obtained in step 2, extract all survey point information of all survey lines with intersections in the overlapping range corresponding to the first survey area.

Survey point retrieval: Repeat steps 1-3 of survey line retrieval.

Step 4: According to the coordinates of the start point of the survey line (x _start , y _start ), the coordinates of the end point of the survey line (x _end , y _end ), the index of the start point of the survey line INDEX _start , and the index of the end point of the survey line INDEX _end obtained in step 2, the The intersection (x ₁ , y ₁ ) (x ₂ , y ₂ ) of a survey line and the overlapping range of the first survey area, calculate two index values INDEX ₁ , INDEX ₂ ;

Step 5: Obtain the entity data storage path information in the survey area, and obtain the corresponding survey point coordinates (X _index1 , Y _index1 ) (X _index2 , Y _index2 ) according to the two indexes of INDEX ₁ and INDEX ₂ in the entity data file, and interpret Whether the measuring point (X _index1 , Y _index1 ) (X _index2 , Y _index2 ) is within the overlapped range corresponding to the corresponding first measuring area; if so, further judge whether the measuring point (X _index1 , Y _index1 ) is far away from the measuring point (X _index2 , Y _index2 ) and whether the consecutive adjacent measuring points close to the starting point or end point of the survey line are within the overlapping range corresponding to the first survey area, until the last measuring point within the overlapping range corresponding to the first survey area is found, Record the index of the measuring point as the critical index THRESHOLD ₁ for the extraction of measuring point data; also judge one by one from the measuring point (X _index2 , Y _index2 ), away from the measuring point (X _index1 , Y _index1 ), close to the starting point or end point of the measuring line The last index of the consecutive adjacent measuring points in the overlapping range corresponding to the first measuring area is used as another critical index THRESHOLD ₂ for extracting the measuring point data.

If not, then judge one by one from (X _index1 , Y _index1 ), whether the adjacent measuring points in the direction of the measuring point (X _index2 , Y _index2 ) are within the overlapping range corresponding to the first measuring area, until the first measuring point is found. For a measuring point in the overlapping range corresponding to the first measuring area, record the index of the measuring point as the critical _{index THRESHOLD 1} for extracting the measuring point data. X _index1 , Y _index1 ) direction continuous adjacent measuring points are judged one by one, find the index of the first measuring point in the overlapping range corresponding to the first measuring area, as another critical index THRESHOLD ₂ for measuring point data extraction .

When the overlapping range corresponding to the survey line and the first survey area has only one intersection point, a critical index is calculated according to the above method, and at the same time, it is further judged whether the start or end point of the survey line is within the overlap range corresponding to the first survey area. The index of the start or end point is then used as another critical index.

Step 6, directly extract the magnetic field value record information of all measuring points in the two critical index ranges in the survey line according to the two critical indexes obtained in step 5;

Step 7, repeat steps 3-6 until the entity data of all survey lines that meet the conditions in the survey area have been extracted;

Step 8, repeat step 1, continue to obtain the survey area code of the next survey area, the overlap range corresponding to the first survey area, and the entity data storage path from the survey area data extraction list List2, and perform steps 2-7;

By analogy, until all the survey area data in List2 is extracted, the data is provided to the user using standard exchange coding data (XYZ format or CSV format), and data graphics can be provided at the same time.

As will be appreciated by those skilled in the art, the embodiments of the present application may be provided as a method, a system, or a computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.

The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the present application. It will be understood that each flow and/or block in the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to the processor of a general purpose computer, special purpose computer, embedded processor or other programmable data processing device to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing device produce Means for implementing the functions specified in a flow or flow of a flowchart and/or a block or blocks of a block diagram.

These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory result in an article of manufacture comprising instruction means, the instructions The apparatus implements the functions specified in the flow or flow of the flowcharts and/or the block or blocks of the block diagrams.

These computer program instructions can also be loaded on a computer or other programmable data processing device to cause a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process such that The instructions provide steps for implementing the functions specified in the flow or blocks of the flowcharts and/or the block or blocks of the block diagrams.

The foregoing descriptions of specific exemplary embodiments of the present invention have been presented for purposes of illustration and description. These descriptions are not intended to limit the invention to the precise form disclosed, and obviously many changes and modifications are possible in light of the above teachings. The exemplary embodiments were chosen and described for the purpose of explaining certain principles of the invention and their practical applications, to thereby enable one skilled in the art to make and utilize various exemplary embodiments and various different aspects of the invention. Choose and change. The scope of the invention is intended to be defined by the claims and their equivalents.

Claims (9)

1. a retrieval and extraction method of aeromagnetic survey data, is characterized in that, described retrieval and extraction method comprises: Obtain user retrieval conditions, wherein the user retrieval conditions include spatial retrieval scope and extraction precision, the spatial retrieval scope included in the user retrieval conditions is the initial spatial retrieval scope, and the extraction precision includes survey area level, survey line level, and survey line level. one of the point classes; Obtain the parameter data of all the survey areas intersecting with the spatial retrieval range, and score all the survey areas intersecting with the spatial retrieval scope according to the preset principle; Among all the survey areas intersecting with the spatial retrieval range, the survey area with the highest score is obtained as the first survey area; Acquiring the spatial range where the spatial range of the first survey area and the spatial retrieval range intersect and overlap, as the overlap range corresponding to the first survey area; Eliminate the overlapped part of the overlapping range corresponding to the first survey area from the spatial retrieval range to obtain the updated spatial retrieval range; Among the other survey areas that intersect with the spatial retrieval range and have not been or have been used as the first survey area, take the survey area with the highest score as the updated first survey area, and repeatedly acquire the space of the first survey area The spatial range where the range intersects and overlaps with the updated spatial retrieval range is taken as the step of the overlapping range corresponding to the updated first survey area, until the initial spatial retrieval range is completely covered or compared with all survey areas intersecting with the spatial retrieval range complete; Data is extracted according to the extraction precision and the overlapping range, wherein the overlapping range is the sum of the overlapping spatial ranges of all the spatial ranges of the first survey area and the spatial retrieval range. 2. The retrieval and extraction method according to claim 1, wherein, according to the extraction precision and the overlapping range, extracting data comprises: When the extraction accuracy is the survey area level, all the entity data of the first survey area corresponding to the overlapping range is taken as the data extraction spatial range of the first survey area; Extract the data according to the sum of the spatial ranges of the data of the first survey area. 3. The retrieval and extraction method according to claim 1, wherein, according to the extraction precision and the overlapping range, extracting data comprises: When the extraction accuracy is at the non-measurement area level, obtain information of all survey lines in a first survey area corresponding to the overlapped range; According to the information of all survey lines in the first survey area, calculate the data extraction space range of the first survey area; Obtain the information of all survey lines in another first survey area corresponding to the overlapping range, and repeat the steps of calculating the data extraction space range of the first survey area until all survey areas corresponding to the overlap range are calculated Finish; Extract the data according to the sum of the spatial ranges of the data of the first survey area. 4. The retrieval and extraction method according to claim 3, characterized in that, when the extraction accuracy is at the survey line level, according to the information of all survey lines in the first survey area, calculate the value of the first survey area. Data extraction spatial extent, including: Acquiring information of all surveying points corresponding to surveying lines whose overlapping ranges corresponding to the first surveying region have intersections among the information of all surveying lines in the first surveying region; The information of all survey points corresponding to all survey lines whose overlapping ranges corresponding to the first survey area have intersections are taken as the data extraction spatial range of the first survey area. 5 . The retrieval and extraction method according to claim 3 , wherein when the extraction accuracy is at the measuring point level, according to all the survey line information, the data extraction space range of the first survey area is calculated, comprising: 5 . : In all survey line information in the first survey area, if there is an intersection point between a survey line and the overlap range corresponding to the first survey area, calculate the overlap range corresponding to the survey line and the first survey area the intersection point; Judging the measuring points on the survey line, obtaining information of all the measuring points in the overlapping range corresponding to the first survey area of the measuring points of the survey line, and using the information of all the measuring points as the The extraction space range of the survey line; The above steps are repeated until the calculation of the extraction space range of the survey line where the overlapping range corresponding to the first survey area in all survey line information has an intersection point is completed; The extraction space range of all survey lines with intersections of overlapping ranges corresponding to the first survey area is taken as the data extraction space range of the first survey area. 6. The retrieval and extraction method according to claim 5, wherein the acquisition of the information of all the measuring points in the overlapping range corresponding to the measuring points of the survey line in the first survey area comprises: Calculate the index marks INDEX ₁ and INDEX ₂ of the survey points near the two intersection points of the survey lines respectively according to the following formulas;

Wherein, INDEX ₁ and INDEX ₂ are integers, (X _cross1 , Y _cross1 ), (X _cross2 , Y _cross2 ) are two intersection coordinates in the overlapping range corresponding to the survey line and the first survey area, (X _start , Y _start ) is the start point coordinate of the survey line, (X _end , Y _end ) is the end point coordinate, INDEX _start is the start point index of the survey line, and INDEX _end is the end point index of the survey line; Obtain the coordinates (X _index1 , Y _index1 ) (X _index2 , Y _index2 ) of the measuring points corresponding to the two index values of the INDEX ₁ and INDEX ₂ on the survey line, and determine the measuring points (X _index1 , Y _index1 ) (X Whether _index2 , Y _index2 ) is within the overlapping range corresponding to the first survey area; If it is, then judge one by one from (X _index1 , Y _index1 ), away from the measuring point (X _index2 , Y _index2 ), and close to the starting point or end point of the measuring line. Whether the adjacent measuring points are in the corresponding overlap of the first measuring area Within the range, until the last measurement point in the overlapping range corresponding to the first measurement area is found, record the index of this measurement point, as the critical index THRESHOLD ₁ of measurement point data extraction; Similarly, judge one by one from the measuring point (X _index2 , Y _index2 ), away from the measuring point (X _index1 , Y _index1 ), and close to the starting point or end point of the measuring line. The last adjacent measuring point is in the overlap corresponding to the first measuring area. The index of measuring points within the range, as another critical index THRESHOLD ₂ for measuring point data extraction; All measuring point data of the survey line between two critical indices THRESHOLD ₁ and THRESHOLD ₂ are used as information of all survey points of the survey line in the overlapping range corresponding to the first survey area. 7. The retrieval and extraction method according to claim 1, wherein the user retrieval conditions further comprise attribute constraints; The parameter data of all the survey areas intersecting with the spatial retrieval range are the parameter data of all survey areas that meet the attribute constraints and intersect with the spatial retrieval scope; The scoring of all survey areas intersecting with the spatial retrieval range according to the preset principle includes: According to the preset principle, the weights are calculated for all the survey areas that meet the attribute constraints and intersect with the spatial retrieval range, the total weights are sorted, and the survey area with the largest total weight value is regarded as the survey area with the highest score. 8. The retrieval and extraction method according to claim 1, characterized in that, before acquiring the parameter data of all survey areas intersecting with the spatial retrieval range, the method further comprises: Obtain aeromagnetic raw encoded data; A unified standard exchange coding data format is formed after the aeronautical geophysical data processing of a survey area is stored as an entity file in binary form, and the entity file of the survey area stores the survey point field value information of the survey area and Index information for each measurement point. 9 . The retrieval and extraction method according to claim 1 , wherein after extracting the data, the method further comprises: feeding back the extracted data to the user through an icon. 10 .

Images