CN110501733A - A kind of adaptive grid VRS is generated and method of servicing - Google Patents
A kind of adaptive grid VRS is generated and method of servicing Download PDFInfo
- Publication number
- CN110501733A CN110501733A CN201910679958.0A CN201910679958A CN110501733A CN 110501733 A CN110501733 A CN 110501733A CN 201910679958 A CN201910679958 A CN 201910679958A CN 110501733 A CN110501733 A CN 110501733A
- Authority
- CN
- China
- Prior art keywords
- grid points
- grid
- double difference
- difference correction
- cor
- 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
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/03—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers
- G01S19/07—Cooperating elements; Interaction or communication between different cooperating elements or between cooperating elements and receivers providing data for correcting measured positioning data, e.g. DGPS [differential GPS] or ionosphere corrections
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
It is generated the invention discloses a kind of adaptive grid VRS and method of servicing, reads whole net reference station coordinates information first, according to equidistantly respectively along east to, the north to the grid points for covering entire Reference network are generated, obtaining each grid points coordinate;Then double difference Correction of Errors number information between the whole net reference station of one epoch of reading, calculates the double difference Correction of Errors number of the satellites in view of all grid points;Then judge that the double difference Correction of Errors number difference of different grid points, marked invalid grid points obtain this epoch Irregular Geogrid point;Judge that grid points user's number of request, the grid points in label service provide this epoch grid points VRS service;It repeats the above steps, until having handled all epoch.Real-time adaptive of the present invention generates grid points, has taken into account the calculation amount and user service precision of CORS system;And the calculation amount of CORS system is significantly reduced come dynamic generation grid points VRS according to user's request.
Description
Technical field
The invention belongs to satellite navigation and positioning data processing techniques, and in particular to a kind of adaptive grid VRS is generated and clothes
Business method.
Background technique
Currently, all kinds of GNSS reference stations that China establishes there are thousands of, save, city-level CORS system has more or less a hundred.But it is more
A independent region CORS system, there are the repeated construction of GNSS reference station resource, servers to share, user service range and
The problems such as limited amount.It is me by the CORS system that the CORS system of dispersion becomes whole nation unified management, maintenance and service
The certainty of state GNSS real-time high-precision oriented development.For nationwide CORS system, the user scope sum number of service
Amount is extremely huge, broadcasts mode using the region common VRS of CORS system, can greatly limit online user's quantity simultaneously, can not
Meets the needs of a wide range of, large user's amount.For this purpose, thering is scholar to propose to use a kind of grid VRS technology, i.e., in service range etc.
Interval divides grid points and establishes grid VRS, is provided services to the user by grid points.System need to only maintain each grid points
VRS reduces the calculation amount of system, relieves the limit value of number of users.
However, the factors such as it is unevenly distributed since reference station is unevenly distributed, surveys area's atmospheric condition, the lattice equidistantly established
The VRS observation precision that site obtains be it is different, it is inconsistent to will cause user service precision.On the other hand, using etc. between
Away from grid points, it if grid points spacing is too small, will increase system-computed amount, waste computing resource;If grid points spacing is excessive, can drop
Low grid points VRS observation precision, to reduce user service precision.
Summary of the invention
A kind of adaptive grid VRS generation and method of servicing in order to overcome the deficiencies in the prior art, the present invention provides,
Grid points spacing is adjusted in real time according to factors such as reference station distribution situation, atmospheric characteristics, combines calculation amount and service precision.
The technical scheme adopted by the invention is that: a kind of adaptive grid VRS is generated and method of servicing, which is characterized in that
The following steps are included:
Step 1: read whole net reference station coordinates information, according to equidistantly respectively along east to, the north to generate covering it is whole
The grid points of a Reference network obtain each grid points coordinate;
Step 2: reading double difference Correction of Errors number information between the whole net reference station of an epoch;
Step 3: calculating the double difference Correction of Errors number of the satellites in view of all grid points;
Step 4: judging that the double difference Correction of Errors number difference of different grid points, marked invalid grid points obtain this epoch not
Regular grid point;
Step 5: judging that grid points user's number of request, the grid points in label service provide this epoch grid points VRS step
Service;
Step 6: step 2- step 5 is repeated, until having handled all epoch.
Preferably, the specific implementation of step 3 includes following sub-step:
Step 3.1: reading the coordinate of a grid points, select away from 1 nearest reference station of the grid points as main reference
It stands, is denoted as main website, obtain the double difference Correction of Errors information at the N number of auxiliary reference station nearest with the main website;The auxiliary reference station,
It is denoted as extension station;N is preset value;
Step 3.2: calculating the double difference Correction of Errors number of all satellites in view of the grid points;
This step further comprises following sub-step:
Step 3.2.1: reading the double difference Correction of Errors number of one satellites in view of main website, and finds and miss with the satellite double difference
The extension station of poor correction;
Step 3.2.2: linear interpolation obtains the double difference Correction of Errors number of the satellites in view at grid points;
Step 3.2.3: repeating step 3.2.1- step 3.2.2, until obtaining all satellites in view of main website at grid points
Double difference Correction of Errors number;
Step 3.3: repeating step 3.1- step 3.2, change until obtaining the double difference error of satellites in view of all grid points
Positive number.
Preferably, the specific implementation of step 3.2.2 includes following sub-step:
(1) the M extension station nearest away from grid points is selected, carries out linear interpolation together with main website, obtains at grid points that this can
Depending on the double difference Correction of Errors number Cor of satellite1;M is preset value;
(2) judge in all triangles for the extension station composition that main website and step 3.2.1 obtain, if this can be surrounded by having
The triangle of grid points;
Such triangle if it exists then further finds qualified smallest triangle, i.e. the three of triangle push up
Point is nearest at a distance from grid points, carries out linear interpolation with the main website and 2 extension stations for constituting the triangle, and obtaining should at grid points
The double difference Correction of Errors number Cor of satellites in view2, compare Cor2With Cor1If meeting:
|Cor2-Cor1|≤0.2m;
Then show that the interpolation result is effective, by Cor2Double difference Correction of Errors number as the satellites in view at grid points;If
It is unsatisfactory for the inequality, then shows that the interpolation result is unreliable, it is believed that double difference error of the satellite at grid points can not be obtained
Correction;
Such triangle if it does not exist, then directly by Cor1Double difference Correction of Errors as the satellites in view at grid points
Number.
Preferably, the specific implementation process of step 4 is: reading the double difference Correction of Errors number of a grid points, and therewith
The double difference Correction of Errors number of grid points afterwards is compared, if difference meets a certain specified threshold, by grid points mark later
It is denoted as invalid grid points;It is finished until to the processing judgement of all grid points.
Preferably, in relatively grid points double difference Correction of Errors number difference, can only compare with it apart from phase in step 4
Close grid points double difference Correction of Errors number;
Specific implementation includes following sub-step:
(1) grid points with behavior master, in more every row;
Specific implementation includes following sub-step:
1) a grid points Grid of a line is read(i,j)Double difference Correction of Errors number information Cor(i,j);
2) next grid points Grid of the row is read(i,k)Double difference Correction of Errors number information Cor(i,k), k=j+ at this time
1;
3) judge whether two grid points satellites in view are identical, if it is different, then executing 4);If they are the same, then successively relatively more every
The double difference Correction of Errors number difference of satellites in view, takes the maximum value of its absolute difference to change as the double difference error between grid points
Positive number discrepancy delta Cor, formula indicate are as follows:
Wherein, subscript 1,2 ..., s indicate satellites in view serial number;If 4) Δ Cor > 0.02m, executes;If Δ Cor≤
0.02m then marks grid points Grid(i,k)For invalid grid points, i.e. bDeleted(i,k)=1;Continue to read next grid points
Grid(i,k)Double difference Correction of Errors number information Cor(i,k), at this time 3) k=k+1, repeats;
4) enable j=k, repeat 2) -3) step, until the row is disposed;
5) next line grid points information is read, is repeated the above steps, until all rows are disposed;
(2) based on column, compare the grid points in each column;
Specific implementation includes following sub-step:
1) a grid points Grid of a column is read(i,j)Double difference Correction of Errors number information Cor(i,j);If Grid(i,j)For
Invalid grid points, i.e. bDeleted(i,j)=1, then continue the double difference Correction of Errors number information for reading next grid points, until
Grid(i,j)For effective grid points;If Grid(i,j)It is always invalid grid points, then executes 5);
2) next grid points Grid of the column is read(k,j)Double difference Correction of Errors number information Cor(k,j), k=i+ at this time
1;If grid points Grid(k,j)For invalid grid points, i.e. bDeleted(k,j)=1, then re-execute 1) -2), i=k+1 at this time;It is no
Then, it executes 3);
3) judge whether two grid points satellites in view are identical, if not identical, execute 4);If they are the same, then successively compare
The double difference Correction of Errors number difference of every satellites in view, takes the maximum value of its absolute difference as the double difference error between grid points
Correction discrepancy delta Cor, formula indicate are as follows:
Wherein, subscript 1,2 ..., s indicate satellites in view serial number;If 4) Δ Cor > 0.02m, executes;If Δ Cor≤
0.02m then marks grid points Grid(k,j)For invalid grid points, i.e. bDeleted(k,j)=1;Continue to read next grid points
Grid(k,j)Double difference Correction of Errors number information Cor(k,j), k=k+1 at this time, if Grid(k,j)For effective grid points, then repeatedly 3),
Otherwise k=k+1 is enabled, is executed 4);
4) enable i=k, repeat 1) -3) step, until the column processing finishes;
5) next column grid points information is read, is repeated the above steps, until all column processing finish.
Preferably, the specific implementation of step 5 includes following sub-step:
Step 5.1: an effective grid points information is read, if the grid points have user's request, according to the grid points
Double difference Correction of Errors number generates the VRS Step Information of the grid points and broadcasts to user;If the grid points no user is requested,
Stop the generation of grid points VRS Step Information;
Step 5.2: step 5.1 is repeated, until having handled all effective grid points.
Preferably, when initially obtaining Irregular Geogrid point, user's number of request is defaulted as 0, later, root in step 5.1
The user received according to each epoch resets user's number of request to the request message of grid points;
Specific implementation process is: receiving a certain epoch message, if user's number of request of the grid points is 0, and the epoch should
Grid points are requested by a user, then user's number of request of the grid points are set as 3600, and generate the VRS information distributing of the grid points
To user;If user's number of request of the grid points is 0, and the grid points are not requested by a user, then do not generate the VRS of the grid points
Information;If user's number of request of the grid points be greater than 0, no matter the epoch grid points whether be requested by a user, generate the lattice
The VRS information of site, and user's number of request of the grid points is subtracted 1.
Compared to the prior art, the invention has the advantages that and the utility model has the advantages that
(1) adaptive generation unequal spacing grid points account for the factors such as reference station distribution situation, atmospheric characteristic, guarantee
The VRS observation precision of whole grid site, has ensured the consistency of user service precision, has taken into account the calculation amount of CORS system
With user service precision.
(2) grid points are generated in real time by epoch, it is contemplated that the validity of reference station and the dynamic change of atmospheric conditions.
(3) grid points VRS information is generated according to user's request, and VRS information is not generated to the grid points of no user request,
Significantly reduce the calculation amount of CORS system.
Detailed description of the invention
Fig. 1 is the flow chart of the embodiment of the present invention.
Specific embodiment
In order to illustrate more clearly of the present invention and/or technical solution in the prior art, Detailed description of the invention sheet will be compareed below
The specific embodiment of invention.It should be evident that drawings in the following description are only some embodiments of the invention, for this
For the those of ordinary skill of field, without creative efforts, it can also be obtained according to these attached drawings others
Attached drawing, and obtain other embodiments.
Referring to Fig.1, a kind of adaptive grid VRS provided by the invention is generated and method of servicing, comprising the following steps:
Step 1: read whole net reference station coordinates information, according to equidistantly respectively along east to, the north to generate covering it is whole
The grid points of a Reference network obtain each grid points coordinate;
This step specifically includes:
Step 1.1: reading whole net reference station coordinates information.
Step 1.2: according to spacing d, usually 5km covers entire Reference network to, the north to generating along east respectively
N row m column grid points, and obtain each grid points coordinate.
Step 2: reading double difference Correction of Errors number information between the whole net reference station of an epoch;
Step 3: calculating the double difference Correction of Errors number of the satellites in view of all grid points;
This step specifically includes:
Step 3.1: reading the coordinate of a grid points, select away from 1 nearest reference station of the grid points as main reference
It stands (main website), obtains the double difference Correction of Errors information of the 10 auxiliary reference stations (extension station) nearest with the main website.
Step 3.2: calculating the double difference Correction of Errors number of all satellites in view of the grid points;
Step 3.2 specifically includes:
Step 3.2.1: reading the double difference Correction of Errors number of one satellites in view of main website, and finds and miss with the satellite double difference
The extension station of poor correction.
Step 3.2.2: linear interpolation obtains the double difference Correction of Errors number of the satellites in view at grid points.
In the present embodiment, the interpolation strategies specifically used are as follows:
The 2-3 extension station nearest away from grid points is selected, carries out linear interpolation together with main website, obtains at grid points this visually
The double difference Correction of Errors number Cor of satellite1;
Judge in all triangles for the extension station composition that main website and 3.2.1 step obtain, if the grid can be surrounded by having
The triangle of point, such triangle, then further find qualified smallest triangle, i.e. the three of triangle if it exists
Vertex is nearest at a distance from grid points, carries out linear interpolation with the main website and 2 extension stations for constituting the triangle, obtains at grid points
The double difference Correction of Errors number Cor of the satellites in view2, compare Cor2With Cor1If meeting
|Cor2-Cor1|≤0.2m
Then show that the interpolation result is effective, by Cor2Double difference Correction of Errors number as the satellites in view at grid points;If
It is unsatisfactory for the inequality, then shows that the interpolation result is unreliable, it is believed that double difference error of the satellite at grid points can not be obtained
Correction;
Such triangle if it does not exist, then directly by Cor1Double difference Correction of Errors as the satellites in view at grid points
Number.
3.2.3 3.2.1-3.2.2 is repeated, until obtaining double difference Correction of Errors of all satellites in view of main website at grid points
Number.
3.3 repeat 3.1-3.2, the double difference Correction of Errors number of the satellites in view until obtaining all grid points.
Step 4: judging that the double difference Correction of Errors number difference of different grid points, marked invalid grid points obtain this epoch not
Regular grid point;
Specific implementation process is: reading the double difference Correction of Errors number of a grid points, and the double difference with grid points later
Correction of Errors number is compared, if difference meets a certain specified threshold, grid points later are labeled as invalid grid points;Directly
It is finished to the processing judgement of all grid points.
More specifically, when each epoch starts to calculate, grid points are set as effectively, i.e. bDeleted=0;Due to grid
Point double difference Correction of Errors number is obtained by its neighbouring reference station interpolation, the grid points double difference Correction of Errors number of entire Reference network
To be obtained by different linear interpolation planes, in relatively grid points double difference Correction of Errors number difference, can only compare and its
The grid points double difference Correction of Errors number of close proximity, specific steps include:
(1) grid points with behavior master, in more every row;
Further comprise:
1) a grid points Grid of a line is read(i,j)Double difference Correction of Errors number information Cor(i,j);
2) next grid points Grid of the row is read(i,k)Double difference Correction of Errors number information Cor(i,k), k=j+ at this time
1;
3) judge whether two grid points satellites in view are identical, if it is different, then executing 4);If they are the same, then successively relatively more every
The double difference Correction of Errors number difference of satellites in view, takes the maximum value of its absolute difference to change as the double difference error between grid points
Positive number discrepancy delta Cor, formula indicate are as follows:
Wherein, subscript 1,2 ..., s indicate satellites in view serial number.If 4) Δ Cor > 0.02m, executes;If Δ Cor≤
0.02m then marks grid points Grid(i,k)For invalid grid points, i.e. bDeleted(i,k)=1;Continue to read next grid points
Grid(i,k)Double difference Correction of Errors number information Cor(i,k), at this time 3) k=k+1, repeats;
4) enable j=k, repeat 2) -3) step, until the row is disposed;
5) next line grid points information is read, is repeated the above steps, until all rows are disposed.
(2) based on column, compare the grid points in each column;
Further comprise:
1) a grid points Grid of a column is read(i,j)Double difference Correction of Errors number information Cor(i,j);If Grid(i,j)For
Invalid grid points, i.e. bDeleted(i,j)=1, then continue the double difference Correction of Errors number information for reading next grid points, until
Grid(i,j)For effective grid points;If Grid(i,j)It is always invalid grid points, then executes 5);
2) next grid points Grid of the column is read(k,j)Double difference Correction of Errors number information Cor(k,j), k=i+ at this time
1;If grid points Grid(k,j)For invalid grid points, i.e. bDeleted(k,j)=1, then re-execute 1) -2), i=k+1 at this time;It is no
Then, it executes 3);
3) judge whether two grid points satellites in view are identical, if not identical, execute 4);If they are the same, then successively compare
The double difference Correction of Errors number difference of every satellites in view, takes the maximum value of its absolute difference as the double difference error between grid points
Correction discrepancy delta Cor, formula indicate are as follows:
Wherein, subscript 1,2 ..., s indicate satellites in view serial number.If 4) Δ Cor > 0.02m, executes;If Δ Cor≤
0.02m then marks grid points Grid(k,j)For invalid grid points, i.e. bDeleted(k,j)=1;Continue to read next grid points
Grid(k,j)Double difference Correction of Errors number information Cor(k,j), k=k+1 at this time, if Grid(k,j)For effective grid points, then repeatedly 3),
Otherwise k=k+1 is enabled, is executed 4);
4) enable i=k, repeat 1) -3) step, until the column processing finishes;
5) next column grid points information is read, is repeated the above steps, until all column processing finish.
Step 5: judging that grid points user's number of request, the grid points in label service provide this epoch grid points VRS clothes
Business;
This step specifically includes:
Step 5.1: an effective grid points information is read, if the grid points have user's request, according to the grid points
Double difference Correction of Errors number generates the VRS information of the grid points and broadcasts to user;If the grid points no user is requested, stop
Grid points VRS information generates.
In this implementation, when initially obtaining Irregular Geogrid point, user's number of request is defaulted as 0, later, according to each epoch
The user received resets user's number of request to the request message of grid points, and specific steps include:
A certain epoch message is received, if user's number of request of the grid points is 0, and the epoch grid points are asked by user
It asks, then user's number of request nDataReq of the grid points is set as 3600, and generate the VRS information distributing of the grid points to use
Family;If user's number of request nDataReq of the grid points is 0, and the grid points are not requested by a user, then do not generate the grid points
VRS information;If user's number of request nDataReq of the grid points be greater than 0, no matter the epoch grid points whether asked by user
It asks, generates the VRS information of the grid points, and enable user's number of request nDataReq=nDataReq-1 of the grid points.
Step 5.2: step 5.1 is repeated, until having handled all effective grid points.
Step 6: step 2- step 5 is repeated, until having handled all epoch.
Be described in above-described embodiment to illustrate the present invention, though text in be illustrated by specific term, not
Can be limited the scope of protection of the present invention with this, be familiar with this technical field personage can understand spirit of the invention with it is right after principle
It changes or modifies and reaches equivalent purpose, and this equivalent change and modification, should all be covered by scope of the claims institute circle
Determine in scope.
Claims (7)
1. a kind of adaptive grid VRS is generated and method of servicing, which comprises the following steps:
Step 1: whole net reference station coordinates information is read, according to equidistantly covering entire ginseng to, the north to generating along east respectively
The grid points for examining station net, obtain each grid points coordinate;
Step 2: reading double difference Correction of Errors number information between the whole net reference station of an epoch;
Step 3: calculating the double difference Correction of Errors number of the satellites in view of all grid points;
Step 4: judging that it is irregular to obtain this epoch for the double difference Correction of Errors number difference of different grid points, marked invalid grid points
Grid points;
Step 5: judging that grid points user's number of request, the grid points in label service provide this epoch grid points VRS step clothes
Business;
Step 6: step 2- step 5 is repeated, until having handled all epoch.
2. adaptive grid VRS according to claim 1 is generated and method of servicing, which is characterized in that the specific reality of step 3
Now include following sub-step:
Step 3.1: reading the coordinate of a grid points, select away from 1 nearest reference station of the grid points as main reference station, note
For main website, the double difference Correction of Errors information at the N number of auxiliary reference station nearest with the main website is obtained;The auxiliary reference station, is denoted as
Extension station;Wherein, N is preset value;
Step 3.2: calculating the double difference Correction of Errors number of all satellites in view of the grid points;
This step further comprises following sub-step:
Step 3.2.1: the double difference Correction of Errors number of one satellites in view of main website is read, and finds and changes with the satellite double difference error
The extension station of positive number;
Step 3.2.2: linear interpolation obtains the double difference Correction of Errors number of the satellites in view at grid points;
Step 3.2.3: repeating step 3.2.1- step 3.2.2, double at grid points until obtaining all satellites in view of main website
Mistake difference correction;
Step 3.3: repeating step 3.1- step 3.2, the double difference Correction of Errors number of the satellites in view until obtaining all grid points.
3. adaptive grid VRS according to claim 2 is generated and method of servicing, which is characterized in that the tool of step 3.2.2
Body realizes to include following sub-step:
(1) the M extension station nearest away from grid points is selected, carries out linear interpolation together with main website, this is obtained at grid points and visually defends
The double difference Correction of Errors number Cor of star1;Wherein, M is preset value;
(2) judge in all triangles for the extension station composition that main website and step 3.2.1 obtain, if the grid can be surrounded by having
The triangle of point;
Such triangle if it exists, then further find qualified smallest triangle, i.e., the three of triangle vertex with
The distance of grid points is nearest, carries out linear interpolation with the main website and 2 extension stations for constituting the triangle, obtains at grid points this visually
The double difference Correction of Errors number Cor of satellite2, compare Cor2With Cor1If meeting:
|Cor2-Cor1|≤0.2m;
Then show that the interpolation result is effective, by Cor2Double difference Correction of Errors number as the satellites in view at grid points;If being unsatisfactory for
The inequality then shows that the interpolation result is unreliable, it is believed that can not obtain double difference Correction of Errors number of the satellite at grid points;
Such triangle if it does not exist, then directly by Cor1Double difference Correction of Errors number as the satellites in view at grid points.
4. adaptive grid VRS according to claim 1 is generated and method of servicing, which is characterized in that the specific reality of step 4
Existing process is: reading the double difference Correction of Errors number of a grid points, and carries out with the double difference Correction of Errors number of grid points later
Compare, if difference meets a certain specified threshold, grid points later are labeled as invalid grid points;Until to all grid points
Processing judgement finishes.
5. adaptive grid VRS according to claim 1 or 4 is generated and method of servicing, it is characterised in that: in step 4, In
When comparing grid points double difference Correction of Errors number difference, the grid points double difference Correction of Errors number with its close proximity can only be compared;
Specific implementation includes following sub-step:
(1) grid points with behavior master, in more every row;
Specific implementation includes following sub-step:
1) a grid points Grid of a line is read(i,j)Double difference Correction of Errors number information Cor(i,j);
2) next grid points Grid of the row is read(i,k)Double difference Correction of Errors number information Cor(i,k), k=j+1 at this time;
3) judge whether two grid points satellites in view are identical, if it is different, then executing 4);If they are the same, then successively more every can
Depending on the double difference Correction of Errors number difference of satellite, take the maximum value of its absolute difference as the double difference Correction of Errors number between grid points
Discrepancy delta Cor, formula indicate are as follows:
Wherein, subscript 1,2 ..., s indicate satellites in view serial number;If 4) Δ Cor > 0.02m, executes;If Δ Cor≤0.02m,
Then mark grid points Grid(i,k)For invalid grid points, i.e. bDeleted(i,k)=1;Continue to read next grid points Grid(i,k)
Double difference Correction of Errors number information Cor(i,k), at this time 3) k=k+1, repeats;
4) enable j=k, repeat 2) -3) step, until the row is disposed;
5) next line grid points information is read, is repeated the above steps, until all rows are disposed;
(2) based on column, compare the grid points in each column;
Specific implementation includes following sub-step:
1) a grid points Grid of a column is read(i,j)Double difference Correction of Errors number information Cor(i,j);If Grid(i,j)It is invalid
Grid points, i.e. bDeleted(i,j)=1, then continue the double difference Correction of Errors number information for reading next grid points, until
Grid(i,j)For effective grid points;If Grid(i,j)It is always invalid grid points, then executes 5);
2) next grid points Grid of the column is read(k,j)Double difference Correction of Errors number information Cor(k,j), k=i+1 at this time;If lattice
Site Grid(k,j)For invalid grid points, i.e. bDeleted(k,j)=1, then re-execute 1) -2), i=k+1 at this time;Otherwise, it holds
Row 3);
3) judge whether two grid points satellites in view are identical, if not identical, execute 4);If they are the same, then successively more every
The double difference Correction of Errors number difference of satellites in view, takes the maximum value of its absolute difference as the double difference Correction of Errors between grid points
Number discrepancy delta Cor, formula indicate are as follows:
Wherein, subscript 1,2 ..., s indicate satellites in view serial number;If 4) Δ Cor > 0.02m, executes;If Δ Cor≤0.02m,
Then mark grid points Grid(k,j)For invalid grid points, i.e. bDeleted(k,j)=1;Continue to read next grid points Grid(k,j)
Double difference Correction of Errors number information Cor(k,j), k=k+1 at this time, if Grid(k,j)Then repeatedly 3) k is otherwise enabled for effective grid points
4)=k+1 is executed;
4) enable i=k, repeat 1) -3) step, until the column processing finishes;
5) next column grid points information is read, is repeated the above steps, until all column processing finish.
6. adaptive grid VRS according to claim 1 is generated and method of servicing, which is characterized in that the specific reality of step 5
Now include following sub-step:
Step 5.1: an effective grid points information is read, if the grid points have user's request, according to the double difference of the grid points
Correction of Errors number generates the VRS Step Information of the grid points and broadcasts to user;If the grid points no user is requested, stop
Grid points VRS Step Information generates;
Step 5.2: step 5.1 is repeated, until having handled all effective grid points.
7. adaptive grid VRS according to claim 1 or 6 is generated and method of servicing, it is characterised in that: in step 5.1,
When initially obtaining Irregular Geogrid point, user's number of request is defaulted as 0, later, the user's plaid matching received according to each epoch
The request message of site resets user's number of request;
Specific implementation process is: a certain epoch message is received, if user's number of request of the grid points is 0, and the epoch grid
Point is requested by a user, then user's number of request of the grid points is set as 3600, and generates the VRS information distributing of the grid points to use
Family;If user's number of request of the grid points is 0, and the grid points are not requested by a user, then do not generate the VRS letter of the grid points
Breath;If user's number of request of the grid points be greater than 0, no matter the epoch grid points whether be requested by a user, generate the grid
The VRS information of point, and user's number of request of the grid points is subtracted 1.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910679958.0A CN110501733B (en) | 2019-07-26 | 2019-07-26 | Self-adaptive grid VRS generation and service method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910679958.0A CN110501733B (en) | 2019-07-26 | 2019-07-26 | Self-adaptive grid VRS generation and service method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110501733A true CN110501733A (en) | 2019-11-26 |
CN110501733B CN110501733B (en) | 2022-11-11 |
Family
ID=68586925
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910679958.0A Active CN110501733B (en) | 2019-07-26 | 2019-07-26 | Self-adaptive grid VRS generation and service method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110501733B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111954226A (en) * | 2020-07-30 | 2020-11-17 | 宁波冶金勘察设计研究股份有限公司 | Triangular network of reference station based on Delaunay heterogeneous CORS system and integration method |
WO2021121062A1 (en) * | 2019-12-20 | 2021-06-24 | 西安空间无线电技术研究所 | "cloud-end" architecture-based ppp-rtk enhancement method and system |
CN113596722A (en) * | 2021-08-04 | 2021-11-02 | 武汉攀达时空科技有限公司 | Method for generating and serving adjacent user data product sharing MCBI |
CN116132987A (en) * | 2023-04-19 | 2023-05-16 | 武汉大学 | RTK data service method, device, electronic equipment and storage medium |
CN117008163A (en) * | 2023-06-07 | 2023-11-07 | 广东省国土资源测绘院 | Method for generating and refining grid VRS by using digital elevation model and related equipment |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2040090A1 (en) * | 2007-09-18 | 2009-03-25 | Leica Geosystems AG | Method for accuracy estimation of network based corrections for a satellite-aided positioning system |
CN106855632A (en) * | 2016-12-30 | 2017-06-16 | 广州市中海达测绘仪器有限公司 | A kind of broadcast type VRS localization methods and system |
EP3343247A1 (en) * | 2016-12-29 | 2018-07-04 | Thales | Method for optimally adjusting give error bounds or for optimally calculating residual variances of igp points of an ionospheric grid for correcting an sbas system and sbas system for carrying out said method |
-
2019
- 2019-07-26 CN CN201910679958.0A patent/CN110501733B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2040090A1 (en) * | 2007-09-18 | 2009-03-25 | Leica Geosystems AG | Method for accuracy estimation of network based corrections for a satellite-aided positioning system |
EP3343247A1 (en) * | 2016-12-29 | 2018-07-04 | Thales | Method for optimally adjusting give error bounds or for optimally calculating residual variances of igp points of an ionospheric grid for correcting an sbas system and sbas system for carrying out said method |
CN106855632A (en) * | 2016-12-30 | 2017-06-16 | 广州市中海达测绘仪器有限公司 | A kind of broadcast type VRS localization methods and system |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021121062A1 (en) * | 2019-12-20 | 2021-06-24 | 西安空间无线电技术研究所 | "cloud-end" architecture-based ppp-rtk enhancement method and system |
CN111954226A (en) * | 2020-07-30 | 2020-11-17 | 宁波冶金勘察设计研究股份有限公司 | Triangular network of reference station based on Delaunay heterogeneous CORS system and integration method |
CN113596722A (en) * | 2021-08-04 | 2021-11-02 | 武汉攀达时空科技有限公司 | Method for generating and serving adjacent user data product sharing MCBI |
CN116132987A (en) * | 2023-04-19 | 2023-05-16 | 武汉大学 | RTK data service method, device, electronic equipment and storage medium |
CN117008163A (en) * | 2023-06-07 | 2023-11-07 | 广东省国土资源测绘院 | Method for generating and refining grid VRS by using digital elevation model and related equipment |
CN117008163B (en) * | 2023-06-07 | 2024-03-26 | 广东省国土资源测绘院 | Method for generating and refining grid VRS by using digital elevation model and related equipment |
Also Published As
Publication number | Publication date |
---|---|
CN110501733B (en) | 2022-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110501733A (en) | A kind of adaptive grid VRS is generated and method of servicing | |
CN105629263B (en) | A kind of troposphere atmosphere delay estimation error correcting method and correction system | |
CN106569239B (en) | A kind of broadcast type network RTK location technology | |
CN106528597B (en) | The mask method and device of point of interest | |
CN107426695A (en) | Gridding differential data broadcasting method, server and storage medium based on broadcast | |
CN106445676A (en) | Distributed data calculation-based task allocation method and task allocation apparatus | |
CN107976702A (en) | A kind of position correcting method based on CORS, positioning terminal and alignment system | |
CN101882175B (en) | Computer graphics method for manufacturing centrifugal rating plate | |
CN110418361B (en) | Multi-engine resolving and high-precision position service method of large-scale CORS network | |
CN105738934B (en) | The quick fixing means of URTK fuzzinesses of additional atmospheric information dynamic constrained | |
CN110954931A (en) | Positioning method, positioning device and computer readable storage medium | |
CN106875047A (en) | Reservoir watershed Runoff Forecast method and system | |
CN111352129A (en) | Method and device for monitoring differential quality and computer medium | |
CN113325448A (en) | Large-altitude-difference CORS network resolving method considering troposphere delay reconstruction | |
CN111208535B (en) | Calculation method based on international reference ionosphere total electron content abnormal value correction | |
CN105608714A (en) | Framing image generation method based on CAD and GIS cooperation and system thereof | |
CN105096062A (en) | CORS (continuously operating reference station) application management system | |
CN109241212B (en) | Retrieval method for historical rainfall based on mesoscale numerical atmospheric mode and high resolution | |
CN114125699A (en) | Network RTK service method for reconstruction by using virtual reference station | |
CN110011719A (en) | It is a kind of to obtain control method using the star upper rail of GPS orbit determination data | |
CN109408041A (en) | A kind of the map coordinates system conversion method and electronic equipment of language based on programming | |
CN111307125B (en) | Inclined-axis cylindrical projection method based on GNSS and ground ranging combined adjustment | |
CN109344439A (en) | A kind of modeling of building block formula and its simulative construction method based on BIM | |
CN105610487A (en) | Satellite observation resource cross-domain coordinated scheduling method | |
CN110456399A (en) | A kind of continuous operational reference station and its operation 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 |