CN103389080A - 一种基于地理信息应用城市独立坐标系参数获取方法 - Google Patents
一种基于地理信息应用城市独立坐标系参数获取方法 Download PDFInfo
- Publication number
- CN103389080A CN103389080A CN2013103617329A CN201310361732A CN103389080A CN 103389080 A CN103389080 A CN 103389080A CN 2013103617329 A CN2013103617329 A CN 2013103617329A CN 201310361732 A CN201310361732 A CN 201310361732A CN 103389080 A CN103389080 A CN 103389080A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- urban
- independent
- formula
- parameter
- 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
Images
Landscapes
- Navigation (AREA)
Abstract
本发明请求保护一种基于地理信息应用城市独立坐标系参数获取方法,涉及利用高斯投影技术,计算该城市独立坐标系近似参数,从而实现地理信息成果相互转换的方法。本发明通过获取三个某城市独立坐标系与国家通用坐标系同名点成果,高斯投影正算公式进行变形,简化成一个关于该独立坐标系中央子午线L0的一个一元高次方程,并可利用迭代法解算出中央子午线L0,进而解算该独坐标系其他参数,包括:参考椭球长半径a,中央子午线L0,北方向偏移量X0和东方向偏移量Y0。这样便可以实现该独立坐标系与国家通用坐标系之间的同名点成果转换与应用。本发明可广泛用于地理空间数据建设,数字城市建设等与地理空间信息坐标转换应用。
Description
技术领域
本发明涉及利用高斯投影技术,通过估算城市独立坐标系参数,从而实现地理信息在各坐标系间相互转换的方法,具体是一种基于地理信息应用城市独立坐标系参数获取方法。
背景技术
随着计算机及信息技术的进步,地理信息产业取得了长足的发展。作为地理信息产业核心的资产──地理信息数据的重要性也日益凸显。我国现有的地理空间信息横向上分布在不同委办局,纵向上又分散在市、区(县)、乡(镇)等多级行政级别,呈现出典型的“多重分布”的复杂特性,在坐标体系上,我国先后建成了北京54坐标系、西安80坐标系、2000国家坐标系,国际上较通用的有WGS84坐标系,而在城市规划、建设、管理等实际工作中,又多采用城市独立坐标系,导致在地理信息应用过程中,需要通过坐标转换来整合不同空间数据。
目前,坐标转换理论很成熟,相关的方法研究也比较多,最常用的有四参数、七参数。其中四参数是平面仿射变换,比较简单,缺点是没有考虑椭球体的变形,不能进行大范围高精度的坐标转换。而七参数转换基于参考椭球,理论严谨,转换精度高,控制范围大,是地理信息应用中常用的一种方法。但七参数法需要知道各城市独立坐标系的参数,但这些数据涉密,在实际应用中很难获取。
建立城市独立坐标系的目的是为了满足投影变形小于2.5cm/km的要求。我国大地测量法和有关测量规范中明确规定,国家大地测量控制网依高斯投影方法按3°带或6°带进行分带和计算,并把观测成果归算到参考椭球面上。这样的规定,不但符合高斯投影的分带原则和计算方法,与国际惯例相一致,而且也便于大地测量成果的统一、使用和互算。但由于我国国土辽阔,地形复杂,我国的北京54坐标系、西安80坐标系的参考椭球长半径比国际椭球半径大,与地球表面存在一定的套合差,所以国家在对原有北京54坐标系、西安80坐标系进行优化的基础上,建立2000国家坐标系。同时根据高斯投影变化规律,离中央子午线越远,变形越大,按照国家坐标分带不能满足变形小于2.5cm/km的要求,因此我国大部分城市和地区布设城市测量控制网都基于北京54坐标系、西安80坐标系或2000国家坐标系参考椭球,采用任意带的独立高斯投影平面直角坐标系,归算测量结果的参考面也是根据需要选定的方法来建立城市独立坐标系。
在计算城市独立坐标系参数中,通常利用高斯投影理论,进行大地坐标和高斯平面直角坐标间的相互转换,其中城市独立坐标系采用高斯平面直角坐标,北京54坐标系、西安80坐标系或2000国家坐标系采用大地坐标,方法如下:
1、大地坐标(B,L,H):大地坐标是采用大地纬度B、经度L和大地高程H来描述空间位置的。纬度是空间的点与参考椭球面的法线与赤道面的夹角;经度是空间的点与参考椭球的自转轴所在的面与参考椭球的起始子午面的夹角λ;大地高程是空间的点沿着参考椭球的法线方向到参考椭球面的距离,参见图1。
2、高斯平面直角坐标(x,y,H):其中,x,y表示平面坐标,H表示大地高程。椭球体中心O在椭圆柱中心轴上,椭球体南北极与椭圆柱相切,并使某一子午线与椭圆柱相切。此子午线称中央子午线。然后将椭球体面上的点、线按正形投影条件投影到椭圆柱上,再沿椭圆柱N、S点母线割开,并展成平面,即成为高斯投影平面,参见图2。在此平面上:
①中央子午线是直线,其长度不变形,离开中央子午线的其他子午线是弧形,凹向中央子午线。离开中央子午线越远,变形越大。
②投影后赤道是一条直线,赤道与中央子午线保持正交。
③离开赤道的纬线是弧线,凸向赤道。
3、空间直角坐标(X,Y,Z)的坐标原点位于参考椭球的中心,Z轴指向参考椭球的北极,X轴指向起始子午面与赤道的交点,Y轴位于赤道面上切按右手系于X轴呈90度夹角,某点中的坐标可用该点在此坐标系的各个坐标轴上的投影来表示。参见图1。
大地坐标向高斯平面坐标转换(高斯投影正算公式):
X=a(1-e2)(A′arcB-B′sin2B+C′sin4B-D′sin6B+E′sin8B-F′sin10B+G′sin12B)
其中:
ρ=180×3600/π
e─第一偏心率
e′─第二偏心率
η2=e'2cos2B
t=tgB
B为纬度,单位弧度
M─子午圈曲率半径
N─卯酉圈曲率半径
在进行城市独立坐标系与国家通用坐标系进行转换时,必须要知道城市独立坐标系的参数,但城市独立坐标系的参数涉密,在实际应用中很难获取。因此要实现城市独立坐标系与国家通用坐标系之间的相互转换,亟需一种能获得城市独立坐标系的参数的方法。
发明内容
本发明的目的在于利用高斯投影技术,计算城市独立坐标系参数,包括参考椭球扁率,参考椭球长半径,中央子午线,北方向偏移量和东方向偏移量,在避免城市独立坐标系参数的公布的前提下,而实现城市独立坐标系与北京54坐标系、西安80坐标系或2000国家坐标系成果之间的相互转换。
通常城市独立坐标系的参数包括参考椭球扁率e,参考椭球长半径a,中央子午线L0,北方向偏移量X0和东方向偏移量Y0,其中参考椭球扁率e已知,与国家通用坐标系参考椭球保持一致,其他几个参数未知。所以确定一个城市独立坐标系参数,就是要计算出参考椭球长半径a,中央子午线L0,北方向偏移量X0和东方向偏移量Y0四个参数。已知某城市独立坐标系和国家通用坐标系两种坐标系下的同名点成果各三个,通过对高斯投影正算公式进行变换,得到一个关于中央子午线L0的一元高次方程,可以通过迭代算法解算出L0,进而可以计算出椭球长半径,北方向偏移量和东方向偏移量,从而获取城市独立坐标系参数。因而本发明的技术方案是:
1)利用GPS联测某城市独立坐标系四等以上的控制点(为了提高坐标计算的精度,所以同名点成果的测量等级尽量高),获得至少三个该城市独立坐标系同名点成果P1、P2、P3,在与测得所述城市独立坐标系同名点成果相同的三点上利用GPS与国家通用坐标系联测四等以上控制点,得到国家通用坐标系同名点成果P1′、P2′、P3′;
2)记
高斯投影正算公式简化为:
x=X+NF(l) (3)
高斯投影正算公式简化为:
y=NF1(l) (4)
其中,经差l=L-L0,以弧度为单位,L0为中央子午线,L为国家通用坐标系同名点成果的经度,B为国家通用坐标系同名点成果的纬度,x,y为城市独立坐标系同名点成果的平面坐标。
3)由于城市独立坐标系参考椭球通常基于国家通用坐标系参考椭球,采用膨胀法建设,所以城市独立坐标系参考椭球的扁率e总是与对应国家通用坐标系参考椭球的扁率相同,但是二者参考椭球的长半径a不同。这这种情况下城市独立坐标系的经度变化△L=0,纬度变化△B≈0,将P1和P1′、P2和P2′、P3和P3′任意两对同名点成果分别带入式(3)、式(4)整理可得:
记任意两对同名点成果为:国家通用坐标系的记为(B1,L1,H1)和(B2,L2,H2),城市独立坐标系的记为(x1,y1,H1)和(x2,y2,H2)。
x1=X0+X1+N1F(L1-L0) (5)
y1=Y0+N1F1(L1-L0) (6)
x2=X0+X2+N2F(L2-L0) (7)
y2=Y0+N2F1(L2-L0) (9)
x1-x2=X1-X2+N1F(L1-L0)-N2F(L2-L0) (10)
式(11)是一个关于L0的一元高次方程,利用高斯迭代算法可以解出L0。
将L0代入式(10)可以求解出参数椭球长半径a
将a,L0分别代入式(5)、式(6)求解出北方向偏移量X0和东方向偏移量Y0。
这样便求得了城市独立坐标系的各个参数。
4)国家通用坐标系参数是已知的,可查阅的,根据步骤3)中获得的城市独立坐标系参数(中央子午线L0、椭球长半径a、北方偏移量X0和东方偏移量Y0),利用P1、P2、P3和P1′、P2′、P3′的两套坐标系同名点成果,基于Bursa七参数转换模型,计算出坐标转换参数,转换参数包括3个平移参数,3个旋转参数和一个比例因子,从而实现坐标转换。
Bursa七参数转换模型:
式中,[△X △Y △Z]T为3个平移参数,[εX εZ εZ]T为3个旋转参数和m为比例因子,[XS YS ZS]T,[XT YT ZT]T分别为城市独立坐标系同名点成果和国家通用坐标系同名点成果。
具体地,上述国家通用坐标系包括北京54坐标系、西安80坐标系和2000国家坐标系。国家通用坐标系用其他已知坐标系替换,也能使用本发明的方法。
本发明的有益效果是:利用高斯投影技术,只需城市独立坐标系和国家通用坐标系三个同名点成果,就可以近似计算出城市独立坐标系参数,实现两个坐标系间成果的高精度转换,完全满足地理信息坐标转换的需要,避免了独立坐标系参数的公布。本发明可广泛用于不同坐标系的地理空间数据整合和坐标换算,以及在数字城市和智慧城市建设中,进行多源数据空间套合和分析等。
附图说明
图1为大地坐标系示意图;
图2为高斯投影示意图。
具体实施方式
由于城市独立坐标系的参数是涉密的,很难获取,所以不能将本发明方法获得的城市独立坐标系参数与实际的城市独立坐标系参数进行比较,来验证方法的有效性。因此将用本方法获得的城市独立坐标系参数,利用七参数转换模型,将国家通用坐标系的同名点成果转换为独立坐标系转换成果,然后将独立坐标系转换成果与实际测得的独立坐标系同名点成果进行比较,来验证本发明方法的有效性。具体参见以下实施例。
实施例1,某城市独立坐标系与北京54坐标系的转换
首先,利用GPS分别与北京54坐标与某城市独立坐标联测,测得同名点成果8个,如下表1所示。
表1某城市独立坐标系与北京54坐标系同名点成果表
其次,利用1、2点成果根据式11,式10,式5及式6,计算出的坐标参数。
其中:
中央子午线:108:45:05.029260
椭球长半径:6379489.37490557
北方偏移量:-6.15235514800548米
东方偏移量:500144.292887685米
第三,利用该坐标系参数,将1、2、3点成果带入七参数转换模型,计算出七参数为:
X偏移:-5.90702772431541
Y偏移:77.6861151843914
Z偏移:50.5046847679187
X旋转:000:00:00.361474
Y旋转:-000:00:00.356650
Z旋转:-000:00:00.976716
比例因子:0.000180875238540795
第四,利用该七参数进行坐标转换后,比较转换结果与测量值,最大差值为0.036米,完全满足地理信息坐标转换的需要,见表2。
表2某城市独立坐标系与北京54坐标系坐标转换成果较差统计表
实施例2,某城市独立坐标系与西安80坐标系的转换
首先,利用GPS分别与西安80坐标系与某城市独立坐标联测,测得同名点成果5个,如下表3所示。
表3某城市独立坐标系与西安80坐标系的转换同名点成果表
其次,利用1、2点成果根据式11,式10,式5及式6,计算出的坐标参数。
其中:
中央子午线:105:53:58.943278
椭球长半径:6378423.47268279
北方偏移量:-3199971.81273551米
东方偏移量:49971.3935643691米
第三,利用该坐标系参数,将1、2、3点成果带入七参数转换模型,计算出七参数为:
X偏移:-2.75906461090835
Y偏移:-7.97653608962355
Z偏移:-7.82526087134284
X旋转:000:00:00.086500
Y旋转:-000:00:00.710559
Z旋转:-000:00:00.215598
比例因子:0.0000459872988582905
第四,利用该七参数进行坐标转换后,转换结果与测量值进行比较,最大差值为0.0272米,完全满足地理信息坐标转换的需要,见表4。
表4某城市独立坐标系与西安80坐标系的转换坐标转换成果较差统计表
实施例3,某城市独立坐标系与2000国家坐标系的转换
首先,利用GPS分别与2000国家坐标与某城市独立坐标联测,测得同名点成果9个,如下表5所示。
表5独立坐标系与2000国家坐标系同名点成果表
其次,利用1、2点成果根据式11,式10,式5及式6,计算出的坐标参数。
其中:
中央子午线:108:22:17.492463
椭球长半径:6378421.52585621
北方偏移量:-1.9293381461545米
东方偏移量:499965.255096663米
第三,利用该坐标系参数,将1、2、3点成果带入七参数转换模型,计算出七参数为:
X偏移:-0.0144736003057915
Y偏移:1.45046010742954
Z偏移:3.63460416544149
X旋转:000:00:00.087127
Y旋转:-000:00:00.012301
Z旋转:-000:00:00.039853
比例因子:0.0000441302554666525
第四,利用该七参数进行坐标转换后,转换结果与测量结果进行比较,最大差值为0.042米,完全满足地理信息坐标转换的需要,见表6。
表6某城市独立坐标系与2000国家坐标系坐标转换成果较差统计表
本发明利用高斯投影理论,在已知至少三个某城市独立坐标系与国家通用坐标系的同名点,未知该独立坐标系参数的情况下,用于计算该城市的近似坐标系参数,实现地理信息成果的高精度转换,避免了城市独立坐标系参数的公布,促进了地理信息系统的应用与推广。
Claims (2)
1.一种基于地理信息应用城市独立坐标系参数获取方法,其特征在于:
1)利用GPS联测某城市独立坐标系四等以上的控制点,获得至少三个该城市独立坐标系同名点成果P1、P2、P3,在与测得所述城市独立坐标系同名点成果相同的三点上利用GPS与国家通用坐标系联测四等以上控制点,得到国家通用坐标系同名点成果P1′、P2′、P3′;
2)将高斯投影正算公式简化为:
x=X+NF(l) (3)
y=NF1(l) (4)
式中:
其中,经差l=L-L0,以弧度为单位,L0为中央子午线,L为国家通用坐标系同名点成果的经度,B为国家通用坐标系同名点成果的纬度,x,y为城市独立坐标系同名点成果的平面坐标;
3)城市独立坐标系参考椭球的扁率e与国家通用坐标系参考椭球的扁率相同,但是两种参考椭球的长半径a不同,在这种情况下城市独立坐标系的经度变化△L=0,纬度变化△B≈0,将P1和P1′、P2和P2′、P3和P3′任意两对分别带入式(3)、式(4)整理可得:
x1=X0+X1+N1F(L1-L0) (5)
y1=Y0+N1F1(L1-L0) (6)
x2=X0+X2+N2F(L2-L0) (7)
y2=Y0+N2F1(L2-L0) (9)
x1-x2=X1-X2+N1F(L1-L0)-N2F(L2-L0) (10)
式(11)是一个关于L0的一元高次方程,利用高斯迭代算法解出中央子午线L0;
将L0代入式(10)求解出椭球长半径a
将a,L0分别代入式(5)、式(6)求解出北方偏移量X0和东方偏移量Y0;
4)根据步骤3)中获得的城市独立坐标系参数,利用P1、P2、P3和P1′、P2′、P3′的两套坐标系同名点成果,利用Bursa七参数转换模型,计算出坐标转换参数,所述转换参数包括3个平移参数,3个旋转参数和1个比例因子;
Bursa七参数转换模型:
式中,[△X △Y △Z]T为3个平移参数,[εX εZ εZ]T为3个旋转参数,m为比例因子,[XS YS ZS]T,[XT YT ZT]T分别为城市独立坐标系和国家通用坐标系同名点成果。
2.根据权利要求1所述一种基于地理信息应用城市独立坐标系参数获取方法,其特征在于:所述国家通用坐标系包括北京54坐标系、西安80坐标系和2000国家坐标系。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310361732.9A CN103389080B (zh) | 2013-08-19 | 2013-08-19 | 一种基于地理信息应用城市独立坐标系参数获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310361732.9A CN103389080B (zh) | 2013-08-19 | 2013-08-19 | 一种基于地理信息应用城市独立坐标系参数获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103389080A true CN103389080A (zh) | 2013-11-13 |
CN103389080B CN103389080B (zh) | 2015-08-05 |
Family
ID=49533426
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310361732.9A Active CN103389080B (zh) | 2013-08-19 | 2013-08-19 | 一种基于地理信息应用城市独立坐标系参数获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103389080B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123695A (zh) * | 2014-06-23 | 2014-10-29 | 唐文兴 | 一种实现坐标转换方法 |
CN104567840A (zh) * | 2015-01-27 | 2015-04-29 | 国家测绘地理信息局大地测量数据处理中心 | 新旧城市独立坐标系融合方法 |
CN104598746A (zh) * | 2015-01-27 | 2015-05-06 | 国家测绘地理信息局大地测量数据处理中心 | 坐标转换模型适用性判别方法 |
CN105869127A (zh) * | 2016-04-14 | 2016-08-17 | 中国电建集团西北勘测设计研究院有限公司 | 基于指定轴线变形最小标准减小任意带高斯投影变形方法 |
CN105869128A (zh) * | 2016-04-14 | 2016-08-17 | 中国电建集团西北勘测设计研究院有限公司 | 基于最大投影变形最小标准减小任意带高斯投影变形方法 |
CN105931199A (zh) * | 2016-04-14 | 2016-09-07 | 中国电建集团西北勘测设计研究院有限公司 | 基于综合投影变形最小标准减小任意带高斯投影变形方法 |
CN105976314A (zh) * | 2016-05-10 | 2016-09-28 | 福州市勘测院 | 顾及相同中央投影经线的不同参考椭球投影平面坐标系统转换方法 |
CN106202000A (zh) * | 2016-07-22 | 2016-12-07 | 武汉大学 | 国家三维坐标系与任意地方平面坐标系间七参数转换方法 |
CN106341936A (zh) * | 2016-08-31 | 2017-01-18 | 大连海事大学 | 一种内河智能航标位置监控方法 |
CN106683404A (zh) * | 2016-12-06 | 2017-05-17 | 华南理工大学 | 一种通过手机定位技术获取公交客流od的方法 |
CN107908808A (zh) * | 2017-08-11 | 2018-04-13 | 山东交通学院 | 基于抵偿高程面或平均高程面的3°分带坐标转换系统 |
CN108731648A (zh) * | 2018-03-15 | 2018-11-02 | 广州市城市规划勘测设计研究院 | 2000独立坐标系参数获取方法、装置与计算机可读存储介质 |
CN110427450A (zh) * | 2019-08-09 | 2019-11-08 | 深圳四维集思技术服务有限公司 | 基于互联网地图的油气管线数据处理方法和系统 |
CN114509051A (zh) * | 2022-01-26 | 2022-05-17 | 中交二公局第三工程有限公司 | 一种道路工程中无人机测量放样方法 |
CN114546941A (zh) * | 2022-04-28 | 2022-05-27 | 长江水利委员会水文局 | 一种数字地形图空间坐标系与投影的全信息无损转换方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030212528A1 (en) * | 2002-01-25 | 2003-11-13 | Chandler Larry S. | Inversion-conforming data sets processing |
CN102223709A (zh) * | 2011-06-10 | 2011-10-19 | 重庆市地理信息中心 | 基于cors系统实时获取特定基准的高斯平面直角坐标的方法 |
CN102288158A (zh) * | 2011-08-02 | 2011-12-21 | 甘肃铁道综合工程勘察院有限公司 | 基于法截面子午线椭球控制高斯投影变形的方法 |
-
2013
- 2013-08-19 CN CN201310361732.9A patent/CN103389080B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030212528A1 (en) * | 2002-01-25 | 2003-11-13 | Chandler Larry S. | Inversion-conforming data sets processing |
CN102223709A (zh) * | 2011-06-10 | 2011-10-19 | 重庆市地理信息中心 | 基于cors系统实时获取特定基准的高斯平面直角坐标的方法 |
CN102288158A (zh) * | 2011-08-02 | 2011-12-21 | 甘肃铁道综合工程勘察院有限公司 | 基于法截面子午线椭球控制高斯投影变形的方法 |
Non-Patent Citations (2)
Title |
---|
张黎等: "CORS系统实时地方坐标测量的保密模式研究", 《城市勘测》 * |
耿晓燕: "地方独立坐标系向2000国家大地坐标系转换研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104123695A (zh) * | 2014-06-23 | 2014-10-29 | 唐文兴 | 一种实现坐标转换方法 |
CN104567840A (zh) * | 2015-01-27 | 2015-04-29 | 国家测绘地理信息局大地测量数据处理中心 | 新旧城市独立坐标系融合方法 |
CN104598746A (zh) * | 2015-01-27 | 2015-05-06 | 国家测绘地理信息局大地测量数据处理中心 | 坐标转换模型适用性判别方法 |
CN105869127B (zh) * | 2016-04-14 | 2018-10-30 | 中国电建集团西北勘测设计研究院有限公司 | 基于指定轴线变形最小标准减小任意带高斯投影变形方法 |
CN105869127A (zh) * | 2016-04-14 | 2016-08-17 | 中国电建集团西北勘测设计研究院有限公司 | 基于指定轴线变形最小标准减小任意带高斯投影变形方法 |
CN105869128A (zh) * | 2016-04-14 | 2016-08-17 | 中国电建集团西北勘测设计研究院有限公司 | 基于最大投影变形最小标准减小任意带高斯投影变形方法 |
CN105931199A (zh) * | 2016-04-14 | 2016-09-07 | 中国电建集团西北勘测设计研究院有限公司 | 基于综合投影变形最小标准减小任意带高斯投影变形方法 |
CN105931199B (zh) * | 2016-04-14 | 2019-01-15 | 中国电建集团西北勘测设计研究院有限公司 | 基于综合投影变形最小标准减小任意带高斯投影变形方法 |
CN105869128B (zh) * | 2016-04-14 | 2018-10-30 | 中国电建集团西北勘测设计研究院有限公司 | 基于最大投影变形最小标准减小任意带高斯投影变形方法 |
CN105976314A (zh) * | 2016-05-10 | 2016-09-28 | 福州市勘测院 | 顾及相同中央投影经线的不同参考椭球投影平面坐标系统转换方法 |
CN105976314B (zh) * | 2016-05-10 | 2019-04-19 | 福州市勘测院 | 顾及相同中央投影经线的不同参考椭球投影平面坐标系统转换方法 |
CN106202000B (zh) * | 2016-07-22 | 2019-01-18 | 武汉大学 | 国家三维坐标系与任意地方平面坐标系间七参数转换方法 |
CN106202000A (zh) * | 2016-07-22 | 2016-12-07 | 武汉大学 | 国家三维坐标系与任意地方平面坐标系间七参数转换方法 |
CN106341936A (zh) * | 2016-08-31 | 2017-01-18 | 大连海事大学 | 一种内河智能航标位置监控方法 |
CN106683404A (zh) * | 2016-12-06 | 2017-05-17 | 华南理工大学 | 一种通过手机定位技术获取公交客流od的方法 |
CN106683404B (zh) * | 2016-12-06 | 2019-10-18 | 华南理工大学 | 一种通过手机定位技术获取公交客流od的方法 |
CN107908808A (zh) * | 2017-08-11 | 2018-04-13 | 山东交通学院 | 基于抵偿高程面或平均高程面的3°分带坐标转换系统 |
CN108731648A (zh) * | 2018-03-15 | 2018-11-02 | 广州市城市规划勘测设计研究院 | 2000独立坐标系参数获取方法、装置与计算机可读存储介质 |
CN108731648B (zh) * | 2018-03-15 | 2020-12-22 | 广州市城市规划勘测设计研究院 | 2000独立坐标系参数获取方法、装置与计算机可读存储介质 |
CN110427450A (zh) * | 2019-08-09 | 2019-11-08 | 深圳四维集思技术服务有限公司 | 基于互联网地图的油气管线数据处理方法和系统 |
CN114509051A (zh) * | 2022-01-26 | 2022-05-17 | 中交二公局第三工程有限公司 | 一种道路工程中无人机测量放样方法 |
CN114546941A (zh) * | 2022-04-28 | 2022-05-27 | 长江水利委员会水文局 | 一种数字地形图空间坐标系与投影的全信息无损转换方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103389080B (zh) | 2015-08-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103389080B (zh) | 一种基于地理信息应用城市独立坐标系参数获取方法 | |
CN104537606B (zh) | 地理坐标投影变化方法 | |
Janssen | Understanding coordinate reference systems, datums and transformations | |
CN105203104B (zh) | 一种适用于高精度惯导系统的重力场建模方法 | |
CN105160197A (zh) | 一种综合性地理空间数据坐标转换方法及系统 | |
CN102223709B (zh) | 基于cors系统实时获取特定基准的高斯平面直角坐标的方法 | |
Okiemute et al. | Practical local geoid model determination for mean sea level heights of surveys and stable building projects | |
CN102567627A (zh) | 基于卫星重力梯度观测数据的圆环面调和分析方法 | |
Soler et al. | On deflection of the vertical components and their transformations | |
CN101959121A (zh) | 一种地理位置服务标识编码技术 | |
Gaspar | Using empirical map projections for modeling early nautical charts | |
Yakubu et al. | Ramification of datum and ellipsoidal parameters on post processed differential global positioning system (DGPS) data–A case study | |
Ollikainen et al. | The Finnish coordinate reference systems | |
Čechurová et al. | Cartometric analysis of the Czechoslovak version of 1: 75 000 scale sheets of the Third Military Survey (1918–1956) | |
Moscovici et al. | Mathematical cartography based on georeferencing maps | |
Lu et al. | Reference ellipsoid and the geodetic coordinate system | |
Deakin | A guide to the mathematics of map projections | |
SOSSA et al. | CARTOMETRIC INVESTIGATION OF THE ACCURACY PLAN OF LVIV IN 1894 | |
Basil et al. | Determination of Deflection of the Vertical Components: Im-plications on Terrestrial Geodetic Measurement | |
Nikolli et al. | The coordinate reference systems of spatial data infrastructures in Albania | |
Lu et al. | Gauss and UTM Conformal Projections and the Plane Rectangular Coordinate System | |
de Lange | Geoobjects and Reference Systems | |
Nurçe | Development of Classical and Modern Geodetic Reference Systems of Albania | |
Sami | Sudan Reference System | |
CN106126682B (zh) | 产业地图区域地斑智能绘制方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: 400020 Jiangbei District, Chongqing electric measuring Village No. 231 Patentee after: Chongqing geographic information and Remote Sensing Application Center (Chongqing surveying and mapping product quality inspection and testing center) Address before: The 401121 northern New District of Chongqing Municipality Neptune technology building D block 4 Building Patentee before: Chongqing Geographical Information Center |
|
CP03 | Change of name, title or address |