CN110415165B - 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 - Google Patents
一种极区等间距格网点的平面坐标与大地坐标相互转换方法 Download PDFInfo
- Publication number
- CN110415165B CN110415165B CN201910658922.4A CN201910658922A CN110415165B CN 110415165 B CN110415165 B CN 110415165B CN 201910658922 A CN201910658922 A CN 201910658922A CN 110415165 B CN110415165 B CN 110415165B
- Authority
- CN
- China
- Prior art keywords
- coordinates
- plane
- polar region
- point
- geodetic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 230000005484 gravity Effects 0.000 description 4
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 1
Images
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),过该点的子午线与赤道的交点为Q,Q和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.257 222 101)
地心引力常数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.257 222 101)
地心引力常数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,过该点的子午线与赤道的交点为Q,Q和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 (1)
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),过该点的子午线与赤道的交点为Q,Q和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的极区平面直角坐标可按下式计算:
通过上式,可将极区平面直角坐标系中等间距各点的平面坐标转换为椭球面上的大地坐标,完成极区平面直角坐标与大地坐标的相互转换。
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 CN110415165A (zh) | 2019-11-05 |
CN110415165B true 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) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110910306B (zh) * | 2019-11-28 | 2022-06-21 | 深圳市信义科技有限公司 | 一种将全景地图中球面坐标转换为直线段的方法 |
CN111366130A (zh) * | 2020-03-03 | 2020-07-03 | 山东交通学院 | 一种高程异常内插推估方法 |
CN112037122A (zh) * | 2020-08-31 | 2020-12-04 | 四川易利数字城市科技有限公司 | 一种在游戏引擎中使用地理坐标的方法 |
CN113569197B (zh) * | 2021-07-28 | 2023-06-30 | 国家海洋信息中心 | 一种基于地球椭球面两点等距离点解算方法 |
CN113917917B (zh) * | 2021-09-24 | 2023-09-15 | 四川启睿克科技有限公司 | 室内仿生多足机器人避障方法、装置及计算机可读介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 | 中国人民解放军海军工程大学 | 用地球椭球模型改进的捷联惯导极地横向导航方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7872647B2 (en) * | 2002-07-19 | 2011-01-18 | Rockwell Collins Simulation And Training Solutions Llc | System and method for modeling a spheroid world database |
-
2019
- 2019-07-22 CN CN201910658922.4A patent/CN110415165B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
---|
基于统一横向坐标系的极区地球椭球模型导航方法;卞鸿巍;林秀秀;王荣颖;马恒;;中国惯性技术学报(第05期);全文 * |
基于贝叶斯正则化BP神经网络的GPS高程转换;宋雷;黄腾;方剑;周旭华;;西南交通大学学报(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN110415165A (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110415165B (zh) | 一种极区等间距格网点的平面坐标与大地坐标相互转换方法 | |
CN105203104B (zh) | 一种适用于高精度惯导系统的重力场建模方法 | |
CN110361683B (zh) | 基于双目标函数粒子群优化的磁力计校正方法 | |
CN106871932B (zh) | 基于金字塔搜索地形匹配的星载激光在轨指向检校方法 | |
WO2016095749A1 (zh) | 电子地图中查询空间对象和建立空间索引的方法及装置 | |
CN108426575B (zh) | 用地球椭球模型改进的捷联惯导极地横向导航方法 | |
CN103389080A (zh) | 一种基于地理信息应用城市独立坐标系参数获取方法 | |
CN111460071B (zh) | 高精地图的偏转方法、装置、设备及可读存储介质 | |
CN111380557A (zh) | 一种无人车全局路径规划方法及装置 | |
CN103335649A (zh) | 一种惯性导航系统极区导航参数解算方法 | |
CN110440753B (zh) | 考虑地球曲率的高精度dem航空重力远区地形改正方法 | |
CN105116470A (zh) | 基于地物tin模型的日照时数计算 | |
CN111060078A (zh) | 基于卫星观测角误差估计的定位方法 | |
CN114564031A (zh) | 一种基于旋翼无人机实现风机巡检拍摄的路径规划方法 | |
JP2018054561A (ja) | ジャイロコンパス測量方法 | |
CN111027137B (zh) | 基于遥测数据的航天器动力学模型高精度动态构建方法 | |
CN113959450B (zh) | 基于法向量位置模型的惯性导航阻尼方法 | |
CN115857036B (zh) | 一种新型球面均匀分布网格在球谐分析中的应用 | |
CN114777730B (zh) | 一种地基太阳高度角的计算方法和装置 | |
CN115683091B (zh) | 一种基于太阳偏振中性面分时反演重构的自主定位方法 | |
CN105486304B (zh) | 太阳同步轨道星敏感器遇月亮的预报方法 | |
CN115951589B (zh) | 基于最大化Kozachenko-Leonenko熵的恒星均匀选取方法 | |
CN114547847B (zh) | 基于平方公里格网系统的遥感卫星数据源解算方法 | |
CN117556845A (zh) | 一种东北半球地理实体二维位置码快速编码方法 | |
Yan et al. | One kind of networking mode ATC system coordinate projected optimization method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |