CN106611004B - Points of interest attribute display methods based on vector regular quadrangle grid - Google Patents
Points of interest attribute display methods based on vector regular quadrangle grid Download PDFInfo
- Publication number
- CN106611004B CN106611004B CN201510702927.4A CN201510702927A CN106611004B CN 106611004 B CN106611004 B CN 106611004B CN 201510702927 A CN201510702927 A CN 201510702927A CN 106611004 B CN106611004 B CN 106611004B
- Authority
- CN
- China
- Prior art keywords
- grid
- interest
- point
- attribute
- points
- 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/29—Geographical information databases
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Processing Or Creating Images (AREA)
- Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
Abstract
The present invention provides a kind of points of interest attribute display methods based on vector regular quadrangle grid, comprising: to map perspective view carries out grid processing;New figure layer is established in the top of map projection's figure, and navigates to grid Q;Grid Q is encoded;And record the corresponding relationship of grid Q Yu its point of interest P for being included;Spatial index is established as a result, which records the mapping relations between each point of interest and its affiliated grid coding;All mapping relations recorded to spatial index are for statistical analysis, obtain the statistical attribute for all points of interest that same grid is included, and using the statistical attribute of the point of interest as grid attribute;Based on default rendering rule, grid attribute progress classification symbolization is shown.Advantage are as follows: achieve the effect that replace point of interest conventional point component attributes using the face element that grid indicates, therefore, more intuitive can easily observe the attribute information of point of interest, the hiding information of point of interest easy to find and the regularity of distribution etc..
Description
Technical field
The invention belongs to survey and draw technical field of geographic information, and in particular to a kind of interest based on vector regular quadrangle grid
Point attribute display method.
Background technique
Point of interest (point of interest, POI) is a term in GIS-Geographic Information System, and referring to all can be with
It is abstracted as geographic object a little, closely related geographical entity of more especially living with people, for example, school, bank, meal
Shop, gas station, hospital, supermarket etc..Point of interest is mainly used for that the address of things or event is described, can be in very great Cheng
Enhance the descriptive power and query capability to things or event location on degree, improves the accuracy and speed of geo-location.
Currently, point of interest in map projection's figure to show method main are as follows: according to the attribute information of point of interest, to interest
Point carry out symbol rendering, the point of interest of symbolism is then shown in the correspondence geographical location of map projection's figure, with this reach by
Point of interest is presented in the effect in map projection's figure.By taking gas station's point of interest as an example, gas station's interest point symbol can be pre-defined,
For example, dot symbol, triangle point symbol or rectangle symbols etc.;Then, it directly shows and adds in map projection's figure diverse geographic location
Petrol station symbol.As shown in Figure 1, showing schematic diagram for tradition of the point of interest in map projection's figure.
Above-mentioned point of interest shows method in map projection's figure, in conjunction with Fig. 1, mainly has the following deficiencies:
In large scale, such as all gas station's points of interest in the whole world are presented in map projection's figure, because emerging
Interest point quantity is excessive, causes shown interest point symbol very mixed and disorderly, it is difficult to know interest by observing map perspective view
The regularity of point distribution, for example, gas station is in which region gather, and it is sparse etc. in which region.
Summary of the invention
In view of the defects existing in the prior art, the present invention provides a kind of points of interest attribute based on vector regular quadrangle grid
Display methods can effectively solve the above problems.
The technical solution adopted by the invention is as follows:
The present invention provides a kind of points of interest attribute display methods based on vector regular quadrangle grid, comprising the following steps:
Step 1, three-dimensional map is projected to plane, obtains map projection's figure;Wherein, include in map projection's figure
Have it is multiple by the form of show interest point identification, also, the actual geographic coordinate (X, Y) of each point of interest be given value;
Step 2, the grid radius r value of predefined regular quadrangle grid, and the map projection is schemed to carry out at grid
Reason;
Step 3, it in map projection's figure after grid, chooses any point and is defined as initial coordinate origin O, this is initial
The actual geographic coordinate value of coordinate origin O is (X0, Y0);If the grid coordinate value of grid where initial coordinate origin O is (R0,
C0), also, R0=0, C0=0, meaning are as follows: R0The line number of grid, C where representing initial coordinate origin O0Represent initial coordinate
The row number of grid where origin O, it may be assumed that grid where initial coordinate origin O is the grid of the 0th row the 0th column;
Step 4, in map projection's figure after grid, for by the form of show each point of interest, remembered
For point of interest P, it is performed both by following operation:
Step 4.1, the actual geographic coordinate value of point of interest P is (XP, YP) and be given value;
Step 4.2, grid coordinate value (R of the point of interest P relative to initial coordinate origin O is calculatedP, CP);Wherein, RP
Represent grid line number of the point of interest P relative to initial coordinate origin O;CPRepresent lattice of the point of interest P relative to initial coordinate origin O
Train number;The affiliated grid of point of interest P is denoted as grid Q;
Step 4.3, according to the grid line number R of the affiliated grid Q of point of interest PPWith grid row number CP, the lattice of grid Q are calculated
Actual geographic coordinate value (the X of net central point O 'O’, YO’);
Step 4.4, according to the actual geographic coordinate value (X of grid central point O 'O’, YO’), four of grid Q are calculated
The actual geographic coordinate value on vertex;
Step 4.5, new figure layer is established in the top of map projection's figure, and in the new figure layer, navigates to four of grid Q
Grid Q is finally drawn out on vertex;
Meanwhile grid Q is encoded, make it have unique encoded radio;And record the interest that grid Q is included with it
The corresponding relationship of point P;
Step 5, spatial index is established as a result, which records each point of interest and its affiliated grid encodes it
Between mapping relations;
All mapping relations recorded to spatial index are for statistical analysis, and obtain that same grid included is all emerging
The statistical attribute of interest point, and using the statistical attribute of the point of interest as grid attribute;
Step 6, based on default rendering rule, grid attribute progress classification symbolization is shown.
Preferably, in step 1, three-dimensional map is projected to plane, specifically:
Using Mercator projection mode, three-dimensional map is projected to plane.
Preferably, in step 4.2, by following formula, lattice of the point of interest P relative to initial coordinate origin O are calculated
Net coordinate value (RP, CP):
RP=│ Y0-YP│÷2r+1;
CP=│ X0-XP│÷2r+1。
Preferably, in step 4.3, by following formula, the actual geographic that the grid central point O ' of grid Q is calculated is sat
Scale value (XO’, YO’):
XO’=2r CP-r+X0;
YO’=Y0-2r RP+r。
Preferably, in step 5, the statistical attribute for all points of interest that same grid is included includes:
The quantity for all points of interest that same grid is included.
Points of interest attribute display methods provided by the invention based on vector regular quadrangle grid has the advantage that
The present invention counts the statistical attribute for the point of interest that each grid is included first by map projection's map grid networking,
Then, the statistical attribute of point of interest is assigned to corresponding grid, the face element that reaching is indicated using grid is passed instead of point of interest
Therefore the effect of system point component attributes more intuitive can easily observe that the attribute information of point of interest, point of interest easy to find are hidden
Information and the regularity of distribution of hiding etc..
Detailed description of the invention
Fig. 1 is that tradition of the point of interest in map projection's figure shows schematic diagram;
Fig. 2 is that the process of the points of interest attribute display methods provided by the invention based on vector regular quadrangle grid is illustrated
Figure;
Fig. 3 is calibration example of the initial coordinate origin O provided by the invention in map projection's figure;
Fig. 4 is calibration example of any point of interest P provided by the invention in map projection's figure;
Fig. 5 is four vertex M1, M2, M3, M4's provided by the invention that grid Q is calculated according to grid central point O '
The schematic diagram of actual geographic coordinate value;
Fig. 6 is the points of interest attribute display schematic diagram that present invention rendering obtains.
Specific embodiment
In order to which the technical problems, technical solutions and beneficial effects solved by the present invention is more clearly understood, below in conjunction with
Accompanying drawings and embodiments, the present invention will be described in further detail.It should be appreciated that specific embodiment described herein only to
It explains the present invention, is not intended to limit the present invention.
The present invention provides a kind of points of interest attribute display methods based on vector regular quadrangle grid, basic ideas are as follows:
In map projection's figure, grid will be carried out by statistical regions, and count the point of interest that each grid is included
Then the statistical attribute of point of interest is assigned to corresponding grid by statistical attribute, reach and replaced using the face element that grid indicates
The effect of point of interest conventional point component attributes, due to showing the statistical attribute of point of interest using face element, can more intuitive side
Just the attribute information for observing point of interest, the hiding information of point of interest easy to find and the regularity of distribution etc..
It is the process of the points of interest attribute display methods provided by the invention based on vector regular quadrangle grid with reference to Fig. 2
Schematic diagram, specifically includes the following steps:
Step 1, three-dimensional map is projected to plane, obtains map projection's figure;Wherein, include in map projection's figure
Have it is multiple by the form of show interest point identification, also, the actual geographic coordinate (X, Y) of each point of interest be given value;
In this step, Mercator projection mode can be used, three-dimensional map is projected to plane.Mercator projection: again
Name " conic projection of isogonism positive axis ", was drafted by Dutch map scholar Mercator (Mercator) in 1569, it is assumed that the earth is enclosed
In a hollow cylinder, equator is in contact with cylinder, and then imaginary earth center has a lamp again, the figure on spherical surface
Shape projects on cylindrical body, then cylindrical body is unfolded, and here it is " Mercator's throwings that a width standard parallel is zero degree (i.e. equator)
The world map that shadow " is drawn out.
Step 2, the grid radius r value of predefined regular quadrangle grid, and the map projection is schemed to carry out at grid
Reason;
Step 3, it in map projection's figure after grid, chooses any point and is defined as initial coordinate origin O, this is initial
The actual geographic coordinate value of coordinate origin O is (X0, Y0);If the grid coordinate value of grid where initial coordinate origin O is (R0,
C0), also, R0=0, C0=0, meaning are as follows: R0The line number of grid, C where representing initial coordinate origin O0Represent initial coordinate
The row number of grid where origin O, it may be assumed that grid where initial coordinate origin O is the grid of the 0th row the 0th column;
It is calibration example of the initial coordinate origin O in map projection's figure with reference to Fig. 3.Herein, initial coordinate origin O is
A bit arbitrarily selected in map projection's figure.
Step 4, in map projection's figure after grid, for by the form of show each point of interest, remembered
It is calibration example of any point of interest P in map projection's figure with reference to Fig. 4 for point of interest P;It is performed both by following operation:
Step 4.1, the actual geographic coordinate value of point of interest P is (XP, YP) and be given value;
Step 4.2, grid coordinate value (R of the point of interest P relative to initial coordinate origin O is calculatedP, CP);Wherein, RP
Represent grid line number of the point of interest P relative to initial coordinate origin O;CPRepresent lattice of the point of interest P relative to initial coordinate origin O
Train number;The affiliated grid of point of interest P is denoted as grid Q;
Specifically, grid coordinate value of the point of interest P relative to initial coordinate origin O can be calculated by following formula
(RP, CP):
RP=│ Y0-YP│÷2r+1;
CP=│ X0-XP│÷2r+1。
Step 4.3, according to the grid line number R of the affiliated grid Q of point of interest PPWith grid row number CP, the lattice of grid Q are calculated
Actual geographic coordinate value (the X of net central point O 'O’, YO’);
Specifically, the actual geographic coordinate value (X of the grid central point O ' of grid Q by following formula, can be calculatedO’,
YO’):
XO’=2r CP-r+X0;
YO’=Y0-2r RP+r。
Step 4.4, according to the actual geographic coordinate value (X of grid central point O 'O’, YO’), four of grid Q are calculated
The actual geographic coordinate value on vertex;
It can refer to Fig. 5, be the reality that four vertex M1, M2, M3, M4 of grid Q are calculated according to grid central point O '
The schematic diagram of geographic coordinate values, it is only necessary to can be calculated according to trigonometric function relationship, this will not be detailed here.
Step 4.5, new figure layer is established in the top of map projection's figure, and in the new figure layer, navigates to four of grid Q
Grid Q is finally drawn out on vertex;
Meanwhile grid Q is encoded, make it have unique encoded radio;And record the interest that grid Q is included with it
The corresponding relationship of point P;
Step 5, spatial index is established as a result, which records each point of interest and its affiliated grid encodes it
Between mapping relations;
All mapping relations recorded to spatial index are for statistical analysis, and obtain that same grid included is all emerging
The statistical attribute of interest point, for example, the quantity for all points of interest that same grid is included, and the statistical attribute of the point of interest is made
For grid attribute;
Step 6, based on default rendering rule, grid attribute progress classification symbolization is shown.
For example, for different grid different colors can be rendered into, substantially according to the point of interest quantity that it is included
Thought can be with are as follows: the point of interest quantity for being included is more, and color is deeper.With reference to Fig. 6, the point of interest category obtained for present invention rendering
Property display schematic diagram.
It can be seen that the points of interest attribute display methods provided by the invention based on vector regular quadrangle grid, have with
Lower advantage:
The present invention counts the statistical attribute for the point of interest that each grid is included first by map projection's map grid networking,
Then, the statistical attribute of point of interest is assigned to corresponding grid, the face element that reaching is indicated using grid is passed instead of point of interest
Therefore the effect of system point component attributes more intuitive can easily observe that the attribute information of point of interest, point of interest easy to find are hidden
Information and the regularity of distribution of hiding etc..
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also answered
Depending on protection scope of the present invention.
Claims (5)
1. a kind of points of interest attribute display methods based on vector regular quadrangle grid, which comprises the following steps:
Step 1, three-dimensional map is projected to plane, obtains map projection's figure;It wherein, include more in map projection's figure
It is a by the form of show interest point identification, also, the actual geographic coordinate (X, Y) of each point of interest be given value;
Step 2, the grid radius r value of predefined regular quadrangle grid, and the map projection is schemed to carry out grid processing;Its
In, grid radius r value is the 1/2 of the side length of regular quadrangle grid;
Step 3, it in map projection's figure after grid, chooses any point and is defined as initial coordinate origin O, the initial coordinate
The actual geographic coordinate value of origin O is (X0, Y0);If the grid coordinate value of grid where initial coordinate origin O is (R0, C0),
Also, R0=0, C0=0, meaning are as follows: R0The line number of grid, C where representing initial coordinate origin O0Represent initial coordinate origin
The row number of grid where O, it may be assumed that grid where initial coordinate origin O is the grid of the 0th row the 0th column;
Step 4, in map projection's figure after grid, for by the form of show each point of interest, be denoted as emerging
Interesting point P, is performed both by following operation:
Step 4.1, the actual geographic coordinate value of point of interest P is (XP, YP) and be given value;
Step 4.2, grid coordinate value (R of the point of interest P relative to initial coordinate origin O is calculatedP, CP);Wherein, RPIt represents
Grid line number of the point of interest P relative to initial coordinate origin O;CPThe grid that point of interest P is represented relative to initial coordinate origin O arranges
Number;The affiliated grid of point of interest P is denoted as grid Q;
Step 4.3, according to the grid line number R of the affiliated grid Q of point of interest PPWith grid row number CP, it is calculated in the grid of grid Q
Actual geographic coordinate value (the X of heart point O 'O’, YO’);
Step 4.4, according to the actual geographic coordinate value (X of grid central point O 'O’, YO’), four vertex of grid Q are calculated
Actual geographic coordinate value;
Step 4.5, new figure layer is established in the top of map projection's figure, and in the new figure layer, navigates to four tops of grid Q
Point finally draws out grid Q;
Meanwhile grid Q is encoded, make it have unique encoded radio;And record the point of interest P that grid Q is included with it
Corresponding relationship;
Step 5, spatial index is established as a result, which records between each point of interest and its affiliated grid coding
Mapping relations;
All mapping relations recorded to spatial index are for statistical analysis, obtain all points of interest that same grid is included
Statistical attribute, and using the statistical attribute of the point of interest as grid attribute;
Step 6, based on default rendering rule, grid attribute progress classification symbolization is shown.
2. the points of interest attribute display methods according to claim 1 based on vector regular quadrangle grid, which is characterized in that
In step 1, three-dimensional map is projected to plane, specifically:
Using Mercator projection mode, three-dimensional map is projected to plane.
3. the points of interest attribute display methods according to claim 1 based on vector regular quadrangle grid, which is characterized in that
In step 4.2, by following formula, grid coordinate value (R of the point of interest P relative to initial coordinate origin O is calculatedP, CP):
RP=│ Y0-YP│÷2r+1;
CP=│ X0-XP│÷2r+1。
4. the points of interest attribute display methods according to claim 1 based on vector regular quadrangle grid, which is characterized in that
In step 4.3, by following formula, the actual geographic coordinate value (X of the grid central point O ' of grid Q is calculatedO’, YO’):
XO’=2r CP-r+X0;
YO’=Y0-2r RP+r。
5. the points of interest attribute display methods according to claim 1 based on vector regular quadrangle grid, which is characterized in that
In step 5, the statistical attribute for all points of interest that same grid is included includes:
The quantity for all points of interest that same grid is included.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510702927.4A CN106611004B (en) | 2015-10-26 | 2015-10-26 | Points of interest attribute display methods based on vector regular quadrangle grid |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510702927.4A CN106611004B (en) | 2015-10-26 | 2015-10-26 | Points of interest attribute display methods based on vector regular quadrangle grid |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106611004A CN106611004A (en) | 2017-05-03 |
CN106611004B true CN106611004B (en) | 2019-04-12 |
Family
ID=58613904
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510702927.4A Active CN106611004B (en) | 2015-10-26 | 2015-10-26 | Points of interest attribute display methods based on vector regular quadrangle grid |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106611004B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110390045B (en) * | 2018-04-12 | 2021-12-17 | 腾讯大地通途(北京)科技有限公司 | Interest point recommendation method and device based on location service |
CN108984729B (en) * | 2018-07-12 | 2021-01-22 | 武大吉奥信息技术有限公司 | Method and device for efficiently displaying dynamic legends in WebGIS environment |
CN111047217B (en) * | 2019-12-27 | 2022-05-24 | 中国科学院城市环境研究所 | Urban functional space compactness measuring method combined with electronic map interest points |
CN114297324A (en) * | 2021-11-08 | 2022-04-08 | 广东省国土资源测绘院 | Coordinate conversion method, equipment, medium and product of geographic data |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101154319A (en) * | 2006-09-29 | 2008-04-02 | 爱信艾达株式会社 | Traffic information generating apparatus and traffic information generating method |
CN102639961A (en) * | 2009-10-08 | 2012-08-15 | 泰为信息科技公司 | Navigation system with map compression and method of operation thereof |
CN103119544A (en) * | 2010-05-16 | 2013-05-22 | 诺基亚公司 | Method and apparatus for presenting location-based content |
CN104025163A (en) * | 2011-05-11 | 2014-09-03 | 吉祥地方全球游戏有限公司 | Systems and methods for processing, storing, and displaying map data |
CN104202547A (en) * | 2014-08-27 | 2014-12-10 | 广东威创视讯科技股份有限公司 | Method for extracting target object in projection picture, projection interaction method and system thereof |
-
2015
- 2015-10-26 CN CN201510702927.4A patent/CN106611004B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101154319A (en) * | 2006-09-29 | 2008-04-02 | 爱信艾达株式会社 | Traffic information generating apparatus and traffic information generating method |
CN102639961A (en) * | 2009-10-08 | 2012-08-15 | 泰为信息科技公司 | Navigation system with map compression and method of operation thereof |
CN103119544A (en) * | 2010-05-16 | 2013-05-22 | 诺基亚公司 | Method and apparatus for presenting location-based content |
CN104025163A (en) * | 2011-05-11 | 2014-09-03 | 吉祥地方全球游戏有限公司 | Systems and methods for processing, storing, and displaying map data |
CN104202547A (en) * | 2014-08-27 | 2014-12-10 | 广东威创视讯科技股份有限公司 | Method for extracting target object in projection picture, projection interaction method and system thereof |
Also Published As
Publication number | Publication date |
---|---|
CN106611004A (en) | 2017-05-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106611004B (en) | Points of interest attribute display methods based on vector regular quadrangle grid | |
Claessen et al. | Flexible linked axes for multivariate data visualization | |
Krempel | Network visualization | |
CN106324002A (en) | Carbonatite pore structure characterization method based on rock classification and multi-scale digital cores | |
CN104750854A (en) | Mass three-dimensional laser point cloud compression storage and rapid loading and displaying method | |
Butkiewicz et al. | Multi-focused geospatial analysis using probes | |
CN102722885A (en) | Method for accelerating three-dimensional graphic display | |
US9251610B2 (en) | Location info-graphics visualizations | |
CN108198231A (en) | Electric power GIS vector graphics real-time drawing method, storage medium | |
CN105653881B (en) | Flow field visualization method based on multi-density level | |
Cox et al. | Case study: 3D displays of Internet traffic | |
CN109859109A (en) | A kind of Serial Scale PDF map seamless organization and display methods | |
CN108961406A (en) | Geographical information visualization method, apparatus and user terminal | |
CN106844519A (en) | Remote sensing special topic product method for visualizing based on three-dimensional digital earth | |
Zhou et al. | Topological simplifications of hypergraphs | |
Brainerd et al. | Interactive map projections and distortion | |
CN104504053A (en) | Location three-dimensional graph display method based on streetscape attribute information | |
CN106873931A (en) | Raster data visualization device, method and computer system | |
CN102547270B (en) | Method and system for obtaining availability information of adjacent blocks of current block | |
CN102364522A (en) | Large-scale weather volume data drawing method based on global division | |
CN105183745A (en) | Method based on data visualization analysis | |
Murray et al. | StretchPlot: Interactive visualization of multi-dimensional trajectory data | |
CN109544688A (en) | A kind of volume drawing fusion method based on opacity transfer function | |
Soares et al. | Depicting more information in enriched squarified treemaps with layered glyphs | |
CN108073657A (en) | More attribute mobile object 3D displays |
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 |