CN107945247A - The method that drawing area settles quantity map - Google Patents
The method that drawing area settles quantity map Download PDFInfo
- Publication number
- CN107945247A CN107945247A CN201711096451.XA CN201711096451A CN107945247A CN 107945247 A CN107945247 A CN 107945247A CN 201711096451 A CN201711096451 A CN 201711096451A CN 107945247 A CN107945247 A CN 107945247A
- Authority
- CN
- China
- Prior art keywords
- sampling point
- base station
- grid sampling
- settling amount
- latitude
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W4/00—Services specially adapted for wireless communication networks; Facilities therefor
- H04W4/02—Services making use of location information
- H04W4/021—Services related to particular areas, e.g. point of interest [POI] services, venue services or geofences
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04W—WIRELESS COMMUNICATION NETWORKS
- H04W64/00—Locating users or terminals or network equipment for network management purposes, e.g. mobility management
- H04W64/003—Locating users or terminals or network equipment for network management purposes, e.g. mobility management locating network equipment
Abstract
The present invention discloses a kind of method of drawing area sedimentation quantity map.This method comprises the following steps：Include the following steps：Obtain the settling amount of ground base station and the latitude and longitude coordinates of base station.Determine the Grid Sampling point in sample area and the region, calculate the longitude and latitude of Grid Sampling point.The latitude and longitude coordinates of latitude and longitude coordinates, base station, the settling amount of base station and weighting algorithm based on Grid Sampling point calculate the settling amount of Grid Sampling point.Settling amount drawing area settling amount based on each Grid Sampling point.The present invention is not required a variety of data only to need the latitude and longitude coordinates of the settling amount of ground base station, latitude and longitude coordinates and each Grid Sampling point, so, required parameter is few and algorithm is simple.
Description
Technical field
The present invention relates to the base station settling amount drawing area settling amount of the ground strengthening system positioned based on aeronautical satellite
The method of figure.
Background technology
With the fast development of Chinese big-dipper satellite and internet, the demand of high accuracy positioning is more strong.It is high-precision fixed
There is very important importance in position to unmanned, track level navigation, Navigation of Pilotless Aircraft, precision agriculture etc..However, including
Civil signal positioning accuracy including GPS is usually all at ten meters or so, it is impossible to meets the needs of above-mentioned scene, so needing ground
Base strengthening system or satellite-based augmentation system are modified the positional information of satellite positioning, it is main eliminate include systematic error,
Error including propagated error etc..Ground strengthens technology by establishing fixed base stations on ground come to needing high accuracy positioning to set
It is standby that differential data is provided, have difference resolving ability equipment position can be modified by the differential data of acquisition so as to
Reach high accuracy positioning.
Since the calculating of differential data is calculated by changeless ground base station, so the displacement of base station is to high-precision
The effect of degree positioning plays the role of important.Wherein, the settling amount of base station is exactly an important index, and settling amount refers to base station
Change on elevation.While base station settling amount is obtained, generally require draw sedimentation quantity map come check neighboring area or
It is the settling amount of area-of-interest, so that it is determined that the influence to locating effect.Therefore, the method for drawing area settling amount, which becomes, needs
Ask, and parameter is less required for this method, algorithm is simple.
The content of the invention
The present invention solves the problems, such as be how to design drawing area settle quantity map method and this method needed for parameter it is few
And algorithm simple question.
The application scenarios of the present invention are illustrated first：The service side for providing differential data periodically obtains the position of ground base station
Change is moved, the geographical location of base station is corrected.The region nearer apart from base station, settling amount should be with the settling amount phases of base station
Closing property is larger, and the area coherence apart from base station farther out is less or even unrelated.It can be assumed that centered on base station, with ground
Positional distance is managed to weigh the parameter of correlation, the settling amount of certain point in zoning.If want the settling amount of drawing area
Map, the method that can use grid matrix, using grid overlay area and calculates the grid lattice point i.e. settling amount of sampled point, i.e.,
The sedimentation quantity map in the region can be drawn out.Based on this, the present invention provides a kind of method of drawing area sedimentation quantity map.The party
Method includes the following steps：Include the following steps：Obtain the settling amount of ground base station and the latitude and longitude coordinates of base station；Determine sample region
Grid Sampling point in domain and the region, calculates the longitude and latitude of Grid Sampling point；Latitude and longitude coordinates, base based on Grid Sampling point
Latitude and longitude coordinates, the settling amount of base station and the weighting algorithm stood calculate the settling amount of Grid Sampling point；Based on each Grid Sampling point
Settling amount drawing area settling amount.
Compared with prior art, the present invention has at least the following advantages：
The present invention is not required a variety of data only to need the settling amount of ground base station, latitude and longitude coordinates and each Grid Sampling point
Latitude and longitude coordinates, so, required parameter is few and algorithm is simple.
Brief description of the drawings
Fig. 1 is the flow chart of the method for drawing area settling amount of the present invention；
Fig. 2 is the schematic diagram that Grid Sampling point is chosen by grid.
Embodiment
For technology contents, construction feature, institute's reached purpose and effect that the present invention will be described in detail, below in conjunction with embodiment
And attached drawing is coordinated to be described in detail.
Referring to Fig. 1, the method for drawing area sedimentation quantity map of the present invention is the ground enhancing based on aeronautical satellite positioning
The method of the base station settling amount drawing area settling amount of system, it is more specifically a kind of to be painted using the calculating of ground base station settling amount
The computational methods of area sampling point settling amount processed, include the following steps：
S1：Obtain the settling amount of ground base station and the latitude and longitude coordinates of base station；The step can use the prior art, herein
Repeat no more.
S2：Determine the Grid Sampling point in sample area and the region, calculate the latitude and longitude coordinates of each Grid Sampling point.
The step is specific as follows：First, the sample area in covering and drawing region is determined to, the sample area is in this implementation
It is rectangular region in mode, as shown in Fig. 2, then, determining that the grid is adopted in the sample area according to equidistant step-length
Sampling point, definite Grid Sampling point are as shown in Figure 2.The number of Grid Sampling point influences the performance of sedimentation quantity map.Grid is adopted
Sampling point is more, and map level sense is stronger.The acquisition of Grid Sampling point is the same as the meshgrid functions of Matlab and the Numpy of Python
The mode that meshgrid functions in math library obtain grid matrix is consistent.
S3：The latitude and longitude coordinates of latitude and longitude coordinates, base station, the settling amount of base station and weighting algorithm based on Grid Sampling point
Calculate the settling amount of Grid Sampling point.The step is specific as follows：
First, Grid Sampling point is calculated to the distance of each base station：Use Haversine formula or Vincenty formula
.Simple in order to state, the present invention repeats no more detailed process of the formula according to calculation of longitude ＆ latitude distance.Present invention assumes that make
With Haversine formula calculating distance, the distance of Grid Sampling point and base station is with formulae express：
Wherein, d is the two point distances that are calculated,Respectively Grid Sampling point and the longitude of base station
With latitude coordinate.
Then, influence of each base station to Grid Sampling point settling amount is weighed by using gaussian radial basis function to weigh
Weight, calculation are as follows：
W=e-dε (2)
Wherein, w is the weight calculated, distances of the d by the Grid Sampling point that formula (1) calculates to base station.ε is desirable
The parameter of adjusting,When the parasang kilometer of calculating,εAdjustable range be usually 0.1~10.According to required drafting effect
Debug εParameter, value see that the figure drawn out seems that intuitively just takes that value according to subjective judgement,εEffect be to away from From zooming in and out.No longer adjust the distance in addition to ε and carry out any Nonlinear Processing, so that the measurement being consistent.Below to calculate one
Exemplified by the settling amount of a sampled point i：
Assuming that one to share N number of base station and corresponding settling amount be s in sample arean(1≤n≤N).So calculate every
One base station is for the weight of Grid Sampling point i：
Wherein, dniRepresent Grid Sampling point i to the distance of nth base station.Weight is normalized operation：
Calculate the settling amount s of Grid Sampling point iiFor：
S4：Settling amount drawing area sedimentation quantity map based on each Grid Sampling point.
In this step, there may be sawtooth effect for the map that the sampled point settling amount calculated in aforementioned manners is drawn
Should, it can be carried out with Gaussian filter smoothly.
Assuming that it is [75.0,135.0] to draw longitude range, latitude scope is the regional subsidence quantity map of [15.0,55.0]
Exemplified by the more specific detail above method it is as follows：
Since selected areas is rectangular area, it is possible to directly as sample area.In sample area, step is set
A length of 0.0625, then 616001 Grid Sampling points are produced altogether, are calculated according to the above method of the present invention and are based on base station settling amount
Grid Sampling point settling amount.
It should be noted that, although present embodiment is stated with step S1, S2, S3 and S4, this statement is not to suitable
The restriction of sequence, simply illustrates that this method includes this four steps, its order can adjust, such as, can first step S2, then step
S1, that is to say, that the latitude and longitude coordinates of calculating Grid Sampling point are first, and obtain the latitude and longitude coordinates of base station rear.
In conclusion those skilled in the art is it is recognised that the present invention at least has the following advantages that：
1st, the settling amount, latitude and longitude coordinates and each Grid Sampling point of the invention that are not required a variety of data only to need ground base station
Latitude and longitude coordinates, so, required parameter is few and algorithm is simple.
2nd, the present invention calculates the settling amount of Grid Sampling point using Gaussian radial basis function, and Gaussian radial basis function can be with
Preferable characterization settling amount and the correlation of distance.
3. the parameter that this method needs are adjusted is less, so can easily be accommodated.Only need to contract by adjusting ε and adjusting the distance
Put, it is possible to reach required effect.
Claims (7)
1. a kind of method of drawing area sedimentation quantity map, it is characterized in that：Include the following steps：
Obtain the settling amount of ground base station and the latitude and longitude coordinates of base station；
Determine the Grid Sampling point in sample area and the region, calculate the longitude and latitude of Grid Sampling point；
The latitude and longitude coordinates of latitude and longitude coordinates, base station, the settling amount of base station and weighting algorithm based on Grid Sampling point calculate net
The settling amount of lattice sampled point；
Settling amount drawing area sedimentation quantity map based on each Grid Sampling point.
2. the method for drawing area sedimentation quantity map as claimed in claim 1, it is characterized in that：The net determined in the region
Lattice sampled point includes：The sample area in covering and drawing region is first determined to, it is true in the sample area according to equidistant step-length
The fixed Grid Sampling point.
3. the method for drawing area sedimentation quantity map as claimed in claim 2, it is characterized in that：The sample area is rectangle.
4. the method for drawing area sedimentation quantity map as claimed in claim 1, it is characterized in that：The calculating Grid Sampling point
Settling amount includes the following steps：
According to the latitude and longitude coordinates of Grid Sampling point and the latitude and longitude coordinates of each base station, each Grid Sampling point is calculated to each base
The distance stood；
Weigh influence i.e. weight of each base station to Grid Sampling point settling amount；
The settling amount of each Grid Sampling point is obtained according to weight and distance.
5. the method for drawing area sedimentation quantity map as claimed in claim 4, it is characterized in that：Using gaussian radial basis function
Influence i.e. weight of each base station to Grid Sampling point settling amount is weighed, calculation is as follows：W=e-dε, wherein, w is is calculated
Weight, d is the distance that Grid Sampling point arrives base station, and ε is the parameter that adjusts of needs.
6. the method for drawing area sedimentation quantity map as described in claim 4 or 5, it is characterized in that：It is described according to weight and away from
Include from the settling amount for obtaining each Grid Sampling point：
Weight of each base station relative to Grid Sampling point i is calculated, calculation formula is
Wherein, N be base station quantity, dniRepresent sampled point i to the distance of nth base station；
Weight is normalized operation：
Calculate the settling amount of Grid Sampling point i
7. the method for drawing area sedimentation quantity map as claimed in claim 1, it is characterized in that：Based on each Grid Sampling point-rendering
Regional subsidence amount includes carrying out each Grid Sampling point using Gaussian filter smooth.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711096451.XA CN107945247A (en) | 2017-11-08 | 2017-11-08 | The method that drawing area settles quantity map |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201711096451.XA CN107945247A (en) | 2017-11-08 | 2017-11-08 | The method that drawing area settles quantity map |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN107945247A true CN107945247A (en) | 2018-04-20 |
Family
ID=61934587
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201711096451.XA Pending CN107945247A (en) | 2017-11-08 | 2017-11-08 | The method that drawing area settles quantity map |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107945247A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104504621A (en) * | 2014-10-29 | 2015-04-08 | 中国电建集团成都勘测设计研究院有限公司 | Method for evaluating utilizability of mountainous area engineering site |
| CN106017483A (en) * | 2016-05-06 | 2016-10-12 | 厦门蓝斯通信股份有限公司 | Map vehicle icon drawing method and system and navigation terminal |
| CN106225766A (en) * | 2016-07-01 | 2016-12-14 | 陕西铁道工程勘察有限公司 | Monitoring and managing method is analyzed in structure sedimentation based on 3S technology monitoring |
-
2017
- 2017-11-08 CN CN201711096451.XA patent/CN107945247A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104504621A (en) * | 2014-10-29 | 2015-04-08 | 中国电建集团成都勘测设计研究院有限公司 | Method for evaluating utilizability of mountainous area engineering site |
| CN106017483A (en) * | 2016-05-06 | 2016-10-12 | 厦门蓝斯通信股份有限公司 | Map vehicle icon drawing method and system and navigation terminal |
| CN106225766A (en) * | 2016-07-01 | 2016-12-14 | 陕西铁道工程勘察有限公司 | Monitoring and managing method is analyzed in structure sedimentation based on 3S technology monitoring |
Non-Patent Citations (1)
| Title |
|---|
| 徐琳: "基于ArcGIS的地面沉降数据空间分析技术与应用", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105607093B (en) | A kind of integrated navigation system and the method for obtaining navigation coordinate | |
| CN103399326B (en) | GNSS (global navigation satellite system) dynamic measurement accuracy test system and method | |
| CN106970398B (en) | Satellite visibility analysis and ephemeris forecasting method considering satellite shielding condition | |
| CN104748722B (en) | Utilize the elevation location method of satellite positioning information real time calibration barometric leveling result | |
| WO2011127332A1 (en) | Probe data processing | |
| Niu et al. | Development and evaluation of GNSS/INS data processing software for position and orientation systems | |
| CN105301601A (en) | Global navigation satellite system (GNSS) ionosphere delayed three-dimensional modeling method suitable for global area | |
| Zhang et al. | Using Allan variance to evaluate the relative accuracy on different time scales of GNSS/INS systems | |
| Mu et al. | A GNSS/INS-integrated system for an arbitrarily mounted land vehicle navigation device | |
| Ayres-Sampaio et al. | A comparison between three IMUs for strapdown airborne gravimetry | |
| Islam et al. | An effective approach to improving low-cost GPS positioning accuracy in real-time navigation | |
| Yang et al. | Precise levelling in crossing river over 5 km using total station and GNSS | |
| CN102879013B (en) | Method for correcting influences of atmosphere inclination on ground star observation values | |
| Zhuang et al. | Autonomous WLAN heading and position for smartphones | |
| CN103697857B (en) | Height measurement method and device | |
| KR100538319B1 (en) | Method for constituting gis of river information by updating river pathway changes to digital map via mobile internet | |
| CN110058319B (en) | Magnetotelluric data acquisition method and device and terminal equipment | |
| KR101060575B1 (en) | Space Object Pointing Device | |
| CN111679297A (en) | Noise point drift removal method for GPS positioning track | |
| CN107945247A (en) | The method that drawing area settles quantity map | |
| CN116166680A (en) | Automatic updating and maintaining method and system for railway Beidou reference station control network | |
| Piras et al. | Comparison of two different mass-market IMU generations: bias analyses and real time applications | |
| Shi et al. | BeiDou/GPS relative kinematic positioning in challenging environments including poor satellite visibility and high receiver velocity | |
| Alothman et al. | Validation of regional geoid models for Saudi Arabia using GPS/levelling data and GOCE models | |
| Xue et al. | Formulas for precisely and efficiently estimating the bias and variance of the length measurements |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180420 |