CN110415165A - 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 - Google Patents
一种极区等间距格网点的平面坐标与大地坐标相互转换方法 Download PDFInfo
- Publication number
- CN110415165A CN110415165A CN201910658922.4A CN201910658922A CN110415165A CN 110415165 A CN110415165 A CN 110415165A CN 201910658922 A CN201910658922 A CN 201910658922A CN 110415165 A CN110415165 A CN 110415165A
- Authority
- CN
- China
- Prior art keywords
- coordinates
- polar region
- geodetic
- plane
- point
- 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
- 238000000034 method Methods 0.000 title claims abstract description 14
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 230000005484 gravity Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/06—Topological mapping of higher dimensional structures onto lower dimensional surfaces
- G06T3/067—Reshaping or unfolding 3D tree structures onto 2D planes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/08—Projecting images onto non-planar surfaces, e.g. geodetic screens
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种极区等间距格网点的平面坐标与大地坐标相互转换方法,首先建立极区平面直角坐标系,将极区椭球面采用正轴等角方位投影方法进行投影,在极区平面直角坐标系中建立等间距格网,将极区等间距格网点的平面直角坐标转换大地坐标,再将椭球面上的大地坐标转换为极区平面直角坐标。完成极区平面直角坐标与大地坐标得相互转换。本发明的有益效果是能够高效的将极区等间距格网点的平面坐标与大地坐标进行相互转换。
Description
技术领域
本发明属于地理技术领域,涉及一种极区等间距格网点的平面坐标与大地坐标相互转换的方法。
背景技术
地球重力场模型通常展开成球谐或椭球谐函数的级数,重力场参数如重力异常、大地水准面起伏和垂线偏差等都可以经纬度为因数利用地球重力场模型进行计算。但是在地球两极区域,由于极区具有纬度高、经线收敛快、存在极点等特点,相同经纬度间距格网对应的实际距离有较大差异,在离散点格网插值计算重力场参数过程中精度不均匀。
发明内容
本发明的目的在于提供一种极区等间距格网点的平面坐标转换和大地坐标相互转换的方法,本发明的有益效果是能够高效的将极区等间距格网点的平面坐标转换和为大地坐标进行相互转换。
本发明所采用的技术方案是首先建立极区平面直角坐标系,以北极区域为例,将极区椭球面采用正轴等角方位投影方法进行投影,如图1,在投影面上以极点为坐标原点O,本初子午线方向作为格网纵方向y轴,垂直于y轴方向为x轴建立极区平面直角坐标系,投影后极点为坐标原点O,经线由极点向外成放射状直线,经线间的夹角等于经度差,地理坐标以大地经度与大地纬度表示,通过地面M点和地球椭球体旋转轴的平面,称M点的大地子午面,M点的大地子午面与起始大地子午面本初子午面间的夹角L称为大地经度,通过M点的地球椭球体的法线与赤道平面的夹角B,称为大地纬度,大地经度与大地纬度合称大地经纬度。
极区平面直角坐标转换大地坐标:对于极区平面直角坐标任意点H,其平面直角坐标为(x,y),该点对应椭球面上的大地坐标为(B,L),大地经度按以下公式进行转换:
式中,sgn(x)为符号函数,当x>0时,sgn(x)=1;当x<0时,sgn(x)=-1。
纬度线在投影面上的投影为圆线,投影面上点的纬度与该点至极点的距离相关,设有子午线上两点E和H,E在赤道上,H的纬度为B,以2000国家大地坐标系采用的地球椭球参数计算,当点H位于极点时,E和H间的子午线弧长Ym的值为10001965.729m,在极区平面直角坐标中,设任意点H至极点O的距离为S,则点E和H间的子午线弧长为(Ym-S),点H的纬度可根据子午线弧长(Ym-S)进行计算,令X=(Ym-S),纬度采用迭代方法计算,计算方式为:
直到Bi-1-Bi,小于某一个指定数值,即可停止迭代,式中:
这样,可将极区平面直角坐标系中等间距各点的平面坐标转换为椭球面上的大地坐标。
极区大地坐标转换平面直角坐标:设椭球面上任意点N的大地坐标为(L,B),过该点的子午线与赤道的交点为E,E和N间的子午线弧长X的计算公式为:
X=a(1-e2)[A′(B/ρ)-B′sin2B+C′sin4B-D′sin6B
+E′sin8B-F′sin10B+G′sin12B] (6)
式中各项按下式计算:
则,点A至极点的距离为(Ym-X),点A的极区平面直角坐标可按下式计算:
通过上式,可将极区平面直角坐标系中等间距各点的平面坐标转换为椭球面上的大地坐标。完成极区平面直角坐标与大地坐标得相互转换。
进一步,地球椭球参数选用2000国家大地坐标系采用的地球椭球参数:
长半轴a=6378137m
地球重力场二阶帯谐系数:J2=1.082629832258×10-3
(对应扁率:1/298.257222101)
地心引力常数GM=3.986004418×1014m3/s2
地球自转角速度ω=7.292115×10-5rad/s
附图说明
图1是极区平面直角坐标系
图2极区平面直角坐标转换大地坐标
图3大地坐标转换极区平面直角坐标
具体实施方式
下面结合具体实施方式对本发明进行详细说明。
在投影计算过程中,本发明选用2000国家大地坐标系(China GeodeticCoordinate System 2000,简称CGCS2000)采用的地球椭球参数:
长半轴a=6378137m
地球重力场二阶帯谐系数:J2=1.082629832258×10-3
(对应扁率:1/298.257222101)
地心引力常数GM=3.986004418×1014m3/s2
地球自转角速度ω=7.292115×10-5rad/s
极区平面直角坐标系:
首先建立极区平面直角坐标系,以北极区域为例,将极区椭球面采用正轴等角方位投影方法进行投影,如图1,在投影面上以极点为坐标原点O,本初子午线方向作为格网纵方向y轴,垂直于y轴方向为x轴建立极区平面直角坐标系,投影后极点为坐标原点O,图1中圆为纬线,经线由极点向外成放射状直线,经线间的夹角等于经度差。
地理坐标以大地经度与大地纬度表示,通过地面M点和地球椭球体旋转轴的平面,称M点的大地子午面,M点的大地子午面与起始大地子午面本初子午面间的夹角L称为大地经度,通过M点的地球椭球体的法线与赤道平面的夹角B,称为大地纬度,大地经度与大地纬度合称大地经纬度。
极区平面直角坐标转换大地坐标:
对于极区平面直角坐标任意点H,如图2,其平面直角坐标为(x,y),该点对应椭球面上的大地坐标为(B,L),大地经度按以下公式进行转换:
式中,sgn(x)为符号函数,当x>0时,sgn(x)=1;当x<0时,sgn(x)=-1。
纬度线在投影面上的投影为圆线,投影面上点的纬度与该点至极点的距离相关,设有子午线上两点E和H,E在赤道上,H的纬度为B,以2000国家大地坐标系采用的地球椭球参数计算,当点H位于极点时,E和H间的子午线弧长Ym的值为10001965.729m,在极区平面直角坐标中,设任意点H至极点O的距离为S,则点E和H间的子午线弧长为(Ym-S),点H的纬度可根据子午线弧长(Ym-S)进行计算,令X=(Ym-S),纬度采用迭代方法计算,计算方式为:
直到Bi-1-Bi,小于某一个指定数值,即可停止迭代,式中:
这样,可将极区平面直角坐标系中等间距各点的平面坐标转换为椭球面上的大地坐标。
在极区平面直角坐标中建立等间距格网,格网间距为5.0km,x与y轴坐标的取值范围都为-1000km~1000km之间,格网为401行列,共计160801个格网结点。按以上算法编制程序,将所有格网结点的极区平面直角坐标转换为大地坐标。
极区大地坐标转换平面直角坐标
设椭球面上任意点N的大地坐标为(L,B),如图3,过该点的子午线与赤道的交点为E,E和N间的子午线弧长X的计算公式为:
X=a(1-e2)[A′(B/ρ)-B′sin2B+C′sin4B-D′sin6B
+E′sin8B-F′sin10B+G′sin12B] (6)
式中各项按下式计算:
则,点A至极点的距离为(Ym-X),点A的极区平面直角坐标可按下式计算:
通过上式,可将极区平面直角坐标系中等间距各点的平面坐标转换为椭球面上的大地坐标。完成极区平面直角坐标与大地坐标得相互转换。
按以上算法编制程序,将所有格网结点的极区大地坐标转换为平面直角坐标,所有的转换结果与原始格网结点坐标进行比较,通过两次转换,转换点坐标差均小于0.001m,证明了转换算法的正确性。
以上所述仅是对本发明的较佳实施方式而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施方式所做的任何简单修改,等同变化与修饰,均属于本发明技术方案的范围内。
Claims (2)
1.一种极区等间距格网点的平面坐标与大地坐标相互转换方法,其特征在于:首先建立极区平面直角坐标系,将极区椭球面采用正轴等角方位投影方法进行投影,在投影面上以极点为坐标原点O,本初子午线方向作为格网纵方向y轴,垂直于y轴方向为x轴建立极区平面直角坐标系,投影后极点为坐标原点O,经线由极点向外成放射状直线,经线间的夹角等于经度差,地理坐标以大地经度与大地纬度表示,通过地面M点和地球椭球体旋转轴的平面,称M点的大地子午面,M点的大地子午面与起始大地子午面本初子午面间的夹角L称为大地经度,通过M点的地球椭球体的法线与赤道平面的夹角B,称为大地纬度,大地经度与大地纬度合称大地经纬度;
极区平面直角坐标转换大地坐标:对于极区平面直角坐标任意点H,其平面直角坐标为(x,y),该点对应椭球面上的大地坐标为(B,L),大地经度按以下公式进行转换:
式中,sgn(x)为符号函数,当x>0时,sgn(x)=1;当x<0时,sgn(x)=-1;
纬度线在投影面上的投影为圆线,投影面上点的纬度与该点至极点的距离相关,设有子午线上两点E和H,E在赤道上,H的纬度为B,以2000国家大地坐标系采用的地球椭球参数计算,当点H位于极点时,E和H间的子午线弧长Ym的值为10001965.729m,在极区平面直角坐标中,设任意点H至极点O的距离为S,则点E和H间的子午线弧长为(Ym-S),点H的纬度可根据子午线弧长(Ym-S)进行计算,令X=(Ym-S),纬度采用迭代方法计算,计算方式为:
直到Bi-1-Bi,小于某一个指定数值,即可停止迭代,式中:
这样,可将极区平面直角坐标系中等间距各点的平面坐标转换为椭球面上的大地坐标;
在极区平面直角坐标中建立等间距格网,格网间距为5.0km,x与y轴坐标的取值范围都为-1000km~1000km之间,格网为401行列,共计160801个格网结点;按以上算法编制程序,将所有格网结点的极区平面直角坐标转换为大地坐标;
极区大地坐标转换平面直角坐标:设椭球面上任意点N的大地坐标为(L,B),过该点的子午线与赤道的交点为E,E和N间的子午线弧长X的计算公式为:
X=a(1-e2)[A′(B/ρ)-B′sin2B+C′sin4B-D′sin6B
+E′sin8B-F′sin10B+G′sin12B] (6)
式中各项按下式计算:
则,点A至极点的距离为(Ym-X),点A的极区平面直角坐标可按下式计算:
通过上式,可将极区平面直角坐标系中等间距各点的平面坐标转换为椭球面上的大地坐标,完成极区平面直角坐标与大地坐标的相互转换。
2.按照权利要求1所述一种极区等间距格网点的平面坐标与大地坐标相互转换方法,其特征在于:所述地球椭球参数选用2000国家大地坐标系采用的地球椭球参数:
长半轴a=6378137m
地球重力场二阶帯谐系数:J2=1.082629832258×10-3
(对应扁率:1/298.257222101)
地心引力常数GM=3.986004418×1014m3/s2
地球自转角速度ω=7.292115×10-5rad/s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910658922.4A CN110415165B (zh) | 2019-07-22 | 2019-07-22 | 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910658922.4A CN110415165B (zh) | 2019-07-22 | 2019-07-22 | 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110415165A true CN110415165A (zh) | 2019-11-05 |
CN110415165B CN110415165B (zh) | 2023-06-13 |
Family
ID=68362281
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910658922.4A Active CN110415165B (zh) | 2019-07-22 | 2019-07-22 | 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110415165B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110910306A (zh) * | 2019-11-28 | 2020-03-24 | 深圳市信义科技有限公司 | 一种将全景地图中球面坐标转换为直线段的方法 |
CN111366130A (zh) * | 2020-03-03 | 2020-07-03 | 山东交通学院 | 一种高程异常内插推估方法 |
CN112037122A (zh) * | 2020-08-31 | 2020-12-04 | 四川易利数字城市科技有限公司 | 一种在游戏引擎中使用地理坐标的方法 |
CN113569197A (zh) * | 2021-07-28 | 2021-10-29 | 国家海洋信息中心 | 一种基于地球椭球面两点等距离点解算方法 |
CN113917917A (zh) * | 2021-09-24 | 2022-01-11 | 四川启睿克科技有限公司 | 室内仿生多足机器人避障方法、装置及计算机可读介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040098236A1 (en) * | 2002-07-19 | 2004-05-20 | Evans & Sutherland Computer Corporation | System and method for modeling a spheroid world database |
CN104215242A (zh) * | 2014-09-30 | 2014-12-17 | 东南大学 | 一种基于横向游移坐标系的极区惯性导航方法 |
CN105930306A (zh) * | 2016-04-14 | 2016-09-07 | 中国电建集团西北勘测设计研究院有限公司 | 一种建立工程区域水准椭球的方法 |
CN107908808A (zh) * | 2017-08-11 | 2018-04-13 | 山东交通学院 | 基于抵偿高程面或平均高程面的3°分带坐标转换*** |
CN108426575A (zh) * | 2018-02-12 | 2018-08-21 | 中国人民解放军海军工程大学 | 用地球椭球模型改进的捷联惯导极地横向导航方法 |
-
2019
- 2019-07-22 CN CN201910658922.4A patent/CN110415165B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040098236A1 (en) * | 2002-07-19 | 2004-05-20 | Evans & Sutherland Computer Corporation | System and method for modeling a spheroid world database |
CN104215242A (zh) * | 2014-09-30 | 2014-12-17 | 东南大学 | 一种基于横向游移坐标系的极区惯性导航方法 |
CN105930306A (zh) * | 2016-04-14 | 2016-09-07 | 中国电建集团西北勘测设计研究院有限公司 | 一种建立工程区域水准椭球的方法 |
CN107908808A (zh) * | 2017-08-11 | 2018-04-13 | 山东交通学院 | 基于抵偿高程面或平均高程面的3°分带坐标转换*** |
CN108426575A (zh) * | 2018-02-12 | 2018-08-21 | 中国人民解放军海军工程大学 | 用地球椭球模型改进的捷联惯导极地横向导航方法 |
Non-Patent Citations (2)
Title |
---|
卞鸿巍;林秀秀;王荣颖;马恒;: "基于统一横向坐标系的极区地球椭球模型导航方法", 中国惯性技术学报 * |
宋雷;黄腾;方剑;周旭华;: "基于贝叶斯正则化BP神经网络的GPS高程转换", 西南交通大学学报 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110910306A (zh) * | 2019-11-28 | 2020-03-24 | 深圳市信义科技有限公司 | 一种将全景地图中球面坐标转换为直线段的方法 |
CN110910306B (zh) * | 2019-11-28 | 2022-06-21 | 深圳市信义科技有限公司 | 一种将全景地图中球面坐标转换为直线段的方法 |
CN111366130A (zh) * | 2020-03-03 | 2020-07-03 | 山东交通学院 | 一种高程异常内插推估方法 |
CN112037122A (zh) * | 2020-08-31 | 2020-12-04 | 四川易利数字城市科技有限公司 | 一种在游戏引擎中使用地理坐标的方法 |
CN113569197A (zh) * | 2021-07-28 | 2021-10-29 | 国家海洋信息中心 | 一种基于地球椭球面两点等距离点解算方法 |
CN113569197B (zh) * | 2021-07-28 | 2023-06-30 | 国家海洋信息中心 | 一种基于地球椭球面两点等距离点解算方法 |
CN113917917A (zh) * | 2021-09-24 | 2022-01-11 | 四川启睿克科技有限公司 | 室内仿生多足机器人避障方法、装置及计算机可读介质 |
CN113917917B (zh) * | 2021-09-24 | 2023-09-15 | 四川启睿克科技有限公司 | 室内仿生多足机器人避障方法、装置及计算机可读介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110415165B (zh) | 2023-06-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110415165A (zh) | 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 | |
CN104048659B (zh) | 地图坐标系的转换方法和*** | |
CN103810286B (zh) | 一种用于二维地图与三维地图匹配的坐标点定位方法 | |
CN103389080B (zh) | 一种基于地理信息应用城市独立坐标系参数获取方法 | |
CN111161179A (zh) | 一种基于法向量的点云平滑光顺滤波方法 | |
CN105160197A (zh) | 一种综合性地理空间数据坐标转换方法及*** | |
CN105203104A (zh) | 一种适用于高精度惯导***的重力场建模方法 | |
CN102393183A (zh) | 基于控制网的海量点云快速配准方法 | |
CN104820984B (zh) | 一种卫星遥感三线阵立体影像处理***及方法 | |
CN112461224B (zh) | 一种基于已知姿态角的磁力计标定方法 | |
CN107833282A (zh) | 一种地形建模和网格生成方法及装置 | |
CN109872350A (zh) | 一种新的点云自动配准方法 | |
CN109671150A (zh) | 基于数字地球的机场土方计算方法 | |
CN107908808A (zh) | 基于抵偿高程面或平均高程面的3°分带坐标转换*** | |
CN105180894B (zh) | 一种工测数据与三维数字地球高程拟合的方法 | |
CN108733711B (zh) | 基于三维gis技术的配电线路空间距离获取方法 | |
CN106767811B (zh) | 大椭圆轨道姿态基准确定方法 | |
CN108981734A (zh) | 电子地图道路拓展方法、装置、电子设备及存储介质 | |
CN105004351A (zh) | 基于自适应upf的sins大方位失准角初始对准方法 | |
CN106885585A (zh) | 一种基于光束法平差的星载摄影测量***一体化检校方法 | |
CN105930306B (zh) | 一种建立工程区域水准椭球的方法 | |
CN108508463A (zh) | 基于Fourier-Hermite正交多项式扩展椭球集员滤波方法 | |
CN107462220A (zh) | 面向月基对地观测影像的投影极坐标几何表达方法 | |
CN107067454A (zh) | 一种基于超图组件开发平台的雷达威力图三维显示方法 | |
CN104050869B (zh) | Gis地图覆盖物布点响应方法和*** |
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 |