CN116755619B - Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station - Google Patents
Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station Download PDFInfo
- Publication number
- CN116755619B CN116755619B CN202310664904.3A CN202310664904A CN116755619B CN 116755619 B CN116755619 B CN 116755619B CN 202310664904 A CN202310664904 A CN 202310664904A CN 116755619 B CN116755619 B CN 116755619B
- Authority
- CN
- China
- Prior art keywords
- data
- measurement
- line
- measuring
- scale
- 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
- 238000005259 measurement Methods 0.000 title claims abstract description 430
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000003860 storage Methods 0.000 claims abstract description 31
- 230000004044 response Effects 0.000 claims abstract description 28
- 238000005070 sampling Methods 0.000 claims description 196
- 238000009877 rendering Methods 0.000 claims description 80
- 230000002776 aggregation Effects 0.000 claims description 46
- 238000004220 aggregation Methods 0.000 claims description 46
- 230000000007 visual effect Effects 0.000 claims description 42
- 238000012545 processing Methods 0.000 claims description 29
- 230000011218 segmentation Effects 0.000 claims description 29
- 230000000875 corresponding effect Effects 0.000 claims description 21
- 238000012800 visualization Methods 0.000 claims description 21
- 238000013075 data extraction Methods 0.000 claims description 19
- 238000004590 computer program Methods 0.000 claims description 17
- 230000015654 memory Effects 0.000 claims description 12
- 238000013500 data storage Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 7
- 230000002596 correlated effect Effects 0.000 claims description 5
- 230000003321 amplification Effects 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 230000008569 process Effects 0.000 description 18
- 230000006870 function Effects 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 239000000284 extract Substances 0.000 description 5
- 239000011591 potassium Substances 0.000 description 5
- 229910052700 potassium Inorganic materials 0.000 description 5
- 238000001228 spectrum Methods 0.000 description 5
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 5
- 230000005540 biological transmission Effects 0.000 description 4
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 4
- 238000005404 magnetometry Methods 0.000 description 4
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 description 3
- 229910052776 Thorium Inorganic materials 0.000 description 3
- 229910052770 Uranium Inorganic materials 0.000 description 3
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 239000002131 composite material Substances 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000001276 controlling effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- ZSLUVFAKFWKJRC-UHFFFAOYSA-N thorium Chemical compound [Th] ZSLUVFAKFWKJRC-UHFFFAOYSA-N 0.000 description 2
- NUXLOJGLOHSJEF-UHFFFAOYSA-N 1,1-bis(4-chlorophenyl)ethanol;(4-chlorophenyl)sulfanyl-(2,4,5-trichlorophenyl)diazene Chemical compound C=1C=C(Cl)C=CC=1C(O)(C)C1=CC=C(Cl)C=C1.C1=CC(Cl)=CC=C1SN=NC1=CC(Cl)=C(Cl)C=C1Cl NUXLOJGLOHSJEF-UHFFFAOYSA-N 0.000 description 1
- 241000287196 Asthenes Species 0.000 description 1
- FGUUSXIOTUKUDN-IBGZPJMESA-N C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 Chemical compound C1(=CC=CC=C1)N1C2=C(NC([C@H](C1)NC=1OC(=NN=1)C1=CC=CC=C1)=O)C=CC=C2 FGUUSXIOTUKUDN-IBGZPJMESA-N 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007596 consolidation process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000013441 quality evaluation Methods 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/06—Digital input from, or digital output to, record carriers, e.g. RAID, emulated record carriers or networked record carriers
- G06F3/0601—Interfaces specially adapted for storage systems
- G06F3/0628—Interfaces specially adapted for storage systems making use of a particular technique
- G06F3/0638—Organizing or formatting or addressing of data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/06—Digital input from, or digital output to, record carriers, e.g. RAID, emulated record carriers or networked record carriers
- G06F3/0601—Interfaces specially adapted for storage systems
- G06F3/0602—Interfaces specially adapted for storage systems specifically adapted to achieve a particular effect
- G06F3/0614—Improving the reliability of storage systems
- G06F3/0619—Improving the reliability of storage systems in relation to data integrity, e.g. data losses, bit errors
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/06—Digital input from, or digital output to, record carriers, e.g. RAID, emulated record carriers or networked record carriers
- G06F3/0601—Interfaces specially adapted for storage systems
- G06F3/0628—Interfaces specially adapted for storage systems making use of a particular technique
- G06F3/0638—Organizing or formatting or addressing of data
- G06F3/0643—Management of files
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Human Computer Interaction (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Measuring Magnetic Variables (AREA)
Abstract
The application discloses a method, a device, equipment and a medium for slicing measurement data of an aviation magnetic release comprehensive station, and belongs to the technical field of aviation geophysical prospecting. According to the embodiment of the application, the dividing interval and the number of the target measuring lines are determined according to the measuring scale, so that dividing lines with almost the interval are determined, and dividing is performed based on the same number of the measuring lines, so that the size of each divided data slice is ensured to be consistent as much as possible, a large number of small files are avoided, the deviation of data volume among different data slices is small, the response time of the data slices can be shortened when the data slices are transmitted to a network and loaded and displayed by a browser, and the loading efficiency is improved. In addition, the data of the whole area is only subjected to one-time data slicing treatment, the data of a plurality of layers are not required to be sliced according to different scale levels in advance, repeated storage of the data is avoided, the data is not required to be thinned in advance according to different displayed scales, and the integrity and the accuracy of the measured data are ensured.
Description
Technical Field
The application relates to the technical field of aviation geophysical prospecting, in particular to a method, a device, equipment and a medium for slicing measurement data of an aviation magnetic release comprehensive station.
Background
Aviation geophysical prospecting is short for aviation geophysical prospecting, and refers to flight operations for measuring changes of various physical fields of the earth from the air by using special detecting instruments and equipment carried by an aircraft to investigate geological structures and search for mineral deposits. Among these, aeronautical measurements include four aspects of aeromagnetic (aeromagnetometry), avionic (aeroradiometry), avionic (avionic measurement) and avionic weight (aerogravity measurement). The aviation geophysical prospecting measurement has the advantages of being free of limitation of ground topography and traffic conditions, high in speed, high in efficiency, low in cost, wide in coverage, large in information quantity, obvious in prospecting effect, extremely small in environmental influence and the like. The aeromagnetic and aeromagnetic comprehensive station measurement can be realized simultaneously. The aviation magnetic release comprehensive station measurement is an important technical means for carrying out mineral resource investigation, and particularly has important significance for searching radioactive strategic mineral resources such as uranium ores, thorium ores, potassium ores and the like.
With the development of aviation magnetic discharge, the data sampling rate is continuously improved, and the data volume acquired by aviation magnetic discharge measurement is increased in geometric scale. This huge amount of data requires slicing the data, and then dividing the data into small pieces of data for storage. In the prior art, a pyramid data tile method is generally used for slicing data of the aerospace geophysical prospecting measurement data, generally, one part of measurement data needs to be repeatedly segmented and stored according to a plurality of scale levels, and the data needs to be thinned and stored in advance according to the plurality of scale levels.
Such a scale level corresponds to a slice of data, severely burdening the data storage. In addition, the traditional data slicing method also causes unbalanced data volume among vector tiles, so that a large number of small files are generated, and the problem that the loading and rendering efficiency of a browser is affected when the obtained data slice is loaded and rendered by the subsequent browser is also solved. And such trivial and many duplicate doclets can also lead to difficulties in later data consolidation.
Disclosure of Invention
The embodiment of the application provides a method, a device, equipment and a medium for slicing measurement data of an aviation magnetic release comprehensive station, which can achieve the effects of reducing storage burden, improving response speed and efficiency of loading and display and improving data integrity and accuracy. The technical scheme is as follows:
in one aspect, a method for slicing measurement data of an aviation magnetic release comprehensive station is provided, and the method comprises the following steps:
acquiring measurement data of an aviation magnetic discharge comprehensive station in a certain area, wherein the measurement data of the aviation magnetic discharge comprehensive station comprises measurement metadata and measurement sampling point entity data, and the measurement metadata at least comprises a measurement scale and a measurement line direction;
determining a segmentation interval and the number of target measuring lines according to the measurement scale in the measurement metadata;
Determining a dividing line direction according to the measuring line direction in the measurement metadata;
dividing the measurement sampling point entity data according to the dividing interval, the number of the target measuring lines and the dividing line direction to obtain a plurality of data slices of the measuring area, wherein each data slice comprises a slice metadata file and a measurement sampling point entity data slice file;
a plurality of data slices of the region are stored.
In some embodiments, the determining a split line direction according to the line direction in the measurement metadata includes:
in response to the measurement metadata, the line direction is greater than or equal to 0 degrees and less than or equal to 45 degrees, or the line direction is greater than or equal to 135 degrees and less than or equal to 180 degrees, determining that the dividing line direction is a warp direction, and the line direction is an angle between the line and the weft direction:
and determining that the dividing line direction is the weft direction in response to the measuring metadata that the measuring line direction is more than 45 degrees and less than 135 degrees.
In some embodiments, the storing the plurality of data slices of the region comprises:
and storing a plurality of data slices corresponding to the aviation measurement entity data according to a predefined data storage path setting rule and a slice file naming rule, wherein the data slices comprise slice metadata files and aviation geophysical prospecting measurement entity data slice files.
In some embodiments, the method further comprises:
and responding to a map visualization rendering instruction, displaying a graph corresponding to the size relation according to the size relation between the current display scale and a critical scale, wherein the critical scale is determined based on the distance between the measuring lines, the graph comprises a plane graph, a line graph and a point graph, and the data quantity of different graphs is different.
In some embodiments, the method further comprises:
responding to a map visual rendering instruction, and displaying a region surface graph when the current display scale is smaller than or equal to a first critical scale, wherein the first critical scale is determined based on the distance between the measuring lines and a first distance threshold value of the measuring lines displayed on a screen;
responding to a map visual rendering instruction, displaying a line graph, wherein the current display scale is larger than a first critical scale and smaller than or equal to a second critical scale, the data volume of the line graph is larger than that of the area graph, and the second critical scale is determined based on the distance between the line and a second distance threshold value of the line displayed on a screen;
and responding to a map visual rendering instruction, and displaying a measurement sampling point graph with the current display scale larger than the second critical scale, wherein the data volume of the measurement sampling point graph is larger than the measuring line graph and the area graph.
In some embodiments, the method further comprises:
in response to a map visualization rendering instruction for any plurality of data slices, determining a visualization rendering priority of any plurality of data slices according to an overlapping coverage relationship between a geospatial range of data in the any plurality of data slices and a geospatial range of the display window, the visualization rendering priority being positively correlated with the overlapping coverage relationship;
and according to the visual rendering priority, sequentially obtaining loading and rendering and displaying the plurality of data slices.
In some embodiments, the method further comprises:
responding to a map visual rendering instruction of any plurality of data slices, calculating a dynamic aggregation threshold value of measurement sampling point entity data in any plurality of data slices according to the current display scale and a third interval threshold value of two measurement sampling points displayed on a screen, and carrying out dynamic aggregation processing on the measurement sampling point entity data in any plurality of data slices by taking a measuring line as a unit according to the dynamic aggregation threshold value to obtain target measurement sampling point entity data for visual rendering, wherein the difference value between the distance between adjacent measurement sampling points on the same measuring line in the target measurement sampling point entity data and the third interval threshold value is smaller than a target value;
Performing visual rendering display on the target measurement sampling point entity data;
wherein the dynamic aggregation threshold is: [ dXM/(1000X 1) ];
where d is a third pitch threshold for two measurement sample points displayed on the screen, and l is the measurement sample point pitch, which is = (average aircraft flight speed x 1000)/(3600 x coordinate sampling rate). Wherein, 1:M is the current display scale of the map. The unit of the third interval threshold value displayed on the screen by the two measurement sampling points is millimeter, the unit of the interval between the measurement sampling points is meter, the unit of the average flight speed of the aircraft is kilometer per hour, and the unit of the coordinate sampling rate is secondary per second.
In some embodiments, the method further comprises:
displaying the number of the aggregated measurement sampling points and/or the measurement average value of the currently displayed measurement sampling points on the target measurement sampling points.
In some embodiments, the method further comprises:
responding to a data slice sharing instruction, extracting shared data corresponding to the data extraction range in the plurality of data slices according to the data extraction range in the data slice sharing instruction, and integrating measurement sampling point data in the shared data according to a measurement sequence by taking a measurement line in the shared data as a reference to obtain a shared data file;
And sharing the shared data file to a target address, and providing distribution downloading.
In some embodiments, the dividing the measurement sampling point entity data according to the dividing interval, the number of target lines and the dividing line direction to obtain a plurality of data slices of the area includes:
determining a dividing line position according to the dividing interval, the number of the target measuring lines and the dividing line direction;
and dividing the measurement sampling point entity data based on the dividing line positions to obtain a plurality of data slices of the measuring area.
In some embodiments, the determining a split line location based on the split interval, the number of target lines, and the split line direction includes:
responding to the dividing line direction as the warp direction, starting from 0 degrees of warp, acquiring the starting point positions and the ending point positions of the measuring lines of the target quantity each time, determining the edge positions in the warp direction in the measuring lines of the target quantity according to the starting point positions and the ending point positions, and determining the dividing line positions according to the edge positions and the dividing intervals;
and responding to the dividing line direction as the weft direction, starting from 0-degree weft, acquiring the starting point positions and the ending point positions of the measuring lines of the target quantity each time, determining the edge positions of the measuring lines of the target quantity in the weft direction according to the starting point positions and the ending point positions, and determining the dividing line positions according to the edge positions and the dividing intervals.
In one aspect, there is provided an aircraft magnetic release integrated station measurement data slicing apparatus, the apparatus comprising:
the acquisition module is used for acquiring aviation magnetic discharge comprehensive station measurement data of a certain area, wherein the aviation magnetic discharge comprehensive station measurement data comprise measurement metadata and measurement sampling point entity data, and the measurement metadata at least comprise measurement scale and measurement line direction;
the determining module is used for determining the dividing interval and the number of the target measuring lines according to the measuring scale in the measuring metadata;
the determining module is used for determining the direction of the dividing line according to the direction of the measuring line in the measurement metadata;
the segmentation module is used for segmenting the measurement sampling point entity data according to the segmentation interval, the number of the target measuring lines and the direction of the segmentation line to obtain a plurality of data slices of the measurement area, wherein each data slice comprises a slice metadata file and a measurement sampling point entity data slice file;
and the storage module is used for storing a plurality of data slices of the measuring area.
In some embodiments, the determining module is to:
determining that the dividing line direction is a warp direction in response to the measuring of the metadata, wherein the measuring line direction is greater than or equal to 0 degrees and less than or equal to 45 degrees, or the measuring line direction is greater than or equal to 135 degrees and less than or equal to 180 degrees, and the measuring line direction is an angle between the measuring line and the weft direction;
And determining that the dividing line direction is the weft direction in response to the measuring metadata that the measuring line direction is more than 45 degrees and less than 135 degrees.
In some embodiments, the storage module is to:
and storing a plurality of data slices corresponding to the aviation measurement entity data according to a predefined data storage path setting rule and a slice file naming rule, wherein the data slices comprise slice metadata files and aviation geophysical prospecting measurement entity data slice files.
In some embodiments, the apparatus further comprises:
the display module is used for responding to the map visual rendering instruction, displaying the graph corresponding to the size relation according to the size relation between the current display scale and the critical scale, wherein the critical scale is determined based on the distance between the measuring lines, the graph comprises a plane graph, a line graph and a point graph, and the data quantity of different graphs is different.
In some embodiments, the display module is to:
responding to a map visual rendering instruction, and displaying a region surface graph when the current display scale is smaller than or equal to a first critical scale, wherein the first critical scale is determined based on the distance between the measuring lines and a first distance threshold value of the measuring lines displayed on a screen;
Responding to a map visual rendering instruction, displaying a line graph, wherein the current display scale is larger than a first critical scale and smaller than or equal to a second critical scale, the data volume of the line graph is larger than that of the area graph, and the second critical scale is determined based on the distance between the line and a second distance threshold value of the line displayed on a screen;
and responding to a map visual rendering instruction, and displaying a measurement sampling point graph with the current display scale larger than the second critical scale, wherein the data volume of the measurement sampling point graph is larger than the measuring line graph and the area graph.
In some embodiments, the display module is further to:
in response to a map visualization rendering instruction for any plurality of data slices, determining a visualization rendering priority of any plurality of data slices according to an overlapping coverage relationship between a geospatial range of data in the any plurality of data slices and a geospatial range of the display window, the visualization rendering priority being positively correlated with the overlapping coverage relationship;
and according to the visual rendering priority, sequentially obtaining loading and rendering and displaying the plurality of data slices.
In some embodiments, the display module is further to:
responding to a map visual rendering instruction of any plurality of data slices, calculating a dynamic aggregation threshold value of measurement sampling point entity data in any plurality of data slices according to the current display scale and a third interval threshold value of two measurement sampling points displayed on a screen, and carrying out dynamic aggregation processing on the measurement sampling point entity data in any plurality of data slices by taking a measuring line as a unit according to the dynamic aggregation threshold value to obtain target measurement sampling point entity data for visual rendering, wherein the difference value between the distance between adjacent measurement sampling points on the same measuring line in the target measurement sampling point entity data and the third interval threshold value is smaller than a target value;
performing visual rendering display on the target measurement sampling point entity data;
wherein the dynamic aggregation threshold is: [ dXM/(1000×l) ];
where d is a third pitch threshold for two measurement sample points displayed on the screen, and l is the measurement sample point pitch, which is = (average aircraft flight speed x 1000)/(3600 x coordinate sampling rate). Wherein, 1:M is the current display scale of the map. The unit of the third interval threshold value displayed on the screen by the two measurement sampling points is millimeter, the unit of the interval between the measurement sampling points is meter, the unit of the average flight speed of the aircraft is kilometer per hour, and the unit of the coordinate sampling rate is secondary per second.
In some embodiments, the display module is further to:
displaying the number of the aggregated measurement sampling points and/or the measurement average value of the currently displayed measurement sampling points on the dynamically aggregated target measurement sampling points of the target measurement sampling point sampling entity data.
In some embodiments, the apparatus further comprises a sharing module, where the sharing module is configured to:
responding to a data slice sharing instruction, extracting shared data corresponding to the data extraction range in the plurality of data slices according to the data extraction range in the data slice sharing instruction, and integrating measurement sampling point data in the shared data according to a measurement sequence by taking a measurement line in the shared data as a reference to obtain a shared data file;
and sharing the shared data file to a target address, and providing distribution downloading.
In some embodiments, the segmentation module is to:
determining a dividing line position according to the dividing interval, the number of the target measuring lines and the dividing line direction;
and dividing the measurement sampling point entity data based on the dividing line positions to obtain a plurality of data slices of the measuring area.
In some embodiments, the segmentation module is to:
Responding to the dividing line direction as the warp direction, starting from 0 degrees of warp, acquiring the starting point positions and the ending point positions of the measuring lines of the target quantity each time, determining the edge positions of the measuring lines of the target quantity in the weft direction according to the starting point positions and the ending point positions, and determining the dividing line positions according to the edge positions and the dividing intervals;
and responding to the dividing line direction as the weft direction, starting from 0-degree weft, acquiring the starting point positions and the ending point positions of the measuring lines of the target quantity at each time, determining the edge positions in the warp direction in the measuring lines of the target quantity according to the starting point positions and the ending point positions, and determining the dividing line positions according to the edge positions and the dividing intervals.
In one aspect, an electronic device is provided that includes one or more processors and one or more memories having at least one computer program stored therein, the at least one computer program loaded and executed by the one or more processors to implement various alternative implementations of the aero-magnetic composite station measurement data slicing method described above.
In one aspect, a computer readable storage medium having at least one computer program stored therein is provided, the at least one computer program being loaded and executed by a processor to implement various alternative implementations of the aerial magnetic composite station measurement data slicing method described above.
In one aspect, a computer program product or computer program is provided, the computer program product or computer program comprising one or more program codes, the one or more program codes being stored in a computer readable storage medium. One or more processors of the electronic device are capable of reading the one or more program codes from the computer readable storage medium, the one or more processors executing the one or more program codes such that the electronic device is capable of performing the airborne magnetic and amplification comprehensive station measurement data slicing method of any one of the possible embodiments described above.
According to the embodiment of the application, the dividing interval and the number of the target measuring lines are determined according to the measuring scale, so that dividing lines with almost the interval are determined, and dividing is performed based on the same number of the measuring lines, so that the size of each divided data slice is ensured to be consistent as much as possible, a large number of small files are avoided, the deviation of data volume among different data slices is small, the response time of the data slices can be shortened when the data slices are transmitted to a network and loaded and displayed by a browser, and the loading efficiency is improved. In addition, the data of the whole area is only subjected to one-time data slicing treatment, the data of a plurality of layers are not required to be sliced according to different scale levels in advance, repeated storage of the data is avoided, the data is not required to be thinned in advance according to different displayed scales, and the integrity and the accuracy of the measured data are ensured.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present application, the drawings that are needed in the description of the embodiments will be briefly introduced below, it being obvious that the drawings in the following description are only some embodiments of the present application, and that other drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of an implementation environment of an aircraft magnetic release comprehensive station measurement data slicing method according to an embodiment of the present application;
FIG. 2 is a flowchart of a method for slicing measurement data of an aircraft magnetic release comprehensive station according to an embodiment of the present application;
FIG. 3 is a schematic diagram of a parting line location determination result according to an embodiment of the present application;
FIG. 4 is a flowchart of a method for slicing measurement data of an aircraft magnetic release comprehensive station according to an embodiment of the present application;
fig. 5 is a schematic structural diagram of an aeronautical magnetic release comprehensive station measurement data slicing device according to an embodiment of the present application;
fig. 6 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present application more apparent, the embodiments of the present application will be described in further detail below with reference to the accompanying drawings.
The terms "first," "second," and the like in this application are used to distinguish between identical or similar items that have substantially the same function and function, and it should be understood that there is no logical or chronological dependency between the "first," "second," and "nth" terms, nor is it limited to the number or order of execution. It will be further understood that, although the following description uses the terms first, second, etc. to describe various elements, these elements should not be limited by the terms. These terms are only used to distinguish one element from another element. For example, a first image can be referred to as a second image, and similarly, a second image can be referred to as a first image, without departing from the scope of the various examples. The first image and the second image can both be images, and in some cases, can be separate and distinct images.
The term "at least one" in this application means one or more, the term "plurality" in this application means two or more, for example, a plurality of data packets means two or more.
It should be understood that the terminology used in the description of the various examples herein is for the purpose of describing particular examples only and is not intended to be limiting. As used in the description of various examples and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
It will also be understood that the term "and/or" as used herein refers to and encompasses any and all possible combinations of one or more of the associated listed items. The term "and/or" is an association relationship describing an associated object, meaning that three relationships can exist, e.g., a and/or B, can be represented: a exists alone, A and B exist together, and B exists alone. In addition, the character "/" in the present application generally indicates that the front-rear association object is an or relationship.
It should also be understood that, in the embodiments of the present application, the sequence number of each process does not mean that the execution sequence of each process should be determined by the function and the internal logic of each process, and should not constitute any limitation on the implementation process of the embodiments of the present application.
It should also be understood that determining B from a does not mean determining B from a alone, but can also determine B from a and/or other information.
It will be further understood that the terms "Comprises" and/or "Comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
It should also be understood that the term "if" may be interpreted to mean "when" ("white" or "upon") or "in response to a determination" or "in response to detection". Similarly, the phrase "if a [ stated condition or event ] is detected" may be interpreted to mean "upon a determination" or "in response to a determination" or "upon a detection of a [ stated condition or event ] or" in response to a detection of a [ stated condition or event ], depending on the context.
The following describes the environment in which the present application is implemented.
Fig. 1 is a schematic diagram of an implementation environment of an airborne magnetic station measurement data slicing method according to an embodiment of the present application. The implementation environment includes a terminal 101 or the implementation environment includes a terminal 101 and an avionics integrated station measurement data processing platform 102. The terminal 101 is connected to the aerial magnetic station measurement data processing platform 102 through a wireless network or a wired network.
The terminal 101 can be at least one of a smart phone, a game console, a desktop computer, a tablet computer, an electronic book reader, an MP3 (Moving Picture Experts Group Audio Layer III, moving picture experts compression standard audio layer 3) player, or an MP4 (Moving Picture Experts Group Audio Layer IV, moving picture experts compression standard audio layer 4) player, a laptop portable computer. The terminal 101 is installed and running an application supporting the aerial magnetic recording station measurement data slicing, which can be, for example, a system application, a browser application, a map software.
The terminal 101 can have data slicing, data storage and data display functions, and can perform data slicing processing on acquired measurement data of the aerial magnetic recording and reproducing integrated station, store the obtained data slices, and perform data display based on the stored data slices. The terminal 101 is capable of performing this task independently and also capable of providing data services to the aerial magnetic station survey data processing platform 102. The embodiments of the present application are not limited in this regard.
The aerial magnetic delivery integrated station measurement data processing 102 includes at least one of a server, a plurality of servers, a cloud computing platform, and a virtualization center. The avionics magnetic station measurement data processing platform 102 is configured to provide background services for applications that slice avionics magnetic station measurement data. Optionally, the aviation magnetic release comprehensive station measurement data processing platform 102 performs primary processing work, and the terminal 101 performs secondary processing work; or, the aviation magnetic release comprehensive station measurement data processing platform 102 bears secondary processing work, and the terminal 101 bears primary processing work; alternatively, the avionics magnetic placement integration station measurement data processing platform 102 or the terminal 101, respectively, can be solely responsible for processing tasks. Alternatively, the aviation magnetic release comprehensive station measurement data processing platform 102 and the terminal 101 are cooperatively calculated by adopting a distributed computing architecture.
Optionally, the avionics integrated station measurement data processing platform 102 includes at least one server 1021 and a database 1022, where the database 1022 is configured to store data, and in this embodiment, the database 1022 is configured to store collected avionics integrated station measurement data, and provide data services for the at least one server 1021.
The server may be an independent physical server, a server cluster or a distributed system formed by a plurality of physical servers, or a cloud server providing cloud services, cloud databases, cloud computing, cloud functions, cloud storage, network services, cloud communication, middleware services, domain name services, security services, content delivery networks (Content Delivery Network, CDN), basic cloud computing services such as big data and artificial intelligence platforms. The terminal can be, but is not limited to, a smart phone, a tablet computer, a notebook computer, a desktop computer, a smart speaker, a smart watch, etc.
Those skilled in the art will appreciate that the number of terminals 101 and servers 1021 can be greater or fewer. For example, the number of the terminals 101 and the servers 1021 can be only one, or the number of the terminals 101 and the servers 1021 can be tens or hundreds, or more, and the number and the device type of the terminals or the servers are not limited in the embodiment of the present application.
Fig. 2 is a flowchart of a method for slicing measurement data of an avionic magnetic placement integrated station, which is provided in an embodiment of the present application, and the method is applied to an electronic device, where the electronic device is a terminal or a server, and referring to fig. 2, the method includes the following steps.
201. The electronic equipment acquires measurement data of an aviation magnetic discharge comprehensive station in a certain area, wherein the measurement data of the aviation magnetic discharge comprehensive station comprises measurement metadata and measurement sampling point entity data, and the measurement metadata at least comprises a measurement scale and a measurement line direction.
The aviation magnetic release comprehensive station measurement data are data obtained by aviation magnetic release comprehensive station measurement, wherein the aviation magnetic release comprehensive station measurement can comprise two kinds of measurement of aviation magnetic release and aviation release, and the aviation magnetic release comprehensive station can support the two kinds of measurement, so the aviation magnetic release comprehensive station is called as a comprehensive station of the two kinds of measurement. In some embodiments, the aerial magnetic discharge integrated station measurement may be referred to simply as an integrated station measurement.
The aeronautical magnetic discharge comprehensive station measurement data of a certain area acquired in the step 201 may include two kinds of data: the measurement data of the aviation magnetic amplifying comprehensive station comprises measurement metadata and measurement sampling point entity data. In some embodiments, the avionics integrated station measurement data may include avionics measurement data and avionics measurement data, each data consisting of measurement metadata and measurement sample point entity data.
For measurement metadata in the aviation magnetic release comprehensive station measurement data, the measurement metadata may include a measurement scale and a measurement line direction, where, for example, the measurement scale may be 1:100 ten thousand, or 1:50 ten thousand, or 1: the measurement scale may be any other value, and the measurement scale may be determined by a relevant technician according to the requirement and stored in measurement metadata, which is not limited in the embodiment of the present application. The survey line direction is the direction of the survey path of the aircraft at the time of survey. The line direction can be measured by the angle between the line and the weft direction.
Of course, the measurement metadata may also include other data of aviation measurement, for example, may include any one or more metadata information of measurement type (aviation magnetometry or aviation gamma spectroscopy), line interval (unit: meter), line direction, measurement time, measurement average flying height, aircraft average flying speed (unit: kilometer/hour), coordinate sampling rate (unit: times/second), and the like, may also include other metadata information according to measurement settings, and may include data of a current measurement area and data of a line in the current measurement area, for example, may include a measurement area code, a line number of the line in the measurement area, a start coordinate, an end coordinate, a start index, an end index, and the like of the line. The embodiments of the present application are not limited in this regard.
The measurement sample point entity data may include measurement type, measurement line, measurement sample point information, and the like. For aviation, for exampleThe magnetometry, aviation magnetometry and sampling point entity data comprises coordinate information and attribute information of the measurement sampling points. The attribute information may include, among others, a line number, a reference number, a year of flight, a longitude, a latitude, a magnetic field value, a normal field value, a magnetic daily variable value, a magnetic compensation value, a Δt magnetic anomaly value, an ellipsoidal height, a ground clearance height (true height), and the like. For aviation gamma energy spectrum measurement, aviation gamma energy spectrum measurement sampling point entity data comprise coordinate information and attribute information of measurement sampling points. Wherein the attribute information may include line number, reference number, year of flight, longitude, latitude, total window (Tc) count rate, potassium window (K) count rate, uranium window (U) count rate, thorium window (Th) count rate, potassium% 40 K) Content of uranium 238 U) content and thorium 232 Th), content, total amount, ellipsoidal height, height above ground (true height), etc.
In the embodiment of the application, after the electronic device performs data slicing on the measurement sampling point entity data according to some parameters in the measurement metadata, the data slicing is stored. Some terms involved in the present application are explained first below. The measurement sampling points may be simply referred to as measurement points, and the measurement sampling point data may also be referred to as measurement point data.
The measurement types supported by the aviation magnetic release comprehensive station can comprise aviation magnetic force measurement, aviation gamma energy spectrum measurement and the like. The measurement sampling point entity data is the data obtained by measuring the measurement sampling point through any measurement type.
In a specific example, an aviation magnetic release comprehensive station is used for measuring a measurement network in a measurement area (simply referred to as a measurement area) according to a measurement scale, an aircraft uses global satellite navigation positioning navigation, and when performing flight measurement along a pre-planned measurement line, an aviation geophysical prospecting data recording system is used for acquiring three-dimensional space position data, geophysical field (magnetic field, gamma energy spectrum field) data and the like of measurement sampling points according to a certain sampling rate. And after the data of the measurement sampling points are acquired, the data of the measurement sampling points are obtained later. Typically, the position coordinate data is sampled 2 times per second, the magnetic field data is sampled 10 times per second, and the gamma energy spectrum field is sampled 1 time per second. Of course, the sampling frequency can be set by the relevant technician as desired.
The space between adjacent measuring lines is different according to different requirements of the measuring scale, the average flying speed of the aircraft is also different, namely the space between measuring sampling points is different, the space between measuring sampling points is generally between 1:100 ten thousand and 1:25 ten thousand, the space between measuring sampling points is about 30 meters and 60 meters, the space between measuring sampling points is between 1:10 ten thousand and 1:5 ten thousand, the space between measuring sampling points is about 10 meters and 30 meters, the space between measuring sampling points is between 1:2.5 ten thousand and 1:1 ten thousand, and the space between measuring sampling points is about 10 meters. Thus, the number of measurement sampling points is typically several hundred thousand to several million, depending on the range of the measurement, for the regions having the same measurement scale. The measurement sampling point spacing is the measuring point spacing.
In order to facilitate data processing and quality evaluation, acquired measurement data of the aviation magnetic-discharge comprehensive station are organized by using measuring lines, and in the flight measurement process, flight crews or recording equipment need to carefully record the line number of each measuring line, the time and reference number of entering the measuring line and exiting the measuring line. The measurement sampling point data (i.e., measurement point data) obtained by the flight measurement is strictly organized by the path line, and has the following three characteristics.
Ordering of: the method is characterized in that measurement sampling point data acquired by the aerospace geophysical prospecting measurement are organized according to the navigation line, and are arranged in the sequence of acquisition, namely, in the sequence from small to large according to the reference number.
Continuity: refers to continuous data sampling at given time intervals during the operation of the aerospace geophysical prospecting. The attribute has the meaning that the line data of the overhead measurement can be constrained by the first and last reference numbers of the line, namely the line index.
Uniqueness: in the process of measuring the overhead flight, each track point (marked by a reference number) only belongs to one measuring line, so that the maximum yaw distance can be set to be half of the designed measuring line distance when the measuring line index is edited.
In some embodiments, the measurement sampling point entity data in the measurement data of the avionic integrated magnetic discharge station acquired in the step 201 may be measured by the measurement method described above. When on-line visual sharing and rapid extraction and distribution sharing are required to be carried out on the measurement data of the aviation magnetic release comprehensive station, the data can be firstly sliced and then stored and shared.
In one possible implementation, the avionic integrated station measurement data may be stored in a database from which it may be extracted when the electronic device needs to slice the data. In another possible implementation manner, the measurement data of the aviation magnetic release comprehensive station is just acquired by the acquisition device, and when the acquisition device and the electronic device are the same device, the electronic device can directly acquire the measurement data of the aviation magnetic release comprehensive station. When the acquisition device and the electronic device are not the same device, the electronic device can receive the aviation magnetic discharge comprehensive station measurement data sent by the acquisition device. The foregoing describes several possible implementations for acquiring measurement data of the integrated station for aviation magnetic discharge, and the embodiments of the present application are not limited to what implementation is specifically adopted.
202. And the electronic equipment determines the dividing interval and the number of the target measuring lines according to the measuring scale in the measuring metadata.
The electronic device obtains measurement metadata, that is, obtains a measurement scale therein, and can determine the dividing interval and the number of target measuring lines according to the measurement scale. The dividing interval is used for controlling the interval between dividing lines, the number of the target measuring lines is used for controlling the grouping of the measuring lines when dividing based on the dividing line positions, so that the number of each group of measuring lines is consistent, the interval of the dividing lines is consistent, and the data quantity deviation between the divided data slices is small. That is, the measurement entity data organized by the line can be divided by equally-spaced dividing lines, thereby obtaining data slices having almost uniform sizes.
It will be appreciated that the measurement scale may be different, and the number of lines in the measurement sample point physical data and the data amount of the measurement data may be different. However, in the embodiment of the present application, the data amounts of slice data of different measurement scale areas may be made consistent by adjusting the dividing interval and the number of target measuring lines.
The number of the target lines can be determined by related technicians according to requirements according to the limitation of network transmission data quantity and the limitation of visual rendering, the number of the target lines can be an empirical value, and the number of the target lines corresponding to the same measurement scale is determined according to requirements and is kept consistent as much as possible. For example, assuming that the separation interval is denoted by L and the number of target lines is denoted by N, this step 202 may be: and the electronic equipment determines a segmentation interval L and the number N of target measuring lines according to the measurement scale in the measurement metadata.
203. And the electronic equipment determines the direction of the parting line according to the direction of the measuring line in the measurement metadata.
The electronic device obtains measurement metadata, that is, obtains the direction of the line therein, and can determine the direction of the dividing line according to the direction of the line. This ensures that the line is split sufficiently to be split into different data slices.
In some embodiments, the line direction is different, and equally spaced warp or weft lines may be selected for segmentation as segmentation lines. Specifically, this step 203 may be implemented by the following steps 1 and 2.
In step 1, the electronic device may determine that the dividing line direction is a meridian direction in response to the measuring direction being greater than or equal to 0 degrees and less than or equal to 45 degrees or the measuring direction being greater than or equal to 135 degrees and less than or equal to 180 degrees in the measurement metadata, where the measuring direction is an angle between the measuring line and the latitude direction.
In this step 1, the dividing line is determined to be the warp direction, that is, the dividing line is warp. The earth was overlooked from above the earth's latitude line by 0 degrees (equator), and the latitude line by 0 degrees (equator) was regarded as the X-axis, the direction from west to east was regarded as the positive X-axis direction, and the direction from east to west was regarded as the negative X-axis direction. The direction of the measuring line is the angle between the measuring line and the positive direction of the X axis. The direction of the dividing line is the same as the direction of the measuring line, and is also the angle between the dividing line and the positive direction of the X axis.
And 2, the electronic equipment responds to the measurement metadata that the direction of the measuring line is greater than 45 degrees and less than 135 degrees, and determines that the direction of the dividing line is the weft direction.
In this step 1, the dividing line is determined to be in the weft direction, that is, the dividing line is determined to be the weft. For the above steps 1 and 2, it can be understood that: assuming that alpha is used for representing the direction of the measuring line, if alpha is more than or equal to 0 degree and less than or equal to 45 degrees or alpha is more than or equal to 135 degrees and less than or equal to 180 degrees, equally spaced warp threads can be used as dividing lines for dividing the measured entity data organized according to the measuring line; if the angle alpha is less than 45 degrees and less than 135 degrees, equally-spaced wefts are used as dividing lines to divide the measured entity data organized by the measuring lines.
It should be noted that, the electronic device may execute step 202 first, and then execute step 203. The electronic device may also execute step 203 first and then execute step 202. The electronic device may also perform step 202 and step 203 simultaneously. The sequence of the step 202 and the step 203 is not limited in this application.
204. The electronic equipment divides the measurement sampling point entity data according to the dividing interval, the number of the target measuring lines and the dividing line direction to obtain a plurality of data slices of the measuring area, wherein each data slice comprises a slice metadata file and a measurement sampling point entity data slice file.
If the electronic device determines the dividing interval, the number of the target measuring lines and the dividing line direction, the measuring lines can be grouped, the position of the dividing line is determined by taking the measuring lines of the number of the target measuring lines as a group, and then the determined dividing line is used for dividing the measured sampling point entity data so as to obtain the data slice.
In some embodiments, with the same measurement scale, different measurement types of zones, the parting line spacing is the same, and the parting line may be calculated starting from 0 ° warp or 0 ° weft. Specifically, this step 204 may be implemented by the following steps one and two.
Step one, the electronic equipment determines the position of the dividing line according to the dividing interval, the number of the target measuring lines and the direction of the dividing line.
The electronic device determining the position of the dividing line can also be understood as determining the dividing line, and determining the position of the dividing line can only determine the position from which to divide the measured sample point entity data.
The dividing interval is used for determining the distance between adjacent dividing lines, so that the dividing line positions are determined through the determined dividing interval, the dividing lines with equal intervals can be determined, data slices with almost the same size are divided, and the situation that a large number of small files are generated can be effectively avoided.
The number of target lines is assumed to be N, that is, the dividing line positions are determined by taking N lines as a group at a time, that is, the dividing is performed by taking N lines as a group.
Specifically, in step 203, different dividing line directions are determined, and dividing line calculation may be performed with respect to the different dividing line directions, respectively, with respect to 0 degree warp and with respect to 0 degree weft. Specifically, the method can comprise the following step (1) and step (2).
And (1) responding to the dividing line direction as the warp direction, starting from 0-degree warp, the electronic equipment acquires the starting point position and the end point position of the measuring lines of the target quantity each time, determining the edge position of the measuring lines of the target quantity in the weft direction according to the starting point position and the end point position, and determining the dividing line position according to the edge position and the dividing interval.
The start position and the end position of the line may be stored in the measurement metadata, and the start position and the end position of the line may be extracted from the measurement metadata when the calculation of the position of the dividing line is determined. The start position and the end position may be represented by means of coordinates. For example, the latitude and longitude of the start point and the latitude and longitude of the end point. The latitude and longitude range of the current measuring line can be determined according to the two positions, if the current measuring line is divided by the meridian dividing line, the latitude and longitude range (edge position) can be determined, and therefore the dividing lines with equal intervals are determined for the areas according to the dividing intervals.
And (2) responding to the direction of the dividing line as the weft direction, starting from 0-degree weft, acquiring the starting point position and the end point position of the measuring lines of the target quantity by the electronic equipment each time, determining the edge position of the measuring lines of the target quantity in the weft direction according to the starting point position and the end point position, and determining the dividing line position according to the edge position and the dividing interval.
The dividing line determination method in the weft direction is the same as the dividing line determination method in the warp direction in the step (1), except that the latitude range is determined according to the starting point position and the end point position.
For example, in one specific example, assuming that the number of target lines is N and the dividing interval is L, if the warp line is taken as the dividing line, the coordinates of the start point and the end point of the 1 st to nth lines are obtained, and the longitude coordinates of the two points at the farthest ends in the X-axis direction are obtained by comparison and are respectively recorded as X 1 And x 2 The longitude of the two outermost dividing lines is calculated asThen m= (lon) 2 -lon 1 ) L+1 dividing lines with longitude values ofWherein L is a division interval, and m is a positive integer. If the weft is used as the dividing line, the coordinates of the start point and the end point of the 1 st to the N th measuring lines are obtained, and the Y coordinates of the two points at the farthest end in the Y-axis direction are obtained by comparison and are respectively marked as Y 1 And y 2 Similarly, the latitude values of all the dividing lines can be calculated.
And secondly, the electronic equipment divides the measurement sampling point entity data based on the dividing line position to obtain a plurality of data slices of the measuring area.
After determining the position of the dividing line in the first step, the electronic device may divide the measurement sampling point entity data according to the dividing line and the number of the target measuring lines, so as to obtain a data slice. The segmentation process can be that the electronic equipment extracts data between two adjacent segmentation lines from the measurement sampling point entity data according to the number of the target measurement lines by grouping according to the segmentation lines.
For example, as shown in fig. 3, taking a dotted line as measurement sampling point entity data organized by a measuring line, a solid line as a dividing line, a measuring line direction being smaller than 45 degrees, and a warp line as the dividing line is determined. Through the first step, the equally-spaced dividing lines are determined to divide the group of measuring lines, and small blocks, namely data slices, can be obtained.
Each of the data slices includes a slice metadata file and a measurement sample point entity data slice file.
The measurement sampling point entity data slice file records the line number and measurement sampling point information contained in the data slice, for example: for aeromagnetometry, the aeromagnetometry data includes coordinate information and attribute information of the measurement sampling points. The attribute information may include, among others, a line number, a reference number, a year of flight, a longitude, a latitude, a magnetic field value, a normal field value, a magnetic daily variable value, a magnetic compensation value, a Δt magnetic anomaly value, an ellipsoidal height, a ground clearance height (true height), and the like. For aviation gamma spectrometry, the aviation gamma spectrometry data comprises coordinate information and attribute information of measurement sampling points. Wherein the attribute information may include line number, reference number, year of flight, longitude, latitude, total window (Tc) count rate, potassium window (K) count rate, uranium window (U) count rate, thorium window (Th) count rate, potassium% 40 K) Content of uranium 238 U) content and thorium 232 Th), content, total amount, ellipsoidal height, height above ground (true height), etc.
The slice metadata file records the original line number of the line contained in the corresponding data slice file, the start point and end point coordinates of the line in the slice, the code of the area, and the like. The start point and end point coordinates of the line are new start point and end point coordinates obtained after the data is sliced, and are not the start point and end point coordinates of the original line obtained in step 201.
205. The electronic device stores a plurality of data slices of the region.
After the electronic equipment is divided to obtain a plurality of data slices, the data slices can be stored, the slicing method can fully ensure that the data volume of each data slice is consistent as much as possible, a large number of small files are avoided, the data volume deviation among different data slices is small, the response time is quick in network transmission and browser loading display, and the loading efficiency is high.
In some embodiments, the electronic device may store a plurality of data slices in its own memory or store the data slices in its own database, or may send the plurality of data slices to a server, where the server stores the data slices in its own database.
In some embodiments, one data slice may be saved as one data slice file. The contents of the data slicing H file may be referred to in the description of step 204, and the specific contents of the data slicing H file may be set by a person skilled in the relevant arts according to the needs, which is not limited in this embodiment of the present application.
In some embodiments, after determining that relevant data of one data slice is stored in the data slice file, when the data slice file is stored, the data slice file can be stored in a multi-layer directory according to a file naming rule, so that when the data slice is extracted subsequently, the corresponding data slice can be retrieved layer by layer for rendering display or sharing.
Specifically, the electronic device may store a plurality of data slices corresponding to the aero-survey entity data according to predefined data storage path setting rules and slice file naming rules, where the data slices include slice metadata files and aero-survey entity data slice files.
The data storage path setting rule and the slice file naming rule may be set by a related technician according to requirements, and a possible implementation manner is provided below, which manner is not limited in the embodiments of the present application.
In some embodiments, the save path of the data slice file may be determined by measuring information such as type, scale encoding, region encoding, and parting line position, and the data slice file may be stored according to the save path.
For example, in one specific example, the save path directory of the data slice file is: the method comprises the steps of encoding a named folder with a measurement type and a scale, encoding a named folder with a measurement area, dividing line coordinates with an H/V direction, dividing block folders with contained line numbers, and a data slice file with V/H direction dividing line coordinates and contained line numbers. Wherein H is used for representing that weft is used as a dividing line, and V is used for representing that warp is used as a dividing line.
For the above coding, the following can be specifically defined:
measurement type and scale coding: the measurement type (also referred to as measurement class) refers to the identification code of different types of measurements, such as aerial magnetometry, aerial gamma spectrometry, etc. The scale mainly comprises 1:100 ten thousand, 1:50 ten thousand, 1:25 ten thousand, 1:20 ten thousand, 1:10 ten thousand, 1:5 ten thousand, 1:2.5 ten thousand, 1:1 ten thousand, different scales can use different letters as scale identification codes.
And (3) coding a measuring region: the method comprises a code of a detection zone and a code of a detection zone division direction, wherein the code is marked as V by dividing warp and is marked as H by dividing weft.
Partition block coding: the divided block comprises all data slices between two adjacent dividing lines, the code of the divided block comprises two parts, the first part is the longitude coordinate value of the dividing line at the left side of the block or the latitude coordinate value of the dividing line at the lower side of the block, three decimal places are reserved by taking radian (°) as a unit, the recording mode is six-bit characters, the first three digits of the six-bit characters are integer parts of the coordinate values, 0 can be supplemented in front when the three digits are less, and 0 can be supplemented behind when the three digits are less than the three digits. The second part is mainly used for marking the line numbers of the first line and the last line in the block, and is marked by L, the line numbers are connected by using "_" and if a certain area is divided by warps, wherein the warps are 110.060 DEG when the left side of the first divided block is divided, 20 lines are counted by the lines 110 to 300, and then the block code is recorded as: v110060_l110_300.
For the data slice file itself, its naming convention may be: the method comprises the steps of reserving three decimal places in radian (°) unit, wherein the recording mode comprises six digits, wherein the first three digits of the six digits are integers, 0 is supplemented before the six digits are less than the three digits, the last three digits of the six digits are decimal parts, and 0 is supplemented after the six digits are less than the three digits, wherein the code comprises a line start-stop number (marked by L) and a dividing line coordinate range in the direction of a non-dividing line where the data slice is located. For example, a certain data slice includes 5 lines 110 to 150, where the number of units of a line number is used to indicate that a line is composed of several segments, and thus, the lines 110 to 150 include five lines, i.e., the line 110, the line 120, the line 130, the line 140, and the line 150. The parting line coordinate range of the non-parting line direction is 45 degrees to 45.12 degrees, and then the slice file is named as: l110_150_h045000_045120.
It should be noted that, the above steps one, two and 205 may be performed simultaneously, for example, in step one, a start position and an end position of the target number of measurement lines are obtained each time, an edge position in a warp or weft direction of the target number of measurement lines is determined according to the start position and the end position, and then, after determining a dividing line position according to the edge position and the dividing interval, step two may be performed for the dividing line position determined this time, the number of measurement lines of the target number of measurement lines in the measurement sampling point entity data are divided, and step 205 is performed after obtaining a data slice corresponding to the group of measurement lines and storing the data slice. And then continuing to execute the step one, obtaining the number of the next group of target measuring lines to further determine the positions of the dividing lines, and then continuing to execute the step two to divide and store the next group of measuring lines, and so on until all the measuring lines are divided.
In a specific possible embodiment, in combination with the above steps 201 to 205, as shown in fig. 4, metadata information of the measurement area may be obtained first, where the metadata information includes a measurement scale, a measurement direction, and coordinates of first and second points of all the traces (i.e. the measurement lines), then, according to the measurement scale of the current measurement area, a data slicing interval L (i.e. a slicing interval) and the number N of each group of measurement lines (i.e. the number of target measurement lines) are determined, a data slicing direction (i.e. a slicing line direction) is determined according to the measurement direction of the current measurement area, then, according to the coordinates of first and second points of the 1 st to nth traces and the determined slicing direction, division lines are calculated, and then, according to the measurement data of the 1 st to nth measurement lines, all measured sample point entity data between two adjacent slicing lines are sequentially extracted. And then storing the data slice files according to the predefined slice file naming rules and the entity file storage paths by the measured sampling point entity data obtained by segmentation. And then, continuously acquiring the measurement data of the next group of N measuring lines, and carrying out data segmentation until the measurement data of all the measuring lines are processed.
The above steps 201 to 205 illustrate a data slicing process and a storage process, and after the data slicing is stored, when the user wants to view the measurement sampling point entity data, the user can reload and render the display. The steps of rendering the display for the data slice are described below.
In some embodiments, when the map is required to be visually rendered and displayed, three levels may be set for rendering, a critical scale may be set, and when the size relationship between the current display scale and the critical scale is different, graphics corresponding to different levels may be selected for rendering and displaying.
It will be appreciated that when the current display scale is relatively small, the area of the field measurement area displayed by the same display window will be relatively large, generally without displaying too much detail. When the current display scale is larger, the area of the field measurement area displayed by the same display window is smaller, so that the content which is generally required to be displayed is particularly detailed, and the data volume to be displayed is larger.
Specifically, in response to a map visualization rendering instruction, displaying a graph corresponding to the size relationship in the plurality of data slices according to the size relationship between the current display scale and the critical scale, wherein the target scale is determined based on the line spacing, the graph comprises a measured area graph, a line graph and a sampling point graph, and the data amounts of different graphs are different.
The area pattern is an area range pattern formed by the head and tail points of the route line. The line pattern is a pattern obtained by thinning the line pattern containing all the navigation lines. The sample point pattern is a measurement sample point pattern comprised by a data slice intersecting the display window.
It will be appreciated that the survey area pattern simply delineates a survey area by connecting the end points of the track lines, excluding specific track line locations and measurement sample point data, etc., and is of relatively small data size. The line pattern is a pattern with all the navigation lines thinned, wherein the line is displayed, the data volume of the line is relatively larger than that of the area pattern, but the measurement sampling point data is not required to be displayed, and the data volume of the line is smaller than that of the measurement sampling point pattern. The measurement sampling point graph needs to display the relevant data of the specific measurement sampling points according to the display window, and the data volume is larger than that of the area graph and the line graph.
In some embodiments, the first critical scale and the second critical scale are set, and in the following, it is specifically described which graphic is displayed when the size relationships among the current display scale, the first critical scale and the second critical scale are which respectively.
Case (1): and responding to the map visual rendering instruction, and displaying the area graph with the current display scale being smaller than or equal to a first critical scale, wherein the first critical scale is determined based on the distance between the measuring lines and a first distance threshold value of the measuring lines displayed on the screen.
Case (2): and responding to the map visual rendering instruction, displaying a line graph, wherein the current display scale is larger than the first critical scale and smaller than or equal to the second critical scale, the data volume of the line graph is larger than the area graph, and the second critical scale is determined based on the line spacing and a second spacing threshold value of the line displayed on a screen.
Case (3): and responding to the map visual rendering instruction, and displaying a measurement sampling point graph with the current display scale larger than the second critical scale, wherein the data volume of the measurement sampling point graph is larger than the line graph and the area graph.
The first spacing threshold and the second spacing threshold can be set and modified by a related technician according to requirements and experience, and the first critical scale and the second critical scale can be set and modified by the related technician according to requirements and experience correspondingly because the first critical scale and the second critical scale are respectively determined by the first spacing threshold and the second spacing threshold.
For example, a first pitch threshold of 1 millimeter (mm) may be set, and the first critical scale may be calculated from the first pitch threshold of 1 millimeter and the line pitch. For another example, the second pitch threshold may be set to 10 millimeters (mm), and then the second critical scale may be calculated based on the second pitch threshold and the line pitch, which is only an exemplary illustration and is not limited in this embodiment.
For example, let d be the first and second pitch thresholds 1 And d2, d 1 <d 2 Then two critical scale bars 1:M can be calculated correspondingly 1 And 1:M 2 . The first critical scale calculation formula is: 1:M 1 Display = d 1/(line spacing x 1000), the second critical scale calculation formula is: 1:M 2 Display = d 2/(line spacing x 1000), where d 1 And d 2 In millimeters and the line spacing in meters.
Thus, when the current display scale is 1:M.ltoreq.1:M 1 Displaying the graph of the area; when 1:M 1 The current display scale is 1:M less than or equal to 1:M 2 Displaying the graph of the measuring line; when 1:M 2 When the current display scale is displayed at 1:M, the measurement sampling point graph is displayed.
That is, when the map is browsed online, three data levels are rendered according to different scales of the current display of the map, the first level directly renders a region range graph formed by the head end point and the tail end point of the route line, the second level renders the thinned line graph, the third level directly obtains measurement sampling points in a data slice intersected with the window, the measurement sampling points are dynamically rendered as required, other data slices do not need to be output, and the workload of data transmission and visual rendering is greatly reduced.
In some embodiments, when loading and rendering a display of multiple data slices, the display speed may be slow due to the large amount of data, resulting in a blank display page. To address this, one possible implementation is provided: in response to a map visualization rendering instruction for any of a plurality of data slices, the electronic device may determine a visualization rendering priority for the any of the plurality of data slices according to an overlapping coverage relationship between a geospatial range of data in the any of the plurality of data slices and a geospatial range of the display window, the visualization rendering priority being positively correlated with the overlapping coverage relationship, i.e., the larger the overlapping coverage area, the higher the visualization rendering priority, the smaller the overlapping coverage area, and the lower the visualization rendering priority. And then the electronic equipment sequentially acquires and loads and renders and displays the plurality of data slices according to the visual rendering priority.
In this embodiment, when data loading and rendering display are performed, for any one data slice, if the display window ranges overlap entirely over the data slice range, the data slice may be loaded preferentially. The data slice may be loaded a second time if the display window range partially overlaps the data slice range. If the display window ranges do not overlap each other to cover the data slice range, then the data slice is not loaded.
It can be understood that: in this embodiment, when data loading is performed, a set of data slices to be rendered are faced, ordered by distance from the center of the screen, and the center slice of the screen is preferentially loaded and rendered. Therefore, the user can be ensured to quickly see the content of the screen center, so that the display requirement of the user is met, the loading according to the requirement is realized, and the user experience is improved.
In some embodiments, when the map is visually rendered, if the current display scale is large, the measurement sampling point entity data in the part of the map may also need to be rendered and displayed. When the measurement sampling points are rendered, the display positions of some measurement sampling points are very small and very close to each other possibly due to the current display scale, so that only the data of the milbex can be seen, and the data can be overlapped, for example, when the current display scale is small, the display positions of a plurality of measurement sampling points on a map are overlapped, and the measurement sampling points are not required to be completely rendered.
In the embodiment of the application, the measurement sampling points can be dynamically aggregated according to the current display scale, namely N coincident measurement sampling points only render one aggregation point or only display one measurement sampling point in the map window, so that the number of the rendered measurement sampling points is greatly reduced, real on-demand rendering is achieved, the visualization efficiency is improved, meanwhile, the situation that local data of a map are too dense to be identified is avoided, and the neatness of the map is ensured.
In some embodiments, the dynamic aggregation process may be implemented by a dynamic aggregation threshold. Specifically, in response to a map visual rendering instruction for any data slice, the electronic device calculates a dynamic aggregation threshold of measurement sampling point entity data in any data slice according to the current display scale and a third spacing threshold of two measurement sampling points displayed on a screen, and performs dynamic aggregation processing on the measurement sampling point entity data in any data slice by taking a test line as a unit according to the dynamic aggregation threshold to obtain target measurement sampling point entity data for visual rendering, wherein the difference between the distance between adjacent measurement sampling points on the same test line in the target measurement sampling point entity data and the third spacing threshold is smaller than a target value, and then performs visual rendering display on the target measurement sampling point entity data.
The third interval threshold may be set by a related technician according to requirements, and the value of the third interval threshold is not limited in the embodiment of the present application. For example, the third pitch threshold may be set to 1mm.
For the target value, the target value may be set by the relevant technician as desired, but in some embodiments is less than or equal to the measurement sample point spacing l×1000/M. l is the measurement sampling point spacing = (average aircraft flight speed x 1000)/(3600 x coordinate sampling rate). Wherein, 1:M is the current display scale of the map. Thus the above difference= (calculated aggregation threshold-rounded aggregation threshold) ×1000×measurement sampling point spacing l/M. The aggregation threshold value after the aggregation threshold value is calculated to be a number greater than or equal to 0 and less than 1, so that the difference value is less than 1000 x the distance l/M between measurement sampling points, the target value can be freely set according to the number relation, and the specific value is not limited in the embodiment of the application.
Wherein the dynamic aggregation threshold is: [ d×M/(1000×l) ].
Where d is a third pitch threshold for two measurement sample points displayed on the screen, and 1 is the measurement sample point pitch, which measurement sample point pitch= (average aircraft flight speed x 1000)/(3600 x coordinate sampling rate). Wherein, 1:M is the current display scale of the map. The unit of the third interval threshold value displayed by the two measurement sampling points on the screen is millimeter, the unit of the interval of the measurement sampling points is meter, the unit of the average flight speed of the aircraft is kilometer per hour, and the unit of the coordinate sampling rate is times per second.
The dynamic aggregation threshold may be set so that the dynamic aggregation processing is not performed when the dynamic aggregation threshold satisfies a condition, and the dynamic aggregation processing may be performed when the dynamic aggregation threshold does not satisfy a condition. Whether the dynamic aggregation threshold meets the condition or not can be set by a relevant technician according to the requirement, for example, when the dynamic aggregation threshold is smaller than or equal to 2, the dynamic aggregation processing is not performed, and when the dynamic aggregation threshold is larger than 2, the dynamic aggregation processing is performed, which is not limited in the embodiment of the present application.
In some embodiments, the number of the aggregated measurement sampling points can be identified on any one of the target measurement sampling points after dynamic aggregation according to the requirement, so that the user can intuitively understand that the currently displayed measurement sampling point is the result after aggregation.
Specifically, the electronic device may display the number of aggregated measurement sampling points and/or the measurement mean of the currently displayed measurement sampling points on the target measurement sampling point. Of course, the electronic device may also display other information on the measurement sampling point, and may perform setting modification according to the requirement, which is not limited in the embodiment of the present application.
When rendering and displaying the data slices, the catalog of the data slices can be determined according to the display requirements of users, and the data slices are extracted from the catalog to be displayed.
The first case, designated the zone.
Specifically, the rendering instruction comprises a region code, a first-level storage catalog and a second-level storage catalog of the data slice are determined according to the region code, a measurement scale and a measurement type, then the position of a display window is obtained, the dividing line position of the data slice intersected with the display window is determined based on the geographic space range of the display window, the measurement sampling point data of the data slice under the catalog corresponding to the dividing line position is extracted based on the dividing line position, and then the measurement sampling point data in the data slice is rendered and displayed.
For example, in one specific example, when the user designates only a region, and does not limit a line, the following steps can be divided:
And a, when the user selects the loaded area code, determining a primary catalog and a secondary catalog stored in the data slice according to the scale information, the measurement type information and the area code information in the current area metadata.
Step b, acquiring a geographical space range from four to four of the screen window, wherein the x-direction range is x 1 To x 2 The y direction range is y 1 To y 2 。
And c, determining slice segmentation interval L and slice segmentation direction H/V information according to the current region metadata information.
Step d, if the dividing direction of the area is V, according to x 1 To x 2 The slice division line coordinates intersecting the current screen window range can be calculated, whereby the three-level directory in which the data slice is stored, that is, the directory including the slice division line coordinates described above in all the three-level directories can be determined irrespective of the line code in the directory. By y 1 To y 2 The coordinate range of the transverse dividing line where the data slice intersected with the current screen is located can be calculated, so that the file name of the data slice can be determined, namely, in the determined three-level directory, all slice files containing the coordinate range of the transverse dividing line determined in the step are not considered in the file name, and the code information of the line in the file name is not considered.
And e, acquiring the position information of the measurement sampling points in the data slice file according to the file catalogue determined in the step a and the step b, and performing visual rendering based on the geographic base map.
During rendering, an aggregation threshold C of the measurement sampling points can be determined according to a current map scale 1:M, the distance between two measurement sampling points 1 and a third distance threshold d of the sampling points displayed on a screen, point aggregation processing is carried out on the measurement sampling points, and then the aggregated points are visualized or the aggregated point connecting lines are visualized, so that the visualization speed is improved.
Wherein the aggregation threshold c= [ (d×m)/(l×1000) ], wherein 1= (aircraft average flight speed×1000)/(3600×coordinate sampling rate). d is in millimeters, 1 unit is in units, the average flight speed of the aircraft is in kilometers per hour, and the coordinate sampling rate is times per second.
In some embodiments, the aggregation process may not be performed when the aggregation threshold satisfies the condition. Specifically, when the aggregation value is less than or equal to 2, point aggregation processing is not performed, and all measurement sampling points are subjected to visual rendering.
Step f, loading attribute information of the measurement sampling points: and extracting corresponding geophysical field information according to the display requirement of a user, loading and displaying by taking the measuring line as a unit in a list form, and establishing a linkage relation between the list attribute record and the measuring point on the graph.
The second case, designated area and line.
When a user designates a loading area and designates a loaded line code, the data loading step is basically consistent with the step of the first condition, but the data loading step is only to screen according to the line code in the file name in the storage directory determined in the step b and the slice file determined in the step c, and the line designated by the user is screened for rendering display.
The above steps 201 to 205 illustrate a data slicing process and a storage process, and after the data slicing is stored, the user may want to share the data slice with other electronic devices, so as to realize data sharing. The following is a description of the steps for data slice sharing.
In some embodiments, when the data sharing service is performed, the user needs measurement sampling point data of aviation magnetic discharge measurement, and intuitively sees the space distribution condition of the data, which is more beneficial to the user to evaluate the data quality and perform data optimization according to the needs of the user. But directly issues the measurement sampling points, not only occupies a large amount of resources of the server, but also has slow browsing and searching speeds and poor user experience.
In some embodiments, the sharing process may share on demand. Specifically, in response to a data slice sharing instruction, the electronic device extracts shared data corresponding to a data extraction range in the plurality of data slices according to the data extraction range in the data slice sharing instruction, integrates measurement sampling point data in the shared data according to a measurement sequence based on a measurement line in the shared data to obtain a shared data file, and then shares the shared data file to a target address to provide distribution and downloading.
Specifically, when data is extracted, according to the data extraction range, relevant measurement data can be extracted, wherein the extraction of the relevant measurement data comprises three parts of data extraction rough measurement and calculation, data accurate measurement and data combination.
Part of, rough measurement: and calculating to obtain the data slice which is intersected with the data extraction x-direction range and the data extraction y-direction range according to the four-to-longitude and latitude values of the data extraction range. The data slice calculated at this time includes two cases: 1. and according to the four-longitude and latitude values of the rectangular range of data extraction, calculating to obtain slices which completely contain the x-direction range and the y-direction range of the data slices, and directly and completely extracting the data slices. 2. The x-direction range and y-direction range of the data extraction range are partially intersected with the x-direction range and y-direction range of the data slice, and for such data slice, further accurate measurement is required, and the following two parts can be referred to. The calculation step may be referred to the steps a, b, and c, and will not be described herein.
Two parts, accurate measurement: aiming at the situation of partial intersection, the start and stop point coordinates of each measuring line in the slicing metadata file are acquired and connected into a line, whether each measuring line is intersected with the data extraction space range or not is calculated, and if so, the measuring sampling points falling into the data extraction range are directly extracted.
Three parts, data combination: all slices extracted from the same area are organized according to the principles of order, continuity and uniqueness, the measuring sampling points in the same measuring line are integrated according to the sequence of acquisition or measurement, and a data file is generated.
It should be noted that the data slicing step, the storing step, and the rendering and displaying step may be performed by different electronic devices. For example, the data slicing step and the storage step may be performed by the electronic device 1, for example, the electronic device 1 is a server, and the server performs data slicing on the collected measurement sampling point entity data and stores the data slicing in its own database or storage device, and the data slicing may be extracted and displayed by any electronic device later, for example, the electronic device 2 may extract a plurality of data slices from the database or storage device of the electronic device 1 for rendering and displaying. For another example, the electronic device 1 may be a terminal, where the terminal may extract the collected measurement sampling point entity data from the server, slice the data to obtain a plurality of data slices, send the plurality of data slices to the server, store the data slices by the server, extract the data slices from the server when rendering is performed subsequently, or store the data slices in its own memory, and load and render the data slices when receiving the rendering instruction subsequently.
The foregoing provides several possible implementations, and the present application does not specifically limit the steps described above to be executed by one or more electronic devices and the type of electronic device.
According to the embodiment of the application, the dividing interval and the number of the target measuring lines are determined according to the measuring scale, so that dividing lines with almost the interval are determined, and dividing is performed based on the same number of the measuring lines, so that the size of each divided data slice is ensured to be consistent as much as possible, a large number of small files are avoided, the deviation of data volume among different data slices is small, the response time of the data slices can be shortened when the data slices are transmitted to a network and loaded and displayed by a browser, and the loading efficiency is improved. In addition, the data of the whole area is only subjected to one-time data slicing treatment, the data of a plurality of layers are not required to be sliced according to different scale levels in advance, repeated storage of the data is avoided, the data is not required to be thinned in advance according to different displayed scales, and the integrity and the accuracy of the measured data are ensured.
All the above optional solutions can be combined to form an optional embodiment of the present application, which is not described in detail herein.
Fig. 5 is a schematic structural diagram of an apparatus for slicing measurement data of an aeronautical magnetic station according to an embodiment of the present application, referring to fig. 5, the apparatus includes:
the acquisition module 501 is configured to acquire measurement data of an airborne magnetic integrated station in a certain area, where the measurement data of the airborne magnetic integrated station includes measurement metadata and measurement sampling point entity data, and the measurement metadata includes at least a measurement scale and a measurement line direction;
a determining module 502, configured to determine a segmentation interval and the number of target measurement lines according to the measurement scale in the measurement metadata;
the determining module 502 is configured to determine a dividing line direction according to the direction of the measuring line in the measurement metadata;
a segmentation module 503, configured to segment the measurement sampling point entity data according to the segmentation interval, the number of target measurement lines, and the direction of the segmentation line, to obtain a plurality of data slices of the region, where each data slice includes a slice metadata file and a measurement sampling point entity data slice file;
a storage module 504 is configured to store a plurality of data slices for the region.
In some embodiments, the determining module 502 is configured to:
responding to the measurement metadata, wherein the direction of the measuring line is greater than or equal to 0 degree and less than or equal to 45 degrees, or the direction of the measuring line is greater than or equal to 135 degrees and less than or equal to 180 degrees, determining the direction of the dividing line as the warp direction, and the direction of the measuring line is the angle between the measuring line and the weft direction;
And determining that the dividing line direction is the weft direction in response to the measuring metadata that the measuring line direction is more than 45 degrees and less than 135 degrees.
In some embodiments, the storage module 504 is configured to:
and storing a plurality of data slices corresponding to the aeronautical measurement entity data according to a predefined data storage path setting rule and a slice file naming rule, wherein the data slices comprise slice metadata files and aeronautical geophysical measurement entity data slice files.
In some embodiments, the apparatus further comprises:
the display module is used for responding to the map visual rendering instruction, displaying a graph corresponding to the size relation according to the size relation between the current display scale and the critical scale, wherein the target scale is determined based on the distance between the measuring lines, the graph comprises a plane graph, a line graph and a point graph, and the data quantity of different graphs is different.
In some embodiments, the display module is to:
responding to a map visual rendering instruction, and displaying a region surface graph when the current display scale is smaller than or equal to a first critical scale, wherein the first critical scale is determined based on the distance between the measuring lines and a first distance threshold value of the measuring lines displayed on a screen;
responding to a map visual rendering instruction, displaying a line graph, wherein the current display scale is larger than a first critical scale and smaller than or equal to a second critical scale, the data volume of the line graph is larger than that of the area graph, and the second critical scale is determined based on the distance between the line and a second distance threshold value of the line displayed on a screen;
And responding to the map visual rendering instruction, and displaying a measurement sampling point graph with the current display scale larger than the second critical scale, wherein the data volume of the measurement sampling point graph is larger than the line graph and the area graph.
In some embodiments, the display module is further to:
in response to a map visualization rendering instruction for any of a plurality of data slices, determining a visualization rendering priority for the any of the plurality of data slices according to an overlapping coverage relationship between a geospatial range of data in the any of the plurality of data slices and a geospatial range of the display window, the visualization rendering priority being positively correlated with the overlapping coverage relationship;
and according to the visual rendering priority, sequentially acquiring loading and rendering and displaying the plurality of data slices.
In some embodiments, the display module is further to:
responding to a map visual rendering instruction of any plurality of data slices, calculating a dynamic aggregation threshold value of the measurement sampling point entity data in any plurality of data slices according to the current display scale and a third interval threshold value of two measurement sampling points displayed on a screen, and carrying out dynamic aggregation processing on the measurement sampling point entity data in any plurality of data slices by taking a measuring line as a unit according to the dynamic aggregation threshold value to obtain target measurement sampling point entity data for visual rendering, wherein the difference value between the distance between adjacent measurement sampling points on the same measuring line in the target measurement sampling point entity data and the third interval threshold value is smaller than a target value;
Performing visual rendering display on the target measurement sampling point entity data;
wherein the dynamic aggregation threshold is: [ dXM/(1000X 1) ];
where d is a third pitch threshold for two measurement sample points displayed on the screen, and l is the measurement sample point pitch, which measurement sample point pitch= (average aircraft flight speed x 1000)/(3600 x coordinate sampling rate). Wherein, 1:M is the current display scale of the map. The unit of the third interval threshold value displayed by the two measurement sampling points on the screen is millimeter, the unit of the interval of the measurement sampling points is meter, the unit of the average flight speed of the aircraft is kilometer per hour, and the unit of the coordinate sampling rate is times per second.
In some embodiments, the display module is further to:
displaying the number of the aggregated measurement sampling points and/or the measurement average value of the currently displayed measurement sampling points on the dynamically aggregated target measurement sampling points of the target measurement sampling point sampling entity data.
In some embodiments, the apparatus further includes a sharing module configured to:
responding to a data slice sharing instruction, extracting shared data corresponding to the data extraction range in the plurality of data slices according to the data extraction range in the data slice sharing instruction, and integrating measurement sampling point data in the shared data according to a measurement sequence by taking a measurement line in the shared data as a reference to obtain a shared data file;
And sharing the shared data file to a target address, and providing distribution and downloading.
In some embodiments, the segmentation module 503 is configured to:
determining a dividing line position according to the dividing interval, the number of the target measuring lines and the dividing line direction;
and dividing the measurement sampling point entity data based on the dividing line position to obtain a plurality of data slices of the measuring area.
In some embodiments, the segmentation module 503 is configured to:
responding to the dividing line direction as the warp direction, starting from 0 degree warp, acquiring the starting point position and the end point position of the measuring lines of the target quantity at each time, determining the edge position of the measuring lines of the target quantity in the weft direction according to the starting point position and the end point position, and determining the dividing line position according to the edge position and the dividing interval;
and responding to the dividing line direction as the weft direction, starting from 0-degree weft, acquiring the starting point position and the ending point position of the measuring lines of the target quantity at each time, determining the edge position in the warp direction in the measuring lines of the target quantity according to the starting point position and the ending point position, and determining the dividing line position according to the edge position and the dividing interval.
According to the device provided by the embodiment of the application, the dividing interval and the number of target measuring lines are determined according to the measuring scale, so that the dividing line with the interval being more or less is determined, and the dividing is performed based on the same number of measuring lines, so that the size of each divided data slice is ensured to be as consistent as possible, a large number of small files are avoided, the data quantity deviation among different data slices is small, the response time of the data slices can be accelerated during network transmission and browser loading display, and the loading efficiency is improved. In addition, the data of the whole area is only subjected to one-time data slicing treatment, the data of a plurality of layers are not required to be sliced according to different scale levels in advance, repeated storage of the data is avoided, the data is not required to be thinned in advance according to different displayed scales, and the integrity and the accuracy of the measured data are ensured.
It should be noted that: when the aeronautical magnetic release comprehensive station measurement data slicing device provided by the embodiment is used for slicing aeronautical measurement data, only the division of the functional modules is used for illustration, and in practical application, the functional allocation can be completed by different functional modules according to the needs, namely, the internal structure of the aeronautical magnetic release comprehensive station measurement data slicing device is divided into different functional modules so as to complete all or part of the functions described above. In addition, the apparatus for slicing measurement data of the integrated station for aviation magnetic discharge and the method for slicing measurement data of the integrated station for aviation magnetic discharge provided in the above embodiments belong to the same concept, and detailed implementation processes of the apparatus are shown in the method embodiments, and are not repeated here.
Fig. 6 is a schematic structural diagram of an electronic device provided in an embodiment of the present application, where the electronic device 600 may have a relatively large difference due to different configurations or performances, and may include one or more processors (Central Processing Units, CPU) 601 and one or more memories 602, where at least one computer program is stored in the memories 602, and the at least one computer program is loaded and executed by the processor 601 to implement the method for slicing measurement data of an aeronautical magnetic amplifying comprehensive station provided in each of the above method embodiments. The electronic device can also include other components for implementing device functions, for example, the electronic device can also have wired or wireless network interfaces, input-output interfaces, and the like for input-output. The embodiments of the present application are not described herein.
In an exemplary embodiment, a computer readable storage medium is also provided, such as a memory comprising at least one computer program executable by a processor to perform the airborne magnetic station survey data slicing method of the above-described embodiment. For example, the computer readable storage medium can be Read-Only Memory (ROM), random access Memory (Random Access Memory, RAM), compact disk Read-Only Memory (Compact Disc Read-Only Memory, CD-ROM), magnetic tape, floppy disk, optical data storage device, etc.
In an exemplary embodiment, a computer program product or a computer program is also provided, the computer program product or the computer program comprising one or more program codes, the one or more program codes being stored in a computer readable storage medium. The one or more processors of the electronic device are capable of reading the one or more program codes from the computer-readable storage medium, the one or more processors executing the one or more program codes such that the electronic device is capable of performing the airborne magnetic and amplification integrated station measurement data slicing method described above.
It should be understood that, in various embodiments of the present application, the sequence numbers of the foregoing processes do not mean the order of execution, and the order of execution of the processes should be determined by the functions and internal logic thereof, and should not constitute any limitation on the implementation process of the embodiments of the present application.
It should be understood that determining B from a does not mean determining B from a alone, but can also determine B from a and/or other information.
Those of ordinary skill in the art will appreciate that all or a portion of the steps implementing the above-described embodiments can be implemented by hardware, or can be implemented by a program instructing the relevant hardware, and the program can be stored in a computer readable storage medium, and the above-mentioned storage medium can be a read-only memory, a magnetic disk or an optical disk, etc.
The foregoing description is only of alternative embodiments of the present application and is not intended to limit the present application, but any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the present application are intended to be included within the scope of the present application.
Claims (9)
1. An aircraft magnetic release integrated station measurement data slicing method, comprising:
acquiring measurement data of an aviation magnetic discharge comprehensive station in a certain area, wherein the measurement data of the aviation magnetic discharge comprehensive station comprises measurement metadata and measurement sampling point entity data, and the measurement metadata at least comprises a measurement scale and a measurement line direction;
Determining a segmentation interval and the number of target measuring lines according to the measurement scale in the measurement metadata, wherein the data volume of slice data of different measurement scale areas is consistent by adjusting the segmentation interval and the number of the target measuring lines;
determining that the dividing line direction is a warp direction in response to the measuring of the metadata, wherein the measuring line direction is greater than or equal to 0 degrees and less than or equal to 45 degrees, or the measuring line direction is greater than or equal to 135 degrees and less than or equal to 180 degrees, and the measuring line direction is an angle between the measuring line and the weft direction;
determining that the dividing line direction is the weft direction in response to the measuring metadata in which the measuring line direction is greater than 45 degrees and less than 135 degrees;
dividing the measurement sampling point entity data according to the dividing interval, the number of the target measuring lines and the dividing line direction to obtain a plurality of data slices of the measuring area, wherein each data slice comprises a slice metadata file and a measurement sampling point entity data slice file;
a plurality of data slices of the region are stored.
2. The method of claim 1, wherein the storing the plurality of data slices of the region comprises:
And storing a plurality of data slices corresponding to the measurement sampling point entity data according to a predefined data storage path setting rule and a slice file naming rule, wherein the data slices comprise slice metadata files and aviation geophysical prospecting measurement entity data slice files.
3. The method according to claim 1, wherein the method further comprises:
responding to a map visual rendering instruction, and displaying a region surface graph when the current display scale is smaller than or equal to a first critical scale, wherein the first critical scale is determined based on the distance between the measuring lines and a first distance threshold value of the measuring lines displayed on a screen;
responding to a map visual rendering instruction, displaying a line graph, wherein the current display scale is larger than a first critical scale and smaller than or equal to a second critical scale, the data volume of the line graph is larger than that of the area graph, and the second critical scale is determined based on the distance between the line and a second distance threshold value of the line displayed on a screen;
and responding to a map visual rendering instruction, and displaying a measurement sampling point graph with the current display scale larger than the second critical scale, wherein the data volume of the measurement sampling point graph is larger than the measuring line graph and the area graph.
4. A method according to claim 3, characterized in that the method further comprises:
in response to a map visualization rendering instruction for any of a plurality of data slices, determining a visualization rendering priority of any of the plurality of data slices according to an overlapping coverage relationship between a geospatial range of data in the any of the plurality of data slices and a geospatial range of a display window, the visualization rendering priority being positively correlated with the overlapping coverage relationship;
and according to the visual rendering priority, sequentially obtaining loading and rendering and displaying the plurality of data slices.
5. The method according to claim 3 or 4, characterized in that the method further comprises:
responding to a map visual rendering instruction of any plurality of data slices, calculating a dynamic aggregation threshold value of measurement sampling point entity data in any plurality of data slices according to the current display scale and a third interval threshold value of two measurement sampling points displayed on a screen, and carrying out dynamic aggregation processing on the measurement sampling point entity data in any plurality of data slices by taking a measuring line as a unit according to the dynamic aggregation threshold value to obtain target measurement sampling point entity data for visual rendering, wherein the difference value between the distance between adjacent measurement sampling points on the same measuring line in the target measurement sampling point entity data and the third interval threshold value is smaller than a target value;
Performing visual rendering display on the target measurement sampling point entity data;
wherein the dynamic aggregation threshold is: [ dXM/(1000×l) ];
wherein d is a third spacing threshold of two measurement sampling points displayed on a screen, l is a measurement sampling point spacing, and the measurement sampling point spacing= (average flight speed of aircraft x 1000)/(3600 x coordinate sampling rate); 1:M is the current display scale of the map; the unit of the third interval threshold value displayed on the screen by the two measurement sampling points is millimeter, the unit of the interval between the measurement sampling points is meter, the unit of the average flight speed of the aircraft is kilometer per hour, and the unit of the coordinate sampling rate is secondary per second.
6. The method according to claim 1, wherein the method further comprises:
responding to a data slice sharing instruction, extracting shared data corresponding to the data extraction range in the plurality of data slices according to the data extraction range in the data slice sharing instruction, and integrating measurement sampling point data in the shared data according to a measurement sequence by taking a measurement line in the shared data as a reference to obtain a shared data file;
and sharing the shared data file to a target address, and providing distribution downloading.
7. An airborne magnetic station survey data slicing apparatus, the apparatus comprising:
the acquisition module is used for acquiring aviation magnetic discharge comprehensive station measurement data of a certain area, wherein the aviation magnetic discharge comprehensive station measurement data comprise measurement metadata and measurement sampling point entity data, and the measurement metadata at least comprise measurement scale and measurement line direction;
the determining module is used for determining a segmentation interval and the number of target measuring lines according to the measurement scale in the measurement metadata, wherein the data volume of slice data of different measurement scale areas is consistent by adjusting the segmentation interval and the number of the target measuring lines;
the determining module is configured to determine, in response to the measurement metadata, that the direction of the line is greater than or equal to 0 degrees and less than or equal to 45 degrees, or that the direction of the line is greater than or equal to 135 degrees and less than or equal to 180 degrees, that the direction of the dividing line is a warp direction, where the direction of the line is an angle between the line and a weft direction; determining that the dividing line direction is the weft direction in response to the measuring metadata in which the measuring line direction is greater than 45 degrees and less than 135 degrees;
the segmentation module is used for segmenting the measurement sampling point entity data according to the segmentation interval, the number of the target measuring lines and the direction of the segmentation line to obtain a plurality of data slices of the measurement area, wherein each data slice comprises a slice metadata file and a measurement sampling point entity data slice file;
And the storage module is used for storing a plurality of data slices of the measuring area.
8. An electronic device comprising one or more processors and one or more memories, the one or more memories having stored therein at least one computer program loaded and executed by the one or more processors to implement the airborne magnetic and amplification integrated station measurement data slicing method of any of claims 1-6.
9. A computer readable storage medium having stored therein at least one computer program loaded and executed by a processor to implement the airborne magnetic and radio integrated station measurement data slicing method of any one of claims 1 to 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310664904.3A CN116755619B (en) | 2023-06-06 | 2023-06-06 | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310664904.3A CN116755619B (en) | 2023-06-06 | 2023-06-06 | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116755619A CN116755619A (en) | 2023-09-15 |
CN116755619B true CN116755619B (en) | 2024-01-05 |
Family
ID=87947037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310664904.3A Active CN116755619B (en) | 2023-06-06 | 2023-06-06 | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116755619B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN118013063A (en) * | 2024-02-05 | 2024-05-10 | 中国自然资源航空物探遥感中心 | Geospatial dynamic partition image data storage and retrieval method, system and equipment |
Citations (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU728105A1 (en) * | 1978-06-02 | 1980-04-15 | Vasilev Roman T | Geomagnetic field measuring method |
JPS6095378A (en) * | 1983-10-31 | 1985-05-28 | Mitsubishi Electric Corp | Remote sensor data obtaining device |
US5555018A (en) * | 1991-04-25 | 1996-09-10 | Von Braun; Heiko S. | Large-scale mapping of parameters of multi-dimensional structures in natural environments |
CN101944132A (en) * | 2010-09-30 | 2011-01-12 | 武汉大学 | Tile map data organization method |
CN102298873A (en) * | 2010-06-24 | 2011-12-28 | 科菱航睿空间信息技术有限公司 | Navigation electronic map data storage method and display method based on the same |
CN102607569A (en) * | 2012-01-17 | 2012-07-25 | 光庭导航数据(武汉)有限公司 | Method for automatically generating data of multiple small scales by using large scale water system in navigation |
RU2592042C1 (en) * | 2015-03-30 | 2016-07-20 | Акционерное общество "Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья" | Method of optimising flight trajectory of movable object during aerogeophysical survey and device for its implementation |
CN106687819A (en) * | 2014-08-29 | 2017-05-17 | 斯布克费舍创新私人有限公司 | An aerial survey image capture system |
CN107229234A (en) * | 2017-05-23 | 2017-10-03 | 深圳大学 | The distributed libray system and method for Aviation electronic data |
CN108227023A (en) * | 2017-12-27 | 2018-06-29 | 彩虹无人机科技有限公司 | CH-3 unmanned plane airborne geophysical prospecting magnetic puts combined station application system |
CN108241712A (en) * | 2016-12-27 | 2018-07-03 | 北京四维图新科技股份有限公司 | A kind of map data processing method and device |
CN109656906A (en) * | 2018-11-21 | 2019-04-19 | 国勘数字地球(北京)科技有限公司 | A kind of processing method of the exploration data based on big data |
CN110308492A (en) * | 2019-07-03 | 2019-10-08 | 核工业航测遥感中心 | Method for extracting aerial geophysical prospecting data in center line of curved railway tunnel |
KR102174524B1 (en) * | 2020-04-07 | 2020-11-04 | 안정훈 | Electronic map service method and system |
CN112069285A (en) * | 2020-11-11 | 2020-12-11 | 湖北亿咖通科技有限公司 | Map generation method and device based on three-dimensional high-precision map slice and electronic equipment |
CN113093296A (en) * | 2021-04-15 | 2021-07-09 | 成都纵横自动化技术股份有限公司 | Aeromagnetic measurement method and related assembly |
CN114528454A (en) * | 2022-02-18 | 2022-05-24 | 中国自然资源航空物探遥感中心 | Method for retrieving and extracting aeromagnetic measurement data |
US11460601B1 (en) * | 2021-10-26 | 2022-10-04 | Beijing Research Institute Of Uranium Geology | Method for determining distribution of uranium deposit |
CN115511999A (en) * | 2022-09-16 | 2022-12-23 | 深圳市赛盈地空技术有限公司 | Line measurement segmentation method, device, equipment and medium |
CN115730175A (en) * | 2021-08-31 | 2023-03-03 | 中国石油天然气集团有限公司 | Ocean and aviation gravity and magnetic data correction method and device |
Family Cites Families (1)
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 |
-
2023
- 2023-06-06 CN CN202310664904.3A patent/CN116755619B/en active Active
Patent Citations (21)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU728105A1 (en) * | 1978-06-02 | 1980-04-15 | Vasilev Roman T | Geomagnetic field measuring method |
JPS6095378A (en) * | 1983-10-31 | 1985-05-28 | Mitsubishi Electric Corp | Remote sensor data obtaining device |
US5555018A (en) * | 1991-04-25 | 1996-09-10 | Von Braun; Heiko S. | Large-scale mapping of parameters of multi-dimensional structures in natural environments |
CN102298873A (en) * | 2010-06-24 | 2011-12-28 | 科菱航睿空间信息技术有限公司 | Navigation electronic map data storage method and display method based on the same |
CN101944132A (en) * | 2010-09-30 | 2011-01-12 | 武汉大学 | Tile map data organization method |
CN102607569A (en) * | 2012-01-17 | 2012-07-25 | 光庭导航数据(武汉)有限公司 | Method for automatically generating data of multiple small scales by using large scale water system in navigation |
CN106687819A (en) * | 2014-08-29 | 2017-05-17 | 斯布克费舍创新私人有限公司 | An aerial survey image capture system |
RU2592042C1 (en) * | 2015-03-30 | 2016-07-20 | Акционерное общество "Сибирский научно-исследовательский институт геологии, геофизики и минерального сырья" | Method of optimising flight trajectory of movable object during aerogeophysical survey and device for its implementation |
CN108241712A (en) * | 2016-12-27 | 2018-07-03 | 北京四维图新科技股份有限公司 | A kind of map data processing method and device |
CN107229234A (en) * | 2017-05-23 | 2017-10-03 | 深圳大学 | The distributed libray system and method for Aviation electronic data |
CN108227023A (en) * | 2017-12-27 | 2018-06-29 | 彩虹无人机科技有限公司 | CH-3 unmanned plane airborne geophysical prospecting magnetic puts combined station application system |
CN109656906A (en) * | 2018-11-21 | 2019-04-19 | 国勘数字地球(北京)科技有限公司 | A kind of processing method of the exploration data based on big data |
CN110308492A (en) * | 2019-07-03 | 2019-10-08 | 核工业航测遥感中心 | Method for extracting aerial geophysical prospecting data in center line of curved railway tunnel |
KR102174524B1 (en) * | 2020-04-07 | 2020-11-04 | 안정훈 | Electronic map service method and system |
CN112069285A (en) * | 2020-11-11 | 2020-12-11 | 湖北亿咖通科技有限公司 | Map generation method and device based on three-dimensional high-precision map slice and electronic equipment |
WO2022100212A1 (en) * | 2020-11-11 | 2022-05-19 | 湖北亿咖通科技有限公司 | Three-dimensional highly precise map tile-based generation method, apparatus, and electronic device |
CN113093296A (en) * | 2021-04-15 | 2021-07-09 | 成都纵横自动化技术股份有限公司 | Aeromagnetic measurement method and related assembly |
CN115730175A (en) * | 2021-08-31 | 2023-03-03 | 中国石油天然气集团有限公司 | Ocean and aviation gravity and magnetic data correction method and device |
US11460601B1 (en) * | 2021-10-26 | 2022-10-04 | Beijing Research Institute Of Uranium Geology | Method for determining distribution of uranium deposit |
CN114528454A (en) * | 2022-02-18 | 2022-05-24 | 中国自然资源航空物探遥感中心 | Method for retrieving and extracting aeromagnetic measurement data |
CN115511999A (en) * | 2022-09-16 | 2022-12-23 | 深圳市赛盈地空技术有限公司 | Line measurement segmentation method, device, equipment and medium |
Non-Patent Citations (4)
Title |
---|
Aeromagnetic Anomaly Detection under Low SNR Conditions Using Multiscale Wavelet Energy Accumulation;Jiaxin Zhou;《2020 IEEE 20th International Conference on Communication Technology》;全文 * |
何昀.基于图像水平集分割的航磁图像信息提取方法研究.《 《中国博士学位论文全文数据库》》.2016,全文. * |
地图切片技术分析与简单实现;赵大龙;孙恒宇;;测绘与空间地理信息(第01期);全文 * |
航空物探测量数据共享技术研究与应用;王林飞;薛典军;李芳;何辉;;物探与化探(第S1期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN116755619A (en) | 2023-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11204896B2 (en) | Scalable space-time density data fusion | |
CN108932303B (en) | Distributed dynamic target detection and analysis system for visible light remote sensing image | |
US10012513B2 (en) | Computing method for ridesharing paths, computing apparatus and recording medium using the same | |
CN116755619B (en) | Method, device, equipment and medium for slicing measurement data of aviation magnetic-release comprehensive station | |
US11954317B2 (en) | Systems and method for a customizable layered map for visualizing and analyzing geospatial data | |
JP2019002913A (en) | Air pollution source recognition based on dynamic pairing and comparison | |
CN110832478A (en) | System and method for on-demand services | |
CN106649653A (en) | Preprocessing method for onboard vector map data | |
Sazhina et al. | Optical analysis of a CMB cosmic string candidate | |
CN110853412B (en) | Method and device for identifying abnormal track point | |
CN113989454B (en) | Fusion method, device and system suitable for geological data and geographic information data | |
CN113139258B (en) | Road data processing method, device, equipment and storage medium | |
CN110321435B (en) | Data source dividing method, device, equipment and storage medium | |
CN114743395A (en) | Signal lamp detection method, device, equipment and medium | |
CN110807954B (en) | Method and device for identifying abnormal track point | |
CN110413662B (en) | Multichannel economic data input system, acquisition system and method | |
CN114078178A (en) | Visualizing weather phenomena along a path in three dimensions | |
CN115048479B (en) | Method, device, system and medium for processing aerograph of airborne mobile terminal | |
CN117349467B (en) | Multi-scale dynamic service method, device and storage medium for remote sensing image map | |
CN113065071B (en) | Product information recommendation method and computer equipment | |
CN113326410B (en) | Method and device for analyzing behavior track of relatives based on space-time association | |
CN116932780B (en) | Astronomical knowledge graph construction method, resource searching method, device and medium | |
CN117951086A (en) | Index construction method, index method and related device | |
CN117831352A (en) | Method and system for predicting take-off airport of airplane, electronic equipment and storage medium | |
Kang et al. | Semi-automatic identification of rock discontinuity orientation based on 3D point clouds and its engineering application |
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 |