CN103984748B - A kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics - Google Patents
A kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics Download PDFInfo
- Publication number
- CN103984748B CN103984748B CN201410225392.1A CN201410225392A CN103984748B CN 103984748 B CN103984748 B CN 103984748B CN 201410225392 A CN201410225392 A CN 201410225392A CN 103984748 B CN103984748 B CN 103984748B
- Authority
- CN
- China
- Prior art keywords
- grid
- subdivision
- pdqg
- coding
- data
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/22—Indexing; Data structures therefor; Storage structures
- G06F16/2228—Indexing structures
- G06F16/2246—Trees, e.g. B+trees
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention provides a kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics, methods described is included：Receive the original sun wind data of three-dimensional solar wind model output；Ecliptic plane or meridian plane are divided with PDQG R grid, so that the PDQG R grid of corresponding resolution grade is built automatically, wherein the PDQG R grid is：When radial direction subdivision level is more than warp-wise subdivision level, continue to carry out independent subdivision in radial direction on the basis of PDQG grid；PDQG R grid is encoded, sub-grid that will be each after subdivision gives unique mark so that the coordinate and coding of grid form one-to-one incidence relation；Calculate the center point coordinate of each PDQG R grids；According to the coordinate for the central point for calculating obtained PDQG R grids, original sun wind data is put into PDQG R grid, after initial data matches a PDQG R grid, then the property value of the original sun wind data grid element center point is assigned to.
Description
Technical field
The present invention relates to the rapid tissue field of Solar-terrestrial space physics mass data, and in particular to a kind of large scale day ground
Solar wind data partition, coding and the access method in space.
Background technology
The development of the development of GIS-Geographic Information System, particularly digital earth technology, having the research of geoscience has by force
The data tissue and management means of power.And in face of fast-developing space science, increase with the magnanimity of space exploration data, such as
Detection data in its solar system is organized and managed for what, how to improve the effect of visualization of data access efficiency and data,
The problem of as putting in the urgent need to address in face of space scientist, with most important theories Research Significance and before being widely applied
Scape.
Due to not yet there is the sun wind data of actual detection at present, then the present invention uses Chinese Space synoptic meteorology country
The interplanetary conservation member scholar who won the first place in provincial imperial examinations (SIP-CESE) the three-dimensional solar wind model output of the sun of SIGMA research groups of key lab exploitation
Data, for the tissue of the actual detection data of solar wind from now on, management, analysis technical support is provided.
Whole solar system space is huge incomparable, big surprising of solar wind data volume, to improve the retrieval of mass data and depositing
Speed is taken, it is necessary to build two-dimentional subdivision model, the data of different level of resolution are generated, so as to improve whole solar system spheroid
The data organizational efficiency in space, so as to reach the purpose that quick even real-time rendering is carried out to mass data.
At present, there are three kinds of typical two-dimensional space subdivision models, be longitude and latitude subdivision model, regular polygon subdivision mould respectively
Type and Voronoi sphere adaptive lattice pessimistic concurrency controls.
It is applied in the subdivision of solar wind, polyhedron subdivision model and adaptive lattice pessimistic concurrency control, polyhedron can be excluded
The coordinate of model calculates relatively complicated, and organizational form is different from existing mapping, remotely-sensed data, integrates difficulty, from
Adapt to the subdivision that subdivision model is single engineer's scale, it is difficult to combine with multiresolution.
The content of the invention
It is an object of the present invention to which to overcome above mentioned problem, the present invention provides a kind of solar wind of large scale solar-terrestrial physics
Data partition, coding and access method.
To achieve the above object, the present invention provides a kind of solar wind data partition of large scale solar-terrestrial physics, encodes and deposit
Method is taken, methods described is included：
Step 101) receive the original sun wind data that three-dimensional solar wind model is exported；
Step 102) ecliptic plane or meridian plane with PDQG-R grid are subjected to subdivision, the LOD for subdivision is set according to demand
Grade, so that the PDQG-R grid of corresponding resolution grade is built automatically, wherein the PDQG-R grid is：When radial direction subdivision layer
It is secondary be more than warp-wise subdivision level when, continue on the basis of PDQG grid radial direction carrying out independent subdivision；
Step 103) PDQG-R grid is encoded, sub-grid that will be each after subdivision gives unique mark, makes
The coordinate and coding for obtaining grid form one-to-one incidence relation；
Step 104) the corresponding coordinate positions of central point of each PDQG-R grids is calculated, the coordinate position is expressed as：
(latitude, longitude, radius)；
Step 105) according to step 104) calculate obtained PDQG-R grids central point coordinate, by original solar wind
Data are put into PDQG-R grid, after initial data matches a PDQG-R grid, then by original sun wind data
Property value be assigned to the grid element center point；
Step 106)
Based on step 103) grid coding retrieve original sun wind data, that is, search the corresponding grid position of the coding
Put, and the contents extraction placed in the grid is come out；Visualization for multiresolution data is shown, provides different resolution
The display directly perceived of the data of grade.
Optionally, above-mentioned steps 102) further include：
Step 102-1) by the subdivision of ecliptic plane progress PDQG grid, comprise the following steps that：
Step 102-1-1) by one point of ecliptic plane it is four, per being partly that longitude range is 0 °~90 ° and radius is about
1AU (mean Earth-Sun distance) a quarter disc；
Step 102-1-2) recursion is carried out to each quadrant face with PDQG grid, the 1st subdivision take quadrant
The midpoint on three sides on face, obtains 3 new points；2 on the waist of quadrant two new points are linked to be a parallel, then by the parallel
Midpoint and another new point be linked to be a radial line, form 1 new sub- quadrant and 2 sub- quadrangles；
Step 102-1-3) the 2nd subdivision is then the recursion of the two types sub-grid to the 1st generation；Antithetical phrase
Triangle carries out subdivision by the method for first time subdivision, and for sub- quadrangle grids, then can use quadtree approach to take four
The midpoint on side obtains 4 new points, and the midpoint and the center of circle on two waists are connected with concentric arc line, the midpoint of two parallels is directly connected
Straight line, then 4 new sub- quadrangles of getting back, and then 1 new sub- triangle and 10 quadrangles are just produced, realize to four points
The subdivision of disc higher resolution；
Step 102-1-4) repeat step 102-1-3), until resolution ratio meets the requirement of setting；
Step 102-1-5) when radial resolving power is more than warp-wise resolution ratio, i.e., radial direction subdivision level is more than warp-wise subdivision layer
When secondary, continue, in radially independent subdivision, to obtain PDQG-R grid on the basis of PDQG grid.
Further alternative, above-mentioned steps 103) specific coding of PDQG-R grid is：
Step 103-1) first by one point for four ecliptic plane according to numbered counter-clockwise 0,1,2,3, determine successively each
Quadrant where quadrant face；
Step 103-2) after the 1st subdivision, for each quadrant face, two, outside quadrangle is according to counterclockwise successively
2,3 are encoded to, Center-triangle is encoded to 0；
Step 103-3) after the 2nd subdivision, Center-triangle coding is identical with the coding method of the 1st subdivision unit, sub
Quadrangle encoded left, upper left are respectively 0,1, and bottom right, upper right are respectively 2,3, and the rest may be inferred, and subdivision level often increases by 1
Layer, the corresponding increase by one of coding；
Step 103-4) diametrically, each PDQG grid will by radially 3 subdivisions, 1 point be 8, then any one
Grid adds extra-code after itself coding, and then identifies the grid particular location after radial direction is segmented；
Wherein, the numbering from the center of circle outward is represented with 0,1 successively, and the 1st subdivision is identified with 0,1, the 2nd subdivision, successively
Represented with 00,01,10,11, the 3rd subdivision is represented with 000,001,010,011,100,101,110,111 successively；
Step 103-5) digit mark radius subdivision position behind cut-off rule, the bits per inch word of remaining Morton code is all not
Quaternary number more than 3, corresponding Morton code bit numbers are longer to be represented by：
Morton=q1q2q3…qnqn+1—p
Wherein, q1It is tetrad quadrant flag, "-" is separator, p is in radially 3 subdivisions on the basis of PDQG grid
Radial direction distinguishing identifier position afterwards.
Optionally, above-mentioned steps 105) further include：
Step 105-1) from initial data, first judge radius, the corresponding radius layer of the grid where finding；Based on
Mesh radius layer, matches longitude range；Specially：
A) difference of the radius of grid element center point and the radius of original sample point is less than the half of the mesh radius span；
B) difference of the longitude of grid element center point and the longitude of original sample point is less than the half of the grid longitude span
If meeting two above condition simultaneously, the corresponding grid of initial data is found；
Step 105-2) match grid after, then the property value of initial data is assigned to the grid element center point, circulates institute
Some initial data, all initial data are mapped in PDQG-R grids one by one.
Optionally, above-mentioned steps 106) further include：
Step 106-1) grid coding of reception is decoded, calculate the corresponding grid lattice coordinate position of coding；
Step 106-2) according to coding, corresponding data are extracted from the grid that subdivision is obtained.
Due to grid coding, each grid has been caused to have unique mark, and each grid has unique coordinate,
Therefore mesh coordinate then can be directly calculated from grid coding, which attribute is housed in the grid while can also retrieve
Data.
The data that the invention described above is used are by three-dimensional solar wind model (Chinese Space synoptic meteorology National Key Laboratory
SIGMA research groups exploitation the interplanetary conservation member scholar who won the first place in provincial imperial examinations (SIP-CESE) the three-dimensional solar wind model of the sun) output data, should
Data are that under polar coordinates, three reference axis are latitude, longitude, radius respectively.The sun wind data has irregular sampling
The characteristics of, in heliocentric spheroid sample space, latitude scope is -90 °~90 °, is sampled 55 times, and longitude range is
0~360 °, sample 80 times, radial direction sample range is about 1AU (astronomical unit, mean Earth-Sun distance, 1AU=
149597870.691km), sample 154 times.
Compared with prior art, the technical advantages of the present invention are that：
Because common longitude and latitude grid is intensive in the two poles of the earth latitude and longitude, the ratio for the grid divided on the two poles of the earth and right ascension does not restrain,
It is serious so as to cause grid to draw effect deformation.To solve this problem, present invention introduces degeneration quaternary tree subdivision thought, for Huang
Two kinds of Typical Planars in road face and meridian plane, according to solar wind data characteristicses, it is proposed that a kind of new PDQG-R grid models, and pin
Corresponding encoding scheme is proposed to the grid model.
The subdivision model not only solves the problem of sun center of mass grid is overstocked well, can also meet and radially divide
Resolution is with passing through (latitude) to the nonsynchronous demand of resolution ratio.And the data of a variety of resolution ratio levels can be provided, considerably improve
Data search efficiency, it is adaptable to the subdivision of solar wind irregular sampling data.
Brief description of the drawings
Fig. 1-a are the solar wind data partition flow charts that the present invention is provided；
Fig. 1-b are the solar wind data visualization flow charts that the present invention is provided；
Fig. 2 is ecliptic plane schematic diagram；
Fig. 3 is meridian plane schematic diagram；
Fig. 4 is three subdivision schematic diagrames of PDQG grid；
Fig. 5 is that PDQG-R grid starves 3 subdivision schematic diagrames；
Fig. 6 is quadrant face code schematic diagram；
Fig. 7 is PDQG grid coding schematic diagrams；
Fig. 8 is 3 subdivision schematic diagrames of radial direction.
Embodiment
The method of the invention is described in detail with reference to the accompanying drawings and examples.
In heliocentric solar-terrestrial physics, in order to study magnetic field, temperature and particle density of dishing out is with away from day heart
The situation of change of distance, and this panel data situation of change, it is necessary to cross sun barycenter two kinds of typical cases cut
Plane (i.e. ecliptic plane, meridian plane) builds plane grid subdivision model, divides data into different resolution ratio levels, and right
Grid is encoded, so as to realize the data tissue of huge circular flat, and then improves the access efficiency of mass data.
The present invention supports magnanimity, the rapid tissue of the sun wind data of large scale solar-terrestrial physics, in solar wind two-dimensional space
Subdivision model on the basis of, set up pyramid model, the data of different resolution grade can be provided, user when in use can be with
The data of adequate resolution are selected according to self-demand.
The major function of this patent includes：
1. the two-dimensional space subdivision model based on solar wind
The distribution characteristics of sun wind data in plane to study solar core, puts down for two kinds of ecliptic plane and meridian plane
The uniformly subdivision in face, for huge circular flat, proposes a kind of new PDQG-R models, i.e. the plane degeneration four of radially subdivision
Tree mesh subdivision model, the subdivision model not only solves the problem of grid is overstocked at the centre of sphere, moreover it is possible to meet radial resolving power
More than the demand of (latitude) to resolution ratio of warp.
2. the corresponding encoding scheme of subdivision model
For PDQG-R grid, encoded with improved degeneration Z curves, propose corresponding encoding scheme,
Morton=q1q2q3…qnqn+1—p
Wherein, q1It is tetrad quadrant flag, "-" is separator, p is in radially 3 subdivisions on the basis of PDQG grid
Radial direction distinguishing identifier position afterwards.
By being counted to the number of grid after 5 subdivisions, contrast as follows：
It can be summarized as, when subdivision level is n, the grid number of a quarter disc subdivision generation is num, then meets
Equation below：
That is num=1+21+…+22n-1n≥1
I.e.
And the grid number of common longitude and latitude model is num=4n(n ＞ 0), then can calculate the compression of PDQG grid number
The limit is
When with traditional longitude and latitude grid subdivision, grid quantity is with 22Speed increases.Quantity is bright at the centre of sphere after degeneration
It is aobvious to reduce, the centre of sphere and the sizing grid scale smaller of periphery.Effectively improve the overstocked shortcoming of grid at day heart so that grid
Size tends to homogenization.
On the basis of PDQG grid, in order to meet radial resolving power ratio through (latitude) Xiang Yao great demand, then radial direction after
Continuous subdivision, forms PDQG-R grid, can retain grid at the advantage of PDQG grid, the centre of sphere and be unlikely to overstocked, and can be maximum
Retain the data characteristicses of radial direction.
Embodiment
Below by taking ecliptic plane as an example, the particle density data instance exported with a certain moment Solar wind model is described in detail
The step of subdivision, coding, retrieval.
Step 101) data of three-dimensional solar wind model output are received, and extract the data of ecliptic plane；
Step 102) ecliptic plane is divided with PDQG-R grid, the LOD grades for subdivision are set according to demand, so that
Automatically the PDQG-R grid of corresponding resolution grade is built, the PDQG-R grid is：
PDQG grid refers to plane degenerate quadtree grid model (Plane Degenerated Quadtree Grid),
It introduces the thought of degeneration quaternary tree, with reference to common longitude and latitude grid, the higher-quality grid model of formation.In order to solve bidimensional
The problem of resolution ratio is inconsistent is spent, when radial direction subdivision level is more than warp-wise subdivision level, is continued on the basis of PDQG grid
Independent subdivision is carried out in radial direction, to solve radial resolving power and the nonsynchronous problem of warp-wise resolution ratio, i.e. PDQG-R grid.
Below exemplified by sun ecliptic plane, the subdivision method of PDQG grid is first introduced, the subdivision method step for intending using is as follows：
1) be first, four by one point of ecliptic plane, per be partly longitude range be 0 °~90 °, radius be about the four of 1AU
/ mono- disc is as shown in Figure 2.
2) recursion is carried out to each quadrant face with PDQG grid, the 1st subdivision take three on quadrant face
The midpoint on side, obtains 3 new points, and 2 new points on the waist of quadrant two are linked to be into a parallel, then by the midpoint of the parallel and separately
New point is linked to be a radial line, forms 1 new sub- quadrant (it is triangle that can be approximately considered) and 2 sub- quadrangles, in detail
See Fig. 4, wherein, what solid line was represented is the result of the 1st subdivision.
3) the 2nd subdivision, then be the recursion of the two types sub-grid to the 1st generation.To sub- triangle by the
The method of subdivision carries out subdivision, and for sub- quadrangle grids, then can use common quadtree approach, takes four edges
Midpoint, obtains 4 new points, the midpoint and the center of circle on two waists is connected with concentric arc line, the midpoint of two parallels is directly connected directly
Line, then 4 new sub- quadrangles of getting back thus produce 1 new sub- triangle and 10 (i.e. 2+4+4) individual quadrangles, real
Now to the subdivision of quadrant face higher resolution.What Fig. 4 chain-dotted lines were represented is the 2nd subdivision.
4) repeat step 3, are required until resolution ratio is met.Wherein, what Fig. 4 dotted lines were represented is the result of the 3rd subdivision.
In order to solve the problem of two dimension resolution ratio are inconsistent, when primarily radial resolution ratio is more than warp-wise resolution ratio, i.e.,
When radial direction subdivision level is more than warp-wise subdivision level, it may be considered that continue on the basis of PDQG grid in radially independent subdivision,
To meet radial resolving power and the nonsynchronous demand of warp-wise resolution ratio.Therefore, PDQG-R grid is herein proposed.
Because quadrant face data dimension is 25*28, it can be seen that radial direction sampling density is the 2 of warp-wise3Times.So,
Continue to can reach resolution requirement in radial direction 3 subdivisions of progress on the basis of PDQG, as shown in Figure 5.
It is the schematic diagram segmented to some grid at radially 3 times on the basis of 3 subdivisions of PDQG grid in Fig. 5, its
In, solid line is radially the 1st subdivision, and chain-dotted line is radially the 2nd subdivision, and dotted line is radially the 3rd subdivision.Thus, one
PDQG grids then become 8 sub-grids.Then, PDQG-R grid is generated.
PDQG-R grid on each level is generated on the basis of the PDQG grid of correspondence level, wherein radially
The number of times of subdivision can also change as needed, be PDQG grid when radial direction subdivision number of times is 0.
The step is the small grid that the huge circular flat of ecliptic plane, meridian plane is divided into different resolution grade, with
It is easy to the organization and management of the sun wind data of the plane.The data or model that can be used for handling actual detection are defeated
The sun wind data gone out.
Step 103) PDQG-R grid is encoded, the coded system is：
The coding is to give unique mark by the sub-grid after each subdivision so that coordinate, attribute number in grid
One-to-one incidence relation is formed according to coding, in order to carry out efficient data directory, the specific coding rule of PDQG-R grid
It is then as follows：
1) first by one point for four ecliptic plane according to numbered counter-clockwise 0,1,2,3, each quadrant face is determined successively
The quadrant at place.
2) after the 1st subdivision, for each quadrant face, two, outside quadrangle according to being encoded to 2 successively counterclockwise,
3, Center-triangle is encoded to 0, it is believed that be that 0,1 merging forms 0, such as Fig. 7 a.
3) after the 2nd subdivision, Center-triangle coding is identical with the coding method of the 1st subdivision unit, and sub- quadrangle is compiled
Code lower-left, upper left are respectively 0,1, and bottom right, upper right are respectively 2,3, and the rest may be inferred.Also, subdivision level often increases by 1 layer, coding
Corresponding increase by one.
4) diametrically, each PDQG grid will be by radially 3 subdivisions, and 1 point is 8, then any one grid can be with
Extra-code is added after itself coding, to identify the grid particular location after radial direction is segmented, as numbering is 010 in Fig. 8
Grid is one of 8 sub-grids being generated after PDQG-R subdivisions.Can be separated with cut-off rule after coding itself, below
Digit mark radius subdivision position, with binary code representation, 1 represents subdivision 1 time, and 2 represent subdivision 2 times, and often pass through one
Secondary subdivision, coding increase by one, the number of times of subdivision is more, and gained subregion is smaller.
Represented successively with 0,1 from the numbering of the center of circle outward, the 1st subdivision is segmented for the 2nd time with 0,1 mark, successively with 00,
01st, 10,11 represent, the 3rd subdivision is represented with 000,001,010,011,100,101,110,111 successively.Can be according to coding
Digit determines the number of times of radially subdivision, and it is dynamic to encode, not only with heredity, and can be fine with multiresolution level
Ground combines.
5) the digit mark radius subdivision position behind cut-off rule, is binary code, the bits per inch word of remaining Morton code is all
3 quaternary number is no more than, corresponding Morton code bit numbers are longer, are represented by：
Morton=q1q2q3…qnqn+1—p
Wherein, q1It is tetrad quadrant flag, "-" is separator, p is in radially 3 subdivisions on the basis of PDQG grid
Radial direction distinguishing identifier position afterwards.
The step is to each grid coding, so that the data in grid are also mapped with coding, so as to support data
Retrieval, supports data tissue, is that data application below or even data analysis both provide support.
Step 104) calculate the corresponding coordinate positions of central point (radius, latitude, longitude) of each PDQG-R grids
(radius, latitude, longitude), here is the calculation of coordinate in the range of [0, pi/2], other quadrants according to
This analogizes.
Radius=maxradius*2^ (- lod) * (i-1)+maxradius*2^ (- lod-1)
Lontitude=0.25*pi/lon_num (i)+(0.5*pi/lon_num (i)) * (j-1)
Latitude=0
Wherein, maxradius is ecliptic plane maximum radius, and lod is the subdivision grade lod values set, and lon_num (i) is
Radially i-th layer of number of grid for being divided, j refers to i-th layer of j-th of grid.
Step 105) initial data is put into PDQG-R grid, that is, calculate the property value of grid：
The centre coordinate contrast of coordinate and PDQG-R grid to each data in initial data, initial data is put into
In the grid nearest from oneself, i.e., each data are only put into a grid nearest from oneself.So may caused result
Have：Some grid only one of which data；Some grids may have multiple data, can finally be calculated according to certain suitable algorithm
Go out the property value of the grid；Some grids do not have data, then are empty, to ensure the integrality and uniformity of data.
1) from initial data, radius is first judged, the mesh radius layer where finding, based on the layer, goes matching
Longitude range, the condition of satisfaction is：
A) difference of the radius of grid element center point and the radius of original sample point is less than the half of the mesh radius span；
B) difference of the longitude of grid element center point and the longitude of original sample point is less than the half of the grid longitude span.
If both meet simultaneously, the corresponding grid of the initial data is found.
2) match after grid, then the property value of initial data is assigned to the grid.All initial data of circulation, will
All initial data are mapped in PDQG-R grids.
With data instance (L, B, R, den) (3.673205e-068.246700e-014.732900e+07
3.585800e-10), unit is respectively (radian, radian, km), with sampled data maximum distance 0.15391E+
09km is basis of the radius as uniformly subdivision, and recurrence subdivision is carried out with the scope.(a solar-terrestrial physics distance is 1AU=
1.469E+8km)。
By above-mentioned steps, can find the corresponding PDQG-R grids of the data is：
002113-1103.673205e-068.246700e-014.732900e+073.585800e-10
It is respectively from left to right：Coding, latitude, longitude, radius, density attributes value
Similarly, the corresponding grid codings of 3.673205e-062.081300e+001.080500e+081.026500e-10
It is：
121221-011
Step 106) data based on the two-dimentional Solar wind model of coding retrieval.
The step is actually decoding process, while the corresponding grid position of the coding is found out, and the grid is put
The contents extraction put comes out.
Exemplified by encoding 033123-100, the corresponding data of the coding are found.
1) decode, analyze the coding, calculate the corresponding grid position of the coding.First 0 expression first quartile, [0, pi/
2] in the range of, 33,123 5 represent that subdivision grade is 5, and radially subdivision coding 100 has 3, represents radially to have segmented 3 times, and be
The 5th layer radially segmented.
33123 corresponding coordinates (latitude, longitude, radius) are calculated below.
1st 3 corresponding grid element center points：
Radius (1)=maxradius*2^ (- 1-1)
Lontitude (1)=0.25*pi/2+ (0.5*pi/2) * (2-1)
2nd 3 corresponding grid element center points：
Radius (2)=Radius (1)+maxradius*2^ (- 2-1)
Lontitude (2)=Lontitude (1)+0.25*pi*2^ (- 2-1)
3rd 1 corresponding grid element center point：
Radius (3)=Radius (2)-maxradius*2^ (- 3-1)
Lontitude (3)=Lontitude (2)+0.25*pi*2^ (- 3-1)
4th 2 corresponding grid element center points：
Radius (4)=Radius (3)+maxradius*2^ (- 4-1)
Lontitude (4)=Lontitude (3) -0.25*pi*2^ (- 4-1)
5th 3 corresponding grid element center points：
Radius (5)=Radius (4)+maxradius*2^ (- 5-1)
Lontitude (5)=Lontitude (4)+0.25*pi*2^ (- 5-1)
Position is segmented plus radial direction, final coordinate value is：
Radius=Radius (5)-maxradius*2^ (- 5-1)+maxradius*2^ (- 5-3-1) * 5
Lontitude=Lontitude (5)
Coordinate value is：3.673205e-061.055379e+001.325670e+08
Due to grid coding, each grid has been caused to have unique mark, and each grid has unique coordinate,
So can then directly retrieve step 104 from grid coding) calculate mesh coordinate, while can also be by encoding this rope
Draw and retrieve have which attribute data in the grid.In a word because of some steps before having, so that the step is most
Simply, projected coordinate is not needed in fact, it is possible to complete data retrieval because coding be it is unique, by coding, coordinate and
Feature data association gets up.But projected coordinate can be used as One function of the invention.
2) according to coding, then corresponding data can be extracted from grid.
It is to have 2 to be encoded to the corresponding initial data of 003031-001：
3.673205e-061.129010e+004.899869e+073.491200e-101.381500e+05
3.673205e-061.129010e+004.899869e+073.649800e-101.391600e+05
It is to have 1 to be encoded to the corresponding initial data of 231121@010：
3.673205E-064.589671E+008.807740E+073.783500E-103.104800E+05
Above-mentioned model and coding are in addition to for retrieval, it may also be used for the visualization of multiresolution data is shown, is provided
The display directly perceived of the data of different resolution grade, is easy to user to get information about data characteristic, so that quick location sense is emerging
The data in interesting region.
It should be noted last that, the above embodiments are merely illustrative of the technical solutions of the present invention and it is unrestricted.Although ginseng
The present invention is described in detail according to embodiment, it will be understood by those within the art that, to the technical side of the present invention
Case is modified or equivalent substitution, and without departure from the spirit and scope of technical solution of the present invention, it all should cover in the present invention
Right among.
Claims (4)
1. a kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics, methods described are included：
Step 101) receive the original sun wind data that three-dimensional solar wind model is exported；
Step 102) ecliptic plane or meridian plane with PDQG-R grid are subjected to subdivision, the LOD grades for subdivision are set according to demand,
So as to build the PDQG-R grid of corresponding resolution grade automatically, wherein the PDQG-R grid is：When radial direction subdivision level is big
When warp-wise subdivision level, continue to carry out independent subdivision in radial direction on the basis of PDQG grid；
Step 103) PDQG-R grid is encoded, sub-grid that will be each after subdivision gives unique mark so that net
The coordinate and coding of lattice form one-to-one incidence relation；
Step 104) the corresponding coordinate positions of central point of each PDQG-R grids is calculated, the coordinate position is expressed as：(latitude
Degree, longitude, radius)；
Step 105) according to step 104) calculate obtained PDQG-R grids central point position coordinates, by original solar wind
Data are put into PDQG-R grid, after initial data matches a PDQG-R grid, then by original sun wind data
Property value be assigned to the grid element center point；
Step 106) be based on step 103) grid coding retrieve original sun wind data, that is, search the corresponding grid of the coding
Position, and the contents extraction placed in the grid is come out；Visualization for multiresolution data is shown, provides different resolutions
The display directly perceived of the data of rate grade；
The step 102) further include：
Step 102-1) by the subdivision of ecliptic plane progress PDQG grid, comprise the following steps that：
Step 102-1-1) by one point of ecliptic plane it is four, per being partly that longitude range is 90 ° and radius is ground average departure on the 1st
From a quarter disc；
Step 102-1-2) recursion is carried out to each quadrant face with PDQG grid, the 1st subdivision is taken on quadrant face
Three sides midpoint, obtain 3 new points；The new point of 2 on the waist of quadrant two is linked to be a parallel, then by the parallel
Point is linked to be a radial line with another new point, forms 1 new sub- quadrant and 2 sub- quadrangles；
Step 102-1-3) the 2nd subdivision is then the recursion of the two types sub-grid to the 1st generation；To sub- triangle
Shape carries out subdivision by the method for first time subdivision, and for sub- quadrangle grids, then can use quadtree approach to take four edges
Midpoint obtains 4 new points, and the midpoint and the center of circle on two waists are connected with concentric arc line, the midpoint of two parallels are connected, then again
4 new sub- quadrangles are obtained, and then just produce 1 new sub- triangle and 10 quadrangles, are realized higher to quadrant face
The subdivision of resolution ratio；
Step 102-1-4) repeat step 102-1-3), until resolution ratio meets the requirement of setting；
Step 102-1-5) when radial resolving power is more than warp-wise resolution ratio, i.e., radial direction subdivision level is more than warp-wise subdivision level
When, continue, in radially independent subdivision, to obtain PDQG-R grid on the basis of PDQG grid.
2. solar wind data partition, coding and the access method of large scale solar-terrestrial physics according to claim 1, its feature
It is, the step 103) specific coding of PDQG-R grid is：
Step 103-1) first by one point for four ecliptic plane according to numbered counter-clockwise 0,1,2,3, each four points are determined successively
Quadrant where disc；
Step 103-2) after the 1st subdivision, for each quadrant face, two, outside quadrangle according to encoding successively counterclockwise
For 2,3, Center-triangle is encoded to 0；
Step 103-3) after the 2nd subdivision, Center-triangle encodes, sub four side identical with the coding method of the 1st subdivision unit
Shape encoded left, upper left are respectively 0,1, and bottom right, upper right are respectively 2,3, and the rest may be inferred, and subdivision level often increases by 1 layer, volume
Code-phase should increase by one；
Step 103-4) diametrically, each PDQG grid will be by radially 3 subdivisions, and 1 point is 8, then any one grid
Extra-code is added after itself coding, and then identifies the grid particular location after radial direction is segmented；
Wherein, the numbering from the center of circle outward is represented with 0,1 successively, and the 1st subdivision is segmented for the 2nd time with 0,1 mark, successively with 00,
01st, 10,11 represent, the 3rd subdivision is represented with 000,001,010,011,100,101,110,111 successively；
Step 103-5) digit mark radius subdivision position behind cut-off rule, the bits per inch word of remaining Morton code is all no more than
3 quaternary number, corresponding Morton code bit numbers are represented by：
Morton=q1q2q3…qnqn+1—p
Wherein, q1Tetrad quadrant flag, "-" is separator, p be on the basis of PDQG grid after radially 3 subdivisions
Radial direction distinguishing identifier position.
3. solar wind data partition, coding and the access method of large scale solar-terrestrial physics according to claim 1, its feature
It is, the step 105) further include：
Step 105-1) from initial data, first judge radius, the corresponding radius layer of the grid where finding；Grid is based on again
Radius layer, matches longitude range；Specially：
A) difference of the radius of grid element center point and the radius of original sample point is less than the half of the mesh radius span；
B) difference of the longitude of grid element center point and the longitude of original sample point is less than the half of the grid longitude span；
If meeting two above condition simultaneously, the corresponding grid of initial data is found；
Step 105-2) match grid after, then the property value of initial data is assigned to the grid element center point, circulated all
Initial data, all initial data are mapped in PDQG-R grid one by one.
4. solar wind data partition, coding and the access method of large scale solar-terrestrial physics according to claim 1, its feature
It is, the step 106) further include：
Step 106-1) grid coding of reception is decoded, calculate the coordinate position of the corresponding grid of coding；
Step 106-2) according to coding, corresponding data are extracted from the grid that subdivision is obtained.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410225392.1A CN103984748B (en) | 2014-05-26 | 2014-05-26 | A kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410225392.1A CN103984748B (en) | 2014-05-26 | 2014-05-26 | A kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103984748A CN103984748A (en) | 2014-08-13 |
| CN103984748B true CN103984748B (en) | 2017-10-03 |
Family
ID=51276721
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410225392.1A Active CN103984748B (en) | 2014-05-26 | 2014-05-26 | A kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103984748B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106021480B (en) * | 2016-05-19 | 2019-09-17 | 福建农林大学 | A kind of parallel spatial division methods and its system based on grid dividing |
| CN110706304B (en) * | 2019-09-23 | 2023-02-03 | 中国科学院海洋研究所 | Visualization algorithm for polar region water vapor transport flux |
| CN111046005B (en) * | 2019-11-01 | 2023-09-01 | 宝略科技(浙江)有限公司 | Gridding coding method for urban three-dimensional live-action data |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102117494A (en) * | 2009-12-30 | 2011-07-06 | 新奥特（北京）视频技术有限公司 | Method and device for delaminating and blocking data in three-dimensional modeling |
| CN102609525A (en) * | 2012-02-10 | 2012-07-25 | 北京大学 | Method for unifying existing longitude and latitude subdividing grids |
-
2014
- 2014-05-26 CN CN201410225392.1A patent/CN103984748B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102117494A (en) * | 2009-12-30 | 2011-07-06 | 新奥特（北京）视频技术有限公司 | Method and device for delaminating and blocking data in three-dimensional modeling |
| CN102609525A (en) * | 2012-02-10 | 2012-07-25 | 北京大学 | Method for unifying existing longitude and latitude subdividing grids |
Non-Patent Citations (2)
| Title |
|---|
| 基于经纬度剖分的空间信息编码研究;李康容，李奇;《四川师范大学学报（自然科学版）》;20100915(第2010年05期);全文 * |
| 网络环境下全球地形数据组织;周勇;《地理空间信息》;20070628(第2007年03期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103984748A (en) | 2014-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102609525B (en) | Method for dissecting and unifying geographic grid based on longitude and latitude | |
| CN106898045B (en) | Large-area true three-dimensional geographic scene self-adaptive construction method based on SGOG tiles | |
| CN103136371A (en) | Subdivision identification generation method and data search method of multi-source space data | |
| CN113010620B (en) | Natural resource data index statistical method and system based on geographic multilevel grid | |
| CN108648271A (en) | A kind of interpolation method generating complicated landform grid model based on GIS data | |
| CN103984748B (en) | A kind of solar wind data partition, coding and the access method of large scale solar-terrestrial physics | |
| Guo et al. | A geographic meshing and coding method based on adaptive Hilbert-Geohash | |
| CN103065361A (en) | Three-dimensional (3d) island sandbox achieving method | |
| CN111475597A (en) | Non-rigid grid coding, spatial object unique identification and query method and device | |
| CN114049462B (en) | Three-dimensional model monomer method and device | |
| Fan et al. | Origin, spatial pattern, and evolution of urban system: Testing a hypothesis of “urban tree” | |
| CN115408382B (en) | Vector space information gridding storage management method for large-scale satellite data | |
| CN101964157B (en) | Sphere rhombic grid recursive subdivision method | |
| CN114090714B (en) | Multi-scale large data query method based on ecological evolution of geographic grid model mining area | |
| CN102074048A (en) | Method for structuring and dispatching digital city model library | |
| CN109885638A (en) | A kind of three-dimensional space indexing means and system | |
| CN109241216B (en) | Method for realizing gridding management by using power grid level storage database architecture | |
| CN105468759B (en) | The frequency spectrum data construction method of space body | |
| CN101763659B (en) | Method for automatically generating large-scale digitization altitude data model based on image set | |
| CN103995877B (en) | Large-scale solar-terrestrial space solar wind data subdivision, encoding and access method | |
| CN116578657A (en) | Geographic position coding and decoding method, device and system | |
| Addae et al. | Enabling geosimulations for global scale: Spherical geographic automata | |
| CN114544286A (en) | Particle random generation method based on particle size distribution | |
| CN114565730A (en) | Method and device for processing map vector sample and computer readable storage medium | |
| Yang et al. | Spatiotemporal simulation of tourist town growth based on the cellular automata model: the case of Sanpo town in Hebei province |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CP01 | Change in the name or title of a patent holder |
Address after: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee after: NATIONAL SPACE SCIENCE CENTER, CAS Address before: 100190 No. two south of Zhongguancun, Haidian District, Beijing 1 Patentee before: Space Science & Applied Research Centre, Chinese Academy of Sciences |
|
| CP01 | Change in the name or title of a patent holder |