适合男性喝的养生茶:投影坐标++地理坐标

投影坐标++地理坐标 1 地理信息中各种坐标系区别和转换总结  
一、北京54坐标到西安80坐标转换小结  

1、北京54和西安80是两种不同的大地基准面，不同的参考椭球体，因而两种地图下，同一个点的坐标是不同的，无论是三度带六度带坐标还是经纬度坐标都是不同的。  
2、数字化后的得到的坐标其实不是WGS84的经纬度坐标，因为54和80的转换参数至今没有公布，一般的软件中都没有54或80投影系的选项，往往会选择WGS84投影。  
3、WGS84、北京54、西安80之间，没有现成的公式来完成转换。  
4、对于54或80坐标，从经纬度到平面坐标（三度带或六度带）的相互转换可以借助软件完成。  
5、54和80间的转换，必须借助现有的点和两种坐标，推算出变换参数，再对待转换坐标进行转换。（均靠软件实现）  
6、在选择参考点时，注意不能选取河流、等高线、地名、高程点，公路尽量不选。这些在两幅地图上变化很大，不能用作参考。而应该选择固定物，如电站，桥梁等。  


二、西安80坐标系与北京54坐标系转换  

  西安80坐标系与北京54坐标系其实是一种椭球参数的转换作为这种转换在同一个椭球里的转换都是严密的，而在不同的椭球之间的转换是不严密，因此不存在一套转换参数可以全国通用的，在每个地方会不一样，因为它们是两个不同的椭球基准。那么，两个椭球间的坐标转换，一般而言比较严密的是用七参数布尔莎模型，即 X 平移， Y 平移， Z 平移， X 旋转（WX）， Y 旋转（WY）， Z 旋转（WZ），尺度变化（DM ）。要求得七参数就需要在一个地区需要 3 个以上的已知点。如果区域范围不大， 最远点间的距离不大于 30Km（ 经验值 ） ，这可以用三参数，即 X 平移， Y 平移， Z 平移，而将 X 旋转， Y 旋转， Z 旋转，尺度变化面DM视为 0 。  


在MAPGIS平台中实现步骤：  


第一步：向地方测绘局（或其它地方）找本区域三个公共点坐标对（即54坐标x，y，z和80坐标x，y，z）；

第二步：将三个点的坐标对全部转换以弧度为单位。（菜单：投影转换/输入单点投影转换，计算出这三个点的弧度值并记录下来）

第三步：求公共点求操作系数（菜单：投影转换/坐标系转换）。如果求出转换系数后，记录下来。

第四步：编辑坐标转换系数。（菜单：投影转换/编辑坐标转换系数。）最后进行投影变换，“当前投影”输入80坐标系参数，“目的投影”输入54坐标系参数。进行转换时系统会自动调用曾编辑过的坐标转换系数。  



三、地理坐标系与投影坐标系的区别  

1、首先理解地理坐标系（Geographic coordinate system），Geographic coordinate system直译为地理坐标系统，是以经纬度为地图的存储单位的。很明显，Geographic coordinate system是球面坐标系统。我们要将地球上的数字化信息存放到球面坐标系统上，如何进行操作呢？地球是一个不规则的椭球，如何将数据信息以科学的方法存放到椭球上？这必然要求我们找到这样的一个椭球体。这样的椭球体具有特点：可以量化计算的。具有长半轴，短  
半轴，偏心率。以下几行便是Krasovsky_1940椭球及其相应参数。  
Spheroid: Krasovsky_1940  
Semimajor Axis: 6378245.000000000000000000  
Semiminor Axis: 6356863.018773047300000000  
Inverse Flattening（扁率）: 298.300000000000010000  
  然而有了这个椭球体以后还不够，还需要一个大地基准面将这个椭球定位。在坐标系统描述中，可以看到有这么一行：  
Datum: D_Beijing_1954  
表示，大地基准面是D_Beijing_1954。  
有了Spheroid和Datum两个基本条件，地理坐标系统便可以使用。  
完整参数：  
Alias:  
Abbreviation:  
Remarks:  
Angular Unit: Degree (0.017453292519943299)  
Prime Meridian（起始经度）: Greenwich (0.000000000000000000)  
Datum（大地基准面）: D_Beijing_1954  
Spheroid（参考椭球体）: Krasovsky_1940  
Semimajor Axis: 6378245.000000000000000000  
2 地理信息中各种坐标系区别和转换总结  
Semiminor Axis: 6356863.018773047300000000  
Inverse Flattening: 298.300000000000010000  


2、接下来便是Projection coordinate system（投影坐标系统），首先看看投影坐标系统中的一些参数。  
Projection: Gauss_Kruger  
Parameters:  
False_Easting: 500000.000000  
False_Northing: 0.000000  
Central_Meridian: 117.000000  
Scale_Factor: 1.000000  
Latitude_Of_Origin: 0.000000  
Linear Unit: Meter (1.000000)  
Geographic Coordinate System:  
Name: GCS_Beijing_1954  
Alias:  
Abbreviation:  
Remarks:  
Angular Unit: Degree (0.017453292519943299)  
Prime Meridian: Greenwich (0.000000000000000000)  
Datum: D_Beijing_1954  
Spheroid: Krasovsky_1940  
Semimajor Axis: 6378245.000000000000000000  
Semiminor Axis: 6356863.018773047300000000  
Inverse Flattening: 298.300000000000010000  
从参数中可以看出，每一个投影坐标系统都必定会有Geographic Coordinate System。  
投影坐标系统，实质上便是平面坐标系统，其地图单位通常为米。  
那么为什么投影坐标系统中要存在坐标系统的参数呢？  
这时候，又要说明一下投影的意义：将球面坐标转化为平面坐标的过程便称为投影。  
好了，投影的条件就出来了：  
a、球面坐标  
b、转化过程（也就是算法）  
也就是说，要得到投影坐标就必须得有一个“拿来”投影的球面坐标，然后才能使用算法去投影！  
即每一个投影坐标系统都必须要求有Geographic Coordinate System参数。  
3、我们现在看到的很多教材上的对坐标系统的称呼很多，都可以归结为上述两种投影。其中包括我们常见的“非地球投影坐标系统”。）：  


大地坐标（Geodetic Coordinate）:大地测量中以参考椭球面为基准面的坐标。地面点P的位置用大地经度L、大地纬度B和大地高H表示。当点在参考椭球面上时，仅用大地经度和大地纬度表示。大地经度是通过该点的大地子午面与起始大地子午面之间的夹角，大地纬度是通过该点的法线与赤道面的夹角，大地高是地面点沿法线到参考椭球面的距离。  


方里网:是由平行于投影坐标轴的两组平行线所构成的方格网。因为是每隔整公里绘出坐标纵线和坐标横线，所以称之为方里网，由于方 里线同时 又是平行于直角坐标轴的坐标网线，故又称直角坐标网。  
在1：1万——1：20万比例尺的地形图上，经纬线只以图廓线的形式直接表现出来，并在图角处注出相应度数。为了在用图时加密成 网，在内外图廓间还绘有加密经纬网的加密分划短线(图式中称“分度带”)，必要时对应短线相连就可以构成加密的经纬线网。1：2 5万地形图上，除内图廓上绘有经纬网的加密分划外，图内还有加密用的十字线。  


我国的1：50万——1：100万地形图，在图面上直接绘出经纬线网，内图廓上也有供加密经纬线网的加密分划短线。  


直角坐标网的坐标系以中央经线投影后的直线为X轴，以赤道投影后的直线为Y轴，它们的交点为坐标原点。这样，坐标系中就出现了四 个象限。纵坐标从赤道算起向北为正、向南为负；横坐标从中央经线算起，向东为正、向西为负。  


虽然我们可以认为方里网是直角坐标，大地坐标就是球面坐标。但是我们在一副地形图上经常见到方里网和经纬度网，我们很习惯的称经 纬度网为大地坐标，这个时候的大地坐标不是球面坐标，她与方里网的投影是一样的（一般为高斯），也是平面坐标  

四、GIS中的坐标系定义与转换  

1. 椭球体、基准面及地图投影  

GIS中的坐标系定义是GIS系统的基础，正确定义GIS系统的坐标系非常重要。GIS中的坐标系定义由基准面和地图投影两组参数确定，而基准面的定义则由特定椭球体及其对应的转换参数确定，因此欲正确定义GIS系统坐标系，首先必须弄清地球椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念及它们之间的关系。  


  
作者：giskaiser2008-4-11 12:10 回复此发言    

--------------------------------------------------------------------------------

3 地理信息中各种坐标系区别和转换总结  

基准面是利用特定椭球体对特定地区地球表面的逼近，因此每个国家或地区均有各自的基准面，我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系，1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国新的大地坐标系--西安80坐标系，目前大地测量基本上仍以北京54坐标系作为参照，北京54与西安80坐标之间的转换可查阅国家测绘局公布的对照表。 WGS1984基准面采用WGS84椭球体，它是一地心坐标系，即以地心作为椭球体中心，目前GPS测量数据多以WGS1984为基准。  

上述3个椭球体参数如下：  

  
椭球体与基准面之间的关系是一对多的关系，也就是基准面是在椭球体基础上建立的，但椭球体不能代表基准面，同样的椭球体能定义不同的基准面，如前苏联的Pulkovo 1942、非洲索马里的Afgooye基准面都采用了Krassovsky椭球体，但它们的基准面显然是不同的。  

地图投影是将地图从球面转换到平面的数学变换，如果有人说：该点北京54坐标值为X=4231898,Y=21655933，实际上指的是北京54基准面下的投影坐标，也就是北京54基准面下的经纬度坐标在直角平面坐标上的投影结果。  

2. GIS中基准面的定义与转换  

虽然现有GIS平台中都预定义有上百个基准面供用户选用，但均没有我们国家的基准面定义。假如精度要求不高，可利用前苏联的Pulkovo 1942基准面(Mapinfo中代号为1001)代替北京54坐标系；假如精度要求较高，如土地利用、海域使用、城市基建等GIS系统，则需要自定义基准面。  

GIS系统中的基准面通过当地基准面向WGS1984的转换7参数来定义，转换通过相似变换方法实现，具体算法可参考科学出版社1999年出版的《城市地理信息系统标准化指南》第76至86页。假设Xg、Yg、Zg表示WGS84地心坐标系的三坐标轴，Xt、Yt、Zt表示当地坐标系的三坐标轴，那么自定义基准面的7参数分别为：三个平移参数ΔX、ΔY、ΔZ表示两坐标原点的平移值；三个旋转参数εx、εy、εz表示当地坐标系旋转至与地心坐标系平行时，分别绕Xt、Yt、Zt的旋转角；最后是比例校正因子，用于调整椭球大小。  


MapX中基准面定义方法如下：  

Datum.Set(Ellipsoid, ShiftX, ShiftY, ShiftZ, RotateX, RotateY, RotateZ, ScaleAdjust, PrimeMeridian)  

其中参数： Ellipsoid为基准面采用的椭球体；  
ShiftX, ShiftY, ShiftZ为平移参数；  
RotateX, RotateY, RotateZ为旋转参数；  
ScaleAdjust为比例校正因子，以百万分之一计；  
PrimeMeridian为本初子午线经度，在我国取0，表示经度从格林威治起算。  

美国国家测绘局(National Imagery and Mapping Agency)公布了世界大多数国家的当地基准面至WGS1984基准面的转换3参数(平移参数)，可从 http://164.214.2.59/GandG/wgs84dt/dtp.html 下载，其中包括有香港Hong Kong 1963基准面、台湾 Hu-Tzu-Shan 基准面的转换3参数，但是没有中国大陆的参数。  

实际工作中一般都根据工作区内已知的北京54坐标控制点计算转换参数，如果工作区内有足够多的已知北京54与WGS84坐标控制点，可直接计算坐标转换的7参数或3参数；当工作区内有3个已知北京54与WGS84坐标控制点时，可用下式计算WGS84到北京54坐标的转换参数(A、B、C、D、E、F)：x54 = AX84 + BY84 + C，y54 = DX84 + EY84 + F，多余一点用作检验；在只有一个已知控制点的情况下(往往如此)，用已知点的北京54与WGS84坐标之差作为平移参数，当工作区范围不大时精度也足够了。  

从Mapinfo中国的URL(http://www.mapinfo.com.cn/download)可下载到包含北京54、西安80坐标系定义的Mapinfow.prj文件，其中定义的北京54基准面参数为：(3,24,-123,-94,-0.02,0.25,0.13,1.1,0)，西安80基准面参数为：(31,24,-123,-94,-0.02,0.25,0.13,1.1,0)，文件中没有注明其参数的来源，我发现它们与Mapinfo参考手册附录G"定义自定义基准面"中的一个例子所列参数相同，因此其可靠性值得怀疑，尤其从西安80与北京54采用相同的7参数来看，至少西安80的基准面定义肯定是不对的。因此，当系统精度要求较高时，一定要对所采用的参数进行检测、验证，确保坐标系定义的正确性。  
WGS84经纬度坐标与北京54坐标或者西安80坐标的关系2008-11-07 17:58一般来讲，GPS直接提供的坐标（B,L,H）是1984年世界大地坐标系(Word Geodetic System 1984即WGS-84)的坐标，其中B为纬度，L为经度，H为大地高即是到WGS-84椭球面的高度。而在实际应用中，我国地图采用的是1954北京坐标系或者1980西安坐标系下的高斯投影坐标（x,y,），不过也有一些电子地图采用1954北京坐标系或者1980西安坐标系下的经纬度坐标（B,L），高程一般为海拔高度h。


GPS的测量结果与我国的54系或80系坐标相差几十米至一百多米，随区域不同，差别也不同，经粗落统计，我国西部相差70米左右，东北部140米左右，南部75米左右，中部45米左右。现就上述几种坐标系进行简单介绍，供大家参阅，并提供各坐标系的基本参数，以便大家在使用过程中自定义坐标系。


1、1984世界大地坐标系

WGS-84坐标系是美国国防部研制确定的大地坐标系，是一种协议地球坐标系。WGS-84坐标系的定义是：原点是地球的质心，空间直角坐标系的Z轴指向BIH（1984.0）定义的地极（CTP）方向，即国际协议原点CIO，它由IAU和IUGG共同推荐。X轴指向BIH定义的零度子午面和CTP赤道的交点，Y轴和Z，X轴构成右手坐标系。WGS-84椭球采用国际大地测量与地球物理联合会第17届大会测量常数推荐值，采用的两个常用基本几何参数：


长半轴a=6378137m；扁率f=1:298.257223563。


2、1954北京坐标系


1954北京坐标系是将我国大地控制网与前苏联1942年普尔科沃大地坐标系相联结后建立的我国过渡性大地坐标系。属于参心大地坐标系，采用了前苏联的克拉索夫斯基椭球体。其长半轴 a=6378245，扁率 f=1/298.3。1954年北京坐标系虽然是苏联1942年坐标系的延伸,但也还不能说它们完全相同。　　


3、1980西安坐标系


1978年，我国决定建立新的国家大地坐标系统,并且在新的大地坐标系统中进行全国天文大地网的整体平差,这个坐标系统定名为1980年西安坐标系。属参心大地坐标系。1980年西安坐标系Xi'an Geodetic Coordinate System 1980 采用1975国际椭球,以JYD 1968.0系统为椭球定向基准,大地原点设在陕西省泾阳县永乐镇,采用多点定位所建立的大地坐标系.其椭球参数采用1975年国际大地测量与地球物理联合会推荐值,它们为:其长半轴a=6378140m; 扁率f=1/298.257。


4 高斯平面直角坐标系和UTM


一般的地图均为平面图,其对应的也是平面坐标.因此,需要将椭球面上各点的大地坐标,按照一定的数学规律投影到平面上成为平面直角坐标.目前世界各国采用最广泛的高斯- 克吕格投影和墨卡托投影(UTM)均是正形投影（等角投影）， 即该投影在小区域范围内使平面图形与椭球面上的图形保持相似。为了限制长度变形,，根据国际测量协会规定,将全球按一定经差分成若干带。我国采用6度带或3度带，6度带是自零度子午线起每隔经度。


高斯平面直角坐标系一般以中央经线（L0）投影为纵轴X, 赤道投影为横轴Y，两轴交点即为各带的坐标原点。为了避免横坐标出现负值，在投影中规定将坐标纵轴西移500公里当作起始轴。为了区别某一坐标系统属于哪一带，通常在横轴坐标前加上带号，如(4231898m,21655933m)，其中21即为带号。 城建坐标多采用三度带的高斯-克吕格投影。同一坐标系下的大地坐标（即经纬度坐标B,L）与其对应的高斯平面直角坐标（x，y）有严格的转换关系。现行的测绘的教科书的一般都有。


5、 地方独立坐标系
在我国许多城市测量与工程测量中，若直接采用国家坐标系下的高斯平面直角坐标，则可能会由于远离中央子午线，或由于测区平均高程较大，而导致长度投影变形较大，难以满足工程上或实用上的精度要求。另一方面，对于一些特殊的测量,如大桥施工测量，水利水坝测量，滑坡变形监测等，采用国家坐标系在实用中也会很不方便。因此，基于限制变形，以及方便实用，科学的目的，在许多城市和工程测量中,常常会建立适合本地区的地方独立坐标系。建立地方独立坐标系，实际上就是通过一些元素的确定来决定地方参考椭球与投影面.地方参考椭球一般选择与当地平均高程相对应的参考椭球，该椭球的中心，轴向和扁率与国家参考椭球相同。其椭球半径α1增大为:α1=α+Δα1,Δα1=Hm+ζ0式中:Hm为当地平均海拔高程,ζ0为该地区的平均高程异常。而地方投影面的确定中,选取过测区中心的经线或某个起算点的经线作为独立中央子午线.以某个特定方便使用的点和方位为地方独立坐标系的起算原点和方位,并选取当地平均高程面Hm为投影面。



既然说到了不同的坐标系，就存在坐标转换的问题。关于坐标转换，首先要搞清楚转换的严密性问题，即在同一个椭球里的坐标转换都是严密的，而在不同的椭球之间的转换这时不严密的。例如，由1954北京坐标系的大地坐标转换到954北京坐标系的高斯平面直角坐标是在同一参考椭球体范畴内的坐标转换，其转换过程是严密的。由1954北京坐标系的大地坐标转换到WGS-84的大地坐标，就属于不同椭球体间的转换。


不同椭球体间的坐标转换在局部地区的采用的常用办法是相似变换法，即利用部分分布相对合理高等级公共点求出相应的转换参数。一般而言，比较严密的是用七参数的相似变换法，即X平移，Y平移，Z平移，X旋转，Y旋转，Z旋转，尺度变化K。要求得七参数就需要在一个地区需要3个以上的已知点，如果区域范围不大，最远点间的距离不大于30Km(经验值)，这可以用三参数，即X平移，Y平移，Z平移，而将X旋转，Y旋转，Z旋转，尺度变化K视为0，所以三参数只是七参数的一种特例。


如果不考虑高程的影响，对于不同椭球体下的高斯平面直角坐标可采用四参数的相似变换法，即四参数（x平移，y平移，尺度变化m，旋转角度α）。如果用户要求的精度低于20米，在一定范围（2＇＊2＇）内，就直接可以用二参数法（ΔB，ΔL）或（Δx，Δy）修正。但在实际操作中，这也取决于选取的公共点是否合理，并保证其足够的精度。