CN106202000B - Seven-parameter transformation method between national three-dimensional system of coordinate and anywhere plane coordinate system - Google Patents
Seven-parameter transformation method between national three-dimensional system of coordinate and anywhere plane coordinate system Download PDFInfo
- Publication number
- CN106202000B CN106202000B CN201610580387.1A CN201610580387A CN106202000B CN 106202000 B CN106202000 B CN 106202000B CN 201610580387 A CN201610580387 A CN 201610580387A CN 106202000 B CN106202000 B CN 106202000B
- Authority
- CN
- China
- Prior art keywords
- coordinate
- coordinate system
- parameter
- dimensional
- local
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Image Processing (AREA)
- Processing Or Creating Images (AREA)
Abstract
A kind of the invention discloses national three-dimensional system of coordinate and anywhere seven-parameter transformation method between plane coordinate system, the present invention can not directly acquire the single seven-parameter transformation parameter of high-precision for the problem that existing algorithm, and then can not obtain in real time local coordinate system achievement based on CORS measurement data；Propose a kind of seven-parameter transformation model of additional additive constant.This method passes through volume of data processing method: additional additive constant significantly reduces the influence of distortion of projection, influence etc. of the geodetic height to ginseng precision is sought is reduced by three Parameter Switch models, significantly improve the precision of conversion parameter, and method, which meets measurement handbook currently on the market, can only receive the requirement that single four/seven parameter inputs, the national coordinate system achievement obtained based on CORS network measure, it can be converted into local coordinate system achievement in real time, this is conducive to expand application space of the provincial, and municipal level CORS in measurement practice process.
Description
Technical field
The invention belongs to technical field of geographic information, it is related to a kind of national three-dimensional system of coordinate and anywhere plane coordinate system
Between seven-parameter transformation method, more particularly to a kind of new high-precision country three-dimensional system of coordinate and anywhere between plane coordinate system
Seven-parameter transformation acquiring method.
Background technique
National Geodetic Datum and local reference data are confidential datas, and each province and city area is in order to protect basic parameter
It is close, when establishing local coordinate system, on the basis of national coordinate system, mostly carry out the projective transformation of central meridian, portion
Region-by-region has further carried out four-parameter conversion (irregular translation, rotation, scaling etc.) to coordinate after transformation, this makes ground
There is more complicated transformational relations between square coordinate system and national coordinate system.And measurement hand currently used in the market
Book can only receive single four/seven parameters input, this makes user ginseng method can not be asked to seek reliably converting using traditional
Parameter, so that coordinates of local systems can not be obtained in real time.Under the background that Provincial CORS is covered comprehensively, this is very not
It is promoted and applied conducive to provincial CORS in area.
Summary of the invention
In order to solve the above-mentioned technical problem, the invention proposes a kind of new high-precision country three-dimensional system of coordinate and arbitrarily
Seven-parameter transformation acquiring method between square plane coordinate system.
The technical scheme adopted by the invention is that: it is a kind of country three-dimensional system of coordinate and anywhere between plane coordinate system seven ginseng
Number conversion method, which comprises the following steps:
Step 1: unified country's three-dimensional system of coordinate and local plane coordinate system central meridian are corrected using one group of additive constant
Local coordinate system achievement；
Step 2: weakening influence of the geodetic height direction to two-dimensional assemblage precision using three Parameter Switch models；
Step 3: seeking correcting the conversion between local plane coordinate system and national three-dimensional system of coordinate using seven parameter models
Relationship；
Step 4: seven parameters of setting and additive constant obtain local coordinate system achievement based on CORS network in real time.
Preferably, the specific implementation of step 1 includes following sub-step:
Step 1.1: by the coordinate series (x of national three-dimensional system of coordinate2000,y2000) pass through projective transformation, by original country
Three-dimensional system of coordinate central meridian L02000It is transformed into local plane coordinate system central meridian L0dlUnder, the state after being corrected
Family three-dimensional system of coordinate coordinate (xgz2000,ygz2000)；
Step 1.2: by the national three-dimensional system of coordinate coordinate (x after correctiongz2000,ygz2000) with corresponding local plane coordinates
It is coordinate (xdl,ydl) ask poor, obtain the difference in the direction x and the direction y；
Step 1.3: seeking average value △ x, △ y, △ x, the △ y of both direction difference respectively as the direction x, y plus normal
Number；
Step 1.4: it is flat that average value △ x, △ y being added separately to the place system after being corrected on Local coordinates plane system
Areal coordinate (xgzdl,ydgzl)。
Preferably, the specific implementation of step 2 includes following sub-step:
Step 2.1: being plane coordinates (x by the place after correctiongzdl,ydgzl) it is converted into longitude and latitude (Bgzdl, Lgzdl)；
Step 2.2: by national three-dimensional system of coordinate (x2000,y2000) it is converted into longitude and latitude (B2000, L2000)；
Step 2.3: obtaining national three-dimensional system of coordinate geodetic height H2000；
Step 2.4: by H2000Instead of local coordinate system geodetic height Hdl, then seek turning between Two coordinate system using three parameters
Relationship is changed, and local coordinate system geodetic height H is updated according to conversion parameterdl。
Preferably, the specific implementation process of step 3 is: by (Bgzdl, Lgzdl) with update after local plane coordinate system
Geodetic height HdlAs target-based coordinate system (Bgzdl, Lgzdl, Hdl), by (B2000, L2000, H2000) it is used as source coordinate system；Seek two coordinates
Transformational relation between system obtains last seven parameter of conversion.
Preferably, the specific implementation process of step 4 is: creating task on handbook, projection is arranged in task attribute
Parameter, i.e., seven parameter of conversion obtained in step 3, and handbook be arranged corresponding additive constant (on handbook default setting plus it is normal
Number is north 0.000, at this moment corresponding additive constant can be added on this default value by east 500000.000)；Eventually by single
The method that seven parameters add additive constant obtains local plane coordinate system in real time.
Compared with the existing method for directly seeking seven parameters, the advantageous effect of present invention is that: it is directed to existing calculation
Method can not directly acquire the single seven-parameter transformation parameter of high-precision, and then can not obtain in real time place based on CORS measurement data
The problem of coordinate system achievement, proposes a kind of seven-parameter transformation model of additional additive constant.This method is by adding additive constant pole
The earth reduces the influence of distortion of projection, reduces influence etc. of the geodetic height to ginseng precision is sought by three Parameter Switch models, greatly
The earth improves the precision of conversion parameter, and method meets measurement handbook currently on the market can only receive single four/seven parameter
The requirement of input can be converted into local coordinate system achievement based on the national coordinate system achievement that CORS network measure obtains in real time,
This is conducive to expand application space of the provincial, and municipal level CORS in measurement practice process.
Detailed description of the invention
Fig. 1: the flow chart of the embodiment of the present invention.
Specific embodiment
Understand for the ease of those of ordinary skill in the art and implement the present invention, with reference to the accompanying drawings and embodiments to this hair
It is bright to be described in further detail, it should be understood that implementation example described herein is merely to illustrate and explain the present invention, not
For limiting the present invention.
The present invention is insufficient for available data processing method, proposes the new high-precision country three-dimensional system of coordinate of one kind and appoints
Seven-parameter transformation acquiring method between the local plane coordinate system of meaning.Pass through unified national coordinate system and local coordinate system center first
Meridian greatly weakens the influence that distortion of projection seeks conversion parameter, also with three Parameter Switch models, reduces
Local coordinate system does not have the influence that geodetic height this defect seeks precision to conversion parameter, finally utilizes seven unified parameter moulds
Type seeks the transformational relation between national coordinate system and local coordinate system.In actual use, the present invention passes through seven ginsengs first
The national coordinate system result transformation that number real-time measurement obtains is then to deduct the influence of additive constant to place, final to obtain in real time
The local coordinate system achievement that user needs.
Referring to Fig.1, it (is data instance, local coordinate system achievement with somewhere place that the present embodiment, which combines specific measured data,
Data have been subjected to processing) further explain is made to the present invention.
Local plane coordinate system and national three-dimensional system of coordinate achievement are obtained first, and obtain each coordinate system central meridian letter
Breath, as shown in table 1,
The national three-dimensional system of coordinate of table 1 and local plane coordinate system achievement
Method point following steps of the invention carry out, and specific steps are described as follows:
Step 1: national three-dimensional system of coordinate and local plane coordinate system central meridian are unified, and reducing distortion of projection influences.
First by national three-dimensional system of coordinate (with (x2000,y2000) for) coordinate series is by projective transformation, by original state
Family three-dimensional system of coordinate central meridian (L02000) it is transformed into local plane coordinate system central meridian (L0dl) under, it obtains correcting it
Country three afterwards maintains coordinate (xgz2000,ygz2000).By this coordinate and corresponding local plane coordinate system coordinate (xdl,ydl) ask
Difference obtains the difference in the direction x and the direction y, and seek respectively both direction difference average value △ x, △ y (△ x, △ y as x,
The additive constant in the direction y).Then this average value is added to local plane coordinates and fastens local plane system coordinate after being corrected
(xgzdl,ydgzl)。
The present embodiment projects to BLH under national three-dimensional system of coordinate under local plane coordinate system central meridian, obtains new
National three-dimensional system of coordinate is corrected under coordinate system as a result, the results are shown in Table 2.By the coordinate acquired and corresponding local plane coordinates
System asks poor, and is averaging respectively in the direction x, y, obtains the additive constant in the direction x, y.
State's family coordinate and additive constant (unit: m) after the projection of table 2
Step 2: further weakening influence of the geodetic height direction to two-dimensional assemblage precision using three parameter model.
By the local plane system coordinate (x after correctiongzdl,ydgzl) it is converted into longitude and latitude (Bgzdl, Lgzdl)；Simultaneously by state
Family three maintains coordinate (x2000,y2000) it is converted into longitude and latitude (B2000, L2000)；It is big to obtain original national three-dimensional system of coordinate simultaneously
The high H in ground2000, first by H2000Instead of local plane coordinate system geodetic height Hdl, then sought between Two coordinate system using three parameters
Transformational relation, and local plane coordinate system geodetic height H is updated according to conversion parameterdl。
In the additive constant correction to coordinates of local systems that the present embodiment will acquire, and it is converted into BLH.Acquired results such as table
Shown in 3,
Local plane system coordinate and national three-dimensional coordinate BLH after the correction of table 3
Local plane coordinate system geodetic height is directly replaced with country's three-dimensional system of coordinate geodetic height in table 3.After correction
Local plane coordinates and primitive state three-dimensional system of coordinate coordinate seek three parameters, further weaken geodetic height inaccuracy to parameter essence
The influence of degree, the results are shown in Table 4.In the present embodiment, local plane coordinate system and national three-dimensional system of coordinate use unified ellipsoid,
So the earth High variation is little, but this step work is also in the case that local plane coordinate system is different with country's three-dimensional system of coordinate
It is meaningful.
The corrected local plane coordinate system coordinate of 4 three parameter of table
Step 3: seeking correcting the conversion between local plane coordinate system and national three-dimensional system of coordinate using seven parameter models
Relationship；
By (Bgzdl, Lgzdl) with update after local plane coordinate system geodetic height HdlAs target-based coordinate system (Bgzdl,
Lgzdl, Hdl), by (B2000, L2000, H2000) it is used as source coordinate system；The transformational relation between Two coordinate system is finally sought, is obtained most
Seven parameter of conversion afterwards；
Using result in table 4 and original national three-dimensional system of coordinate as a result, seeking conversion parameter.In actual use process
In, the national three-dimensional system of coordinate result of real-time measurement is transformed into local plane coordinate system using the conversion parameter sought is existing
Under, the local plane coordinate system of demand is finally then obtained in real time plus additive constant shown in table 2 as a result, for details, reference can be made to table 5,
BLH xy is converted in table to show.
Table 5 utilizes seven-parameter transformation result and conversion front and rear coordinate residual error (unit: m)
Step 4: seven parameters and additive constant is arranged in handbook, obtains local coordinate system achievement in real time based on CORS network.
Task is created on handbook, projective parameter is set in task attribute, i.e. seven parameter of conversion that acquires of previous step, and
In the additive constant that handbook setting responds, (in order to maintain secrecy to actual conversion parameter, seven parameters being arranged herein are by adding
Conversion parameter after close, encryption technology can be found in related existing patent), the additive constant of default setting is north 0.000 on handbook,
Corresponding additive constant △ x, △ y, at this moment can be added on this default value by east 500000.000.Eventually by single seven parameter
The method of additional additive constant realizes the purpose for obtaining local coordinate system in real time.
For the present invention by a series of measurement adjustment data processing means, the method for adding additive constant by single seven parameter is real
Show national three-dimensional geodetic datum coordinate system and anywhere the high-precision between coordinate system is converted, measures handbook meeting
Under the premise of only receiving single four/seven parameters input requirements, the demand of Local Independent Coordinate System is obtained in real time.
The present invention is the coordinate conversion acquiring method for having broad prospect of application in National primary standard course of modernization, is had
Help be widely applied regional coordinate system on the basis of existing provincial CORS network all-round construction is completed and obtain skill in real time
Art.
It should be understood that the part that this specification does not elaborate belongs to the prior art.
It should be understood that the above-mentioned description for preferred embodiment is more detailed, can not therefore be considered to this
The limitation of invention patent protection range, those skilled in the art under the inspiration of the present invention, are not departing from power of the present invention
Benefit requires to make replacement or deformation under protected ambit, fall within the scope of protection of the present invention, this hair
It is bright range is claimed to be determined by the appended claims.
Claims (6)
1. a kind of country's three-dimensional system of coordinate and anywhere seven-parameter transformation method between plane coordinate system, which is characterized in that including
Following steps:
Step 1: unified country's three-dimensional system of coordinate and local plane coordinate system central meridian correct place using one group of additive constant
Coordinate system achievement；
Step 2: weakening influence of the geodetic height direction to two-dimensional assemblage precision using three Parameter Switch models；
Step 3: seeking correcting the transformational relation between local plane coordinate system and national three-dimensional system of coordinate using seven parameter models；
Step 4: seven parameters of setting and additive constant obtain local coordinate system achievement based on CORS network in real time.
2. country's three-dimensional system of coordinate according to claim 1 and anywhere seven-parameter transformation method between plane coordinate system,
It is characterized in that, the specific implementation of step 1 includes following sub-step:
Step 1.1: by the coordinate series (x of national three-dimensional system of coordinate2000,y2000) pass through projective transformation, it is three-dimensional by original country
Coordinate system central meridian L02000It is transformed into local plane coordinate system central meridian L0dlUnder, the country three after being corrected
Tie up coordinate system coordinate (xgz2000,ygz2000)；
Step 1.2: by the national three-dimensional system of coordinate coordinate (x after correctiongz2000,ygz2000) sat with corresponding local plane coordinate system
Mark (xdl,ydl) ask poor, obtain the difference in the direction x and the direction y；
Step 1.3: seeking the average value Δ x of both direction difference, Δ y, Δ x, additive constant of the Δ y as the direction x, y respectively；
Step 1.4: it is that plane is sat that average value Δ x, Δ y are added separately to the place after being corrected on Local coordinates plane system
Mark (xgzdl,ydgzl)。
3. country's three-dimensional system of coordinate according to claim 2 and anywhere seven-parameter transformation method between plane coordinate system,
It is characterized in that, the specific implementation of step 2 includes following sub-step:
Step 2.1: being plane coordinates (x by the place after correctiongzdl,ydgzl) it is converted into longitude and latitude (Bgzdl, Lgzdl)；
Step 2.2: by national three-dimensional system of coordinate (x2000,y2000) it is converted into longitude and latitude (B2000, L2000)；
Step 2.3: obtaining national three-dimensional system of coordinate geodetic height H2000；
Step 2.4: by H2000Instead of local coordinate system geodetic height Hdl, then seek converting between Two coordinate system using three parameters and close
System, and local coordinate system geodetic height H is updated according to conversion parameterdl。
4. country's three-dimensional system of coordinate according to claim 3 and anywhere seven-parameter transformation method between plane coordinate system,
It is characterized in that, the specific implementation process of step 3 is: by (Bgzdl, Lgzdl) with update after local plane coordinate system geodetic height
HdlAs target-based coordinate system (Bgzdl, Lgzdl, Hdl), by (B2000, L2000, H2000) it is used as source coordinate system；It seeks between Two coordinate system
Transformational relation, obtain last seven parameter of conversion.
5. country's three-dimensional system of coordinate according to claim 4 and anywhere seven-parameter transformation method between plane coordinate system,
It is characterized in that, the specific implementation process of step 4 is: creating task on handbook, projective parameter is arranged in task attribute, i.e.,
Seven parameter of conversion obtained in step 3, and corresponding additive constant is set in handbook；Additive constant is added eventually by single seven parameter
Method obtain local plane coordinate system in real time.
6. country's three-dimensional system of coordinate according to claim 5 and anywhere seven-parameter transformation method between plane coordinate system,
It is characterized by: described be arranged corresponding additive constant in handbook, it is that Δ x, Δ y are added to the additive constant of handbook default setting
On.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610580387.1A CN106202000B (en) | 2016-07-22 | 2016-07-22 | Seven-parameter transformation method between national three-dimensional system of coordinate and anywhere plane coordinate system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610580387.1A CN106202000B (en) | 2016-07-22 | 2016-07-22 | Seven-parameter transformation method between national three-dimensional system of coordinate and anywhere plane coordinate system |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106202000A CN106202000A (en) | 2016-12-07 |
| CN106202000B true CN106202000B (en) | 2019-01-18 |
Family
ID=57491249
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610580387.1A Active CN106202000B (en) | 2016-07-22 | 2016-07-22 | Seven-parameter transformation method between national three-dimensional system of coordinate and anywhere plane coordinate system |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106202000B (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106845837B (en) * | 2017-01-23 | 2020-05-22 | 国网山东省电力公司电力科学研究院 | Power transmission and transformation project environment sensitive area monitoring system and method based on big data technology |
| CN107845060A (en) * | 2017-10-31 | 2018-03-27 | 广东中星电子有限公司 | Geographical position and corresponding image position coordinates conversion method and system |
| CN108731648B (en) * | 2018-03-15 | 2020-12-22 | 广州市城市规划勘测设计研究院 | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium |
| CN109947877B (en) * | 2018-07-25 | 2022-04-29 | 武汉光庭信息技术股份有限公司 | Method and system for improving map positioning precision of GIS mobile terminal |
| CN114001650B (en) * | 2021-09-16 | 2023-09-29 | 北京市测绘设计研究院 | Encryption method for conversion parameters of local coordinate system and arbitrary plane coordinate system |
| CN117708960B (en) * | 2024-02-04 | 2024-05-03 | 武汉大学 | Real-time conversion method, device, equipment and medium for plane coordinates and normal height |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103389080A (en) * | 2013-08-19 | 2013-11-13 | 重庆市地理信息中心 | Method for obtaining independent city coordinate system parameter based on geographic information application |
| CN104007457A (en) * | 2014-05-28 | 2014-08-27 | 武汉市测绘研究院 | Parameter-free real-time locating method based on CORS |
| CN105160197A (en) * | 2015-09-23 | 2015-12-16 | 湖北省基础地理信息中心 | Comprehensive geographical spatial data coordinate conversion method and system |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7649540B2 (en) * | 2006-03-24 | 2010-01-19 | Virgil Stanger | Coordinate transformations system and method thereof |
-
2016
- 2016-07-22 CN CN201610580387.1A patent/CN106202000B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103389080A (en) * | 2013-08-19 | 2013-11-13 | 重庆市地理信息中心 | Method for obtaining independent city coordinate system parameter based on geographic information application |
| CN104007457A (en) * | 2014-05-28 | 2014-08-27 | 武汉市测绘研究院 | Parameter-free real-time locating method based on CORS |
| CN105160197A (en) * | 2015-09-23 | 2015-12-16 | 湖北省基础地理信息中心 | Comprehensive geographical spatial data coordinate conversion method and system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106202000A (en) | 2016-12-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106202000B (en) | Seven-parameter transformation method between national three-dimensional system of coordinate and anywhere plane coordinate system | |
| CN102095425B (en) | Method for generating electronic channel chart based on Changjiang River standards | |
| TWI328678B (en) | ||
| CN104463954B (en) | A kind of analogy method and system of 3-D view surface details | |
| CN103810286A (en) | Coordinate point positioning method for matching two-dimensional map with three-dimensional map | |
| CN104007457A (en) | Parameter-free real-time locating method based on CORS | |
| CN109801371B (en) | Network three-dimensional electronic map construction method based on Cesium | |
| CN106447774A (en) | Mapping method from GPS (global positioning system) map to three-dimensional map | |
| AU2003253785A1 (en) | Method of determining the long axis of an object | |
| CN106897445A (en) | It is a kind of that the novel method that operation delimited on border is improved based on ArcGIS | |
| CN103198494A (en) | Building elevation extraction method based on building outer contour lines and laser-point cloud | |
| CN104537031B (en) | The amending method and device of a kind of map datum | |
| CN107066489A (en) | A kind of method of Fast Drawing farmland electronic map | |
| CN106815823A (en) | Bearing calibration and its device are demarcated in a kind of lens distortion | |
| CN107339938A (en) | A kind of special-shaped calibrating block and scaling method for single eye stereo vision self-calibration | |
| CN109558643B (en) | Traffic scene modeling and model unitization method | |
| CN102750734B (en) | The method and system that a kind of virtual three-dimensional earth system is shown | |
| CN108897937B (en) | Method for automatically converting CAD (computer aided design) data of civil aviation airport into DEM (digital elevation model) data | |
| Simon et al. | CAD and GIS techniques in georeferencing maps for the identification and mapping of meadows in Arad county. | |
| CN110298778A (en) | One kind being based on raster data and the modified traffic zone demographics method of area ratio | |
| CN103646027A (en) | Geographic spatial database updating method and system | |
| CN105865413B (en) | A kind of acquisition methods and device of depth of building | |
| CN105488309B (en) | Cutter track track restores display methods, system and engraving machine | |
| CN106875330A (en) | A kind of method that areal model turns Sphere Measurement Model | |
| CN106705848A (en) | Ball-link steel structure grid frame reverse modeling method |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |