CN114528454A - Method for retrieving and extracting aeromagnetic measurement data - Google Patents
Method for retrieving and extracting aeromagnetic measurement data Download PDFInfo
- Publication number
- CN114528454A CN114528454A CN202210150948.XA CN202210150948A CN114528454A CN 114528454 A CN114528454 A CN 114528454A CN 202210150948 A CN202210150948 A CN 202210150948A CN 114528454 A CN114528454 A CN 114528454A
- Authority
- CN
- China
- Prior art keywords
- measuring
- range
- area
- data
- measurement
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/907—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually
- G06F16/909—Retrieval characterised by using metadata, e.g. metadata not derived from the content or metadata generated manually using geographical or spatial information, e.g. location
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/90—Details of database functions independent of the retrieved data types
- G06F16/901—Indexing; Data structures therefor; Storage structures
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Software Systems (AREA)
- Library & Information Science (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
The invention discloses a retrieval and extraction method of aeromagnetic measurement data, which comprises the following steps: acquiring user retrieval conditions and parameter data of all measurement areas intersected with the space retrieval range, and scoring all measurement areas intersected with the space retrieval range according to a preset principle; acquiring the measuring region with the highest score from all the measuring regions intersected with the space retrieval range to serve as a first measuring region; acquiring a space range of the first measuring area, which is intersected and overlapped with the space retrieval range, as an overlapping range corresponding to the first measuring area; removing a part which is intersected and covered with the overlapping range corresponding to the first measuring area from the space retrieval range to obtain an updated space retrieval range; and extracting data according to the extraction precision and the overlapping range. The method for retrieving and extracting the aeromagnetic measurement data can realize full-automatic data retrieval and extraction, and has high extraction efficiency and accurate extraction.
Description
Technical Field
The invention relates to the field of aeromagnetic measurement data application and service, in particular to a retrieval and extraction method of aeromagnetic measurement data.
Background
Aeromagnetic measurement, also called aeromagnetic measurement or aeromagnetic exploration, called aeromagnetic measurement for short, is a geophysical method for loading an aeromagnetic magnetometer and auxiliary equipment matched with the aeromagnetic magnetometer on an aircraft and measuring the intensity or gradient of a geomagnetic field according to a preset measuring line and height above a measuring area (measuring area).
In recent years, with the increasing demand for aeromagnetic survey, new aeromagnetic survey data in China continuously increase. At present, the land area full coverage of China is basically realized, and aeromagnetic data repeated coverage conditions of different scales, heights, measurement accuracy, recording modes and the like exist in the same area for many times. These data are usually stored in units of measurement areas, and in the application of future data, a user is required to screen information such as a measurement scale, a height, measurement accuracy, a recording mode, and a research area range, so as to obtain optimal measurement data meeting the requirements.
In addition, when data is provided, because a user cannot directly identify original measurement data, corresponding gridding data is often generated from the original measurement data, when a user demand range covers a plurality of measurement areas, gridding data of different measurement area units are respectively cut according to user demand range information when the data is provided for the user, and data overlapping areas exist at boundaries among the different measurement areas.
Based on this, the inventor of the application finds that the current aeromagnetic measurement data acquisition process mainly adopts manual screening, and the process is complex and tedious; the data is provided mainly by gridding data, is not original measured magnetic field value data obtained based on aeromagnetic measurement, and is greatly influenced by data processing; when a plurality of measuring areas need to be provided at the same time, the repeated workload is large, and the repeated coverage data at the boundaries of different measuring 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 skilled in the art.
Disclosure of Invention
The invention aims to provide a method for retrieving and extracting aeromagnetic measurement data, which can automatically realize retrieval-level extraction of the aeromagnetic measurement data based on original measurement data.
In order to achieve the above object, the present invention provides a method for retrieving and extracting aeromagnetic measurement data, wherein the method for retrieving and extracting aeromagnetic measurement data comprises: acquiring a user retrieval condition, wherein the user retrieval condition comprises a space retrieval range and extraction precision, the space retrieval range included by the user retrieval condition is an initial space retrieval range, and the extraction precision comprises one of a measurement area level, a measurement line level and a measurement point level; acquiring parameter data of all measurement areas intersected with the space retrieval range, and scoring all measurement areas intersected with the space retrieval range according to a preset principle; acquiring the measuring region with the highest score from all the measuring regions intersected with the space retrieval range as a first measuring region; acquiring a space range of the first measuring area, which is intersected and overlapped with the space retrieval range, as an overlapping range corresponding to the first measuring area; removing a part which is intersected and covered with the overlapping range corresponding to the first measuring area from the space retrieval range to obtain an updated space retrieval range; in the other measurement areas which are intersected with the space retrieval range and are not used as or are used as the first measurement area, taking the measurement area with the highest score as an updated first measurement area, repeatedly obtaining the space range of the first measurement area, which is intersected and overlapped with the updated space retrieval range, and taking the space range as the overlapped range corresponding to the updated first measurement area until the initial space retrieval range is completely covered or all the measurement areas intersected with the space retrieval range are completely compared; and extracting data according to the extraction precision and the overlapping range, wherein the overlapping range is the sum of the spatial ranges of all the first measuring areas and the spatial retrieval range which are intersected and overlapped.
In one embodiment of the present invention, extracting data based on the extraction accuracy and the overlap range includes: when the extraction precision is in a measurement area level, taking all entity data of the first measurement area corresponding to the overlapping range as a data extraction space range of the first measurement area; and extracting data according to the sum of the data extraction space ranges of the first measuring area.
In one embodiment of the present invention, extracting data based on the extraction accuracy and the overlap range includes: when the extraction precision is in a non-measurement-area level, acquiring information of all measurement lines in a first measurement area corresponding to the overlapping range; calculating the data extraction spatial range of the first measuring area according to the information of all measuring lines; acquiring information of all measuring lines in another first measuring area corresponding to the overlapping range, and repeating the step of calculating the data extraction spatial range of the first measuring area until the calculation of all measuring areas corresponding to the overlapping range is completed; and extracting data according to the sum of the data extraction space ranges of the first measuring area.
In an embodiment of the present invention, when the extraction accuracy is in a line measurement level, calculating a data extraction spatial range of the first measurement area according to information of all the lines includes: acquiring information of all measuring points corresponding to measuring lines with intersection points in the overlapping range corresponding to the first measuring area in all measuring line information in the first measuring area; and taking the information of all measuring points corresponding to the measuring lines with the intersection points in the overlapping range corresponding to the first measuring area as the data extraction space range of the first measuring area.
In an embodiment of the present invention, when the extraction accuracy is at a measurement point level, calculating a data extraction spatial range of the first measurement area from the all-line information includes: in all the measuring line information in the first measuring area, if an intersection point exists in the overlapping range corresponding to one measuring line and the first measuring area, calculating the intersection point of the overlapping range corresponding to the measuring line and the first measuring area; judging the measuring points on the measuring line, acquiring the information of all measuring points of the measuring line in the overlapping range corresponding to the first measuring area, and taking the information of all measuring points as the extraction space range of the measuring line; repeating the steps until the calculation of the extraction space range of the measuring line with the intersection point in the overlapping range corresponding to the first measuring area in all the measuring line information is completed; and taking the extraction space range of all the measuring lines with intersection points in the overlapping range corresponding to the first measuring area as the data extraction space range of the first measuring area.
In an embodiment of the present invention, the acquiring information of all the measurement points of the measurement line in the overlapping range corresponding to the first measurement area includes: respectively calculating the INDEX of the measuring points near the two intersection points of the measuring lines according to the following formula1、INDEX2；
Wherein, INDEX1、INDEX2Taking an integer, (X)cross1，Ycross1)、(Xcross2，Ycross2) The overlapping range of the measuring line and the first measuring area has two intersection point coordinates (X)start，Ystatr) Is the coordinate of the starting point of the survey line, (X)end，Yend) Is the end point coordinate, INDEXstartIs the INDEX of the starting point of the survey line, INDEXendIs the index of the end point of the survey line; obtaining the on-line and the INDEX1、INDEX2Measuring point coordinate (X) corresponding to two index valuesindex1，Yindex1)(Xindex2，Yindex2) Judgment measuring point (X)index1，Yindex1)(Xindex2，Yindex2) Whether the first measurement area is in the corresponding overlapping range of the first measurement area or not; if yes, then judge one by one from (X)index1，Yindex1) From, away from the measuring point (X)index2，Yindex2) Whether the continuous adjacent measuring points close to the starting point or the end point of the measuring line are in the overlapping range corresponding to the first measuring area until the last measuring point in the overlapping range corresponding to the first measuring area is found, recording the index of the measuring point, and taking the index as the critical index THRESHOLD extracted from the measuring point data1(ii) a Also judge one by oneBreak self-test point (X)index2，Yindex2) From, away from the measuring point (X)index1，Yindex1) The last measuring point index in the overlapping range corresponding to the first measuring area in the continuous adjacent measuring points close to the starting point or the end point of the measuring line is used as another critical index THRESHOLD extracted from the measuring point data2(ii) a Placing the line under two critical indices THRESHOLD1、THRESHOLD2And all the measuring point data in between are used as the information of all the measuring points of the measuring line in the corresponding overlapping range of the first measuring area.
In an embodiment of the present invention, the user search condition further includes an attribute constraint condition; the parameter data of all the measurement areas intersected with the space retrieval range are the parameter data of all the measurement areas which meet the attribute constraint condition and are intersected with the space retrieval range; the scoring of all measurement areas intersecting the spatial retrieval range according to a preset principle comprises: and calculating the weights of all the measurement areas which meet the attribute constraint conditions and are intersected with the space retrieval range according to a preset principle, sequencing the total weights, and taking the measurement area with the maximum total weight as the measurement area with the highest score.
In an embodiment of the present invention, before acquiring parameter data of all measurement areas intersecting with the spatial search range, the method further includes: acquiring aeromagnetic original coding data; standard exchange coding data formed by processing aviation geophysical data of a measuring area are stored as an entity file in a binary form, and measuring point field value information of the measuring area and index information of each measuring point are stored in the entity file of the measuring area.
In an embodiment of the present invention, after the extracting data, the method further includes: the extracted data is fed back to the user by means of a graphical representation.
Compared with the prior art, the method for retrieving and extracting the aeromagnetic measurement data can realize full-flow automatic screening and extract the optimal aeromagnetic measurement data which is covered by the specified area and meets the user requirements once, and has high extraction efficiency and accurate extraction.
When data is extracted, data can be acquired according to three data extraction precision levels of a measuring area, a measuring line and a measuring point, the requirements of different data scenes are met, particularly, when data at the measuring point level is extracted precisely, the optimal measuring data covered once in a user specified area range can be accurately acquired, the problem that data are repeatedly covered at the boundary of the measuring data from different measuring areas is solved, and the difficulty in later-stage data processing of a user is reduced.
In addition, the invention defines uniform data content and data organization mode aiming at the original aeromagnetic measurement field value and related measurement parameter value, forms a uniform standard exchange coding data format to carry out organization and storage of entity data, and different later-stage measurement area data can be merged and stored into the same data file based on the format, thereby reducing the difficulty of processing a plurality of measurement area files of a user, simultaneously providing data service for the user based on the original aeromagnetic measurement field value, and solving the problem that false magnetic anomaly exists in the current data service based on gridded data.
Drawings
FIG. 1 is a flow chart of a method for retrieving and extracting aeromagnetic measurement data according to an embodiment of the invention;
FIG. 2 is a schematic diagram of an entity data extraction spatial range using region search extraction according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a spatial range of entity data extraction using line search extraction according to an embodiment of the present invention;
FIG. 4 is a schematic diagram of an entity data extraction space range extracted by using a survey point search according to an embodiment of the present invention.
Detailed Description
The following detailed description of the present invention is provided in conjunction with the accompanying drawings, but it should be understood that the scope of the present invention is not limited to the specific embodiments.
Throughout the specification and claims, unless explicitly stated otherwise, the word "comprise", or variations such as "comprises" or "comprising", will be understood to imply the inclusion of a stated element or component but not the exclusion of any other element or component.
Fig. 1 is a flowchart of a method for retrieving and extracting aeromagnetic measurement data according to a preferred embodiment of the present invention. The method for retrieving and extracting aeromagnetic measurement data provided by the embodiment comprises the following steps:
step 1, obtaining a user retrieval condition, wherein the user retrieval condition comprises a space retrieval range and extraction precision, the space retrieval range comprised by the user retrieval condition is an initial space retrieval range, and the extraction precision comprises one of a measurement area level, a measurement line level and a measurement point level.
And 2, acquiring parameter data of all measurement areas intersected with the space retrieval range, and scoring all measurement areas intersected with the space retrieval range according to a preset principle. The preset principle can be customized according to the user requirement and can be expanded.
In one implementation, the user search condition further includes an attribute constraint condition. For example, constraints on attribute information such as a measurement scale, measurement time, a data recording method, an average flying height, and measurement total accuracy.
And the parameter data of all the measurement areas intersected with the space retrieval range are the parameter data of all the measurement areas which meet the attribute constraint condition and are intersected with the space retrieval range.
Correspondingly, the scoring of all the measurement areas intersecting the spatial search range according to the preset principle in the step 2 may include: and calculating the weight of all the measurement areas which meet the attribute constraint condition and are intersected with the space retrieval range according to a preset principle, and taking the measurement area with the maximum weight as the measurement area with the highest score.
In this embodiment, the attribute constraints may include, but are not limited to, attribute information such as a measurement scale, measurement time, a data recording manner, an average flying height, and total measurement accuracy, and the attribute constraints may be customized and extensible; the spatial range information in the search condition is based on the 2000 national geodetic coordinate system. One measuring area is a data storage unit and corresponds to an entity storage file, data in the file are organized according to measuring lines, the measuring lines are stored one by measuring points according to spatial positions, measuring time and a measuring point acquisition sequence, and measuring point indexes are established.
Performing space-attribute integrated retrieval in the region-measuring level metadata information according to retrieval conditions (including a space retrieval range and attribute information retrieval conditions), and acquiring region-measuring coding information meeting the conditions and attribute metadata information related in a preset scoring principle;
and according to the optimal retrieval principle, scoring and sequencing are carried out according to attribute information such as a measurement scale, average flight height, total magnetic measurement precision, a data recording mode and the like selected by a user from the metadata of the measurement area as scoring items and respectively given weights.
Step 3, acquiring the measuring area with the highest score from all the measuring areas intersected with the space retrieval range as a first measuring area;
and 4, acquiring a space range of the first measuring area, which is intersected and overlapped with the space retrieval range, as an overlapping range corresponding to the first measuring area. And recording the measuring region code, the measuring region priority level and the corresponding overlapping range of the first measuring region of the measuring region.
And 5, eliminating the part which is intersected and covered with the overlapping range corresponding to the first measuring area from the space retrieval range to obtain the updated space retrieval range.
And 6, in the other measurement areas which are intersected with the space retrieval range and are not used as or are used as the first measurement area, taking the measurement area with the highest score as an updated first measurement area, repeatedly acquiring the space range of the first measurement area, which is intersected and overlapped with the updated space retrieval range, and taking the space range as the overlapped range corresponding to the updated first measurement area until the initial space retrieval range is completely covered or the comparison of all the measurement areas intersected with the space retrieval range is completed.
Specifically, the area code and the overlap range corresponding to the first area and the entity data storage path information list, which satisfy the search condition and cover the optimal area of the search area once, may be obtained.
And 7, extracting data according to the extraction precision and the overlapping range, wherein the overlapping range is the sum of the spatial ranges of all the first measurement areas and the spatial search range which are intersected and overlapped.
The space search range includes an initial space search range and a space construction range updated each time.
Therefore, by the method for retrieving and extracting the aeromagnetic measurement data, the optimal aeromagnetic measurement data covered by the designated area for one time and meeting the user requirements can be automatically screened and extracted in the whole process, and during data extraction, data can be acquired according to three data extraction precision levels of a measuring area, a measuring line and a measuring point, so that different data scene requirements are met.
In one implementation, step 7 may include:
step 7011, when the extraction precision is a measurement area level, using all the entity data of the first measurement area corresponding to the overlap range as a data extraction spatial range of the first measurement area.
Step 7012, extracting data according to the sum of the spatial ranges of the data extraction of the first measurement area.
It should be noted that, if a user wants to obtain all entity data of the optimal region code, that is, the extraction precision is the region level, the entity data in the region unit and the corresponding region metadata information may be directly extracted according to the storage location information in the region attribute information, and the extracted data information refers to fig. 2.
In yet another implementation, step 7 may include:
step 7021, when the extraction precision is at a non-measurement-area level, obtaining information of all measurement lines in a first measurement area corresponding to the overlap range; the information of the measuring line can include corresponding measuring line codes, measuring line starting point coordinates, measuring line end point coordinates, measuring line starting point indexes, measuring line end point indexes and corresponding measuring area code lists.
Step 7022, calculating a data extraction spatial range of the first measurement area according to the information of all the measurement lines.
Step 7023, obtaining information of all measurement lines in another first measurement area corresponding to the overlapping range, and repeating the step of calculating the data extraction spatial range of the first measurement area until the calculation of all measurement areas corresponding to the overlapping range is completed.
Step 7024, extracting data according to the sum of the spatial ranges of the data extraction of the first measurement area.
In one implementation, when the extraction precision is in a line level, step 7022 may include:
obtaining, from the information of all the survey lines, survey line information in which an intersection exists in an overlap range corresponding to the first survey area, specifically, the obtaining may include: obtaining a measuring line code, a belonging measuring area code, a measuring line starting point coordinate, a measuring line end point coordinate, a measuring line starting point index and a measuring line end point index which accord with conditions; and taking the information of all measuring points corresponding to the measuring lines with the intersection points in the overlapping range corresponding to the first measuring area as the data extraction space range of the first measuring area.
It should be noted that, if a user wishes to obtain the line measurement entity data intersecting the overlapping range of the first measurement area under all the optimal measurement areas, that is, the extraction accuracy is the line measurement level, the entity data in units of line measurement are directly extracted according to the entity data storage location information in the measurement area attribute metadata information and the line measurement coding information, the line measurement starting point index and the line measurement end point index information in the line measurement metadata information, and the extracted data information refers to fig. 3.
In one implementation, when the extraction precision is at a measurement point level, step 7022 may include:
in the total survey line information, if there is an intersection in the overlapping range corresponding to one survey line and the first survey area, calculating an intersection (X) of the overlapping range corresponding to the survey line and the first survey areacrosst，Ycrosst)；
When the overlapping range of the measuring line and the first measuring area has two intersection point coordinates (X)cross1，Ycross1)、(Xcross2，Ycross2) Obtaining the coordinate (X) of the starting point of the survey linestart，Ystart) Endpoint coordinate (X)end，Yend) INDEX of starting point of survey line INDEXstartEndpoint INDEX INDEXendThen, thenRespectively estimating the INDEX of the measuring point near the two intersection points according to the following formula and recording the INDEX as INDEX1、INDEX2。
Wherein, INDEX1、INDEX2And taking an integer.
According to INDEX1、INDEX2Two index values obtain the coordinate (X) of the corresponding index of the measuring point on the measuring lineindex1，Yindex1)(Xindex2，Yindex2) Judgment measuring Point (X)index1，Yindex1)(Xindex2，Yindex2) Whether the first measurement area is in the corresponding overlapping range of the corresponding first measurement area;
if yes, then judge one by one from (X)index1，Yindex1) From a far measuring point (X)index2，Yindex2) Whether the continuous adjacent measuring points close to the starting point or the end point of the measuring line are in the overlapping range corresponding to the first measuring area until the last measuring point in the overlapping range corresponding to the first measuring area is found, recording the index of the measuring point, and taking the index as the critical index THRESHOLD extracted from the measuring point data1(ii) a The measured points (X) are also judged one by oneindex2，Yindex2) From, away from the measuring point (X)index1，Yindex1) The last measuring point index in the overlapping range corresponding to the first measuring area in the continuous adjacent measuring points close to the starting point or the end point of the measuring line is used as another critical index THRESHOLD extracted from the measuring point data2；
If not, then one by one is judged from (X)index1，Yindex1) From, near to the measurement point (X)index2，Yindex2) Whether the measuring points which are continuously adjacent in the direction are in the overlapping range corresponding to the first measuring area or not is judged until the first measuring point which is in the overlapping range corresponding to the first measuring area is foundMeasuring point, recording index of the measuring point as critical index THRESHOLD extracted from measuring point data1The same for the self-test point (X)index2，Yindex2) Start close to measuring point (X)index1，Yindex1) The continuous adjacent measuring points in the direction are judged one by one, the index of the first measuring point in the overlapping range corresponding to the first measuring area is found and is used as another critical index THRESHOLD extracted from the measuring point data2。
It should be noted that since the actually obtained measurement line is only an approximate straight line, only the measurement point index near the intersection can be calculated. The intersection points of the overlapping ranges of the survey lines and the first survey area are calculated using the survey line spatial information in the survey line metadata. The measuring line space information consists of measuring points after thinning.
When the measuring line and the overlapping range corresponding to the first measuring area only have one intersection (X)cross1，Ycross1) Then, a critical index THRESHOLD is obtained by calculation according to the method1And simultaneously further judging whether the starting point or the end point of the measuring line is in the overlapping range corresponding to the first measuring area, if so, taking the index of the starting point or the end point as another critical index THRESHOLD2。
THRESHOLD according to two critical indexes1、THRESHOLD2Extracting all measuring point data of the measuring line in two index ranges, and taking the information of all measuring points as an extraction space range of the measuring line; specifically, a plurality of combined values or all values of the information including latitude and longitude, ellipsoidal height, height above ground (true height), originally measured magnetic field value, compensated magnetic field value, normal magnetic field value, magnetic daily variation value, magnetic compensation value, delta T magnetic abnormal value, and the like can be extracted. Therefore, whether all the measuring points on the measuring line are in the overlapping range corresponding to the first measuring area or not is avoided, and the speed is high.
Repeating the steps until the calculation of the extraction space range of the measuring line with the intersection point in the overlapping range corresponding to the first measuring area in all the measuring line information is completed; and taking the extraction space range of all the measuring lines with intersection points in the overlapping range corresponding to the first measuring area as the data extraction space range of the first measuring area. The extracted data information is shown in fig. 4.
Therefore, when the measuring point-level data is extracted accurately, the optimal measuring data covered once in the range of the specified area of the user can be accurately obtained, the problem that the data are repeatedly covered at the boundary of the measuring data from different measuring areas is solved, and the difficulty in post-processing data of the user is reduced.
Furthermore, when data screening and extraction are performed aiming at the application of large-scale data such as nationwide data, batch and business processing can be realized based on the method, and the advantages are more obvious.
In one implementation, prior to step 2, the method further comprises:
acquiring aeromagnetic original coding data; the aeromagnetic original coding data of a measuring area are processed to form standard exchange coding data to be used as an entity file for storage, and magnetic field value information and measuring point indexes of measuring points of the measuring area are stored in the entity file of the measuring area. The survey line index is composed of a start point index and an end point index of the survey line. The data of the line can be extracted according to the start point index and the end point index. The measuring points on the same measuring line are stored in sequence, and a uniform measuring point index is established for all measuring points in the measuring area.
Specifically, the steps can be realized by the following steps:
acquiring aeromagnetic data of a survey area (survey area); and collecting related metadata information by taking the measuring area as a basic unit, wherein the metadata information comprises two parts of space information and attribute information, and the aeromagnetic data of one measuring 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 covered by the aeromagnetic survey, or called as the range information of the survey area, and is expressed by a surface element by uniformly adopting a 2000 national geodetic coordinate system.
The aeromagnetic data refers to aeromagnetic original encoding data obtained by an aeromagnetic survey system, and standard exchange encoding data is formed after the aeromagnetic original encoding data is processed by aeromagnetic geophysical data to be used as a final data storage file. And storing one by one according to the measuring region- > measuring line- > measuring point mode. And the entity data of one measuring area is used as a storage file, and is organized and stored according to measuring lines and measuring points. Each measuring point comprises information such as latitude and longitude, ellipsoid height, height above ground (true height), originally measured magnetic field value, compensated magnetic field value, normal magnetic field value, magnetic daily variation value, magnetic compensation value, delta T magnetic abnormal value and the like.
The survey area attribute information comprises but is not limited to survey area codes, survey area names, survey time, survey scale, average flight height, total magnetic survey precision, survey line workload, survey instrument names, magnetometer sensitivity, data recording mode, data storage path information and the like, and the attribute information can be expanded according to actual needs; the data storage path information comprises an original coded data storage path and a standard exchange coded data storage path, and the subsequent entity data extraction is based on the standard exchange coded data. Collecting metadata information of each measuring line, wherein the metadata information comprises two parts of space information and attribute information, and one measuring line corresponds to one piece of metadata information; the survey line space information is expressed by line elements, a 2000 national geodetic coordinate system is uniformly adopted, and the survey line space information is formed by connecting the coordinates of a survey line initial survey point, a survey line final survey point and the coordinates of the survey points after rarefaction; the attribute information comprises a survey line code, a region code, a frame number, a survey line direction, a start survey point coordinate, an end survey point coordinate, a start survey point index and an end survey point index, wherein the start point index and the end survey point index form index information of the survey line in the entity file. And storing the measurement region metadata information and the measurement line metadata information based on a spatial database. And storing the aeromagnetic data processed by the standard exchange coding, and updating data storage path information in the attribute of the survey area.
Furthermore, in the embodiment, a unified data content and a data organization mode are defined for the original aeromagnetic measurement field value and the related measurement parameter value, a standard exchange coding data format is formed to organize and store the entity data, and different later-period measurement area data can be uniformly stored in the same data file, so that the difficulty in processing a plurality of measurement area files of a user is reduced, meanwhile, data service is provided for the user based on the original aeromagnetic measurement field value, and the problem of false magnetic anomaly in the current data service based on the gridded data is solved.
In one implementation, the standard exchange encoding data is stored in a binary encoding form during the data storage, and the data obtained by data extraction is XYZ or CSV format data which is stored by adopting plain codes in the standard exchange encoding data.
Therefore, compared with the situation that the original data are directly acquired by the airplane measurement system and can be identified only by professional software, the user cannot directly read the data, the data are distributed in the format of XYZ or CSV and the like, the user can conveniently identify the data by using common data processing software, and the data index is conveniently established.
In one implementation, after the extracting data, the method further includes: the feedback of the extracted data to the user through the diagram may specifically include:
data diagram: and (3) dropping the inquired measuring points (reaction magnetic field data) on a drawing according to a scale, representing by circular points rendered according to colors, wherein the central points of the points are longitude and latitude positions of the measuring points, magnetic fields (or other values) are distributed according to input color bar colors (default), and can be distributed linearly, different colors correspond to different magnetic field values, and magnetic field values (marking magnetic field abnormal values and ground clearance (true height) and decimal point numbers of minimum 0.01nT and 0.1m) are marked on the circular points.
Working degree diagram: expressing by using the surface elements, forming a working degree graph by distinguishing according to different colors according to the finally preferred measuring area, using the working degree graph as an attached drawing, and simultaneously marking a legend by using the metadata information of the measuring area; and simultaneously drawing a query retrieval range on the query retrieval range, wherein the query retrieval range is represented by a surface without filling.
Thus, a magnetic field data map expressed only by actual magnetic field values is output without any change in data.
The retrieval section is explained in detail by the following example.
The user retrieval conditions comprise a specified space retrieval range and attribute constraint conditions (such as one or more of measuring region attribute metadata including a measuring scale, a data recording mode, an average flying height, measuring time, a measuring instrument model and the like); determining corresponding scoring weight value by user for each constraint, and defining the total weight function as M ═ Σ f (x)i)。
Searching the survey area:
step 1, according to a space retrieval range and attribute constraint conditions specified by a user, performing space-attribute integrated retrieval based on measured area metadata, and retrieving all measured areas meeting the user retrieval constraint conditions, wherein the total number of the measured areas is recorded as N;
step 2, according to preset scoring condition information, scoring the data quality of the measuring area, wherein the scoring condition can be defined according to user requirements, and the scaling scale can be defined as 1: 5 ten thousand weight plus 40 points, 1: 10 ten thousand weight plus 35 points, 1: 20 ten thousand weight plus 30 points, 1: 50 ten thousand weight plus 25 points, 1: 100 ten thousand weight plus 20 points; the data recording mode is digital recording, the weight is 30 points, and the analog recording is 15 points; measuring the highest weight of the flying height by 30 points; and finally, counting the weight accumulated value of each condition to obtain a total weight value.
For example, M ═ f (x) may be set1)+f(x2)+f(x3)；
Wherein x is1Weight value, x, of the scale2Weight value, x, indicating the manner in which the data is recorded3When the weight values obtained by the three items are the same, the overlap ratio of the space range of the measurement area and the search range is further compared, the overlap ratio is large, and the additional weight is high.
Acquiring all the measurement areas meeting the requirement of user retrieval constraint conditions and corresponding total weights, sequencing the measurement areas from large to small according to the total weights, determining priority sequencing, wherein the higher the total weight is, the higher the priority is, and obtaining a final measurement area priority List1 which comprises a measurement area number, a measurement area space range, the total weight and the priority;
step 3, obtaining the measurement area with the highest priority from the measurement area priority List1 as a first measurement area, obtaining the measurement area code and the measurement area spatial range information of the measurement area, and deleting the measurement area information from the List 1;
step 4, calculating to obtain a spatial range where the measurement area spatial range of the first measurement area and the spatial retrieval range are intersected and overlapped, recording as an overlapping range corresponding to the first measurement area, adding the spatial range to a measurement area data extraction List2, wherein the spatial range includes a measurement area code, a priority, an overlapping range corresponding to the first measurement area, and an entity data storage path;
step 5, calculating a part which is removed from the space retrieval range and is intersected and covered with the measurement area range to obtain a new space retrieval range;
step 6, repeating the step 3 to obtain the measuring area with the highest priority in the List1, comparing the measuring area with the new space retrieval range obtained in the step 5, and repeating the step 4-5 if the measuring area and the new space retrieval range have overlapping areas; if the two areas do not have the overlapping area, repeating the step 3 to continuously select the next measurement area with the highest priority from the List1 as the updated first measurement area for comparison;
and 7, repeating the steps until the space retrieval range initially set in the user constraint condition is completely covered or all the measured areas in the List1 are compared, and finally obtaining an optimal measured area data extraction List List2 meeting the user constraint condition, wherein the number of the measured areas is recorded as M, and M is less than or equal to N.
And (3) line measurement retrieval:
step 1, obtaining a measuring area code of a measuring area, an overlapping range corresponding to the measuring area and an entity data storage path from a measuring area data extraction List List 2;
step 2, performing space-attribute integrated retrieval in the line measurement metadata information according to the line measurement region codes and the overlapping range information corresponding to the line measurement regions to obtain all line measurement information meeting conditions in the line measurement regions, wherein the information comprises the information of the line measurement codes, the codes of the affiliated line measurement regions, line measurement space information, line measurement starting point coordinates, line measurement end point coordinates, line measurement starting point indexes, line measurement end point indexes and the like;
step 3, obtaining the space information of a measuring line and calculating two intersection points (X) of the overlapping range corresponding to the first measuring areacross1，Ycross1)、(Xcross2，Ycross2)；
And 4, extracting all measuring point information of all measuring lines with intersection points in the overlapping range corresponding to the first measuring area according to the entity data storage path information and the start-stop index information of the measuring lines in the entity file acquired in the step 2.
And (3) measuring point retrieval: repeat steps 1-3 of the inline search.
Step 4, according to the coordinate (x) of the starting point of the survey line obtained in the step 2start，ystart) Line measurement end point coordinate (x)end，yend) INDEX of starting point of measuring line INDEXstartLine end point INDEX INDEXendAnd 3, the intersection point (x) of the overlapping range of the measuring line and the first measuring area obtained in the step 31，y1)(x2，y2) Calculating two INDEX values INDEX1、INDEX2；
Step 5, acquiring entity data storage path information in the measuring area, and storing the information in the entity data file according to INDEX1、INDEX2Two indexes obtain corresponding measuring point coordinates (X)index1，Yindex1)(Xindex2，Yindex2) Interpretation of the measurement points (X)index1，Yindex1)(Xindex2，Yindex2) Whether the first measurement area is in the corresponding overlapping range of the corresponding first measurement area; if yes, further judging the self-test point (X)index1，Yindex1) From, away from the measuring point (X)index2，Yindex2) And whether the continuous adjacent measuring points close to the starting point or the end point direction of the measuring line are in the overlapping range corresponding to the first measuring area or not until the last measuring point in the overlapping range corresponding to the first measuring area is found, recording the index of the measuring point as a critical index THRESHOLD extracted from the measuring point data1(ii) a The self-measured points (X) are also judged one by oneindex2，Yindex2) From, away from the measuring point (X)index1，Yindex1) The last measuring point index in the overlapping range corresponding to the first measuring area in the continuous adjacent measuring points close to the starting point or the end point of the measuring line is used as another critical index THRESHOLD extracted from the measuring point data2。
If not, then judge one by oneFrom (X)index1，Yindex1) From, near to the measurement point (X)index2，Yindex2) Whether the measuring points which are continuously adjacent in the direction are in the overlapping range corresponding to the first measuring area or not is judged until the first measuring point which is in the overlapping range corresponding to the first measuring area is found, the index of the measuring point is recorded and is used as the critical index THRESHOLD extracted from the measuring point data1The same for the self-test point (X)index2，Yindex2) Start close to measuring point (X)index1，Yindex1) The continuous adjacent measuring points in the direction are judged one by one, the index of the first measuring point in the overlapping range corresponding to the first measuring area is found and is used as another critical index THRESHOLD extracted from the measuring point data2。
When the overlap range of the measuring line and the first measuring area has only one intersection point, a critical index is obtained by calculation according to the method, and whether the starting point or the end point of the measuring line is in the overlap range corresponding to the first measuring area is further judged, if so, the index of the starting point or the end point is used as another critical index.
Step 6, directly extracting magnetic field value recording information of all measuring points in the range of the two critical indexes in the measuring line according to the two critical indexes obtained in the step 5;
7, repeating the steps 3-6 until the entity data of all the measuring lines meeting the conditions in the measuring area are extracted;
step 8, repeating the step 1, continuously obtaining the measuring area code, the overlapping range corresponding to the first measuring area and the entity data storage path of the next measuring area from the measuring area data extraction List List2, and executing the steps 2-7;
by analogy, until all the measurement area data in the List2 are extracted, the data are provided to the user by using standard exchange coding data (XYZ format or CSV format), and a data graph can be provided at the same time.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or 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, and the like) 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 application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing descriptions of specific exemplary embodiments of the present invention have been presented for purposes of illustration and description. It is not intended to limit the invention to the precise form disclosed, and obviously many modifications and variations are possible in light of the above teaching. The exemplary embodiments were chosen and described in order to explain certain principles of the invention and its practical application to enable one skilled in the art to make and use various exemplary embodiments of the invention and various alternatives and modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims and their equivalents.
Claims (9)
1. A retrieval and extraction method of aeromagnetic measurement data is characterized by comprising the following steps:
acquiring a user retrieval condition, wherein the user retrieval condition comprises a space retrieval range and extraction precision, the space retrieval range included by the user retrieval condition is an initial space retrieval range, and the extraction precision comprises one of a survey area level, a survey line level and a survey point level;
acquiring parameter data of all measurement areas intersected with the space retrieval range, and scoring all measurement areas intersected with the space retrieval range according to a preset principle;
acquiring the measuring region with the highest score from all the measuring regions intersected with the space retrieval range as a first measuring region;
acquiring a spatial range of the first measurement area, which is intersected and overlapped with the spatial retrieval range, as an overlapping range corresponding to the first measurement area;
removing a part which is intersected and covered with the overlapping range corresponding to the first measuring area from the space retrieval range to obtain an updated space retrieval range;
in the other measurement areas which are intersected with the space retrieval range and are not used as or are used as the first measurement area, taking the measurement area with the highest score as an updated first measurement area, repeatedly obtaining the space range of the first measurement area, which is intersected and overlapped with the updated space retrieval range, and taking the space range as the overlapped range corresponding to the updated first measurement area until the initial space retrieval range is completely covered or all the measurement areas intersected with the space retrieval range are completely compared;
and extracting data according to the extraction precision and the overlapping range, wherein the overlapping range is the sum of the spatial ranges of all the first measuring areas and the spatial retrieval range which are intersected and overlapped.
2. The retrieval extraction method according to claim 1, wherein extracting data based on the extraction accuracy and the overlap range comprises:
when the extraction precision is in a measurement area level, taking all entity data of the first measurement area corresponding to the overlapping range as a data extraction space range of the first measurement area;
and extracting data according to the sum of the data extraction space ranges of the first measuring area.
3. The retrieval extraction method according to claim 1, wherein extracting data based on the extraction accuracy and the overlap range comprises:
when the extraction precision is in a non-measurement-area level, acquiring information of all measurement lines in a first measurement area corresponding to the overlapping range;
calculating the data extraction spatial range of the first measuring area according to the information of all measuring lines in the first measuring area;
acquiring information of all measuring lines in another first measuring area corresponding to the overlapping range, and repeating the step of calculating the data extraction spatial range of the first measuring area until the calculation of all measuring areas corresponding to the overlapping range is completed;
and extracting data according to the sum of the data extraction space ranges of the first measuring area.
4. The retrieval extraction method of claim 3, wherein when the extraction accuracy is in a line measurement level, calculating a data extraction spatial range of the first measurement area based on information of all lines in the first measurement area comprises:
acquiring information of all measuring points corresponding to measuring lines with intersection points in the overlapping range corresponding to the first measuring area from information of all measuring lines in the first measuring area;
and taking the information of all measuring points corresponding to the measuring lines with the intersection points in the overlapping range corresponding to the first measuring area as the data extraction space range of the first measuring area.
5. The retrieval extraction method according to claim 3, wherein when the extraction accuracy is a measurement point level, calculating a data extraction spatial range of the first measurement area from the all-line information includes:
in all the measuring line information in the first measuring area, if an intersection point exists in the overlapping range corresponding to one measuring line and the first measuring area, calculating the intersection point of the overlapping range corresponding to the measuring line and the first measuring area;
judging the measuring points on the measuring line, acquiring the information of all measuring points of the measuring line in the overlapping range corresponding to the first measuring area, and taking the information of all measuring points as the extraction space range of the measuring line;
repeating the steps until the calculation of the extraction space range of the measuring line with the intersection point in the overlapping range corresponding to the first measuring area in all the measuring line information is completed;
and taking the extraction space range of all the measuring lines with intersection points in the overlapping range corresponding to the first measuring area as the data extraction space range of the first measuring area.
6. The retrieval extraction method of claim 5, wherein the acquiring information of all the measuring points of the measuring line in the overlapping range corresponding to the first measuring region comprises:
respectively calculating the INDEX of the measuring points near the two intersection points of the measuring lines according to the following formula1、INDEX2；
Wherein, INDEX1、INDEX2Taking an integer, (X)cross1，Ycross1)、(Xcross2，Ycross2) The overlapping range of the measuring line and the first measuring area has two intersection point coordinates (X)start，Ystart) Is the coordinate of the starting point of the survey line, (X)end，Yend) Is the end point coordinate, INDEXstartIs the INDEX of the starting point of the survey line, INDEXendIs the index of the end point of the survey line;
obtaining the on-line and the INDEX1、INDEX2Measuring point coordinate (X) corresponding to two index valuesindex1，Yindex1)(Xindex2，Yindex2) Judgment measuring point (X)index1，Yindex1)(Xindex2，Yindex2) Whether the first measurement area is in the corresponding overlapping range of the first measurement area or not;
if yes, then judge one by one from (X)index1，Yindex1) From, away from the measuring point (X)index2，Yindex2) Whether the continuous adjacent measuring points close to the starting point or the end point of the measuring line are in the overlapping range corresponding to the first measuring area until the last measuring point in the overlapping range corresponding to the first measuring area is found, recording the index of the measuring point, and taking the index as the critical index THRESHOLD extracted from the measuring point data1；
The measured points (X) are also judged one by oneindex2，Yindex2) From, away from the measuring point (X)index1，Yindex1) The last measuring point index in the overlapping range corresponding to the first measuring area in the continuous adjacent measuring points close to the starting point or the end point of the measuring line is used as another critical index THRESHOLD extracted from the measuring point data2；
Placing the line under two critical indices THRESHOLD1、THRESHOLD2And all the measuring point data in between are used as the information of all the measuring points of the measuring line in the corresponding overlapping range of the first measuring area.
7. The retrieval extraction method of claim 1, wherein the user retrieval condition further includes an attribute constraint condition;
the parameter data of all the measurement areas intersected with the space retrieval range are the parameter data of all the measurement areas which meet the attribute constraint condition and are intersected with the space retrieval range;
the scoring of all measurement areas intersected with the space retrieval range according to a preset principle comprises the following steps:
and calculating the weights of all the measurement areas which meet the attribute constraint condition and are intersected with the space retrieval range according to a preset principle, sequencing the total weights, and taking the measurement area with the maximum total weight as the measurement area with the highest score.
8. The search extraction method of claim 1, wherein prior to obtaining parameter data for all regions intersecting the spatial search range, the method further comprises:
acquiring aeromagnetic original coding data;
the method comprises the steps of processing aviation geophysical data of a measuring area to form a uniform standard exchange coding data format, storing the uniform standard exchange coding data format as an entity file in a binary mode, wherein measuring point field value information of the measuring area and index information of each measuring point are stored in the entity file of the measuring area.
9. The retrieval extraction method of claim 1, wherein after the extracting data, further comprising: the extracted data is fed back to the user by means of a graphical representation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210150948.XA CN114528454B (en) | 2022-02-18 | 2022-02-18 | Method for retrieving and extracting aeromagnetic measurement data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210150948.XA CN114528454B (en) | 2022-02-18 | 2022-02-18 | Method for retrieving and extracting aeromagnetic measurement data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114528454A true CN114528454A (en) | 2022-05-24 |
| CN114528454B CN114528454B (en) | 2022-08-19 |
Family
ID=81623351
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210150948.XA Active CN114528454B (en) | 2022-02-18 | 2022-02-18 | Method for retrieving and extracting aeromagnetic measurement data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114528454B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116755619A (en) * | 2023-06-06 | 2023-09-15 | 中国自然资源航空物探遥感中心 | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104122597A (en) * | 2013-08-22 | 2014-10-29 | 中国科学院遥感与数字地球研究所 | Unmanned aerial vehicle aeromagnetic detecting system and method |
| US20160171010A1 (en) * | 2014-04-10 | 2016-06-16 | Korea Resources Corporation | Aeromagnetic pre-processing system based on graphic user interface and pre-processing method using the thereof |
| CN109541704A (en) * | 2018-12-05 | 2019-03-29 | 加泰科(深圳)科技有限公司 | Three axis fluxgate magnetic airborne survey systems of one kind and correction compensation method |
| CN111666370A (en) * | 2020-07-28 | 2020-09-15 | 中国人民解放军国防科技大学 | Semantic indexing method and device for multi-source heterogeneous space data |
| CN112485837A (en) * | 2020-12-22 | 2021-03-12 | 中国地质科学院地球物理地球化学勘查研究所 | Unmanned aerial vehicle aeromagnetic horizontal gradient measurement system and method |
| CN113282584A (en) * | 2021-05-28 | 2021-08-20 | 福州大学 | Method and system for searching space-time cube data of earth observation image |
| CN113447993A (en) * | 2021-07-26 | 2021-09-28 | 中国人民解放军61540部队 | Compensation flight method and system for magnetic vector measurement and magnetic compensation method and system |
| CN113514895A (en) * | 2021-07-12 | 2021-10-19 | 两山资源科技(北京)有限公司 | Multi-rotor unmanned aerial vehicle aeromagnetic detection system and method |
-
2022
- 2022-02-18 CN CN202210150948.XA patent/CN114528454B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104122597A (en) * | 2013-08-22 | 2014-10-29 | 中国科学院遥感与数字地球研究所 | Unmanned aerial vehicle aeromagnetic detecting system and method |
| US20160171010A1 (en) * | 2014-04-10 | 2016-06-16 | Korea Resources Corporation | Aeromagnetic pre-processing system based on graphic user interface and pre-processing method using the thereof |
| CN109541704A (en) * | 2018-12-05 | 2019-03-29 | 加泰科(深圳)科技有限公司 | Three axis fluxgate magnetic airborne survey systems of one kind and correction compensation method |
| CN111666370A (en) * | 2020-07-28 | 2020-09-15 | 中国人民解放军国防科技大学 | Semantic indexing method and device for multi-source heterogeneous space data |
| CN112485837A (en) * | 2020-12-22 | 2021-03-12 | 中国地质科学院地球物理地球化学勘查研究所 | Unmanned aerial vehicle aeromagnetic horizontal gradient measurement system and method |
| CN113282584A (en) * | 2021-05-28 | 2021-08-20 | 福州大学 | Method and system for searching space-time cube data of earth observation image |
| CN113514895A (en) * | 2021-07-12 | 2021-10-19 | 两山资源科技(北京)有限公司 | Multi-rotor unmanned aerial vehicle aeromagnetic detection system and method |
| CN113447993A (en) * | 2021-07-26 | 2021-09-28 | 中国人民解放军61540部队 | Compensation flight method and system for magnetic vector measurement and magnetic compensation method and system |
Non-Patent Citations (5)
| Title |
|---|
| DANDAN JIANG ET AL.: "Integration of an Aeromagnetic Measurement System Based on an Unmanned Aerial Vehicle Platform and Its Application in the Exploration of the Ma’anshan Magnetite Deposit", 《IEEE ACCESS》 * |
| HU HUI ET AL.: "The design of the pilotless airborne aeromagnetic instrument with compensation algorithm", 《2017 13TH IEEE INTERNATIONAL CONFERENCE ON ELECTRONIC MEASUREMENT & INSTRUMENTS》 * |
| 乔鹏: "地面高精度磁测在航磁异常查证中的应用研究", 《中国优秀硕士学位论文全文数据库 (基础科学辑)》 * |
| 刘强 等: "海洋船载和航空磁力测量数据融合方法研究", 《第二十一届海洋测绘综合性学术研讨会论文集》 * |
| 郭华 等: "吊舱式高温超导全张量磁梯度测量系统研发与应用研究", 《地球物理学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116755619A (en) * | 2023-06-06 | 2023-09-15 | 中国自然资源航空物探遥感中心 | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station |
| CN116755619B (en) * | 2023-06-06 | 2024-01-05 | 中国自然资源航空物探遥感中心 | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114528454B (en) | 2022-08-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103268358B (en) | Multi-source control point image database builds and update method | |
| US8949246B2 (en) | System and method for importing NetCDF data | |
| CN102693266B (en) | Search for method, the navigation equipment and method of generation index structure of database | |
| US9678628B2 (en) | Method for generating control-code by a control-code-diagram | |
| JP3538147B2 (en) | Representative color specification method using reliability | |
| Younis et al. | Detection and annotation of plant organs from digitised herbarium scans using deep learning | |
| CN113434623A (en) | Fusion method based on multi-source heterogeneous space planning data | |
| KR102182654B1 (en) | Method of inspecting digital topography using augmented reality | |
| CN114528454B (en) | Method for retrieving and extracting aeromagnetic measurement data | |
| CN110197035A (en) | A kind of navigation channel change of Underwater Topography analysis system and method | |
| KR20140015812A (en) | Methodology of gis based automatic schematic diagram creation system for supporting input data generation of 1-dimensional water quality model | |
| CN113936214A (en) | Karst wetland vegetation community classification method based on fusion of aerospace remote sensing images | |
| CN117636160A (en) | Automatic high-resolution remote sensing cultivated land block updating method based on semi-supervised learning | |
| JP6174469B2 (en) | Accident analysis utilization support apparatus and method | |
| CN112199886A (en) | Processing method of PRB data deep learning geological map prediction model | |
| CN111427976B (en) | Road freshness obtaining method and device | |
| US20130231897A1 (en) | Systems and methods for efficient analysis of topographical models | |
| CN112634447A (en) | Outcrop rock stratum layering method, device, equipment and storage medium | |
| JP4782610B2 (en) | Satellite image processing system | |
| CN116721227A (en) | Automatic modeling method for three-dimensional geologic model of complex geologic body | |
| CN113496551B (en) | Terrain section line drawing method based on geological outcrop three-dimensional model | |
| CN111667552B (en) | S57 electronic chart depth range rapid judging and filling method and equipment | |
| CN104408316B (en) | Sino-U.S.'s Geotechnical Engineering classification of soils standard handovers device and method | |
| CN102467573B (en) | Digital mapping simplicity compiling method and internal sort method are auxiliary bee-line ranking method | |
| KR101008019B1 (en) | System for fabricating numerical map using method for automatic making independent polygon of contiguity structure |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |