CN110473251A - Custom field spatial data area statistics method based on grid spatial index - Google Patents
Custom field spatial data area statistics method based on grid spatial index Download PDFInfo
- Publication number
- CN110473251A CN110473251A CN201910749369.5A CN201910749369A CN110473251A CN 110473251 A CN110473251 A CN 110473251A CN 201910749369 A CN201910749369 A CN 201910749369A CN 110473251 A CN110473251 A CN 110473251A
- Authority
- CN
- China
- Prior art keywords
- grid
- spatial
- index
- custom field
- space
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
Abstract
The custom field spatial data area statistics method based on grid spatial index that the invention discloses a kind of, cartesian coordinate system dichotomy rule, which is based on, including server end establishes Global Grid subdivision model, and encoded after being iterated two points to spatial data, establish space lattice index;Custom field is drawn on map in mobile terminal, and by custom field coordinate and spatial data figure layer ID synchronized upload to server end, the spatial index grid set of the corresponding spatial index level of Custom Space range is calculated according to space lattice index for server end;Server end calculates the spatial data area of the custom field according to spatial index grid set and spatial data figure layer ID, and calculated result is passed back to mobile terminal and show.Its remarkable result is: improving computational efficiency, computational accuracy is high, takes up less resources.
Description
Technical field
The present invention relates to spatial data handling technical fields, and in particular to a kind of based on the customized of grid spatial index
Ranged space data area statistical method.
Background technique
With the development of Mobile GIS, the demand that target figure layer Space Elements are counted after range is drawn in mobile terminal is got over
Come more universal.Current way is uploaded onto the server after the custom field of drafting is converted to the files such as Json, then
It is obtained after counting area after carrying out space intersection operation with target figure layer, calculated result is finally passed back into mobile terminal and is shown
Show.
However, although this method comparison of computational results is accurate, since polygon intersects operation speed with polygon space
The restriction of degree, whole Statistical Speed are slow.For this problem, needing one kind can be towards the custom field sky of mobile application
Between data area faster statistical approach.
Summary of the invention
In view of the deficiencies of the prior art, the object of the present invention is to provide a kind of custom fields based on grid spatial index
Spatial data area statistics method, this method is by establishing grid spatial index to spatial data, by polygon and polygon
Space intersection operation be converted to central point and polygon comprising operation, to improve custom field statistical space data area
Computational efficiency.
In order to achieve the above objectives, The technical solution adopted by the invention is as follows:
A kind of custom field spatial data area statistics method based on grid spatial index, key be include with
Lower step:
Step 1: server end is based on cartesian coordinate system dichotomy rule and establishes Global Grid subdivision model, and to space
Data are encoded after being iterated two points, establish space lattice index;
Step 2: custom field is drawn in mobile terminal on map, and by custom field coordinate and spatial data figure layer ID
The corresponding sky of Custom Space range is calculated according to space lattice index in synchronized upload to server end, server end
Between index level spatial index grid set;
Step 3: server end calculates institute according to the resulting spatial index grid set of step 2 and spatial data figure layer ID
The spatial data area of custom field is stated, and calculated result is passed back into mobile terminal and is shown.
Further, in step 1, the establishment process of the space lattice index is as follows:
Step 1.1: based on cartesian coordinate system dichotomy rule, establish after Global Grid subdivision model to the earth carry out through
Latitude iteration two is divided, and carries out binary coding to each grid;
Step 1.2: being text code by the binary coding transcoding of each level grid, obtain the net of required spatial data
Trellis coding;
Step 1.3: face data being encoded using grid coding, establishes the space lattice index of face data.
Further, the rule of cartesian coordinate system dichotomy described in step 1.1 are as follows: by longitude in (- 180,0) section
Grid coding be 0, by grid coding of the longitude in (0,180) section be 1, by grid of the latitude in (- 90,0) section
It is encoded to 0, is 1 by grid coding of the longitude in (0,90) section.
It further, is the method for text code by binary coding transcoding in step 1.2 are as follows:
Step A: it for the i-th level grid, is combined according to its code length every five, then uses 0 when less than five
Carry out complement code;
Step B: it is ten's digit by each five bit combinations code conversion, and records the length of complement code;
Step C: according to the conversion table of number and text, ten's digit is converted into text, and combine and to form final sky
Between data text code.
Further, in step 1.3 face data coding mode are as follows: first with grid coding to each angle of face data
Point is encoded, and the central point of face data is then calculated using the angular coordinate of face data, is encoded later to center point, shape
It is indexed at the space lattice of face data.
Further, the acquisition process of spatial index grid set described in step 2 is as follows:
Step 2.1: custom field is drawn in mobile terminal on map, and calculates custom field coordinate;
Step 2.2: be by custom field coordinate transformation upload onto the server end after JSON file, while upload need into
The spatial data figure layer ID of row area statistics;
Step 2.3: server end calculates custom field area S according to custom field coordinate0, and utilize formula Si<
S0Spatial index level i is calculated in/P, wherein SiIt is less than S for the i-th level grid area0The maximum mesh area of/P, P are face
Product computational accuracy controlling elements;
Step 2.4: corresponding space is calculated according to Custom Space range and spatial index level i in server end
It indexes grid set Grid [i].
Further, the value of the areal calculation precision controlling factor P is set as 1000.
Further, the spatial data areal calculation process of custom field described in step 3 is as follows:
Step 3.1: according to spatial index grid set Grid [i] and spatial data figure layer ID, correlation space element is indexed,
The spatial index grid IntersectGrid [i] for obtaining correlation space element Feture [i] and intersecting with Space Elements,
Middle i is spatial index level;
Step 3.2: according to the spatial index grid intersected with Space Elements, obtaining the spatial index grid of next level
IntersectGrid[i+1];
Step 3.3: grid element center point is in Feture [i] in statistical space index grid IntersectGrid [i+1]
Spatial index number of grid N;
Step 3.4: according to formula S=N*Si+1The spatial data area that custom field is calculated is S.
Further, next level space lattice indexes if it does not exist in step 3.2, then directly utilizes spatial index grid
IntersectGrid [i] reference area, specific as follows:
Statistical space indexes spatial index net of the grid element center point in Feture [i] in grid IntersectGrid [i]
Lattice quantity N;
According to formula S=N*SiThe spatial data area that custom field is calculated is S.
Remarkable result of the invention is:
1, efficiency is improved:
The present invention can greatly improve the efficiency in mobile terminal custom field statistical space data area, utilize space
Grid index by the space intersection operation of polygon and polygon be converted to central point and polygon comprising operation, improve meter
Calculate efficiency.
2, computational accuracy is high:
The present invention selects index trellis stage while improving computational efficiency, through custom field area, and uses
Small level-one indexes grid reference area, has taken into account the computational accuracy of area statistics.
3, it takes up less resources:
The present invention directly carries out areal calculation using the spatial index that spatial data carries, and does not need individually to establish calculating net
Lattice.
Detailed description of the invention
Fig. 1 is flow chart of the method for the present invention;
Fig. 2 be iteration two in three times after effect picture;
Fig. 3 is the encoding efficiency figure of face adornment data;
Fig. 4 is space lattice index precision schematic diagram;
Fig. 5 is the schematic diagram of spatial index grid set;
Fig. 6 is the schematic diagram of the spatial index grid intersected with Space Elements;
Fig. 7 is the schematic diagram of the spatial index grid of next level;
Fig. 8 is the schematic diagram for calculating index number of grid in space in gained correlation space element.
Specific embodiment
Specific embodiment and working principle of the present invention will be described in further detail with reference to the accompanying drawing.
As shown in Figure 1, a kind of custom field spatial data area statistics method based on grid spatial index, specific to walk
It is rapid as follows:
Step 1: server end is based on cartesian coordinate system dichotomy rule and establishes Global Grid subdivision model, and to space
Data are encoded after being iterated two points, establish space lattice index.
The establishment process of the Global Grid subdivision model is as follows:
Step 1.1: based on cartesian coordinate system dichotomy rule, establish after Global Grid subdivision model to the earth carry out through
Latitude iteration two is divided, and carries out binary coding to each grid;
The cartesian coordinate system dichotomy rule that is to say: if this plane represents the earth, carry out longitude two first
Point, it that is to say that the grid coding by longitude in (- 180,0) section is 0, by grid coding of the longitude in (0,180) section
It is 1;Then it carries out latitude two to divide, that is to say that the grid coding by latitude in (- 90,0) section is 0, by longitude in (0,90)
Grid coding in section is 1, and pair warp and weft degree is iterated two points later, iteration three times after effect it is as shown in Figure 2.
Step 1.2: being text code by the binary coding transcoding of each level grid, obtain the net of required spatial data
Trellis coding:
It is that 0 and 1 binary system constituted is compiled each grid coding of the earth by above-mentioned Global Grid subdivision model
Code.But this application is inconvenient, and first, encode too long, wasting space;Second is exactly coding lookup hell to pay.Therefore it builds
A Global Grid coding rule is found, binary coding is become our literal codes, it is each existing for binary coding to overcome
Kind defect.Detailed process is as follows:
Step A: assuming that some spatial point (i.e. i=19) of the 19th level grid, then it passes through Global Grid subdivision model
The binary coding length formed after being iterated two points is 38, is combined according to its code length according to every five, no
Foot five when then use 0 carry out complement code, that is to say by two 0 by 38 coding covers be 40 encode (it is understood that
Combining digit may be three, four, six, ten etc., and subsequent step is turned using the number and text adaptable with it
Change table);
Step B: it is ten's digit by each five bit combinations code conversion, and records the length of complement code;
Step C: according to the conversion table of number and text, according to the number and the corresponding relationship of character between 0-31, by ten
Binary digits are converted to text, and combine the text code to form final spatial data.The conversion table is as shown in table 1 below:
1 conversion table of table
Step 1.3: face data being encoded using grid coding, establishes the space lattice index of face data: the face
The coding mode of data are as follows: encoded first with each angle point of the grid coding to face data, then utilize face data
Angular coordinate calculates the central point of face data, encodes later to center point, forms the space lattice index of face data, coding
Effect is as shown in Figure 3.
By above step, the space lattice index of face data can be established, and any level, arbitrarily large may be implemented
Small space lattice index, it is as shown in Figure 4 that space lattice indexes precision.10th its error of level is 20km, indicates this grid
There is 40km on maximum side, and the 15th level is 0.61km, and grid is up to 1.22km, and the 20th level error is 19m, the 25th level error
For 0.6m, it is 1cm that the 30th level, which is 0.01m,.
Step 2: custom field is drawn in mobile terminal on map, and by custom field coordinate and spatial data figure layer ID
The corresponding spatial index layer of Custom Space range is calculated according to space lattice index to server end in synchronized upload
The spatial index grid set of grade.
Wherein, the acquisition process of the spatial index grid set is as follows:
Step 2.1: custom field is drawn in mobile terminal on map, and calculates custom field coordinate;
Step 2.2: be by custom field coordinate transformation upload onto the server end after JSON file, while upload need into
The spatial data figure layer ID of row area statistics;
Step 2.3: server end calculates custom field area S according to custom field coordinate0, and utilize formula Si<
S0Spatial index level i is calculated in/P, wherein SiIt is less than S for the i-th level grid area0The maximum mesh area of/P, P are face
Product computational accuracy controlling elements calculate progress, default setting 1000 for control area;If i is greater than maximum mesh level
Lmax, then i=L is enabledmax。
Step 2.4: corresponding space is calculated according to Custom Space range and spatial index level i in server end
It indexes grid set Grid [i], as shown in Figure 5.
Step 3: server end calculates institute according to the resulting spatial index grid set of step 2 and spatial data figure layer ID
The spatial data area of custom field is stated, and calculated result is passed back into mobile terminal and is shown.
The spatial data areal calculation process of the custom field is as follows:
Step 3.1: according to spatial index grid set Grid [i] and spatial data figure layer ID, correlation space element is indexed,
The spatial index grid IntersectGrid [i] for obtaining correlation space element Feture [i] and intersecting with Space Elements, such as
Shown in Fig. 6, wherein i is spatial index level;
Step 3.2: according to the spatial index grid intersected with Space Elements, obtaining the spatial index grid of next level
IntersectGrid [i+1], as shown in Figure 7;
Step 3.3: grid element center point is in Feture [i] in statistical space index grid IntersectGrid [i+1]
Spatial index number of grid N, as shown in Figure 8;
Step 3.4: according to formula S=N*Si+1The spatial data area that custom field is calculated is S.
Further, next level space lattice indexes if it does not exist in step 3.2, then directly utilizes spatial index grid
IntersectGrid [i] reference area, specific as follows:
Statistical space indexes spatial index net of the grid element center point in Feture [i] in grid IntersectGrid [i]
Lattice quantity N;
According to formula S=N*SiThe spatial data area that custom field is calculated is S.
Scheme of the present invention, by the way that the space intersection operation of polygon and polygon is converted to central point and polygon
Include operation, and directly using spatial data carry spatial index carry out areal calculation, do not need individually establish calculating net
Lattice, to improve computational efficiency.
Technical solution provided by the present invention is described in detail above.Specific case used herein is to this hair
Bright principle and embodiment is expounded, method of the invention that the above embodiments are only used to help understand and its
Core concept.It should be pointed out that for those skilled in the art, in the premise for not departing from the principle of the invention
Under, it can be with several improvements and modifications are made to the present invention, these improvement and modification also fall into the protection of the claims in the present invention
In range.
Claims (9)
1. a kind of custom field spatial data area statistics method based on grid spatial index, it is characterised in that including following
Step:
Step 1: server end is based on cartesian coordinate system dichotomy rule and establishes Global Grid subdivision model, and to spatial data
It is encoded after being iterated two points, establishes space lattice index;
Step 2: custom field is drawn in mobile terminal on map, and custom field coordinate is synchronous with spatial data figure layer ID
It is uploaded to server end, the corresponding Spatial Cable of Custom Space range is calculated according to space lattice index in server end
Draw the spatial index grid set of level;
Step 3: for server end according to the resulting spatial index grid set of step 2 and spatial data figure layer ID, calculating is described certainly
The spatial data area of the range of definition, and calculated result is passed back into mobile terminal and is shown.
2. the custom field spatial data area statistics method according to claim 1 based on grid spatial index,
Be characterized in that: in step 1, the establishment process of the space lattice index is as follows:
Step 1.1: based on cartesian coordinate system dichotomy rule, establishing after Global Grid subdivision model and longitude and latitude is carried out to the earth
Iteration two is divided, and carries out binary coding to each grid;
Step 1.2: being text code by the binary coding transcoding of each level grid, the grid for obtaining required spatial data is compiled
Code;
Step 1.3: face data being encoded using grid coding, establishes the space lattice index of face data.
3. the custom field spatial data area statistics method according to claim 2 based on grid spatial index,
It is characterized in that: the rule of cartesian coordinate system dichotomy described in step 1.1 are as follows: compile grid of the longitude in (- 180,0) section
Code is 0, is 1 by grid coding of the longitude in (0,180) section, is 0 by grid coding of the latitude in (- 90,0) section,
It is 1 by grid coding of the longitude in (0,90) section.
4. the custom field spatial data area statistics method according to claim 2 based on grid spatial index,
It is characterized in that: the method in step 1.2 by binary coding transcoding for text code are as follows:
Step A: it for the i-th level grid, is combined according to its code length every five, 0 progress is then used when less than five
Complement code;
Step B: it is ten's digit by each five bit combinations code conversion, and records the length of complement code;
Step C: according to the conversion table of number and text, ten's digit is converted into text, and combine and to form final space number
According to text code.
5. the custom field spatial data area statistics method according to claim 2 or 4 based on grid spatial index,
It is characterized by: in step 1.3 face data coding mode are as follows: each angle point of face data is carried out first with grid coding
Then coding is calculated the central point of face data using the angular coordinate of face data, encoded later to center point, forming face number
According to space lattice index.
6. the custom field spatial data area statistics method according to claim 1 based on grid spatial index,
Be characterized in that: the acquisition process of spatial index grid set described in step 2 is as follows:
Step 2.1: custom field is drawn in mobile terminal on map, and calculates custom field coordinate;
Step 2.2: being end of uploading onto the server after JSON file by custom field coordinate transformation, while uploading and needing the face of progress
The spatial data figure layer ID of product statistics;
Step 2.3: server end calculates custom field area S according to custom field coordinate0, and utilize formula Si<S0/ P meter
Calculation obtains spatial index level i, wherein SiIt is less than S for the i-th level grid area0The maximum mesh area of/P, P are areal calculation
The precision controlling factor;
Step 2.4: corresponding spatial index is calculated according to Custom Space range and spatial index level i in server end
Grid set Grid [i].
7. the custom field spatial data area statistics method according to claim 6 based on grid spatial index,
Be characterized in that: the value of the areal calculation precision controlling factor P is set as 1000.
8. the custom field spatial data area statistics method according to claim 1 based on grid spatial index,
Be characterized in that: the spatial data areal calculation process of custom field described in step 3 is as follows:
Step 3.1: according to spatial index grid set Grid [i] and spatial data figure layer ID, indexing correlation space element, obtain
Correlation space element Feture [i] and the spatial index grid IntersectGrid [i] intersected with Space Elements, wherein i is
Spatial index level;
Step 3.2: according to the spatial index grid intersected with Space Elements, obtaining the spatial index grid of next level
IntersectGrid[i+1];
Step 3.3: statistical space indexes space of the grid element center point in Feture [i] in grid IntersectGrid [i+1]
Index number of grid N;
Step 3.4: according to formula S=N*Si+1The spatial data area that custom field is calculated is S.
9. the custom field spatial data area statistics method according to claim 8 based on grid spatial index,
Be characterized in that: next level space lattice indexes if it does not exist in step 3.2, then directly utilizes spatial index grid
IntersectGrid [i] reference area, specific as follows:
Statistical space indexes spatial index grid number of the grid element center point in Feture [i] in grid IntersectGrid [i]
Measure N;
According to formula S=N*SiThe spatial data area that custom field is calculated is S.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910749369.5A CN110473251B (en) | 2019-08-14 | 2019-08-14 | Self-defined range spatial data area statistical method based on grid spatial index |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910749369.5A CN110473251B (en) | 2019-08-14 | 2019-08-14 | Self-defined range spatial data area statistical method based on grid spatial index |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110473251A true CN110473251A (en) | 2019-11-19 |
CN110473251B CN110473251B (en) | 2022-03-25 |
Family
ID=68510754
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910749369.5A Active CN110473251B (en) | 2019-08-14 | 2019-08-14 | Self-defined range spatial data area statistical method based on grid spatial index |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110473251B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110989886A (en) * | 2019-11-29 | 2020-04-10 | 广州海格星航信息科技有限公司 | Three-dimensional space grid selection method and device based on space map |
CN111143500A (en) * | 2019-12-27 | 2020-05-12 | 中国联合网络通信集团有限公司 | Visualized area calculation method, terminal, control device and storage medium |
CN112069371A (en) * | 2020-08-11 | 2020-12-11 | 中国水产科学研究院东海水产研究所 | Marine ship effective rescue searching method based on gridding neighborhood |
CN112395794A (en) * | 2020-11-17 | 2021-02-23 | 重庆市地理信息和遥感应用中心 | Automatic parameterized slope model construction method based on subdivision technology |
CN115687480A (en) * | 2022-10-31 | 2023-02-03 | 朱俊丰 | Assembled wisdom garden one-picture system |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107886444A (en) * | 2017-10-11 | 2018-04-06 | 南京大学 | A kind of distribution Runoff calculation method in region of no relief under urbanization |
CN109241217A (en) * | 2018-08-14 | 2019-01-18 | 民政部国家减灾中心 | A kind of the condition of a disaster data space statistics, inquiry and methods of exhibiting based on geographic grid |
CN109446943A (en) * | 2019-01-23 | 2019-03-08 | 国家海洋环境监测中心 | A kind of marine reclamation land regionl development overall process remote-sensing monitoring method |
CN109522382A (en) * | 2018-11-14 | 2019-03-26 | 国家基础地理信息中心 | Spatial data gridding statistical method and device |
US20190155865A1 (en) * | 2015-11-12 | 2019-05-23 | International Business Machines Corporation | Grid-based geofence data indexing |
CN109816022A (en) * | 2019-01-29 | 2019-05-28 | 重庆市地理信息中心 | A kind of image-recognizing method based on three decisions and CNN |
CN109977190A (en) * | 2019-04-01 | 2019-07-05 | 中国测绘科学研究院 | Area inquiry processing method and device towards scale vector map data |
CN109992636A (en) * | 2019-03-22 | 2019-07-09 | 中国人民解放军战略支援部队信息工程大学 | Space-time code method, temporal index and querying method and device |
-
2019
- 2019-08-14 CN CN201910749369.5A patent/CN110473251B/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20190155865A1 (en) * | 2015-11-12 | 2019-05-23 | International Business Machines Corporation | Grid-based geofence data indexing |
CN107886444A (en) * | 2017-10-11 | 2018-04-06 | 南京大学 | A kind of distribution Runoff calculation method in region of no relief under urbanization |
CN109241217A (en) * | 2018-08-14 | 2019-01-18 | 民政部国家减灾中心 | A kind of the condition of a disaster data space statistics, inquiry and methods of exhibiting based on geographic grid |
CN109522382A (en) * | 2018-11-14 | 2019-03-26 | 国家基础地理信息中心 | Spatial data gridding statistical method and device |
CN109446943A (en) * | 2019-01-23 | 2019-03-08 | 国家海洋环境监测中心 | A kind of marine reclamation land regionl development overall process remote-sensing monitoring method |
CN109816022A (en) * | 2019-01-29 | 2019-05-28 | 重庆市地理信息中心 | A kind of image-recognizing method based on three decisions and CNN |
CN109992636A (en) * | 2019-03-22 | 2019-07-09 | 中国人民解放军战略支援部队信息工程大学 | Space-time code method, temporal index and querying method and device |
CN109977190A (en) * | 2019-04-01 | 2019-07-05 | 中国测绘科学研究院 | Area inquiry processing method and device towards scale vector map data |
Non-Patent Citations (2)
Title |
---|
CHERYL DESHA 等: "Pivotal principles for digital earth development in the twenty-first century", 《INTERNATIONAL JOURNAL OF DIGITAL EARTH》 * |
赵学胜 等: "地球剖分格网研究进展综述", 《测绘学报》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110989886A (en) * | 2019-11-29 | 2020-04-10 | 广州海格星航信息科技有限公司 | Three-dimensional space grid selection method and device based on space map |
CN110989886B (en) * | 2019-11-29 | 2022-05-27 | 广州海格星航信息科技有限公司 | Three-dimensional space grid selection method and device based on space map |
CN111143500A (en) * | 2019-12-27 | 2020-05-12 | 中国联合网络通信集团有限公司 | Visualized area calculation method, terminal, control device and storage medium |
CN111143500B (en) * | 2019-12-27 | 2023-07-18 | 中国联合网络通信集团有限公司 | Visual area calculation method, terminal, control device and storage medium |
CN112069371A (en) * | 2020-08-11 | 2020-12-11 | 中国水产科学研究院东海水产研究所 | Marine ship effective rescue searching method based on gridding neighborhood |
CN112395794A (en) * | 2020-11-17 | 2021-02-23 | 重庆市地理信息和遥感应用中心 | Automatic parameterized slope model construction method based on subdivision technology |
CN115687480A (en) * | 2022-10-31 | 2023-02-03 | 朱俊丰 | Assembled wisdom garden one-picture system |
Also Published As
Publication number | Publication date |
---|---|
CN110473251B (en) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110473251A (en) | Custom field spatial data area statistics method based on grid spatial index | |
CN103425772B (en) | A kind of mass data inquiry method with multidimensional information | |
CN105786822B (en) | Search space object and the method and device for establishing spatial index in electronic map | |
US9819964B2 (en) | Limited error raster compression | |
CN105138722A (en) | Digital river-lake network based method for dividing water collection unit of river basin of plain river network region | |
CN107193888B (en) | A kind of urban road network model towards lane grade navigator fix | |
CN109524993A (en) | The typical week power output scene generating method of wind-powered electricity generation photovoltaic for Mid-long Term Optimized Scheduling | |
CN107798054A (en) | A kind of range query method and device based on Trie | |
CN106250579A (en) | The method and device that a kind of power circuit path intelligence generates | |
CN101515335A (en) | Color two dimension bar code with high compression ratio Chinese character coding capability and its coding and decoding method | |
CN106909788A (en) | Line of rabbet joint automatic drafting method based on geopotential unit data in computer meteorology software | |
WO2017211150A1 (en) | Processing method and device for storing fingerprint data in library | |
CN110457317A (en) | A kind of Hilbert curve encoding and coding/decoding method based on Status view | |
CN113010620B (en) | Natural resource data index statistical method and system based on geographic multilevel grid | |
CN107273466B (en) | The discrete grid of spherical triangle encodes the fast conversion method to geographic latitude and longitude coordinate | |
CN109325510A (en) | A kind of image characteristic point matching method based on lattice statistical | |
CN114048204A (en) | Beidou grid space indexing method and device based on database inverted index | |
CN110265996A (en) | It is a kind of suitable for photovoltaic/wind power prediction time characteristic scale modeling method | |
CN111314841A (en) | WSN positioning method based on compressed sensing and improved genetic algorithm | |
CN106204719A (en) | Magnanimity model real-time scheduling method in three-dimensional scenic based on two-dimensional neighbourhood retrieval | |
CN108495252A (en) | Indoor positioning network element Optimal Deployment Method based on genetic algorithm and simulated annealing | |
CN103337084B (en) | A kind of line map spot automatic generation method based on man-made features feature | |
CN102932012B (en) | A kind of deletion RS code blind identification method for coding parameters of error-tolerant code | |
CN109547160A (en) | A kind of cyclic shift network code construction method | |
CN109255097B (en) | 3.3471 power function-shaped open channel water delivery section and optimal section solving method |
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 |