CN114528454B - Method for retrieving and extracting aeromagnetic measurement data - Google Patents
Method for retrieving and extracting aeromagnetic measurement data Download PDFInfo
- Publication number
- CN114528454B CN114528454B CN202210150948.XA CN202210150948A CN114528454B CN 114528454 B CN114528454 B CN 114528454B CN 202210150948 A CN202210150948 A CN 202210150948A CN 114528454 B CN114528454 B CN 114528454B
- 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.)
- Active
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
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
The 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 strength 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 the repeated coverage conditions of aeromagnetic data such as different scales, heights, measurement precision, 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 supply is mainly gridded data, is not based on original measured magnetic field value data obtained by 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 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 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 other measurement areas which are intersected with the space retrieval range and are not used as or used as the first measurement area, taking the measurement area with the highest score as an updated first measurement area, repeatedly acquiring a space range of the first measurement area, which is intersected and overlapped with the updated space retrieval range, and taking the space range as an overlapped range corresponding to the updated first measurement area until the initial space retrieval range is completely covered or all measurement areas intersected with the space retrieval range are 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 region level, all entity data of the first measurement region corresponding to the overlapping range are used as a data extraction space range of the first measurement region; 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 a data extraction space 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 one 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 total 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 measuring point index marks near the two intersection points of the measuring lines according to the following formulaINDEX 1 、INDEX 2 ；
Wherein, INDEX 1 、INDEX 2 Taking an integer, (X) cross1 ，Y cross1 )、(X cross2 ，Y cross2 ) The overlapping range of the measuring line and the first measuring area has two intersection point coordinates (X) start， Y statr ) Is the coordinate of the starting point of the survey line, (X) end ，Y end ) Is the end point coordinate, INDEX start Is INDEX of starting point of survey line, INDEX end Is the index of the end point of the survey line; obtaining the on-line and the INDEX 1 、INDEX 2 Measuring point coordinate (X) corresponding to two index values index1 ，Y index1 )(X index2 ，Y index2 ) Judgment measuring Point (X) index1 ，Y index1 )(X index2 ，Y index2 ) 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 ，Y index1 ) From, away from the measuring point (X) index2 ，Y index2 ) 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 data 1 (ii) a The self-measured points (X) are also judged one by one index2 ，Y index2 ) From, away from the measuring point (X) index1 ，Y index1 ) The last measuring point index which is in the overlapping range corresponding to the first measuring area and is close to the continuously adjacent measuring points in the starting point or end point direction of the measuring line is taken as another critical index THRESHOLD extracted from the measuring point data 2 (ii) a Placing the line under two critical indices THRESHOLD 1 、THRESHOLD 2 In betweenAnd all measuring point data are used as information of all measuring points of the measuring lines in the overlapping range corresponding to 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 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.
In an embodiment of the present invention, before acquiring the parameter data of all the survey areas intersecting the spatial search range, the method further includes: acquiring aeromagnetic original coded 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 aeromagnetic measurement data can realize full-flow automatic screening and extract the optimal aeromagnetic measurement data which can be covered by a specified area for a single time and meet the user requirements, and has high extraction efficiency and accurate extraction.
When data are 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 are extracted precisely, the optimal measuring data covered by a single time 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 post data processing of the 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 values, forms a uniform standard exchange coding data format to organize and store entity data, and different later-period 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 of false magnetic anomaly existing 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 the spatial range of entity data extraction 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 term "comprise" or variations such as "comprises" or "comprising", etc., 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 survey area level, a survey line level and a survey 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 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, 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 mode, 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, obtaining 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 crossed and covered part of the overlapping range corresponding to the first measuring area from the space retrieval range to obtain an 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, a region code that satisfies the search condition and covers the optimal region of the search region once, an overlap range corresponding to the first region, and an entity data storage path information list 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 measuring areas and the spatial retrieval range which are intersected and overlapped.
The space search range comprises an initial space search range and a space construction range updated every 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.
And 7012, extracting data according to the sum of the data extraction spatial ranges of the first measurement areas.
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.
And 7024, extracting data according to the sum of the data extraction spatial ranges of the first measurement areas.
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 of one survey line and the first survey area, calculating an intersection (X) of the overlapping range of the survey line and the first survey area crosst ，Y crosst )；
When the overlapping range of the measuring line and the first measuring area has two intersection point coordinates (X) cross1 ，Y cross1 )、(X cross2 ，Y cross2 ) Obtaining the coordinate (X) of the starting point of the survey line start， Y start ) Endpoint coordinate (X) end ，Y end ) INDEX of starting point of survey line INDEX start Endpoint INDEX INDEX end Then, the INDEX of the measuring point near the two intersection points is estimated and recorded as INDEX according to the following formula 1 、INDEX 2 。
Wherein, INDEX 1 、INDEX 2 Taking an integer.
According to INDEX 1 、INDEX 2 Two index values obtain the coordinate (X) of the measuring point of the corresponding index on the measuring line index1 ，Y index1 )(X index2 ，Y index2 ) Judgment measuring point (X) index1 ，Y index1 )(X index2 ，Y index2 ) 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 ，Y index1 ) From, away from the measuring point (X) index2 ，Y index2 ) 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 data 1 (ii) a The self-measured points (X) are also judged one by one index2 ，Y index2 ) From, away from the measuring point (X) index1 ，Y index1 ) The last measuring point index which is in the overlapping range corresponding to the first measuring area and is close to the continuously adjacent measuring points in the starting point or end point direction of the measuring line is taken as another critical index THRESHOLD extracted from the measuring point data 2 ；
If not, then one by one is judged from (X) index1 ，Y index1 ) From, near to the measurement point (X) index2 ，Y index2 ) 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 data 1 The same for the self-test point (X) index2 ，Y index2 ) Start close to measuring point (X) index1 ，Y index1 ) 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 data 2 。
It should be noted that since the actually obtained survey line is only an approximate straight line, only the survey 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 ，Y cross1 ) Then, a critical index THRESHOLD is obtained by calculation according to the method 1 And 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 THRESHOLD 2 。
According to two critical indexes, THRESHOLD 1 、THRESHOLD 2 Extracting 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 the longitude and latitude, the ellipsoid height, the height from the ground (true height), the originally measured magnetic field value, the compensated magnetic field value, the normal earth magnetic field value, the magnetic daily variation value, the magnetic compensation value, the Δ T magnetic abnormal value, and the like may 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; processing aeromagnetic original encoding data of a measuring area to form standard exchange encoding data serving as an entity file for storage, wherein 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 survey area spatial information is mainly boundary coordinate information of a coverage area covered by the aviation magnetic survey, or called survey area range information, and is expressed by a surface element by uniformly adopting a 2000 national geodetic coordinate system.
The aeromagnetic data refers to aeromagnetic original coding data acquired by an aeromagnetic survey system, and standard exchange coding data is formed after the aeromagnetic original coding 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 altitude, total magnetic survey precision, survey line workload, survey instrument names, magnetometer sensitivity, data recording modes, 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 measuring line space information is expressed by line elements, a 2000 national geodetic coordinate system is uniformly adopted, and the measuring line space information is formed by connecting the coordinates of an initial measuring point and a final measuring point of the measuring line and the coordinates of the measuring point after rarefaction; the attribute information comprises a survey line code, a belonged survey area 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 subjected to standard exchange coding processing, and updating data storage path information in the survey area attribute.
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 graph specifically includes:
data diagram: the inquired measuring points (reaction magnetic field data) are dropped on a picture according to a scale, a dot rendered according to color is used for representing, the center point of the dot is the longitude and latitude position of the measuring point, a magnetic field (or other values) is distributed according to the input color of a color bar (default), the magnetic field can be distributed linearly, different colors correspond to different magnetic field values, and magnetic field numerical values (marking magnetic field abnormal values and ground clearance (true height) and decimal point numerical values of 0.01nT and 0.1m) are marked on the dot.
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 )。
And (3) searching the measuring 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 accumulated weight values of all the conditions to obtain a total weight value.
For example, M ═ f (x) may be set 1 )+f(x 2 )+f(x 3 )；
Wherein x is 1 Weight value, x, of the scale 2 Weight value, x, indicating a data recording mode 3 When 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.
Obtaining all the test areas meeting the requirements of user retrieval constraint conditions and corresponding total weights, sequencing the test areas and the corresponding total weights from large to small, determining priority sequencing, wherein the higher the total weight is, the higher the priority is, and obtaining a final test area priority List List1 which comprises test area numbers, test area space ranges, total weights and priorities;
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 a measurement area spatial range of the first measurement area intersects and overlaps with a spatial retrieval range, recording the spatial range 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 the 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, an overlapping range corresponding to a measuring area and an entity data storage path of the measuring area 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 code and the overlapping range information corresponding to the line measurement region to obtain all line measurement information meeting conditions in the line measurement region, wherein the line measurement information comprises information such as the line measurement code, the affiliated line measurement code, the line measurement space information, the line measurement starting point coordinate, the line measurement terminal point coordinate, the line measurement starting point index, the line measurement terminal point index 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 area cross1 ，Y cross1 )、(X cross2 ，Y cross2 )；
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: and repeating the steps 1-3 of the line measurement retrieval.
Step 4, according to the coordinate (x) of the starting point of the measuring line obtained in the step 2 start ，y start ) Line measurement end point coordinate (x) end ，y end ) INDEX of starting point of measuring line INDEX start Line end point INDEX INDEX end And 3, the intersection point (x) of the overlapping range of the measuring line and the first measuring area obtained in the step 3 1 ，y 1 )(x 2 ，y 2 ) Calculating two INDEX values INDEX 1 、INDEX 2 ；
Step 5, acquiring entity data storage path information in the measuring area, and storing the information in the entity data file according to INDEX 1 、INDEX 2 Two indexes obtain corresponding measuring point coordinates (X) index1 ，Y index1 )(X index2 ，Y index2 ) Interpretation of the measurement points (X) index1 ，Y index1 )(X index2 ，Y index2 ) 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 ，Y index1 ) From, away from the measuring point (X) index2 ，Y index2 ) 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 data 1 (ii) a The self-measured points (X) are also judged one by one index2 ，Y index2 ) From, away from the measuring point (X) index1 ，Y index1 ) The last measuring point index which is in the overlapping range corresponding to the first measuring area and is close to the continuously adjacent measuring points in the starting point or end point direction of the measuring line is taken as another critical index THRESHOLD extracted from the measuring point data 2 。
If not, then one by one is judged from (X) index1 ，Y index1 ) From, near to the measurement point (X) index2 ，Y index2 ) 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 a critical index THRESHOLD extracted from the measuring point data 1 Same for self-test point (X) index2 ，Y index2 ) Get upNear the measurement point (X) index1 ，Y index1 ) 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 data 2 。
When the overlap range of the measuring line and the first measuring area has only one intersection, 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 taken 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 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 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, wherein the method specifically comprises the following steps of: the user retrieval conditions also comprise attribute constraint conditions, weight calculation is carried out on all measurement areas which meet the attribute constraint conditions and are intersected with the space retrieval range according to a preset principle, and the measurement area with the maximum weight is used as the measurement area with the highest score;
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 other measurement areas which are intersected with the space retrieval range and are not used as or used as the first measurement area, taking the measurement area with the highest score as an updated first measurement area, repeatedly acquiring a space range of the first measurement area, which is intersected and overlapped with the updated space retrieval range, and taking the space range as an overlapped range corresponding to the updated first measurement area until the initial space retrieval range is completely covered or all measurement areas intersected with the space retrieval range are 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 of claim 1, wherein extracting data based on the extraction precision 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 a 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 all the measurement line information in the first measurement area includes:
in all the survey line information in the first survey area, if an intersection point exists in the overlapping range corresponding to one survey line and the first survey area, calculating the intersection point of the overlapping range corresponding to the survey line and the first survey 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 measuring point INDEX marks INDEX near two intersection points of the measuring lines according to the following formula 1 、INDEX 2 ；
Wherein, INDEX 1 、INDEX 2 Taking an integer, (X) cross1 ,Y cross1 )、(X cross2 ，Y cross2 ) The overlapping range of the measuring line and the first measuring area has two intersection point coordinates (X) start ,Y start ) Is the coordinate of the starting point of the survey line, (X) end ，Y end ) Is the end point coordinate, INDEX start Is the INDEX of the starting point of the survey line, INDEX end Is the index of the end point of the survey line;
obtaining the data on the survey line and the INDEX 1 、INDEX 2 Measuring point coordinate (X) corresponding to two index values index1 ,Y index1 )(X index2 ,Y index2 ) Judgment measuring point (X) index1 ,Y index1 )(X index2 ,Y index2 ) Whether the measurement area is in the overlapping range corresponding to the first measurement area or not;
if yes, then judge one by one from (X) index1 ,Y index1 ) From, away from the measuring point (X) index2 ,Y index2 ) Whether continuous adjacent measuring points close to the starting point or end point direction 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 a critical index THRESHOLD extracted from the measuring point data 1 ；
The measured points (X) are also judged one by one index2 ,Y index2 ) From, away from the measuring point (X) index1 ,Y index1 ) 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 data 2 ；
Placing the line under two critical indices THRESHOLD 1 、TSHOLD 2 And 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 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.
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 CN114528454A (en) | 2022-05-24 |
| CN114528454B true 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) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116755619B (en) * | 2023-06-06 | 2024-01-05 | 中国自然资源航空物探遥感中心 | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104122597A (en) * | 2013-08-22 | 2014-10-29 | 中国科学院遥感与数字地球研究所 | Unmanned aerial vehicle aeromagnetic detecting system and method |
| CN109541704A (en) * | 2018-12-05 | 2019-03-29 | 加泰科(深圳)科技有限公司 | Three axis fluxgate magnetic airborne survey systems of one kind and correction compensation 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 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR101462252B1 (en) * | 2014-04-10 | 2014-11-14 | 한국광물자원공사 | Aeromagnetic pre-processing method based on graphic user interface using the aeromagnetic pre-processing system |
| CN111666370B (en) * | 2020-07-28 | 2022-04-22 | 中国人民解放军国防科技大学 | 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 |
| 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 (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104122597A (en) * | 2013-08-22 | 2014-10-29 | 中国科学院遥感与数字地球研究所 | Unmanned aerial vehicle aeromagnetic detecting system and method |
| CN109541704A (en) * | 2018-12-05 | 2019-03-29 | 加泰科(深圳)科技有限公司 | Three axis fluxgate magnetic airborne survey systems of one kind and correction compensation 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 |
Non-Patent Citations (3)
| Title |
|---|
| 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;Dandan Jiang et al.;《IEEE Access》;20201015;第8卷;1-5 * |
| 吊舱式高温超导全张量磁梯度测量系统研发与应用研究;郭华 等;《地球物理学报》;20220110;第65卷(第1期);360-369 * |
| 地面高精度磁测在航磁异常查证中的应用研究;乔鹏;《中国优秀硕士学位论文全文数据库 (基础科学辑)》;20130215(第2期);A011-235 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114528454A (en) | 2022-05-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP3848509A1 (en) | Crack analysis device, crack analysis method, and crack analysis program | |
| JP3747589B2 (en) | Image feature amount comparison device and recording medium storing image feature amount comparison program | |
| CN103268358B (en) | Multi-source control point image database builds and update method | |
| JP5086157B2 (en) | Image interpretation support method | |
| US8949246B2 (en) | System and method for importing NetCDF data | |
| US20210286104A1 (en) | Systems, apparatus, and methods of watershed modeling and assessment | |
| WO2008105611A1 (en) | Database auto-building method for link of search data in gis system using cad drawings | |
| US4685068A (en) | Generic database generator system and method | |
| CN114528454B (en) | Method for retrieving and extracting aeromagnetic measurement data | |
| US20130156307A1 (en) | Systems and methods for efficiently and accurately detecting changes in spatial feature data | |
| KR20140015812A (en) | Methodology of gis based automatic schematic diagram creation system for supporting input data generation of 1-dimensional water quality model | |
| KR102182654B1 (en) | Method of inspecting digital topography using augmented reality | |
| CN110136174A (en) | A kind of target object tracking and device | |
| CN110197035A (en) | A kind of navigation channel change of Underwater Topography analysis system and method | |
| CN113989414A (en) | Method and system for drawing cross-section diagram and computer storage medium | |
| CN112199886A (en) | Processing method of PRB data deep learning geological map prediction model | |
| US20130231897A1 (en) | Systems and methods for efficient analysis of topographical models | |
| JP6174469B2 (en) | Accident analysis utilization support apparatus and method | |
| CN113936214A (en) | Karst wetland vegetation community classification method based on fusion of aerospace remote sensing images | |
| JP4782610B2 (en) | Satellite image processing system | |
| KR20190063062A (en) | Method and Device for visualizing the change in topography | |
| CN115357675A (en) | Method and system for establishing image control point database through image control point standardization processing | |
| CN112634447B (en) | Outcrop stratum layering method, device, equipment and storage medium | |
| CN104408316B (en) | Sino-U.S.'s Geotechnical Engineering classification of soils standard handovers device and method | |
| CN111667552B (en) | S57 electronic chart depth range rapid judging and filling method and equipment |
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 |