怎么把空间直角坐标换算成大地经纬度_百度作业帮
怎么把空间直角坐标换算成大地经纬度
怎么把空间直角坐标换算成大地经纬度
是投影平面直角坐标系吧?首先你要有转换参数,一般网上都能下到相应的转换软件,比如54-80-WGS84等,输入相应数值就可得出.
这个需专业的软件，如果用手算公式太多，太麻烦。有个导线测量平差的测量用CausCoor.exe计算80坐标， 可以利用一些公开的GPS解算软件里面有自带
别手算了，找软件，百度，或在cass软件中有转换功能的，有很多种，当然不在一个坐标系中的转换是需要求取转换参数的
你知道纬度的推算公式吗
适合编程的公式后使用快捷导航没有帐号？
扫一扫，访问微社区
只需一步, 快速开始
只需一步，快速开始
查看: 2915|回复: 11
求教ncl中把经纬度坐标转换成距离坐标（单位是km）？
在线时间 小时
登录后查看更多精彩内容~
才可以下载或查看，没有帐号？
最近在用ncl画wrfout数据的水平和垂直剖面图，默认的坐标为经纬度坐标，现在我想把坐标转换成距离坐标（km）。请问该怎么操作呢？麻烦大家给点指导。非常感谢l
给自己设置一个个性的签名吧~
在线时间 小时
WRF网格一般是投影面上的等面积网格，要是换距离也是投影面上的，不如经纬度坐标吧。
给自己设置一个个性的签名吧~
在线时间 小时
坐等答案，我要需要
给自己设置一个个性的签名吧~
在线时间 小时
同求！！！
给自己设置一个个性的签名吧~
在线时间 小时
同求！！！
给自己设置一个个性的签名吧~
在线时间 小时
Please try to generate plots by RIP
给自己设置一个个性的签名吧~
在线时间 小时
好好学习，天天向上
给自己设置一个个性的签名吧~
在线时间 小时
量值转换一下
给自己设置一个个性的签名吧~
在线时间 小时
dy& && &= f@DY/1000.& && && && && &&&; dx (km)
&&west_east = ispan(0,mx-1,1)*dy& && & ; calculate x values
&&west_east@long_name = &west_east&
&&west_east@units = &km&
&&w&west_east& &&&= west_east& && && & ; associate &x& values with w
给自己设置一个个性的签名吧~
在线时间 小时
支持一下子
给自己设置一个个性的签名吧~1258人阅读
转自:http://blog.csdn.net/lijie45655/article/details/6771524
手持GPS，正在学习当中，晚上我在谷歌地球上找到端午节要去的目的地，找出了它的经纬度坐标，然后把数值输入GPS，结果一看地图，纳闷了，相差了好几十公里。我又以重庆长江大桥北桥头（石板坡）作为参照物，纬度是N29°33′09.79″，经度是E106°33′29.63″，然后我直接在GPS上输入了N29.33097°，E106.33296°，结果一看航点位置，居然跑到江津去了，吓我一大跳，这个差别也太大了撒。赶紧看看是怎么回事，结果是我自己不懂，以为经纬度标识方式就一种，结果差之毫厘谬以千里，下面就是正确的经纬度格式及转换方式，贴出来供我等菜鸟学习学习，哈哈哈哈，又学到了一点知识。
经纬度格式分为三种：度、度-分、度-份-秒,
R% O$ C- h: ^! Z5 b9 f& u6 U/ z
1.） ddd.ddddd °【度 . 度 格式】的十进制小数部分（5位）
2.） ddd°mm.mmm’ 【度 . 分 . 分 格式】的十进制小数部分（3位）
& q' Z8 a1 d- M( \) L. u( Y1 E9 Y1 t: R
3.)& &ddd°mm’ss’’ 【度 . 分 . 秒 格式】
9 Q$ r1 y% v&&U3 v/ F
Google 使用的是第三种格式&&度。分’秒’’
度分转换：&$
u* X! }! u% I$ P* N0 y
将度分单位数据转换为度单位数据&
度=度+分/60&&
F3 Y4 d/ ^. ^/ A8 _( J/ K
W; n3 ~; t1 v7 r
经度 = 116°20.12’&#
\0 p9 P% y2 g( {% A. F$ N& e: }. d1 Q
纬度 = 39°12.34’&
经度 = 116 + 20.12 / 60 = 116.33533°&
纬度 = 39 + 12.34 / 60 = 39.20567°&
D9 ]- j5 W# \. O& S
& A! {: |3 I$ y2 T8 ?
度分秒转换：&+
O/ F4 U0 F/ f( z
将度分秒单位数据转换为度单位数据&-
w4 B2 _8 T3 F
度 = 度 + 分 / 60 + 秒 / 60 / 60&:
H$ d$ [) @9 i3 b- b6 E) @; C& E
经度 = 116°20’43”&
纬度 = 39°12’37”&
经度 = 116 + 20 / 60 + 43 / 60 / 60 = 116.34528°&,
E* o&&O/ k9 [0 h5 w# B+ @
纬度 = 39 + 12 / 60 + 37 / 60 / 60 = 39.21028°&
0 x4 p7 W4 h$ g
其格式的经纬值先转换度，再采用上面的运算。,
经过从新计算，石板坡长江大桥北桥头换算过来应该是这样表达：7
m4 q* V1 g7 ]4 U* Z&&I
N29.552719°，E106.55823°
把这组数据输入GPS，哈哈哈哈哈，不偏不倚，正好在北桥头4
i% m* ]: R6 ]3 {& p2 l# o
回头把端午节的目的地的坐标重新计算了下输入GPS，嘿嘿，也正确！
这下开心啦
一直在做遥感图像处理的东西，多是一些数字图像处理的分类算法，忘了自己竟是学地学这块的，现在连基本的经纬线都忘了，真是惭愧！现补充如下，以备后用。
关于经纬度的相关定义：
&&&&&&& 经线也称子午线，是人类为度量方便而假设出来的辅助线，定义为地球表面连接南北两极的大圆线上的半圆弧。任两根经线的长度相等，相交于南北两极点。每一根经线都有其相对应的数值，称为经度。经线指示南北方向。
本初子午线 ：
&&&&&& 本初子午线又称“首子午线”或“零子午线”，是计算东西经度的起点，以经过英国伦敦东南格林尼治的经线为本初子午线，作为计算地理的起点和世界标准“时区”的起点。本初子午线以东为东经，以西为西经，全球经度测量均以本初子午线与赤道的交点E点作为经度原点。
&子午线命名的由来：
&&&&& 某一天体视运动轨迹中，同一子午线上的各点该天体在上中天（午）与下中天（子）出现的时刻相同。
&&& 纬线是人类为度量方便而假设出来的辅助线，定义为地球表面某点随地球自转所形成的轨迹。纬线平行于赤道，经线垂直于赤道，任何一根纬线都是圆形而且两两平行。纬线的长度是赤道的周长乘以纬线的纬度的余弦，所以赤道最长，离赤道越远的纬线，周长越短，到了两极就缩为0。纬线指示东西方向。
北回归线（23°26' 22& N）
&&&&&&太阳在北半球能够直射到的离赤道最远的位置，其纬度值为黄赤交角，是一条纬线，大约在北纬23度26分的地方。
赤道（0°N）
&&&& 通过地球中心划一个与地轴成直角相交的平面，在地球表面相应出现一个和地球的极距离相等的假想圆圈。赤道的纬度是0°。是地球表面的点随地球自转产生的轨迹中周长最长的圆周线，赤道半径 Km ；两极半径 Km；平均半径 Km ；赤道周长 40075.7Km。如果把地球看做一个绝对的球体的话，赤道距离南北两极相等，是一个大圆。它把地球分为南北两半球，其以北是北半球，以南是南半球，是划分纬度的基线，是南北纬线的起点。
南回归线（23° 26' 22& S）
&&&&&& 太对应于北回归线，阳在地球上的直射点在一年内到达的最南点所在的纬线。
　&&& 是地球上一个地点离一根被称为本初子午线的南北方向走线以东或以西的度数。它是通过某地的经线面与本初子午面所成的二面角，在本初子午线以东的经度叫东经，在本初子午线以西的叫西经。东经用“E”表示，西经用“W”表示。比如：E116.33533°,W116.33533°.经度的每一度被分为60角分，每一分被分为60秒。也可以用小数点位表示。有时西经也被写成负数,如：-116.33533°。
&&&&& 是指某点与地球球心的连线和地球赤道面所成的线面角，其数值在0至90度之间。位于赤道以北的点的纬度叫北纬，记为N；位于赤道以南的点的纬度称南纬，记为S。
纬度数值在0至30度之间的地区称为低纬度地区；纬度数值在30至60度之间的地区称为中纬度地区；纬度数值在60至90度之间的地区称为高纬度地区。赤道是0°纬线，北纬度的最大值为90°，即北极点；南纬度的最大值为90°，即南极点。
&&&& 垂直于某地参考椭球体表面的线。
经纬度坐标系统：
&&&&&&& 是经度与纬度的合称组成一个坐标系统。又称为地理坐标系统，它是一种利用三度空间的球面来定义地球上的空间的球面坐标系统，能够标示地球上的任何一个位置。
一个经度和一个纬度一起确定地球上一个地点的精确位置。一个完整的图面地理坐标系统由一个完整的参考椭球体系和一个完整的投影体系组成。
天文地理坐标系统和大地坐标系统
&&&&&& 由经线和纬线构成两组互相正交的曲线坐标网叫地理坐标网。由经纬度构成的地理坐标系统又叫地理坐标系。地理坐标分为天文地理坐标和大地地理坐标。天文地理坐标是用天文测量方法确定的，大地地理坐标是用大地测量方法确定的。我们在地球椭球面上所用的地理坐标系属于大地地理坐标系，简称大地坐标系。
&&&&&& 大地坐标系是以以参考椭球面为基准面的坐标，地面点P的位置用大地经度L、大地纬度B和大地高H表示。当点在参考椭球面上时，仅用大地经度和大地纬度表示。虽然如此，但是，大地坐标系与经纬度网并不一致，大地坐标系以球面的形式描述了真实的地球坐标，而经纬度网则是在制图时，为满足平面作图的需要（即地图投影），参照一定的投影公式制作成平面地图，并在图上套上经纬线。常用作世界地图的投影有墨卡托投影、高尔投影、摩尔威特投影、等差分纬线多圆锥投影、格灵顿投影、桑森投影、乌尔马耶夫投影等。而横轴墨卡托投影即为高斯-克吕格投影，一般用于我国的地理制图。
&我国的高斯-克吕格投影体系
&&&&&& 它是一种横轴、椭圆柱面、等角投影。其投影过程可简述如下：椭圆柱面与地球椭球在某一子午圈上相切，这条子午圈叫做投影的中央子午线，又称轴子午线，它也是高斯投影后的平面直角坐标系的纵轴（一般定义为x轴）；地球的赤道面与椭圆柱面相交成一条直线，这条直线与中央子午线正交，它是平面直角坐标系的横轴（y轴）；把椭圆柱面展开，就得出以（x，y）为坐标的平面直角坐标系。按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差，又要使带数不致过多以减少换带计算工作，据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。通常按经差6度或3度分为六度带或三度带。六度带自0度子午线起每隔经差6度自西向东分带，带号依次编为第1、2…60带。三度带是在六度带的基础上分成的，它的中央子午线与六度带的中央子午线和分带子午线重合，即自1.5度子午线起每隔经差3度自西向东分带，带号依次编为三度带第
1、2…120带。我国的经度范围西起73°东至135°，可分成六度带十一个，各带中央经线依次为75°、81°、87°、……、117°、123°、129°、135°，或三度带二十二个。
经纬度格式
经纬度格式分为三种：度、度-分、度-份-秒
1.&& ddd.ddddd °【度 . 度 格式】的十进制小数部分（5位）
2.&& ddd°mm.mmm’ 【度 . 分 . 分 格式】的十进制小数部分（3位）
3.& ddd°mm’ss’’ 【度 . 分 . 秒 格式】
度分转换：&
将度分单位数据转换为度单位数据&
度=度+分/60&
经度 = 116°20.12’&
纬度 = 39°12.34’&
经度 = 116 + 20.12 / 60 = 116.33533°&
纬度 = 39 + 12.34 / 60 = 39.20567°
度分秒转换：&
将度分秒单位数据转换为度单位数据&
度 = 度 + 分 / 60 + 秒 / 60 / 60&
经度 = 116°20’43”&
纬度 = 39°12’37”&
经度 = 116 + 20 / 60 + 43 / 60 / 60 = 116.34528°&
纬度 = 39 + 12 / 60 + 37 / 60 / 60 = 39.21028°
在实际的处理中，计算相应像元点的地理坐标：
GDAL数据集有两种模式描述栅格位置（用点/线坐标系）以及地理参考坐标系之间的关系：首要的也是最普遍的是使用仿射转换，另一种则是GCPs（多控制点定位方式）
1用GDAL获取四个角点的地理坐标
GDALDataset::：：GetGeoTransform()
声明一个存储变量：
double&&&&&&& geoTransform[6];
geoTransform[0]是左上角像元的东坐标；
geoTransform[3]是左上角像元的北坐标；
geoTransform[1]是影像宽度上的分辨率；
geoTransform[5]是影像高度上的分辨率；
geoTransform[2]是旋转, 0表示上面为北方；
geoTransform[4]是旋转, 0表示上面为北方；
2相应的放射变换公式：
Xp = geoTransform [0] +Xpixel*geoTransform [1]+Yline*geoTransform [2];
Yp = geoTransform [3] +&Xpixel*geoTransform [4] +&YlineL*geoTransform [5];
其中Xpixel、YlineL分别地物表示单纯地影像上图像遍历时（x,y）坐标表示，点/线坐标系是从左上角(0,0)点到右下角，也就是坐标轴从左到右增长，从上到下增长的坐标系。
2 GCPs控制点坐标变换：
int nGCPs=poDataset-&GetGCPCount();&&&//获得控制点数目
const GDAL_GCP* pGCPs = poDataset-&GetGCPs(); //获得GCP控制点，是一个字段
double& *geoTransform& ;//与前面的仿射变换的相似
GDALGCPsToGeoTransform( nGCPs, pGCPs, geoTransform, TRUE ); //由GCPs获得仿射变换参数
int x,y;//任意点影像行列数
double xx, //任意点影像对应的地理坐标
double ex[nGCPs];//不同控制点算出的地理坐标X
double ey[nGCPs];//不同控制点算出的地理坐标Y
double pRadio[nGCPs]; //任意点相对与控制点的权值
double total=0; //权值和
double sx=0; //pRadio[i]*ex[i]
double sy=0; //pRadio[i]*ey[i]
for(int num=0;num&nGCPs;num++)
&&&&&&&&double gX,gY,tX,tY;
&&&&&&&&tX=pGCPs[num].dfGCPP //GCP位置X
&&&&&&&&tY=pGCPs[num].dfGCPL&&//GCP位置Y
&&&&&&&&gX=pGCPs[num].dfGCPX;&&&&&//地理位置X
&&&&&&&&gY=pGCPs[num].dfGCPY;&&&&&//地理位置Y
&&&&&&&&pRadio[num]=1/(double)(sqrt((tX-x)*(tX-x)+(tY-y)*(tY-y)));
&&&&&&&&total+=pRadio[num];
&&&&&&&&ex[num]=gX+(x-tX)*geoTransform[1];
&&&&&&&&ey[num]=gY+(y-tY)*geoTransform[5];
&&&&&&&&sx+=pRadio[num]*ex[num];
&&&&&&&&sy+=pRadio[num]*ey[num];
对于遥感图像而言，如果想显示在屏幕上，由于图像过大或者屏幕的关系，常常做成金字塔影像，或者拉伸图像数据，使图像产生变形以此来满足显示的需求。
假如已经将经纬度转换成小数点位，即以度为单位，已知道屏幕的高(y)和宽(h),地理坐标区域的范围(maxLon,minLon,maxLat,minLat)..这里我们知道了这些已知的参数
我们可以算出每像素所代表的经度和纬度(有人称这个为比例因子):
&&&&&&&&公式:scaleX = h/((maxLon-minLon)*3600)& ----------X轴上每像素代表的经度秒数;
&&&&&&& 公式:scaleY = y/((maxLat-minLat)*3600)&&-----------Y轴上每像素代表的纬度秒数;
算出该地理坐标区域中的任何一点(lon,lat)在屏幕上的坐标
公式:screenX = lon*3600/scaleX;& ---------屏幕坐标X轴坐标
公式:screenY = lat*3600/scaleY; ----------屏幕坐标Y轴坐标,&&
lon和lat为任意地点的地理坐标，前面的仿射变换已经算出
假如需要占满整个屏幕：
&公式:minX = minLon*3600/scaleX;&&& 区域左边置最左端
公式:minY = minLat*3600/scaleY;&&&& 区域上面置最上端
当地地理范围区域占满整个屏幕时,我们需要用到第三步计算出来的&screenX和screenY两个参数
该区域中的任何一点的公式如下:&&&
&公式:X&= screenX - minX = (lon - minLon)*3600/scaleX;&&&&&
由于纬度的方向和屏幕Y轴是相反的,&
&公式:screenMaxLat = (maxLat - minLat)*3600/scaleY;
&公式:screenLat = (lat - minLat)*3600/scaleY;
公式:Y = screenMaxLat - screenLat = (maxLat - lat)*3600/scaleY;&
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：222753次
积分：3831
积分：3831
排名：第3922名
原创：158篇
转载：104篇
评论：46条
(1)(1)(1)(3)(1)(7)(12)(1)(2)(1)(1)(4)(2)(1)(3)(2)(4)(1)(1)(1)(6)(2)(1)(1)(1)(1)(1)(2)(4)(2)(10)(1)(1)(1)(21)(1)(3)(2)(1)(3)(4)(5)(2)(4)(9)(3)(8)(2)(3)(3)(1)(2)(2)(6)(11)(34)(15)(8)(2)(8)(1)(16)GPS坐标转换经纬度及换算方法
GPS坐标和经纬度的算法和概率不太一样，但是我们可能会将他们互通起来用，下面先贴上我做的转换工具：。里面实现了gps到谷歌地图百度地图经纬度的转换。不含糊，下面将他们之间的联系。
GPS坐标系我本身不太了解它跟谷歌地图经纬度有多大区别，于是搜了一下，看看他们的区别：
地形图坐标系：我国的地形图采用高斯－克吕格平面直角坐标系。在该坐标系中，横轴：赤道，用Ｙ表示；纵轴：中央经线，用Ｘ表示；坐标原点：中央经线与赤道的交点，用0表示。赤道以南为负，以北为正；中央经线以东为正，以西为负。我国位于北半球，故纵坐标均为正值，但为避免中央经度线以西为负值的情况，将坐标纵轴西移500公里。
北京54坐标系：1954年我国在北京设立了大地坐标原点，采用克拉索夫斯基椭球体，依此计算出来的各大地控制点的坐标，称为北京54坐标系。
GS84坐标系：即世界通用的经纬度坐标系。
经纬度表示：1884年国际经度会议规定，以通过英国伦敦格林威治天文台子午仪中心的经线为0°经线。从0°经线往东叫东经，往西叫西经，东、西各分180°。习惯上以西经20°和东经160°为分界把地球分为东西两个半球。假如从地轴的正中间将地球切成南北两半，上边的一半叫北半球，下边的一半叫南半球。被切的这个平面，叫赤道面。赤道面与地球表面相交的线叫赤道。纬线从赤道往两极越来越短，到了两极就缩小成一个点了。科学家们把赤道定为0°纬线，从赤道向两极各分为90°，赤道以南叫南纬，赤道以北叫北纬。在计算机或GPS上经纬度经常用度、分、秒和度.度、分.分、秒.秒的混合方式进行表示，度、分、秒间的进制是60进制，度.度、分.分、秒.秒的进制是100进制，换算时一定要注意。可以近似地认为每个纬度之间的距离是不变的111KM,每分间
1.85KM，每秒间31.8M。经度间的距离随纬度增高逐渐减小，可按以下公式计算：
经度1°长度=111.413cosφ一0.094cos3φ公里(纬度φ处)。
一般从GPS得到的数据是经纬度。经纬度有多种表示方法。
1.)& ddd.ddddd， 度 .
度的十进制小数部分（5位）例如：31.12035&
2.) ddd.mm.mmm，度 . 分 . 分的十进制小数部分（3位）例如 31&10.335&
3.) ddd.mm.ss, 度 . 分 . 秒 例如 31&12’42&P
地球上任何一个固定的点都可以用确定的经纬度表示出来。
关于经纬度坐标转换的方法
一、十进制转换成经纬度
把经纬度转换成十进制的方法很简单
如下就可以了
Decimal Degrees = Degrees + minutes/60 + seconds/3600
例：57°55’56.6&P =57+55/60+56.6/
114°65’24.6&P=114+65/60+24.6/3600=结果自己算!
如把经纬度 (longitude，latitude)
(205.，57.8)
转换据成坐标(Degrees,minutes,seconds)(205°23’44.1&P，57°55’56.6&P)。
步骤如下：
1， 直接读取”度”：205
2，(205.-205)*60=23. 得到”分”：23
3，(23.-23)*60=44. 得到”秒”：44.1
采用同样的方法可以得到纬度坐标：57°55’56.6&P
好大一篇，都是别的地方copy过来的！
&本人使用的gps模块返回的数据格式是0183gps协议，返回的十进制经纬度，在地图显示需以上转换
以上网友发言只代表其个人观点，不代表新浪网的观点或立场。1919人阅读
转自：/gis/9e/338.html
现在好多人在使用ArcObject的时候都可能需要作经纬度和平面坐标的相互转换。由于经纬度是球面坐标，平面坐标是X-Y的笛卡尔坐标系统，所以这是一个看起来比较难的问题。
& & & 好多人一上来就搬出地图学、地图投影学或者测绘学中的投影公式（如高斯投影 或 墨卡托投影），又是基准坐标又是角度的搞的人头都大了。
&&&&& 实际上要想实现这个功能非常easy。 ARC Engine中的IPoint就可以进行投影和反投影运算了。
投影过程（C#）：
///& flatref 投影的坐标系统，这里的54013是世界投影& ，世界投影所有经纬度都可以转换为平面坐标，但是由于投影面积大失真也会比较大（相当于把整个地球劈成一片片的，然后拉伸最后贴到平面上，失真当然大了）。当然也可以选择精度更高的平面如：esriSRProjCS_Beijing1954GK_23N&&&&&&对应数&#& 仅仅把北京附近的地球平面拉伸铺在平面上，由于投影面积变小，所以投影经度提高。但是由于面积变小，所以有些经度纬度不能转换，比如所美国的精度纬度用北京投影就投不了。
&&&&&&&&flatref = pfactory.CreateProjectedCoordinateSystem(54013);
//没什么说的，标准大地经纬度，可以将X-Y逆投影为经度和纬度
& earthref = pfactory.CreateGeographicCoordinateSystem((int)esriSRGeoCSType.esriSRGeoCS_NAD1983);
&&&&&&& /// 将经纬度点转换为平面坐标。
&&&&&&& private IPoint GetProject(double x, double y)
&&&&&&&&&&& IPoint pt = new PointClass();
&&&&&&&&&&
&&&&&&&&&&& pt.PutCoords(x, y);
&&&&&&&&&&& IGeometry geo = (IGeometry)
&&&&&&&&&&& geo.SpatialReference =
&&&&&&&&&&& geo.Project(flatref);
&&&&&&&&&&&
&&&&&&&&&& /// 将平面坐标转换为经纬度。
&&&&&&&& private IPoint GetGeo(double x, double y)
&&&&&&&&&&& IPoint pt = new PointClass();
&&&&&&&&&&& pt.PutCoords(x, y);
&&&&&&&&&&& IGeometry geo = (IGeometry)
&&&&&&&&&&& geo.SpatialReference =
&&&&&&&&&&& geo.Project(earthref);
&&&&&&&&&&& double xx = pt.X;
&&&&&&&&&&&&
&实际上IPoint的投影和任何地图都没什么大关系，完全可以不用地图，直接调用IPoint进行投影的转换
* 以上用户言论只代表其个人观点，不代表CSDN网站的观点或立场
访问：126673次
积分：1912
积分：1912
排名：第10865名
原创：60篇
转载：36篇
评论：25条
(2)(3)(1)(1)(3)(2)(4)(2)(2)(4)(3)(6)(9)(16)(27)(11)}