CN108731648A - 2000 coordinate system parameter acquiring methods, device and computer readable storage medium - Google Patents
2000 coordinate system parameter acquiring methods, device and computer readable storage medium Download PDFInfo
- Publication number
- CN108731648A CN108731648A CN201810215665.2A CN201810215665A CN108731648A CN 108731648 A CN108731648 A CN 108731648A CN 201810215665 A CN201810215665 A CN 201810215665A CN 108731648 A CN108731648 A CN 108731648A
- Authority
- CN
- China
- Prior art keywords
- parameter
- control point
- coordinate system
- projection
- plane
- 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
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C15/00—Surveying instruments or accessories not provided for in groups G01C1/00 - G01C13/00
Abstract
The present invention provides a kind of 2000 coordinate system parameter acquiring methods, device and computer readable storage medium, this method includes the central meridian and perspective plane elevation according to target cities coordinate system, and the first projective parameter is arranged；According to the first projective parameter, the CGCS2000 geodetic coordinates to overlapping control point in multiple control nets of the target cities obtained in advance carries out gauss projection and obtains the first gauss projection coordinate；The existing plane coordinates at control point and preset Parameter Switch model are overlapped according to the first gauss projection coordinate, the target cities obtained in advance, calculates the first plane conversion parameter；According to the rotation angle parameter and scale parameter of the first plane conversion parameter, best translation parameters and best projection parameter are determined.The projective parameter and translation parameters of 2000 coordinate system of target cities can be accurately obtained by this method, it realizes that 2000 coordinate system of target cities and existing coordinate system are steadily excessive, improves coordinate consistency between 2000 coordinate system of target cities and the existing coordinate system in city.
Description
Technical field
The present invention relates to coordinate system parameters acquiring technology fields, and in particular to a kind of 2000 coordinate system parameter acquiring sides
Method, device and computer readable storage medium.
Background technology
Currently, State Bureau of Surveying and Mapping establish the requirement of 2000 coordinate system of city under, domestic a part of city is established
2000 coordinate systems, such as：2000 Urban Independent Coordinate System of Wuhan, 2000 Urban Independent Coordinate System of Ningbo and Nanchang
2000 Urban Independent Coordinate Systems；It is most typical but there are many problems during establishing 2000 Urban Independent Coordinate System
Be can not ensure the consistency of new and old coordinate, such as：2000 Urban Independent Coordinate System of Wuhan is due to City East, western border district
And to be more than the regional distortion of projection of 145m transfinite geodetic height, and new old information connection problem is not accounted for, new and old coordinate system is sat
Indicate certain difference；2000 Urban Independent Coordinate System of Ningbo is unsatisfactory for being less than 1,/40,000 due to City East, western both sides distortion of projection
The requirement of (2.5cm/km) only needs to translate between the old and new's coordinate, but translation parameters be 100km big number A add it is one small
Number B, but different regions B is inconsistent, is not convenient to use；2000 Urban Independent Coordinate System of Nanchang is thrown due to east, western border district
Shadow deformation is unsatisfactory for requiring, and coordinate differs greatly between new and old coordinate system, planning region 1:500 topographic map data need to carry out coordinate and turn
Changing can just be continuing with, and need to put into a large amount of manpower, material resources and financial resources and carry out graph conversion work.Therefore, CGCS2000 is based on reference to ellipse
During ball creates 2000 coordinate system of city, how to realize and both met the specification that shadow boundary is no more than 1,/40,000
It is required that and meet and historical summary seamless connection so that achievement can continue use under existing coordinate system, avoid Two coordinate system not
Consistency leads to huge conversion work amount and expense, becomes the newly-built independent seat in the city 2000 based on CGCS2000 reference ellipsoids
Mark system urgent problem to be solved.
Invention content
The object of the present invention is to provide a kind of 2000 coordinate system parameter acquiring methods, device and computer-readable storages
Medium can accurately obtain the projective parameter and translation parameters of target cities 2000 coordinate system of the city based on 2000 spheroids,
It realizes that 2000 coordinate system of target cities and existing coordinate system are steadily excessive, improves 2000 coordinate system of target cities and city
Coordinate consistency between the existing coordinate system in city, so that new old information seamless connection, while achievement can continue under existing coordinate system
It uses, reduces conversion work amount and expense.
An embodiment of the present invention provides a kind of 2000 coordinate system parameter acquiring methods, including：
Obtain the CGCS2000 geodetic coordinates that control point is overlapped in multiple control nets of target cities and coincidence control
The plane coordinates of the corresponding target cities coordinate system of system point；
According to the central meridian and perspective plane elevation of the target cities coordinate system, the first projective parameter is set；
According to first projective parameter, gauss projection is carried out to the CGCS2000 geodetic coordinates for overlapping control point,
Calculate the first gauss projection coordinate for overlapping control point；
According to the first gauss projection coordinate, the plane coordinates and the preset Parameter Switch mould for overlapping control point
Type calculates the first plane conversion parameter between the first gauss projection coordinate and the plane coordinates for overlapping control point；
According to the rotation angle parameter and scale parameter of the first plane conversion parameter, determine best translation parameters and
Best projection parameter.
Preferably, the rotation angle parameter and scale parameter according to the first plane conversion parameter determines best
Translation parameters and best projection parameter, specifically include；
According to the first plane conversion parameter, the second projective parameter is set；
According to second projective parameter, gauss projection is carried out to the CGCS2000 geodetic coordinates for overlapping control point,
Calculate the second gauss projection coordinate for overlapping control point；
Turned according to the second gauss projection coordinate, the plane coordinates and the preset parameter for overlapping control point
Mold changing type calculates the between the second gauss projection coordinate and the plane coordinates for overlapping control point for overlapping control point
Two plane conversion parameters；
According to the rotation angle parameter and scale parameter of the second plane conversion parameter, the best translation parameters is determined
And the best projection parameter.
Preferably, the central meridian and perspective plane elevation according to the target cities coordinate system, setting first
Projective parameter specifically includes：
First projective parameter includes the first central meridian and the first perspective plane elevation；
According to the central meridian of the target cities coordinate system, the numerical intervals of first central meridian are determined,
And first central meridian is set according to the first interval of setting；
According to the projection elevation face of the target cities coordinate system, the numerical intervals in first projection elevation face are determined,
And according to the second interval setting first projection elevation face of setting.
Preferably, described according to the first plane conversion parameter, the second projective parameter is set, is specifically included：
Second projective parameter includes the second central meridian and the second perspective plane elevation；
According to the first plane conversion parameter, the numerical intervals and described second of second central meridian are determined
The numerical intervals of perspective plane elevation；
Second central meridian is set according to the third interval of setting and according to the 4th interval setting institute of setting
State the second projection elevation face.
Preferably, the third interval according to setting is arranged second central meridian and according to the 4th of setting
Interval setting second projection elevation face, specifically includes；
When second central meridian is L, according to the 4th interval setting the first projection elevation of the setting
Face obtains the parameter combination that N groups second central meridian projects elevation face with described second；Wherein, L is to belong to described the
In the numerical intervals of two central meridians and according to the spaced central meridian of the third of setting；
When second projection elevation face is H, second central meridian is set according to the third interval of setting, is obtained
Obtain the parameter combination that M groups second central meridian projects elevation face with described second；H is to belong to the second projection elevation
In the numerical intervals in face and according to the 4th spaced projection elevation face of setting.
Preferably, the rotation angle parameter and scale parameter according to the second plane conversion parameter determines best
Translation parameters and best projection parameter, specifically include：
The rotation angle parameter for meeting the second plane conversion parameter is obtained from M+N the second plane conversion parameters
I-th of second plane conversion parameters of α ≈ 0 and the scale parameter m ≈ 1 of the second plane conversion parameter；Wherein, according to N groups
The parameter combination and M groups second central meridian and institute that second central meridian projects elevation face with described second
The parameter combination in the second projection elevation face is stated, M+N the second plane conversion parameters are can get；
Obtain the translation parameters of described i-th second plane conversion parameters；Wherein, described i-th second plane conversion ginsengs
Several translation parameters is the best translation parameters；
The best projection parameter includes best central meridian and best projection elevation face；
According to the rotation angle parameter α and the second central meridian L of the M+N the second plane conversion parameters,
Fit the first fitting function of the rotation angle parameter and second central meridian；
According to first fitting function, the corresponding best central meridian as rotation angle parameter α=0 is calculated；
Elevation face H is projected according to the scale parameter m of the M+N the second plane conversion parameters and described second, is intended
Close out the second fitting function that the scale parameter projects elevation face with described second；
According to second fitting function, the corresponding best projection elevation face as scale parameter m=1 is calculated.
Preferably, the 2000 coordinate system parameter acquiring method further includes：
According to the first gauss projection coordinate, the plane coordinates and the preset Parameter Switch mould for overlapping control point
Type calculates the internal accuracy for overlapping control point；
Obtain the external check point of the multiple control net；Wherein, the external check point is in the multiple control net
The misaligned control point in control point is overlapped with described；
According to the coincidence control point and the external check point, the external accuracy for overlapping control point is calculated；
According to the internal accuracy and external accuracy for overlapping control point, the abnormal control in the coincidence control point is determined
Point processed；
The exception control point in the coincidence control point is rejected, and recalculates the coincidence after rejecting the exception control point
The corresponding first plane conversion parameter in control point.
The embodiment of the present invention additionally provides a kind of 2000 coordinate system parameter obtaining devices, including：
Coordinate obtaining module, the CGCS2000 the earth that control point is overlapped in multiple control nets for obtaining target cities are sat
It is marked with and the plane coordinates for overlapping the corresponding target cities coordinate system in control point；
First projective parameter setup module is used for the central meridian according to the target cities coordinate system and perspective plane
The first projective parameter is arranged in elevation；
First gauss projection coordinate calculation module is used for according to first projective parameter, to the coincidence control point
CGCS2000 geodetic coordinates carries out gauss projection, calculates the first gauss projection coordinate for overlapping control point；
First plane conversion parameter calculating module, for according to the first gauss projection coordinate for overlapping control point, institute
Plane coordinates and preset Parameter Switch model are stated, the first gauss projection coordinate for overlapping control point is calculated and is put down with described
The first plane conversion parameter between areal coordinate；
Coordinate system parameters acquisition module, for being joined according to the rotation angle parameter and scale of the first plane conversion parameter
Number, determines best translation parameters and best projection parameter.
The embodiment of the present invention additionally provides a kind of 2000 coordinate system parameter obtaining devices, including processor, memory
And it is stored in the memory and is configured as the computer program executed by the processor, the processor executes institute
2000 above-mentioned coordinate system parameter acquiring methods are realized when stating computer program.
The embodiment of the present invention additionally provides a kind of computer readable storage medium, which is characterized in that described computer-readable
Storage medium includes the computer program of storage, wherein described computer-readable deposit is controlled when the computer program is run
Equipment where storage media executes 2000 above-mentioned coordinate system parameter acquiring methods.
Compared with the existing technology, a kind of beneficial effect for 2000 coordinate system parameter acquiring method that inventive embodiments provide
Fruit is：The 2000 coordinate system parameter acquiring method includes：It obtains and overlaps control in multiple control nets of target cities
The CGCS2000 geodetic coordinates of point and the plane coordinates for overlapping the corresponding target cities coordinate system in control point；According to described
The first projective parameter is arranged in the central meridian and perspective plane elevation of target cities coordinate system；According to the first projection ginseng
Number carries out gauss projection to the CGCS2000 geodetic coordinates for overlapping control point, and calculate the coincidence control point first is high
This projection coordinate；Turned according to the first gauss projection coordinate, the plane coordinates and the preset parameter for overlapping control point
Mold changing type calculates the first plane conversion ginseng between the first gauss projection coordinate and the plane coordinates for overlapping control point
Number；According to the rotation angle parameter and scale parameter of the first plane conversion parameter, best translation parameters and best is determined
Projective parameter.2000 coordinate system of target cities can be accurately obtained by the 2000 coordinate system parameter acquiring method
Projective parameter, to improve coordinate consistency between 2000 coordinate system of target cities and old coordinate system, so that new and old
Data seamless connection, while achievement can continue use under existing coordinate system, reduce conversion work amount and expense.The present invention is implemented
Example also provides 2000 coordinate system parameter obtaining devices of one kind and computer readable storage medium.
Description of the drawings
Fig. 1 is a kind of flow chart of 2000 coordinate system parameter acquiring method provided in an embodiment of the present invention；
Fig. 2 is a kind of schematic diagram of 2000 coordinate system parameter obtaining device provided in an embodiment of the present invention.
Specific implementation mode
Following will be combined with the drawings in the embodiments of the present invention, and technical solution in the embodiment of the present invention carries out clear, complete
Site preparation describes, it is clear that described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.It is based on
Embodiment in the present invention, those of ordinary skill in the art are obtained every other without creative efforts
Embodiment shall fall within the protection scope of the present invention.
Referring to Fig. 1, it is a kind of flow of 2000 coordinate system parameter acquiring method provided in an embodiment of the present invention
Figure；The 2000 coordinate system parameter acquiring method, including：
S100：It obtains and overlaps the CGCS2000 geodetic coordinates at control point and described heavy in multiple control nets of target cities
Close the plane coordinates of the corresponding target cities coordinate system in control point；
S200：According to the central meridian and perspective plane elevation of the target cities coordinate system, setting the first projection ginseng
Number；
S300：According to first projective parameter, Gauss is carried out to the CGCS2000 geodetic coordinates for overlapping control point
Projection calculates the first gauss projection coordinate for overlapping control point；
S400：According to the first gauss projection coordinate, the plane coordinates and the preset parameter for overlapping control point
Transformation model calculates the first plane conversion between the first gauss projection coordinate and the plane coordinates for overlapping control point
Parameter；
S500：According to the rotation angle parameter and scale parameter of the first plane conversion parameter, best translation ginseng is determined
Number and best projection parameter.
In the present embodiment, the multiple control net includes Guangzhou CORS base stations major network, Guangzhou quasigeoid
Frame net, Precision Quasi-geoid Determining in Guangzhou City high precise control net, Guangzhou GPS chopped-off heads (second-class) plane in 2000 of refining are surveyed
Amount control net；The target cities coordinate system is Guangzhou Coordinate System system, and the CGCS2000 geodetic coordinates is that CGCS2000 the earth is sat
Mark.
Before S100 steps, the CGCS2000 geodetic coordinates for precalculating the control point of multiple control nets is needed, including：
On the basis of the control point of province's CORS systems of target cities, compensating computation is carried out to the control points the CORS major network of target cities,
Obtain the corresponding CGCS2000 geodetic coordinates in control point of the control points CORS major network；Using the control points CORS major network as base
Standard carries out compensating computation processing to the first control net of the target cities, the second control net and third control net respectively, point
Do not show that the corresponding CGCS2000 in control point of the first control net, the second control net and third control net is big
Ground coordinate.Wherein, the first control net, the second control net and the third control net respectively above-mentioned Guangzhou seemingly
Level Geoid frame net, Precision Quasi-geoid Determining in Guangzhou City high precise control net, Guangzhou GPS chopped-off heads in 2000
(second-class) plane survey controls net；The control points CORS major network is Guangzhou CORS base station major networks.
Specifically, the control point for obtaining the existing CGCS2000 geodetic coordinates of Guangzhou CORS base station major networks is corresponding
CGCS2000 geodetic coordinates, and according to the CGCS2000 geodetic coordinates at the control point, using the Guangdong Province control points CORS as start of calculation
Point is carried out compatible row analysis to Fixed Initial Point, is rejected incompatible control point using the POWERNET adjustment Softwares of Wuhan University
Afterwards, the CGCS2000 geodetic coordinates at the remaining control point of compensating computation Guangzhou CORS base station major networks.Then, with Guangzhou CORS bases
The control point of quasi- station owner net is Fixed Initial Point, using identical method, successively to Precision Quasi-geoid Determining in Guangzhou City frame net, wide
State city Precision Quasi-geoid Determining high precise control net, Guangzhou GPS chopped-off heads (second-class) plane survey in 2000 control net carry out
Adjustment calculates Precision Quasi-geoid Determining in Guangzhou City frame net, Precision Quasi-geoid Determining in Guangzhou City high precise control net, Guangzhou
The CGCS2000 geodetic coordinates of city's GPS chopped-off heads (second-class) plane survey in 2000 control net.
In S300 steps, according to first projective parameter, to the coincidence control point under 2000 reference ellipsoids
CGCS2000 geodetic coordinates carries out gauss projection, calculates the first gauss projection coordinate for overlapping control point；
Wherein, gauss projection model is as follows：
Wherein, a2The major semiaxis after (diminution) is expanded for ellipsoid；e2Square of the first eccentricity of ellipsoid；N1It is expanded for ellipsoid
Preceding radius of curvature in prime vertical；B2Latitude after being expanded for ellipsoid；B1Latitude before being expanded for ellipsoid, B1≠B2；I=L-L0, L is
Central meridian；H is perspective plane elevation；ρ is a constant；η=e ' cos B；N=a/W；X indicates to entreat meridianal arc length；B
Indicate the latitude at certain control point.
The preset Parameter Switch model is:
Wherein, x1, y1For the plane coordinates of the corresponding target cities coordinate system in the coincidence control point；x2, y2It is described heavy
Close the gauss projection coordinate at control point；Δ x, Δ y are translation parameters；α is rotation parameter；M is scale parameter.
The projection of target cities 2000 coordinate system of the city based on 2000 spheroids can be accurately obtained by the above method
Parameter and translation parameters realize that 2000 coordinate system of target cities and existing coordinate system are steadily excessive, improve target cities
Coordinate consistency between the existing coordinate system of 2000 coordinate systems and city, so that new old information seamless connection, has simultaneously
Achievement can continue use under coordinate system, reduce conversion work amount and expense.
In a kind of optional embodiment, S500：According to the rotation angle parameter and ruler of the first plane conversion parameter
Parameter is spent, best translation parameters and best projection parameter is determined, specifically includes；
According to the first plane conversion parameter, the second projective parameter is set；
According to second projective parameter, gauss projection is carried out to the CGCS2000 geodetic coordinates for overlapping control point,
Calculate the second gauss projection coordinate for overlapping control point；
Turned according to the second gauss projection coordinate, the plane coordinates and the preset parameter for overlapping control point
Mold changing type calculates the between the second gauss projection coordinate and the plane coordinates for overlapping control point for overlapping control point
Two plane conversion parameters；
According to the rotation angle parameter and scale parameter of the second plane conversion parameter, the best translation parameters is determined
And the best projection parameter.
In a kind of optional embodiment, the central meridian and perspective plane according to the target cities coordinate system
Elevation is arranged the first projective parameter, specifically includes：
First projective parameter includes the first central meridian and the first perspective plane elevation；
According to the central meridian of the target cities coordinate system, the numerical intervals of first central meridian are determined,
And first central meridian is set according to the first interval of setting；
According to the projection elevation face of the target cities coordinate system, the numerical intervals in first projection elevation face are determined,
And according to the second interval setting first projection elevation face of setting.
It is described according to the first plane conversion parameter in a kind of optional embodiment, the second projective parameter, tool are set
Body includes：
Second projective parameter includes the second central meridian and the second perspective plane elevation；
According to the first plane conversion parameter, the numerical intervals and described second of second central meridian are determined
The numerical intervals of perspective plane elevation；
Second central meridian is set according to the third interval of setting and according to the 4th interval setting institute of setting
State the second projection elevation face.
In a kind of optional embodiment, the third interval according to setting be arranged second central meridian and
According to the 4th interval setting second projection elevation face of setting, specifically include；
When second central meridian is L, according to the 4th interval setting the first projection elevation of the setting
Face obtains the parameter combination that N groups second central meridian projects elevation face with described second；Wherein, L is to belong to described the
In the numerical intervals of two central meridians and according to the spaced central meridian of the third of setting；
When second projection elevation face is H, second central meridian is set according to the third interval of setting, is obtained
Obtain the parameter combination that M groups second central meridian projects elevation face with described second；H is to belong to the second projection elevation
In the numerical intervals in face and according to the 4th spaced projection elevation face of setting.
In a kind of optional embodiment, the rotation angle parameter and scale according to the second plane conversion parameter
Parameter determines best translation parameters and best projection parameter, specifically includes：
The rotation angle parameter for meeting the second plane conversion parameter is obtained from M+N the second plane conversion parameters
I-th of second plane conversions ginseng when the scale parameter m ≈ second thresholds of α ≈ first thresholds and the second plane conversion parameter
Number；Wherein, the parameter combination and M groups that elevation face are projected with described second according to N groups second central meridian described the
The parameter combination that two central meridians project elevation face with described second can get M+N the second plane conversion parameters；
Obtain the translation parameters of described i-th second plane conversion parameters；Wherein, described i-th second plane conversion ginsengs
Several translation parameters is the best translation parameters；
The best projection parameter includes best central meridian and best projection elevation face；
According to the rotation angle parameter α and the second central meridian L of the M+N the second plane conversion parameters,
Fit the first fitting function of the rotation angle parameter and second central meridian；
According to first fitting function, calculate when rotation angle parameter α is equal to the first threshold it is corresponding it is best in
Entreat meridian；
Elevation face H is projected according to the scale parameter m of the M+N the second plane conversion parameters and described second, is intended
Close out the second fitting function that the scale parameter projects elevation face with described second；
According to second fitting function, the corresponding best projection when scale parameter m is equal to the second threshold is calculated
Elevation face.
Wherein, the first threshold is 0, and the second threshold is 1.It is calculated and described first by above-mentioned Gauss model
The one-to-one first gauss projection coordinate of parameter combination that central meridian projects elevation face with described first, then by upper
The Parameter Switch model stated, calculate with the one-to-one first plane conversion parameter of the first gauss projection coordinate, according to described
Rotation angle parameter α and scale parameter m in first plane conversion parameter, determination are likely to occur corresponding projection when m ≈ 1 and α ≈ 0
The numberical range of elevation face and central meridian, i.e., the numberical range of described second central meridian and second projection
The numberical range in elevation face.Specifically can by calculate m ≈ 1 and α ≈ 0 in the numberical range of first central meridian and
The probability of occurrence of the numberical range in first projection elevation face, obtains according to the probability of occurrence of m ≈ 1 and α ≈ 0 in described second
Entreat meridianal numberical range and the numberical range in second projection elevation face.
In the same manner, it is calculated by above-mentioned Gauss model and projects elevation face with second central meridian and described second
The one-to-one second gauss projection coordinate of parameter combination calculated with second then by above-mentioned Parameter Switch model
The one-to-one second plane conversion parameter of gauss projection coordinate, according to rotation angle parameter α in the second plane conversion parameter
And scale parameter m, search for corresponding i-th of second plane conversion parameters when m ≈ 1 and α ≈ 0, wherein i-th of second planes
Translation parameters in conversion parameter is the best translation parameters.
First fitting function is：α=a1*L-b1；Second fitting function is：M=a2*H+b2；Wherein, a1、b1
For two constants of first fitting function；a2、b2For two constants of second fitting function.
In a kind of optional embodiment, the 2000 coordinate system parameter acquiring method further includes：
According to the first gauss projection coordinate, the plane coordinates and the preset Parameter Switch mould for overlapping control point
Type calculates the internal accuracy for overlapping control point；
Obtain the external check point of the multiple control net；Wherein, the external check point is in the multiple control net
The misaligned control point in control point is overlapped with described；
According to the coincidence control point and the external check point, the external accuracy for overlapping control point is calculated；
According to the internal accuracy and external accuracy for overlapping control point, the abnormal control in the coincidence control point is determined
Point processed；
The exception control point in the coincidence control point is rejected, and recalculates the coincidence after rejecting the exception control point
The corresponding first plane conversion parameter in control point.
Referring to Fig. 2, it, which is the embodiment of the present invention, additionally provides a kind of showing for 2000 coordinate system parameter obtaining device
It is intended to, the 2000 coordinate system parameter obtaining device includes：
Coordinate obtaining module 1 overlaps the CGCS2000 the earth at control point in multiple control nets for obtaining target cities
Coordinate and the plane coordinates for overlapping the corresponding target cities coordinate system in control point；
First projective parameter setup module 2 is used for the central meridian according to the target cities coordinate system and projection
The first projective parameter is arranged in face elevation；
First gauss projection coordinate calculation module 3 is used for according to first projective parameter, to the coincidence control point
CGCS2000 geodetic coordinates carries out gauss projection, calculates the first gauss projection coordinate for overlapping control point；
First plane conversion parameter calculating module 4, for according to the first gauss projection coordinate for overlapping control point, institute
Plane coordinates and preset Parameter Switch model are stated, the first gauss projection coordinate for overlapping control point is calculated and is put down with described
The first plane conversion parameter between areal coordinate；
Coordinate system parameters acquisition module 5, for the rotation angle parameter and scale according to the first plane conversion parameter
Parameter determines best translation parameters and best projection parameter.
In the present embodiment, the multiple control net includes Guangzhou CORS base stations major network, Guangzhou quasigeoid
Frame net, Precision Quasi-geoid Determining in Guangzhou City high precise control net, Guangzhou GPS chopped-off heads (second-class) plane in 2000 of refining are surveyed
Amount control net；The target cities coordinate system is Guangzhou Coordinate System system, and the CGCS2000 geodetic coordinates is that CGCS2000 the earth is sat
Mark.
The 2000 coordinate system parameter obtaining device further includes CGCS2000 Geodetic Coordinate Calculation modules, for mesh
On the basis of the control point of province's CORS systems for marking city, compensating computation is carried out to the control points the CORS major network of target cities, is obtained
The corresponding CGCS2000 geodetic coordinates in control point of the control points CORS major network；On the basis of the major network of the control points CORS,
Compensating computation processing is carried out to the first control net of the target cities, the second control net and third control net respectively, respectively
Obtain the corresponding CGCS2000 the earth in control point of the first control net, the second control net and third control net
Coordinate.Wherein, the first control net, the second control net and third control net are respectively above-mentioned Guangzhou like big
Ground-level is refined frame net, Precision Quasi-geoid Determining in Guangzhou City high precise control net, Guangzhou GPS chopped-off heads (two in 2000
Deng) plane survey control net；The control points CORS major network is Guangzhou CORS base station major networks.
Specifically, the control point for obtaining the existing CGCS2000 geodetic coordinates of Guangzhou CORS base station major networks is corresponding
CGCS2000 geodetic coordinates, and according to the CGCS2000 geodetic coordinates at the control point, using the Guangdong Province control points CORS as start of calculation
Point is carried out compatible row analysis to Fixed Initial Point, is rejected incompatible control point using the POWERNET adjustment Softwares of Wuhan University
Afterwards, the CGCS2000 geodetic coordinates at the remaining control point of compensating computation Guangzhou CORS base station major networks.Then, with Guangzhou CORS bases
The control point of quasi- station owner net is Fixed Initial Point, using identical method, successively to Precision Quasi-geoid Determining in Guangzhou City frame net, wide
State city Precision Quasi-geoid Determining high precise control net, Guangzhou GPS chopped-off heads (second-class) plane survey in 2000 control net carry out
Adjustment calculates Precision Quasi-geoid Determining in Guangzhou City frame net, Precision Quasi-geoid Determining in Guangzhou City high precise control net, Guangzhou
The CGCS2000 geodetic coordinates of city's GPS chopped-off heads (second-class) plane survey in 2000 control net.
The first gauss projection coordinate calculation module 3 is used for according to first projective parameter, in 2000 reference ellipsoids
Under to it is described overlap control point CGCS2000 geodetic coordinates carry out gauss projection, calculate it is described overlap control point the first Gauss
Projection coordinate；
Wherein, gauss projection model is as follows：
Wherein, a2The major semiaxis after (diminution) is expanded for ellipsoid；e2Square of the first eccentricity of ellipsoid；N1It is expanded for ellipsoid
Preceding radius of curvature in prime vertical；B2Latitude after being expanded for ellipsoid；B1Latitude before being expanded for ellipsoid, B1≠B2；I=L-L0, L is
Central meridian；H is perspective plane elevation；ρ is a constant；η=e ' cosB；N=a/W；X indicates to entreat meridianal arc length；B tables
Show the latitude at certain control point.
The preset Parameter Switch model is:
Wherein, x1, y1For the plane coordinates of the corresponding target cities coordinate system in the coincidence control point；x2, y2It is described heavy
Close the gauss projection coordinate at control point；Δ x, Δ y are translation parameters；α is rotation parameter；M is scale parameter.
The projection of target cities 2000 coordinate system of the city based on 2000 spheroids can be accurately obtained by above-mentioned apparatus
Parameter and translation parameters realize that 2000 coordinate system of target cities and existing coordinate system are steadily excessive, improve target cities
Coordinate consistency between the existing coordinate system of 2000 coordinate systems and city, so that new old information seamless connection, has simultaneously
Achievement can continue use under coordinate system, reduce conversion work amount and expense.
In a kind of optional embodiment, the coordinate system parameters acquisition module 5 includes；
Second projective parameter setup module, for according to the first plane conversion parameter, the second projective parameter to be arranged；
Second gauss projection coordinate calculation module is used for according to second projective parameter, to the coincidence control point
CGCS2000 geodetic coordinates carries out gauss projection, calculates the second gauss projection coordinate for overlapping control point；
Second plane conversion parameter calculating module, for according to the second gauss projection coordinate for overlapping control point, institute
Plane coordinates and the preset Parameter Switch model are stated, the second gauss projection coordinate for overlapping control point and institute are calculated
State the second plane conversion parameter between the plane coordinates for overlapping control point；
Coordinate system parameters determining module, for being joined according to the rotation angle parameter and scale of the second plane conversion parameter
Number, determines the best translation parameters and the best projection parameter.
In a kind of optional embodiment, first projective parameter includes the first central meridian and the first perspective plane
Elevation；
The first projective parameter setup module is additionally operable to the central meridian according to the target cities coordinate system, really
The numerical intervals of fixed first central meridian, and first central meridian is set according to the first interval of setting；
The first projective parameter setup module is additionally operable to the projection elevation face according to the target cities coordinate system, really
The numerical intervals in fixed first projection elevation face, and according to the second interval setting first projection elevation face of setting.
In a kind of optional embodiment, second projective parameter includes the second central meridian and the second perspective plane
Elevation；
The second projective parameter setup module is additionally operable to, according to the first plane conversion parameter, determine described second
The numerical intervals of the numerical intervals of central meridian and second perspective plane elevation；
The second projective parameter setup module is additionally operable to the described second central meridian of third interval setting according to setting
Line and the 4th interval setting second projection elevation face according to setting.
In a kind of optional embodiment, the second projective parameter setup module is additionally operable to when second center
When noon line is L, according to the 4th interval setting first projection elevation face of the setting, N groups the second center is obtained
The parameter combination that noon line projects elevation face with described second；Wherein, L is to belong in the numerical intervals of second central meridian
And according to the spaced central meridian of the third of setting；
The second projective parameter setup module is additionally operable to when second projection elevation face is H, according to the of setting
Second central meridian is arranged in three intervals, obtains M groups second central meridian and projects elevation face with described second
Parameter combination；H is in the numerical intervals for belong to second projection elevation face and according to the 4th spaced projection of setting
Elevation face.
In a kind of optional embodiment, the coordinate system parameters determining module includes：
Second plane conversion parameter acquiring unit meets institute for being obtained from M+N the second plane conversion parameters
State the rotation angle parameter α ≈ first thresholds of the second plane conversion parameter and the scale parameter m ≈ of the second plane conversion parameter
I-th of second plane conversion parameters when second threshold；Wherein, it is thrown according to N groups second central meridian and described second
The parameter combination that parameter combination and M groups second central meridian in shadow elevation face project elevation face with described second, can
Obtain M+N the second plane conversion parameters；
Translation parameters acquiring unit, the translation parameters for obtaining described i-th second plane conversion parameters；Wherein, institute
The translation parameters for stating i-th of second plane conversion parameters is the best translation parameters；
The best projection parameter includes best central meridian and best projection elevation face；
First function fitting unit, for according to the rotation angle parameter α of the M+N the second plane conversion parameters with
The second central meridian L fits the first fitting function of the rotation angle parameter and second central meridian；
Central meridian acquiring unit, for according to first fitting function, calculating described in being equal to as rotation angle parameter α
Corresponding best central meridian when first threshold；
Second function fitting unit is used for the scale parameter m according to a second plane conversion parameters of the M+N and institute
The second projection elevation face H is stated, the second fitting function that the scale parameter projects elevation face with described second is fitted；
Elevation face acquiring unit is projected, for according to second fitting function, calculating when scale parameter m is equal to described the
Corresponding best projection elevation face when two threshold values.
Wherein, the first threshold is 0, and the second threshold is 1.It is calculated and described first by above-mentioned Gauss model
The one-to-one first gauss projection coordinate of parameter combination that central meridian projects elevation face with described first, then by upper
The Parameter Switch model stated, calculate with the one-to-one first plane conversion parameter of the first gauss projection coordinate, according to described
Rotation angle parameter α and scale parameter m in first plane conversion parameter, determination are likely to occur corresponding projection when m ≈ 1 and α ≈ 0
The numberical range of elevation face and central meridian, i.e., the numberical range of described second central meridian and second projection
The numberical range in elevation face.Specifically can by calculate m ≈ 1 and α ≈ 0 in the numberical range of first central meridian and
The probability of occurrence of the numberical range in first projection elevation face, obtains according to the probability of occurrence of m ≈ 1 and α ≈ 0 in described second
Entreat meridianal numberical range and the numberical range in second projection elevation face.
In the same manner, it is calculated by above-mentioned Gauss model and projects elevation face with second central meridian and described second
The one-to-one second gauss projection coordinate of parameter combination calculated with second then by above-mentioned Parameter Switch model
The one-to-one second plane conversion parameter of gauss projection coordinate, according to rotation angle parameter α in the second plane conversion parameter
And scale parameter m, search for corresponding i-th of second plane conversion parameters when m ≈ 1 and α ≈ 0, wherein i-th of second planes
Translation parameters in conversion parameter is the best translation parameters.
First fitting function is：α=a1*L-b1；Second fitting function is：M=a2*H+b2；Wherein, a1、b1
For two constants of first fitting function；a2、b2For two constants of second fitting function.
In a kind of optional embodiment, the 2000 coordinate system parameter obtaining device further includes：
Internal accuracy computing module, for being sat according to the first gauss projection coordinate for overlapping control point, the plane
It is marked with and preset Parameter Switch model, calculates the internal accuracy for overlapping control point；
External check point acquisition module, the external check point for obtaining the multiple control net；Wherein, the external inspection
Epipole is to overlap the misaligned control point in control point with described in the multiple control net；
External accuracy computing unit, for according to the coincidence control point and the external check point, calculating described heavy
Close the external accuracy at control point；
Exception control point confirmation module, for according to the internal accuracy and external accuracy for overlapping control point, determining
Exception control point in the coincidence control point；
Exception control point rejects module, for rejecting the exception control point in the coincidence control point, and recalculates and picks
Except the corresponding first plane conversion parameter in coincidence control point after the exception control point.
The embodiment of the present invention additionally provides a kind of 2000 coordinate system parameter obtaining devices, including processor, memory
And it is stored in the memory and is configured as the computer program executed by the processor, the processor executes institute
2000 above-mentioned coordinate system parameter acquiring methods are realized when stating computer program.
Illustratively, the computer program can be divided into one or more module/units, one or more
A module/unit is stored in the memory, and is executed by the processor, to complete the present invention.It is one or more
A module/unit can be the series of computation machine program instruction section that can complete specific function, and the instruction segment is for describing institute
State implementation procedure of the computer program in the 2000 coordinate system parameter obtaining device.For example, the computer program
Coordinate obtaining module can be divided into, the CGCS2000 that control point is overlapped in multiple control nets for obtaining target cities is big
Ground coordinate and the plane coordinates for overlapping the corresponding target cities coordinate system in control point；First projective parameter setup module,
For the central meridian and perspective plane elevation according to the target cities coordinate system, the first projective parameter is set；First is high
This projection coordinate's computing module, for according to first projective parameter, being sat to the CGCS2000 the earth for overlapping control point
Mark carries out gauss projection, calculates the first gauss projection coordinate for overlapping control point；First plane conversion parameter calculating module,
For according to it is described overlap control point the first gauss projection coordinate, the plane coordinates and preset Parameter Switch model,
Calculate the first plane conversion parameter between the first gauss projection coordinate and the plane coordinates for overlapping control point；Coordinate
It is parameter acquisition module, for the rotation angle parameter and scale parameter according to the first plane conversion parameter, determines best
Translation parameters and best projection parameter.
The 2000 coordinate system parameter obtaining device can be desktop PC, notebook, palm PC and cloud
Hold the computing devices such as server.The 2000 coordinate system parameter obtaining device may include, but be not limited only to, and processor is deposited
Reservoir.It will be understood by those skilled in the art that the schematic diagram is only the example of 2000 coordinate system parameter obtaining devices,
The restriction to 2000 coordinate system parameter obtaining devices is not constituted, may include components more more or fewer than diagram, or
Person combines certain components or different components, such as the 2000 coordinate system parameter obtaining device can also include defeated
Enter output equipment, network access equipment, bus etc..
Alleged processor can be central processing unit (Central Processing Unit, CPU), can also be it
His general processor, digital signal processor (Digital Signal Processor, DSP), application-specific integrated circuit
(Application Specific Integrated Circuit, ASIC), ready-made programmable gate array (Field-
Programmable Gate Array, FPGA) either other programmable logic device, discrete gate or transistor logic,
Discrete hardware components etc..General processor can be microprocessor or the processor can also be any conventional processor
Deng the processor is the control centre of the 2000 coordinate system parameter obtaining device, and various interfaces and circuit is utilized to connect
Connect the various pieces of entire 2000 coordinate system parameter obtaining device.
The memory can be used for storing the computer program and/or module, and the processor is by running or executing
Computer program in the memory and/or module are stored, and calls the data being stored in memory, described in realization
The various functions of 2000 coordinate system parameter obtaining devices.The memory can include mainly storing program area and storage data
Area, wherein storing program area can storage program area, needed at least one function application program (such as sound-playing function,
Image player function etc.) etc.；Storage data field can be stored uses created data (such as audio data, electricity according to mobile phone
Script for story-telling etc.) etc..In addition, memory may include high-speed random access memory, can also include nonvolatile memory, such as
Hard disk, memory, plug-in type hard disk, intelligent memory card (Smart Media Card, SMC), secure digital (Secure
Digital, SD) card, flash card (Flash Card), at least one disk memory, flush memory device or other volatibility are solid
State memory device.
Wherein, if the integrated module/unit of the 2000 coordinate system parameter obtaining device is with SFU software functional unit
Form realize and when sold or used as an independent product, can be stored in a computer read/write memory medium.
Based on this understanding, the present invention realizes all or part of flow in above-described embodiment method, can also pass through computer journey
Sequence is completed to instruct relevant hardware, and the computer program can be stored in a computer readable storage medium, the meter
Calculation machine program is when being executed by processor, it can be achieved that the step of above-mentioned each embodiment of the method.Wherein, the computer program packet
Include computer program code, the computer program code can be source code form, object identification code form, executable file or
Certain intermediate forms etc..The computer-readable medium may include：Any reality of the computer program code can be carried
Body or device, recording medium, USB flash disk, mobile hard disk, magnetic disc, CD, computer storage, read-only memory (ROM, Read-
Only Memory), random access memory (RAM, Random Access Memory), electric carrier signal, telecommunication signal and
Software distribution medium etc..It should be noted that the content that the computer-readable medium includes can be according in jurisdiction
Legislation and the requirement of patent practice carry out increase and decrease appropriate, such as in certain jurisdictions, according to legislation and patent practice, meter
Calculation machine readable medium does not include electric carrier signal and telecommunication signal.
The embodiment of the present invention additionally provides a kind of computer readable storage medium, which is characterized in that described computer-readable
Storage medium includes the computer program of storage, wherein described computer-readable deposit is controlled when the computer program is run
Equipment where storage media executes 2000 above-mentioned coordinate system parameter acquiring methods.
Compared with the existing technology, a kind of beneficial effect for 2000 coordinate system parameter acquiring method that inventive embodiments provide
Fruit is：The 2000 coordinate system parameter acquiring method includes：It obtains and overlaps control in multiple control nets of target cities
The CGCS2000 geodetic coordinates of point and the plane coordinates for overlapping the corresponding target cities coordinate system in control point；According to described
The first projective parameter is arranged in the central meridian and perspective plane elevation of target cities coordinate system；According to the first projection ginseng
Number carries out gauss projection to the CGCS2000 geodetic coordinates for overlapping control point, and calculate the coincidence control point first is high
This projection coordinate；Turned according to the first gauss projection coordinate, the plane coordinates and the preset parameter for overlapping control point
Mold changing type calculates the first plane conversion ginseng between the first gauss projection coordinate and the plane coordinates for overlapping control point
Number；According to the rotation angle parameter and scale parameter of the first plane conversion parameter, best translation parameters and best is determined
Projective parameter.2000 coordinate system of target cities can be accurately obtained by the 2000 coordinate system parameter acquiring method
Projective parameter, to improve coordinate consistency between 2000 coordinate system of target cities and old coordinate system, so that new and old
Data seamless connection, while achievement can continue use under existing coordinate system, reduce conversion work amount and expense.The present invention is implemented
Example also provides 2000 coordinate system parameter obtaining devices of one kind and computer readable storage medium.
The above is the preferred embodiment of the present invention, it is noted that for those skilled in the art
For, various improvements and modifications may be made without departing from the principle of the present invention, these improvements and modifications are also considered as
Protection scope of the present invention.
Claims (10)
1. a kind of 2000 coordinate system parameter acquiring methods, which is characterized in that including：
It obtains and overlaps the CGCS2000 geodetic coordinates at control point and the coincidence control point in multiple control nets of target cities
The plane coordinates of corresponding target cities coordinate system；
According to the central meridian and perspective plane elevation of the target cities coordinate system, the first projective parameter is set；
According to first projective parameter, gauss projection is carried out to the CGCS2000 geodetic coordinates for overlapping control point, is calculated
The first gauss projection coordinate for overlapping control point；
According to it is described overlap control point the first gauss projection coordinate, the plane coordinates and preset Parameter Switch model,
Calculate the first plane conversion parameter between the first gauss projection coordinate and the plane coordinates for overlapping control point；
According to the rotation angle parameter and scale parameter of the first plane conversion parameter, best translation parameters and best is determined
Projective parameter.
2. 2000 coordinate system parameter acquiring method as described in claim 1, which is characterized in that described according to described first
The rotation angle parameter and scale parameter of plane conversion parameter determine best translation parameters and best projection parameter, specific to wrap
It includes；
According to the first plane conversion parameter, the second projective parameter is set；
According to second projective parameter, gauss projection is carried out to the CGCS2000 geodetic coordinates for overlapping control point, is calculated
The second gauss projection coordinate for overlapping control point；
According to the second gauss projection coordinate, the plane coordinates and the preset Parameter Switch mould for overlapping control point
Type, calculates the second gauss projection coordinate for overlapping control point and second between the plane coordinates for overlapping control point is flat
Face conversion parameter；
According to the rotation angle parameter and scale parameter of the second plane conversion parameter, determine the best translation parameters and
The best projection parameter.
3. 2000 coordinate system parameter acquiring method as described in claim 1, which is characterized in that described according to the target
The central meridian and perspective plane elevation of city coordinate system are arranged the first projective parameter, specifically include：
First projective parameter includes the first central meridian and the first perspective plane elevation；
According to the central meridian of the target cities coordinate system, the numerical intervals of first central meridian are determined, and press
According to the first interval of setting, first central meridian is set；
According to the projection elevation face of the target cities coordinate system, the numerical intervals in first projection elevation face are determined, and press
According to the second interval setting first projection elevation face of setting.
4. 2000 coordinate system parameter acquiring method as claimed in claim 2, which is characterized in that described according to described first
Plane conversion parameter is arranged the second projective parameter, specifically includes：
Second projective parameter includes the second central meridian and the second perspective plane elevation；
According to the first plane conversion parameter, the numerical intervals of second central meridian and second projection are determined
The numerical intervals of face elevation；
Second central meridian is set according to the third interval of setting and according to the 4th interval setting described the of setting
Two projection elevation faces.
5. 2000 coordinate system parameter acquiring method as claimed in claim 4, which is characterized in that described according to the of setting
Second central meridian and the 4th interval setting second projection elevation face according to setting is arranged in three intervals, specifically
Including；
When second central meridian is L, according to the 4th interval setting first projection elevation face of the setting, obtain
Obtain the parameter combination that N groups second central meridian projects elevation face with described second；Wherein, L is to belong in described second
It entreats in meridianal numerical intervals and according to the spaced central meridian of the third of setting；
When second projection elevation face is H, second central meridian is set according to the third interval of setting, obtains M
The parameter combination that group second central meridian projects elevation face with described second；H is to belong to second projection elevation face
Numerical intervals in and according to the 4th spaced projection elevation face of setting.
6. 2000 coordinate system parameter acquiring method as claimed in claim 5, which is characterized in that described according to described second
The rotation angle parameter and scale parameter of plane conversion parameter determine best translation parameters and best projection parameter, specific to wrap
It includes：
The rotation angle parameter α ≈ for meeting the second plane conversion parameter are obtained from M+N the second plane conversion parameters
I-th of second plane conversion parameters of the scale parameter m ≈ second thresholds of first threshold and the second plane conversion parameter；
Wherein, the parameter combination and M groups that elevation face are projected with described second according to N groups second central meridian described the
The parameter combination that two central meridians project elevation face with described second can get M+N the second plane conversion parameters；
Obtain the translation parameters of described i-th second plane conversion parameters；Wherein, described i-th second plane conversion parameters
Translation parameters is the best translation parameters；
The best projection parameter includes best central meridian and best projection elevation face；
According to the rotation angle parameter α and the second central meridian L of the M+N the second plane conversion parameters, fitting
Go out the first fitting function of the rotation angle parameter and second central meridian；
According to first fitting function, corresponding best center when rotation angle parameter α is equal to the first threshold is calculated
Noon line；
Elevation face H is projected according to the scale parameter m of the M+N the second plane conversion parameters and described second, is fitted
The second fitting function that the scale parameter projects elevation face with described second；
According to second fitting function, the corresponding best projection elevation face when scale parameter m is equal to second threshold is calculated.
7. 2000 coordinate system parameter acquiring method as described in claim 1, which is characterized in that 2000 independent coordinate
It is that parameter acquiring method further includes：
According to it is described overlap control point the first gauss projection coordinate, the plane coordinates and preset Parameter Switch model,
Calculate the internal accuracy for overlapping control point；
Obtain the external check point of the multiple control net；Wherein, the external check point be during the multiple control is netted with institute
It states and overlaps the misaligned control point in control point；
According to the coincidence control point and the external check point, the external accuracy for overlapping control point is calculated；
According to the internal accuracy and external accuracy for overlapping control point, the exception control in the coincidence control point is determined
Point；
The exception control point in the coincidence control point is rejected, and recalculates the coincidence control after rejecting the exception control point
The corresponding first plane conversion parameter of point.
8. a kind of 2000 coordinate system parameter obtaining devices, which is characterized in that including：
Coordinate obtaining module, overlap in multiple control nets for obtaining target cities the CGCS2000 geodetic coordinates at control point with
And the plane coordinates for overlapping the corresponding target cities coordinate system in control point；
First projective parameter setup module, for high according to the central meridian of the target cities coordinate system and perspective plane
The first projective parameter is arranged in journey；
First gauss projection coordinate calculation module is used for according to first projective parameter, to the coincidence control point
CGCS2000 geodetic coordinates carries out gauss projection, calculates the first gauss projection coordinate for overlapping control point；
First plane conversion parameter calculating module, for according to the first gauss projection coordinate for overlapping control point, described flat
Areal coordinate and preset Parameter Switch model calculate the first gauss projection coordinate for overlapping control point and are sat with the plane
The first plane conversion parameter between mark；
Coordinate system parameters acquisition module is used for rotation angle parameter and scale parameter according to the first plane conversion parameter,
Determine best translation parameters and best projection parameter.
9. a kind of 2000 coordinate system parameter obtaining devices, including processor, memory and it is stored in the memory
And it is configured as the computer program executed by the processor, the processor is realized when executing the computer program as weighed
Profit requires 2000 coordinate system parameter acquiring methods described in any one of 1 to 7.
10. a kind of computer readable storage medium, which is characterized in that the computer readable storage medium includes the calculating of storage
Machine program, wherein equipment where controlling the computer readable storage medium when the computer program is run is executed as weighed
Profit requires 2000 coordinate system parameter acquiring methods described in any one of 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810215665.2A CN108731648B (en) | 2018-03-15 | 2018-03-15 | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810215665.2A CN108731648B (en) | 2018-03-15 | 2018-03-15 | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108731648A true CN108731648A (en) | 2018-11-02 |
| CN108731648B CN108731648B (en) | 2020-12-22 |
Family
ID=63940492
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810215665.2A Active CN108731648B (en) | 2018-03-15 | 2018-03-15 | 2000 independent coordinate system parameter obtaining method, device and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108731648B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114001650A (en) * | 2021-09-16 | 2022-02-01 | 北京市测绘设计研究院 | Method for encrypting conversion parameters of earth coordinate system and arbitrary plane coordinate system |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3595451B2 (en) * | 1998-07-24 | 2004-12-02 | ペンタックス株式会社 | Laser surveying equipment |
| CN103389080A (en) * | 2013-08-19 | 2013-11-13 | 重庆市地理信息中心 | Method for obtaining independent city coordinate system parameter based on geographic information application |
| CN103413352A (en) * | 2013-07-29 | 2013-11-27 | 西北工业大学 | Scene three-dimensional reconstruction method based on RGBD multi-sensor fusion |
| CN103827631A (en) * | 2011-09-27 | 2014-05-28 | 莱卡地球系统公开股份有限公司 | Measuring system and method for marking a known target point in a coordinate system |
| CN104123695A (en) * | 2014-06-23 | 2014-10-29 | 唐文兴 | Method for realizing coordinate conversion |
| CN104567840A (en) * | 2015-01-27 | 2015-04-29 | 国家测绘地理信息局大地测量数据处理中心 | New and old urban independent coordinate system fusion method |
| CN104598746A (en) * | 2015-01-27 | 2015-05-06 | 国家测绘地理信息局大地测量数据处理中心 | Applicability discrimination method for coordinate transformation model |
| CN104766302A (en) * | 2015-02-05 | 2015-07-08 | 武汉大势智慧科技有限公司 | Method and system for optimizing laser scanning point cloud data by means of unmanned aerial vehicle images |
| CN105388507A (en) * | 2015-11-24 | 2016-03-09 | 中国电建集团西北勘测设计研究院有限公司 | Method for determining ellipsoid parameters based on area GNSS and precise distance measurement scale ratio estimation method |
| CN106202000A (en) * | 2016-07-22 | 2016-12-07 | 武汉大学 | Seven-parameter transformation method between country's three-dimensional system of coordinate and anywhere plane coordinate system |
| CN106289195A (en) * | 2016-08-29 | 2017-01-04 | 长江空间信息技术工程有限公司(武汉) | The method for building up of plateau distance heavy construction Measurement and Control System |
| CN107167119A (en) * | 2017-03-29 | 2017-09-15 | 中交烟台环保疏浚有限公司 | The data processing method of distortion of projection |
-
2018
- 2018-03-15 CN CN201810215665.2A patent/CN108731648B/en active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP3595451B2 (en) * | 1998-07-24 | 2004-12-02 | ペンタックス株式会社 | Laser surveying equipment |
| CN103827631A (en) * | 2011-09-27 | 2014-05-28 | 莱卡地球系统公开股份有限公司 | Measuring system and method for marking a known target point in a coordinate system |
| CN103413352A (en) * | 2013-07-29 | 2013-11-27 | 西北工业大学 | Scene three-dimensional reconstruction method based on RGBD multi-sensor fusion |
| CN103389080A (en) * | 2013-08-19 | 2013-11-13 | 重庆市地理信息中心 | Method for obtaining independent city coordinate system parameter based on geographic information application |
| CN104123695A (en) * | 2014-06-23 | 2014-10-29 | 唐文兴 | Method for realizing coordinate conversion |
| CN104567840A (en) * | 2015-01-27 | 2015-04-29 | 国家测绘地理信息局大地测量数据处理中心 | New and old urban independent coordinate system fusion method |
| CN104598746A (en) * | 2015-01-27 | 2015-05-06 | 国家测绘地理信息局大地测量数据处理中心 | Applicability discrimination method for coordinate transformation model |
| CN104766302A (en) * | 2015-02-05 | 2015-07-08 | 武汉大势智慧科技有限公司 | Method and system for optimizing laser scanning point cloud data by means of unmanned aerial vehicle images |
| CN105388507A (en) * | 2015-11-24 | 2016-03-09 | 中国电建集团西北勘测设计研究院有限公司 | Method for determining ellipsoid parameters based on area GNSS and precise distance measurement scale ratio estimation method |
| CN106202000A (en) * | 2016-07-22 | 2016-12-07 | 武汉大学 | Seven-parameter transformation method between country's three-dimensional system of coordinate and anywhere plane coordinate system |
| CN106289195A (en) * | 2016-08-29 | 2017-01-04 | 长江空间信息技术工程有限公司(武汉) | The method for building up of plateau distance heavy construction Measurement and Control System |
| CN107167119A (en) * | 2017-03-29 | 2017-09-15 | 中交烟台环保疏浚有限公司 | The data processing method of distortion of projection |
Non-Patent Citations (4)
| Title |
|---|
| YANG Y 等: "Combined Adjustment Project of National Astronomical Geodetic Networks and 2000 National GPS Control Network", 《PROGRESS IN NATIONAL SCIENCE》 * |
| 王维 等: "参心独立坐标系向地心独立坐标系转换地块面积变化分析和计算方法研究", 《测绘与空间地理信息》 * |
| 耿晓燕: "地方独立坐标系向2000国家大地坐标系转换研究", 《中国优秀硕士学位论文全文数据库基础科学辑》 * |
| 蔡元波 等: "宁波市2000坐标系建立及推广应用研究", 《城市勘探》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114001650A (en) * | 2021-09-16 | 2022-02-01 | 北京市测绘设计研究院 | Method for encrypting conversion parameters of earth coordinate system and arbitrary plane coordinate system |
| CN114001650B (en) * | 2021-09-16 | 2023-09-29 | 北京市测绘设计研究院 | Encryption method for conversion parameters of local coordinate system and arbitrary plane coordinate system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108731648B (en) | 2020-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105785410B (en) | It is a kind of that the method and system of location information is provided | |
| CN105163337B (en) | A method of the mobile network data geography mapping based on coverage prediction emulation | |
| CN109451532B (en) | Method and device for checking position of base station | |
| CN108513305B (en) | Base station address selection method, device and equipment for target area | |
| CN107277772A (en) | A kind of wireless access independent positioning method and device, computer-readable recording medium | |
| CN109688532A (en) | A kind of method and device dividing city function region | |
| CN108267154A (en) | A kind of map-indication method and device | |
| CN103822638B (en) | The treating method and apparatus of customer position information | |
| CN109670005A (en) | A kind of location matching method of BIM model and three-dimensional geographic scenes | |
| CN107976702A (en) | A kind of position correcting method based on CORS, positioning terminal and alignment system | |
| CN108375985A (en) | A kind of soil three-dimensional planning and designing platform and its design method | |
| CN110351665B (en) | Method, apparatus and computer-readable storage medium for habitual identification of a user | |
| CN111935733B (en) | Method, device and equipment for checking base station position information and computer storage medium | |
| CN109949692A (en) | Road network method, apparatus, computer equipment and storage medium | |
| CN108897937B (en) | Method for automatically converting CAD (computer aided design) data of civil aviation airport into DEM (digital elevation model) data | |
| CN108731649B (en) | 2000 mapping reference frame unifying method, device and computer readable storage medium | |
| CN108416514B (en) | Neighbor cell planning method and device based on geographic information system | |
| CN105844031A (en) | Mobile phone positioning data based urban traffic corridor identification method | |
| CN107169080A (en) | A kind of geospatial analysis system being combined based on GIS and spatial database | |
| CN108731648A (en) | 2000 coordinate system parameter acquiring methods, device and computer readable storage medium | |
| CN109163708A (en) | Optimization method, device and the medium of unmanned plane large scale topographical map house line | |
| CN107784408A (en) | Wind resource assessment method, device and system based on terrain classification | |
| CN106255126A (en) | A kind of method and device identifying to arrange flowers in base station location district | |
| CN115588086A (en) | Map dividing method, map dividing device, computer readable storage medium and processor | |
| CN108156657A (en) | A kind of method for positioning LTE network high speed user |
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 |