GIS的基础知识

论坛

一些坐标系统WGS1984、UTM的知识
发布时间: 2006-8-07 11:20    作者: 御天六龙    信息来源: 100GPS-中国最大的GPS网|导航地图软件下载|报价团购手机PDA车载蓝牙方案GPS价格|模块定位卫星系统论文
字体:        打印
广告位招租
看到坛子里,对SMARTCOMGPS的地图校准大多一知半解,为此,特地搜集了一些关于的知识,供大家了解、参考。

自已也跟着大家一起学习!!!

MapInfo地图投影的添加
——吐鲁番地区坐标系

新疆大学资源环境学院地信98-1班 李风云
指导老师:安放舟 教授

【摘要】GIS中的坐标系定义是GIS系统的基础,正确定义GIS系统的坐标系非常重要。在MapInfo地图投影中缺少针对中国很多地区大比例尺
和中比例尺的投影方式,通过对MapInfo中Mapinfow.prj文件的研究,以吐鲁番地区为例,添加了高斯-克吕格(Gauss-Krüger)投影,能解
决一些大比例尺对变形小的要求。
【关键字】地图投影 坐标系 MapInfo

一、前言
地球高低不平、极其复杂的自然表面。为研究和工作方便,常将地球近似为为一个旋转椭球体,称为地球椭球体。地球椭球体的表面是一个
不可展的曲面,地图以平面方式表示地球表面(全部或一部分)。将地球椭球体上的点的坐标投影到平面坐标的方法称为地图投影。地图投影
的种类繁多,不同的投影方式具有不同的形态和变形特征。根据不同的使用目的,可以采用不同的投影方式,一种投影对一种目的是有用的,
而对另一种目的则可能不适合。
地图投影的选取决定于地图的应用及其比例尺大小,在桌面数字环境下,用户可对每一个新创建文件进行投影设置和选择。MapInfo以两级目
录菜单的形式提供了300多个预定义坐标系,当用户要使用其他坐标系或创建新的坐标系时,还可以通过修改投影参数文件(Mapinfow.prj)
来实现。这个数据文件以分行形式记录每一个预定义坐标系的参数表,如坐标系名称、投影代码、基准面代码、坐标单位、原点经度、原点
纬度、标准纬线1、标准纬线2、方位角、比例系数等。MapInfo系统定义了50个基准面代码,它们包括了全世界和一些国家采用的基准面[1]。
如果用户要采用其他的基准面,并且知道该基准面的数学参数,则可以使用代码在该投影参数文件中定义这个基准面。
虽然现有GIS平台中都预定义有上百个基准面供用户选用,但均没有我们国家的基准面定义。假如精度要求不高,可利用前苏联的
Pulkovo 1942基准面(Mapinfo中代号为1001)代替北京54坐标系;假如精度要求较高,如土地利用、海域使用、城市基建等GIS系统,则需要自
定义基准面。

二、MAPINFOW.PRJ文件的参数分析
MAPINFOW.PRJ在MapInfo软件的安装目录下,可以用记事本打开文件,它的格式如下:
"--- Longitude / Latitude ---", 0, 0, 0, 0., 0., 0., 0., 0., 0., 0., 0.
"Longitude / Latitude", 1, 0
"Longitude / Latitude (Adindan)p4201", 1, 1
"Longitude / Latitude (Afgooye)p4205", 1, 2
"Longitude / Latitude (AGD 66)p4202", 1, 12
"Longitude / Latitude (AGD 84)p4203", 1, 13
"Longitude / Latitude (Ain el Abd 1970)p4204", 1, 3
"Longitude / Latitude (Anna 1 Astro 1965)", 1, 4
"Longitude / Latitude (Arc 1950)p4209", 1, 5
第一行是投影方式的一级菜单,后面是二级菜单。每行的第一部分是位于引号内的坐标系名称。每一行的第二部分是指定投影的代号。行中其
余部分是有关该特定坐标系的参数值。下面一以中国区域等积投影为例进行分析,详细分析各参数的含义。Regional Equal-Area Projections
在MAPINFOW.PRJ文件中的代码如下:
"Equal-Area Projection (China)", 9, 0, 0, 110, 10, 25, 40, 0, 0
表 1
投影 正轴等面积割圆锥
基准面 GRS 80基准面
单位 英里
原点 110°E, 10°N
标准纬线 25°N、 40°N
东伪偏移 0
北伪偏移 0
第一个参数 3 代表投影方式是Albers Equal Area Conic(正轴等面积割圆锥投影,也叫亚尔勃斯投影);第二个参数 0 代表GRS 80基准面
(Datum);第三个参数 0 代表单位(Unit)是英里;第四和第五个是代表原点坐标是 110°E, 10°N;第六和第七个代表两条标准纬线分别为
25°N、 40°N;最后两个参数是东伪偏移(FalseEasting)和北伪偏移(FalseNorthing) [2]。Equal-Area Projection (China)的列表(表1)
。由于参数的数目繁多,每种参数在《MapInfo参考手册》中都有详细的表格。
高斯-克吕格(Gauss-Krüger)投影坐标系在MAPINFOW.PRJ文件需要定义的坐标系参数序列如下:
高斯-克吕格:投影代号(Type),基准面(Datum),单位(Unit),
中央经度(OriginLongitude),原点纬度(OriginLatitude),
比例系数(ScaleFactor),
东伪偏移(FalseEasting),北纬偏移(FalseNorthing)

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

椭球体 代号 年代 长半轴 短半轴 1/扁率
Krassovsky 椭球 3 1940 6378245 6356863 298.3
IAG 75椭球 31 1975 6378140 6356755 298.25722101
WGS 84 椭球 28 1984 6378137.000 6356752.314 298.257223563
表 2
椭球体与基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的
基准面。
虽然现有GIS平台中都预定义有上百个基准面供用户选用,但均没有我们国家的基准面定义。假如精度要求不高,可利用前苏联的Pulkovo
1942基准面(Mapinfo中代号为1001)代替北京54坐标系;假如精度要求较高,如土地利用、城市基建等GIS系统,则需要自定义基准面。
从Mapinfo中国的URL( http://www.mapinfo.com.cn)可下载到包含北京54、西安80坐标系,其中定义:
北京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中两个特殊基准面代号:0和9999,是留给用户自定义的参数(在Mapinfo参考手册的附录G中有详细说明)。这样就可以自定义一
个基准面参数为“9999”的基准面[3]。这样就可以以我国常用的基准面建立坐标系,也可以选择不同的投影方式,以满足系统对精度的要求。

WGS1984基准面采用WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心。WGS1984在Mapinfow.prj文件中以代号104作为基准面参数,
当测数据为WGS84坐标数据时WGS1984基准面参数为:104

3.2高斯-克吕格(Gauss-Krüger)投影
我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯-克吕格投影
(Gauss-Kruger),又叫横轴墨卡托投影(Transverse Mercator);小于50万的地形图采用正轴等角割园锥投影,又叫兰勃特投影(Lambert
Conformal Conic),我国的GIS系统中应该采用与我国基本比例尺地形图系列一致的地图投影系统[5]。
高斯-克吕格(Gauss-Krüger)投影又叫横轴墨卡托(Transverse Mercator)投影,依椭圆柱作为投影面,并与椭球体相切于一条经线上,该经
线即为投影带的中央经线,按等角条件将中央经线东西一定范围内的区域投影到椭圆柱表面上,再展成平面,便构成了横轴等角切椭圆柱投影
[6]。
在GIS系统中均采用6°或3°分带的高斯-克吕格投影,因为一般坐标采用的是6°或3°分带的高斯-克吕格投影坐标。高斯-克吕格投影以6°
或3°分带,每一个分带构成一个独立的平面直角坐标网,投影带中央经线投影后的直线为X轴(纵轴,纬度方向),赤道投影后为Y轴(横轴,
经度方向),为了防止经度方向的坐标出现负值,规定每带的中央经线西移500公里,即东伪偏移值为500公里,由于高斯-克吕格投影每一个
投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,因此规定在横轴坐标前加上带号,如(4231898,21655933)其中21即
为带号,同样所定义的东伪偏移值也需要加上带号,如21带的东伪偏移值为21500000米。假如你的工作区位于21带,即经度在120°至126°
范围,该带的中央经度为123°。

3.3吐鲁番地区的分析
吐鲁番地区位于中国西北部的新疆境内,经度差为4.7°,经度位于87°20’ E 和92°E之间。可以采用3°分带的高斯-克吕格投影,其中央
经线为90°。

3.4 创建新的坐标系
前面以对各参数做了分析,针对吐鲁番地区建立坐标系,要对MAPINFOW.PRJ文件增加坐标系的下列参数(表3):
投影 基准面 单位 中央经度 原点纬度 比例系数 东伪偏移 北伪偏移
8 9999 7 90 0 1 30500000 0
表 3
自定义基准面9999为:
北京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)
WGS84基准面参数为:104
需要向MAPINFOW.PRJ文件增加一个列,要有相应要素的新条目。该过程描述如下:
1 在一个文本编辑器或字处理器中打开MAPINFOW.PRJ文件。
2输入代表新坐标系的下列参数:
"--- 吐鲁番坐标系统---"
"Gauss-Kruger 30 (西安1980 3°分带)", 8, 9999, 31, 24, -123, -94, -0.02, 0.25, 0.13, 1.1, 0, 7, 90, 0, 1, 30500000, 0
"Gauss-Kruger 30 (北京1954 3°分带)", 8, 9999, 3, 24, -123, -94, -0.02, 0.25, 0.13, 1.1, 0, 7, 90, 0, 1, 30500000, 0
"Gauss-Kruger 30 (WGS1984 3°分带)", 8, 104,7, 90, 0, 1, 30500000, 0
3 如果有必要把该条目移到相似坐标系之间的合适位置。
4 保存编辑过的MAPINFOW.PRJ文件。
这样就可以使用自定义的坐标系了,就像MapInfo中提供的坐标系一样。如(图1)所示:

图1

四、讨论
从Mapinfo中国提供的北京54、西安80坐标系,定义的北京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),在这两组数据中,有八个参数是相同的,只有参考椭球体的参数是不同的。北京54、西安80
坐标系的大地原点也是不同的,参考椭球体与基准面之间是一对多的关系,显然这两组参数中至少有一组是错的。网站中没有注明其参数的来
源,也没有各参数的含义,也就无法纠正坐标系在大地原点的上的错误。因此,当系统精度要求较高时,一定要对所采用的参数进行检测、验
证,确保坐标系定义的正确性。
如果能够得知各参数的含义,那么也就可以轻松的把北京54、西安80坐标系的数据进行转换,也可以将其他坐标转换成北京54、西安80坐标系
的数据,只需要通过一个方程进行转换,不许要通过GIS的专业软件。

五、结束语
地图投影在GIS领域中占有重要的地位,高斯-克吕格(Gauss-Krüger)投影在我国又比较通用,通过吐鲁番地区高斯-克吕格
(Gauss-Krüger)投影坐标系的定义,也可以定义一些大比例尺的UTM投影(Universal Transverse Mercatol Projection,
即通用横轴墨卡托投影)、等积割圆锥投影等等,提高系统的投影精度。在研究某一区域时,可以依照系统对地图投影精度的要求,
依该区域的范围,设计独特的坐标系,坐标系是GIS系统的基石,正确设置坐标系是系统成败的关键。






【参考文献】
[1] 蔡孟裔等 新编地图学教程 高等教育出版社 2000
[2] MapInfo 参考手册
[3] MapInfo China 网站 http://www.mapinfo.com.cn
[4] MapX中的坐标系定义与转换 国家基础地理信息中心网站 http://www.gischina.com/bbs
[5] 地图投影和图像到地图的几何纠正 http://www.gisres.com/article-pdf/map_prj.pdf
[6] 胡旒钜,龚建文 地图投影 测绘出版社 1982
[7] 李英奎 地图制图中地图投影的应用与计算机实现 mapinfo应用集锦
[8] Map Projections Poster http://mac.usgs.gov/mac/isb/pubs/mapprojections/projections.html
[9] 张力果等 地图学 高等教育出版社 1998

[ 本帖最后由 御天六龙 于 2006-8-7 11:33 AM 编辑 ]

相关回复

御天六龙 (2006-8-07 11:21:47)
IKONOS影像等数据采用的是WGS84 UTM坐标系统

WGS84指的是椭球系统,而UTM指的是坐标系统,坐标系统以椭球系统为基础
把同是WGS84椭球系统下的不同坐标系统数据互相转换.
假设这么一个前提:你的原始数据和你所想要的数据同是WGS84下的数据
你的原始数据是采用地理坐标系统的数据(也就是经纬度),可以参考以下VB脚本:

UserInput_Lon = InputBox("请输入你所需地点的经度值"&Chr(13) & Chr(10)&Chr(13) & Chr(10)&"格式:dd.dddddd")
UserInput_Lat = InputBox("请输入你所需地点的纬度值"&Chr(13) & Chr(10)&Chr(13) & Chr(10)&"格式:dd.dddddd")
'计算用户UTM区域号
UTM_Zone=((Int(UserInput_Lon)-(Int(UserInput_Lon) Mod 6)+180)/6+1)
UTM_Zone_Center_Lon=UTM_Zone*6-180-3
'以下定义了WGS84椭球的参数
PI=3.1415926
WGS84_a=6378137
WGS84_1e2=0.00669437999013
WGS84_2e2=0.006739496742227
WGS84_N=WGS84_a/((1-WGS84_1e2*(sin(UserInput_Lat*PI/180))^2)^(1/2))
WGS84_m=cos(UserInput_Lat*PI/180)*(UserInput_Lon-UTM_Zone_Center_Lon)*PI/180
WGS84_nita2=(WGS84_2e2)*(cos(UserInput_Lat*PI/180)^2)
WGS84_t=tan(UserInput_Lat*PI/180)
'采用IUGG1975椭球计算子午线长度,如有WGS84参数则将下式替换为WGS84数值将会更精确
IUGG1975_PRO1=111133.0047*UserInput_Lat
IUGG1975_PRO2=16038.5282*sin(2*UserInput_Lat*PI/180)
IUGG1975_PRO3=16.8326*sin(4*UserInput_Lat*PI/180)
IUGG1975_PRO4=0.0220*sin(6*UserInput_Lat*PI/180)
IUGG1975_X=IUGG1975_PRO1-IUGG1975_PRO2+IUGG1975_PRO3-IUGG1975_PRO4
'计算WGS84与IUGG1975混合椭球下的UTM坐标
POR1=(WGS84_m^2)/2
PRO2=(5-WGS84_t^2+9*WGS84_nita2+4*WGS84_nita2*WGS84_nita2)*(WGS84_m^4)/24
PRO3=(6*(UserInput_Lon-UTM_Zone_Center_Lon)-58*WGS84_t^2+WGS84_t^4)*WGS84_m^6/720
PRO4=(1-WGS84_t^2+WGS84_nita2)*WGS84_m^3/6
PRO5=(5-18*WGS84_t^2+WGS84_t^4+14*WGS84_nita2-58*WGS84_nita2*WGS84_t^2)*WGS84_m^5/120
'计算出的用户输入点的UTM坐标
UserInput_Y=IUGG1975_X+WGS84_N*WGS84_t*(PRO1+PRO2+PRO3)
UserInput_X=WGS84_N*(WGS84_m+PRO4+PRO5)+500000
MyVar=MsgBox("UTM_Zone="&UTM_Zone&" UTM_X="&UserInput_X&" UTM_Y="&UserInput_Y,0,"计算结果:")

以上代码是我以前进行大致计算用的,可能有错误。
把以上代码粘贴到文本文件,并将文本文件改名为.vbs双击运行


另外一种方法是把你的点存为格式化文本(例如.csv)文件,
然后使用globalmapper导入,(file->open generic ascii text file)
导入后设置想要的坐标系统(tools->config->projection)
软件会自动帮你计算出新的坐标系统下的坐标值
完成之后把数据另存为你想要的格式就可以了(例如.csv)(file->export)

当然,在globalmapper的file菜单下有一个直接进行转换的选项
file->batch covert/reproject,用这个也许更方便
御天六龙 (2006-8-07 11:28:33)
转自: http://www.cnblogs.com/Tangf/archive/2006/04/17/377638.html
地理坐标系统简介
地理坐标系,也可称为真实世界的坐标系,是用于确定地物在地球上位置的坐标系。一个特定的地理坐标系是由一个特定的椭球体和一种特定的地图投影构成,其中椭球体是一种对地球形状的数学描述,而地图投影是将球面坐标转换成平面坐标的数学方法。绝大多数的地图都是遵照一种已知的地理坐标系来显示坐标数据。
1.地球椭球体
地球是一个表面很复杂的球体,人们以假想的平均静止的海水面形成的“大地体”为参照,推求出近似的椭球体,理论和实践证明,该椭球体近似一个以地球短轴为轴的椭园而旋转的椭球面,这个椭球面可用数学公式表达,将自然表面上的点归化到这个椭球面上,就可以计算了。下面列举了一些常用的一些椭球及参数:
1)海福特椭球(1910)                   我国52年以前采用的椭球
  a=6378388m b=6356911.9461279m α=0.33670033670
2)克拉索夫斯基椭球(1940 Krassovsky)           北京54坐标系采用的椭球
  a=6378245m b=6356863.018773m α=0.33523298692
3)1975年I.U.G.G推荐椭球(国际大地测量协会1975)      西安80坐标系采用的椭球
  a=6378140m b=6356755.2881575m α=0.0033528131778
4)WGS-84椭球(GPS全球定位系统椭球、17届国际大地测量协会) WGS-84坐标系椭球
  a=6378137m b=6356752.3142451m α=0.00335281006247
最常用的地理坐标系是经纬度坐标系,这个坐标系可以确定地球上任何一点的位置,如果我们将地球看作一个椭球体,而经纬网就是加在地球表面的地理坐标参照系格网,经度和纬度是从地球中心对地球表面给定点量测得到的角度,经度是东西方向,而纬度是南北方向,经线从地球南北极穿过,而纬线是平行于赤道的环线。地理坐标可分为天文地理坐标和大地地理坐标:天文地理坐标是用天文测量方法确定的,大地地理坐标是用大地测量方法确定的。我们在地球椭球面上所用的地理坐标系属于大地地理坐标系,简称大地坐标系。
确定椭球的大小后,还要进行椭球定向,即把旋转椭球面套在地球的一个适当的位置,这一位置就是该地理坐标系的“坐标原点”,是全部大地坐标计算的起算点,俗称“大地原点”。
需要说明的是经纬度坐标系不是一种平面坐标系,因为度不是标准的长度单位,不可用其量测面积长度;平面坐标系(又称笛卡儿坐标系),因其具有以下特性:可量测水平X方向和竖直Y方向的距离,可进行长度、角度和面积的量测,可用不同的数学公式将地球球体表面投影到二维平面上而得到广泛的应用。而每一个平面坐标系都有一特定的地图投影方法。
2.地图投影
是为解决由不可展的椭球面描绘到平面上的矛盾,用几何透视方法或数学分析的方法,将地球上的点和线投影到可展的曲面(平面、园柱面或圆锥面)上,将此可展曲面展成平面,建立该平面上的点、线和地球椭球面上的点、线的对应关系。
地图投影的过程是可以想象用一张足够大的纸去包裹地球,将地球上的地物投射到这张纸上。地球表面投影到平面上、圆锥面或者圆柱面上,然后把圆锥面、圆柱面沿母线切开后展成平面。根据这张纸包裹的方式,地图投影又可以分成:方位投影、圆锥投影和圆柱投影。根据这张纸与地球相交的方式,地图投影又可以分成切投影和割投影,在切线或者割线上的地物是没有变形的,而距离切线或者割线越远变形越大。
还有不少投影直接用解析法得到。根据所借助的几何面不同可分为伪方位投影、伪圆锥投影、伪圆柱投影等。
地图投影会存在两种误差,形状变化(也称角度变化)或者面积变化。投影以后能保持形状不变化的投影,称为等角投影 (Conformal mapping),它的优点除了地物形状保持不变以外,在地图上测量两个地物之间的角度也能和实地保持一致,这非常重要,当在两地间航行必须保持航向的准确;或者另外一个例子是无论长距离发射导弹还是短距离发射炮弹,发射角度必须准确测量出来。因此等角投影是最常被使用的投影。等角投影的缺点是高纬度地区地物的面积会被放大。投影以后能保持形状不变化的投影,称为等面积投影 (Equivalent mapping),在有按面积分析需要的应用中很重要,显示出来的地物相对面积比例准确,但是形状会有变化,假设地球上有个圆,投影后绘制出来即变成个椭圆了。还有第三种投影,非等角等面积投影,意思是既有形状变化也有面积变化,这类投影既不等角也不等积,长度、角度、面积都有变形。其中有些投影在某个主方向上保持长度比例等于1,称为等距投影。
每一种投影都有其各自的适用方面。例如,墨卡托投影适用于海图,其面积变形随着纬度的增高而加大,但其方向变形很小;横轴墨卡托投影的面积变形随着距中央经线的距离的加大而增大,适用于制作不同的国家地图。等角投影常用于航海图、风向图、洋流图等。现在世界各国地形图采用此类投影比较多。等积投影用于绘制经济地区图和某些自然地图。对于大多数数学地图和小比例尺普通地图来说,应优先考虑等积的要求。地理区域,诸如国家、水域和地理分类地区(植被、人口、气候等)相对分布范围,显然是十分重要的内容。 任意投影常用作数学地图,以及要求沿某一主方向保持距离正确的地图。常用作世界地图的投影有墨卡托投影、高尔投影、摩尔威特投影、等差分纬线多圆锥投影、格灵顿投影、桑森投影、乌尔马耶夫投影等。下面对我国地形图所采用的高斯克吕格投影进行简单的介绍。
2.1高斯-克吕格直角坐标
  高斯-克吕格投影(Gauss_Krivger)属于等角横切椭圆柱投影,是设想用一个椭圆柱横套在地球椭球的外面,并与设定的中央经线相切。其经纬线互相垂直,变形最大位于赤道与投影带最外一条经线的交点上,常用于纬度较高地区。
  高斯-克吕格投影分带规定:该投影是我国国家基本比例尺地形图的数学基础,为控制变形,采用分带投影的方法,在比例尺 1:2.5万-1:50万图上采用6°分带,对比例尺为 1:1万及大于1:1万的图采用3°分带。
  6°分带法:从格林威治零度经线起,每6°分为一个投影带,全球共分为60个投影带,东半球从东经0°-6°为第一带,中央经线为3°,依此类推,投影带号为1-30。其投影代号n和中央经线经度L0的计算公式为:L0=(6n-3)°;西半球投影带从180°回算到0°,编号为31-60,投影代号n和中央经线经度L0的计算公式为L0=360-(6n-3)°。
  3°分带法:从东经1°30′起,每3°为一带,将全球划分为120个投影带,东经1°30′-4°30′,...178°30′-西经178°30′,...1°30′-东经1°30′。
  东半球有60个投影带,编号1-60,各带中央经线计算公式:L0=3°n ,中央经线为3°、6°...180°。
  西半球有60个投影带,编号1-60,各带中央经线计算公式:L0=360°-3°n ,中央经线为西经177°、...3°、0°。
我国规定将各带纵坐标轴西移500公里,即将所有y值加上500公里,坐标值前再加各带带号。以18带为例,原坐标值为y=243353.5m,西移后为y=743353.5,加带号通用坐标为y=18743353.5 。
为了方便大家对不同比例尺的地形图检索,最后对我国地形图的分幅与编号规则进行简单的介绍。
3.我国地形图分幅与编号
  我国基本比例尺地形图分幅与编号,以1:100万地形图为基础,延伸出1:50万、1:25万、1:10万,再以1:10万为基础,延伸出1:5万、1:2.5万及1:1万三种比例尺。
  1:100万从赤道起向两极每纬差4°为一行,至88°,南北半球各分为22横列,依次编号A、B、... V;由精度180°西向东每6°一列,全球60列,以1-60表示,如海南所在1:100万图在第5行,第49列,其编号为 E-49 。
  在1:100万图上,按经差3°纬差2°分成四幅1:50万地形图,编为A、B、C、D,如 E-49-A。按经差1°30′纬差1°分成16幅1:25万地形图,编为[1]、...[16],如 E-49-[1]。按经差30′纬差20′分成144幅1:10万地形图,编为1、...144,如 E-49-1。即后三种比例尺各自独立地与1:100万地图的图号联系。
  1:10万图上每经差15′纬差10′分成四幅1:5万地形图,编为A、B、C、D,如 E-49-1-A。 
  1:5万图上每经差7′30″纬差5′分成四幅1:2.5万,编为1、2、3、4,如 E-49-1-A-1。
  1:10万图上每经差3′45″纬差2′30″分成64幅1:1万地形图,编为(1)、...(64),如E-49-1-A-(1)。
  1:1万图上每经差1′52″纬差1′15″分成四幅1:5000地形图,编为a、b、c、d,如E-49-1-A-(1)-a。

[ 本帖最后由 御天六龙 于 2006-8-7 11:31 AM 编辑 ]
御天六龙 (2006-8-07 11:37:01)
坐标转换请参考下面资料:


       我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),这是一个等角横切椭圆柱投影,又叫横轴墨卡托投影(Transverse Mercator);小于50万的地形图采用等角正轴割园锥投影,又叫兰勃特投影(Lambert Conformal Conic);海上小于50万的地形图多用等角正轴圆柱投影,又叫墨卡托投影(Mercator)。一般应该采用与我国基本比例尺地形图系列一致的地图投影系统。

大地基准面的选择


地图坐标系由大地基准面和地图投影确定,大地基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的大地基准面,我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的IAG 75地球椭球体建立了我国新的大地坐标系--西安80坐标系, 目前GPS定位所得出的结果都属于WGS84坐标系统,WGS84基准面采用WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心的坐标系。因此相对同一地理位置,不同的大地基准面,它们的经纬度坐标是有差异的。


本程序中采用的3个椭球体参数如下(源自“全球定位系统测量规范 GB/T 18314-2001”):


椭球体


长半轴


短半轴


Krassovsky


6378245


6356863.0188


IAG 75


6378140


6356755.2882


WGS 84


6378137


6356752.3142


椭球体与大地基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面,如前苏联的Pulkovo 1942、非洲索马里的Afgooye基准面都采用了Krassovsky椭球体,但它们的大地基准面显然是不同的。在目前的GIS商用软件中,大地基准面都通过当地基准面向WGS84的转换7参数来定义,即三个平移参数ΔX、ΔY、ΔZ表示两坐标原点的平移值;三个旋转参数εx、εy、εz表示当地坐标系旋转至与地心坐标系平行时,分别绕Xt、Yt、Zt的旋转角;最后是比例校正因子,用于调整椭球大小。北京54、西安80相对WGS84的转换参数至今没有公开,实际工作中可利用工作区内已知的北京54或西安80坐标控制点进行与WGS84坐标值的转换,在只有一个已知控制点的情况下(往往如此),用已知点的北京54与WGS84坐标之差作为平移参数,当工作区范围不大时,如青岛市,精度也足够了。


以(32°,121°)的25°,40°双标准纬度,原点经度110°,原点纬度10°的兰勃托投影结果为例,北京54及WGS84基准面,两者投影结果在东西方向差距约18米,南北方向差距约45米(见下表),对于几十或几百万的地图来说,这一误差无足轻重,但在工程地图中还是应该加以考虑的。



 


输入坐标(度)


北京54 兰勃托投影(米)


WGS84 兰勃投影(米)


纬度值(X)


32


2531235


2531190


经度值(Y)


121


1028805


1028787



兰勃托投影

(1)兰勃托投影性质


兰勃托(Lambert)投影,又名"等角正割圆锥投影”,由德国数学家兰勃特(J.H.Lambert)在1772年拟定。设想用一个正圆锥切于或割于球面,应用等角条件将地球面投影到圆锥面上,然后沿一母线展开,即为兰勃托投影平面。投影后纬线为同心圆弧,经线为同心圆半径。


兰勃托投影采用双标准纬线相割,与采用单标准纬线相切比较,其投影变形小而均匀,兰勃托投影的变形分布规律是:


a) 角度没有变形,即投影前后对应的微分面积保持图形相似,故亦可称为正形投影;


b) 等变形线和纬线一致,即同一条纬线上的变形处处相等;


c) 两条标准纬线上没有任何变形;


d) 在同一经线上,两标准纬线外侧为正变形(长度比大于1),而两标准纬线之间为负变形(长度比小于1)。因此,变形比较 均匀,变形绝对值也比较小;


e) 同一纬线上等经差的线段长度相等,两条纬线间的经纬线长度处处相等。


我国1:100万地形图采用了兰勃托投影,其分幅原则与国际地理学会规定的全球统一使用的国际百万分之一地图投影一致。纬度按纬差4°分带,从南到北共分成15个投影带,每个投影带单独计算坐标,每带两条标准纬线,第一标准纬线为图幅南端纬度加30′的纬线,第二标准纬线为图幅北端纬度减30′的纬线,这样处于同一投影带中的各图幅的坐标成果完全相同,不同带的图幅变形值接近相等,因此每投影带只需计算其中一幅图(纬差4°,经差6°)的投影成果即可。由于是纬差4°分带投影的,所以当沿着纬线方向拼接地图时,不论多少图幅,均不会产生裂隙;但是,当沿着经线方向拼接时,因拼接线分别处于上下不同的投影带,投影后的曲率不同,致使拼接时产生裂隙。


(2)兰勃托投影坐标


以图幅的原点经线(一般是中央经线)作纵坐标x轴,原点经线与原点纬线(一般是最南端纬线)的交点作为原点,过此点的切线作为横坐标y轴,构成兰勃托平面直角坐标系,此投影两标准纬线无变形。


在我国1:100万兰勃托投影地形图中,由于经纬网图形是以中央经线为轴左右对称的(原点经线选择图幅的中央经线,原点纬线选择最南端纬线)。因此,只要计算右方经差为1°、2°、3°的经纬线交点的坐标,左方的经纬线交点的坐标,只需y值为负即可。一幅图的直角坐标成果可以在同一纬度带中通用。




单点转换


单点转换步骤如下:


(1)选择是兰勃托正转换还是反转换,缺省为经纬度转换到兰勃托投影坐标,投影坐标单位为米。


(2)选择大地基准面,缺省北京54,如果是GPS定位数据别忘了切换为WGS84。


(3)输入第一标准纬度及第二标准纬度,单位度。


(4)输入原点纬度及原点经度,单位度,缺省为零。


(5)如正向投影,选择经纬度输入数据格式,有三个选项,缺省为十进制度格式。具体输入方式如下例:


格 式


原始纬度值


原始经度值


输入纬度值


输入经度值


十进制度


35.445901°


122.997344°


35.445901


122.997344


度分


35°26.7541′


122°59.8406′


3526.7541


12259.8406


度分秒


35°26′45.245″


122°59′50.438″


352645.245


1225950.438


(6)正投影按选定格式在“输入”栏输入经纬度值,反投影输入以米为单位的X、Y坐标值。


(7)单击“单点转换”按钮。


(8)在“输出”栏查看计算结果。





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

--  作者:dongf333
--  发布时间:2005-6-1 8:07:00

--  
批量转换  

    批量转换步骤如下:   


   (1)准备好需要转换的输入数据文件,要求是文本文件,分两列,第一列纬度值或纵向坐标值,第二列经度值或横向坐标值,两列之间用空格分开。正向投影时,纬度值及经度值格式可以有三种选择(见表),缺省当作十进制度处理;反向投影时,纵向及横向坐标值必须以米为单位。


         下例为度分秒格式(WGS84)的兰勃托正投影输入数据文件 testdata.txt


         350000.000  1220000.000
         351600.519  1225959.506
         345800.101  1225959.8
         343600.336  1230000.26
         341400.018  1225959.897
         335159.17   1225959.46
         333000.08   1230000.28
         332300.674  1232000.103


   (2)选择是兰勃托正转换还是反转换,缺省为经纬度转换到兰勃托投影坐标,投影坐标单位为米。


   (3)选择大地基准面,缺省北京54,如果是GPS定位数据别忘了切换为WGS84。


   (4)输入第一标准纬度及第二标准纬度,单位度。


   (5)输入原点纬度及原点经度,单位度,缺省为零。


   (6)如正向投影,选择输入数据文件中的经纬度输入数据格式,有三个选项,缺省为十进制度格式。


   (7)单击“批量转换”按钮。弹出打开文件对话框,输入你的数据文件名。


   (8)输入转换结果文件名,单击“保存”后,程序开始进行计算。


   (9)打开输出文件查看计算结果,结果分五列,第一序号,第二列输入纬度值或纵向坐标值,第三列输入经度值或横向坐标值,第四列转换后纬度值或纵向坐标值,第五列转换后经度值或横向坐标值。   


    下例为度分秒格式(WGS84)的25°,40°双标准纬度,原点经度110°,原点纬度10°的正投影转换结果数据文件 result.txt  


      1  350000      1220000      2869187.6   1084769.1
      2  351600.519  1225959.506  2908969.7   1171137.8
      3  345800.101  1225959.8    2876173.4   1175175.2
      4  343600.336  1230000.26   2836123.8   1180108.2
      5  341400.018  1225959.897  2796066.3   1185021.2
      6  335159.17   1225959.46   2756002.2   1189933.1
      7  333000.08   1230000.28   2716003.3   1194869.2
      8  332300.674  1232000.103  2707084.4   1226946




GIS中的坐标系定义与转换
  自 "Mapinfo上的GIS系统开发"一文在计算机世界网上刊登后,有好几位网友向我询问坐标系定义与转换方面的问题,问题可归结为 (1) 地图在Mapinfo上显示得很好,但在MapX中却显示不出来或显示得不对;(2) GPS定位得到的WGS84坐标怎么往北京54坐标地图上转。这些问题也是曾经困惑我的问题,在此我谈谈我个人的一些认识及经验,供各位参考,也希望相关方面的专业人士能给予纠正及补充。
1. 椭球体、基准面及地图投影
  GIS中的坐标系定义是GIS系统的基础,正确定义GIS系统的坐标系非常重要。GIS中的坐标系定义由基准面和地图投影两组参数确定,而基准面的定义则由特定椭球体及其对应的转换参数确定,因此欲正确定义GIS系统坐标系,首先必须弄清地球椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念及它们之间的关系。
  基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的基准面,我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国新的大地坐标系--西安80坐标系,目前大地测量基本上仍以北京54坐标系作为参照,北京54与西安80坐标之间的转换可查阅国家测绘局公布的对照表。 WGS1984基准面采用WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心,目前GPS测量数据多以WGS1984为基准。
  上述3个椭球体参数如下:
  椭球体   Mapinfo中代号 年代  长半轴     短半轴     1/扁率
  Krassovsky    3      1940   6378245     6356863     298.3
  IAG 75       31     1975  6378140     6356755     298.25722101
  WGS 84      28     1984   6378137.000  6356752.314 298.257223563
  椭球体与基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面,如前苏联的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的基准面定义肯定是不对的。因此,当系统精度要求较高时,一定要对所采用的参数进行检测、验证,确保坐标系定义的正确性。
3. GIS中地图投影的定义
  我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),又叫横轴墨卡托投影(Transverse Mercator);小于50万的地形图采用正轴等角割园锥投影,又叫兰勃特投影(Lambert Conformal Conic);海上小于50万的地形图多用正轴等角园柱投影,又叫墨卡托投影(Mercator),我国的GIS系统中应该采用与我国基本比例尺地形图系列一致的地图投影系统。 
  在MapX中坐标系定义由基准面、投影两部分参数组成,方法如下:



CoordSys.Set(Type,[Datum],[Units],[OriginLongitude],[OriginLatitude],[StandardParallelOne],[StandardParallelTwo],[Azimuth],[ScaleFactor],[FalseEasting],[FalseNorthing],[Range],[Bounds], [AffineTransform]) 


其中参数:Type表示投影类型,Type为1时地图坐标以经纬度表示,它是必选参数,它后面的参数都为可选参数;
     Datum为大地基准面对象,如果采用非地球坐标(NonEarth)无需定义该参数;
     Units为坐标单位,如Units为7表示以米为单位;
     OriginLongitude、OriginLatitude分别为原点经度和纬度;
     StandardParallelOne、StandardParallelTwo为第一、第二标准纬线;
     Azimuth为方位角,斜轴投影需要定义该参数;
     ScaleFactor为比例系数;
     FalseEasting, FalseNorthing为东伪偏移、北伪偏移值;
     Range为地图可见纬度范围;
     Bounds为地图坐标范围,是一矩形对象,非地球坐标(NonEarth)必须定义该参数;   
     AffineTransform为坐标系变换对象。
   相应高斯-克吕格投影、兰勃特投影、墨卡托投影需要定义的坐标系参数序列如下:
   高斯-克吕格:投影代号(Type),基准面(Datum),单位(Unit),
         中央经度(OriginLongitude),原点纬度(OriginLatitude),
         比例系数(ScaleFactor),
         东伪偏移(FalseEasting),北纬偏移(FalseNorthing)
   兰勃特: 投影代号(Type),基准面(Datum),单位(Unit),
       中央经度(OriginLongitude),原点纬度(OriginLatitude),
       标准纬度1(StandardParallelOne),标准纬度2(StandardParallelTwo),
       东伪偏移(FalseEasting),北纬偏移(FalseNorthing)
   墨卡托: 投影代号(Type),基准面(Datum),单位(Unit),
       原点经度(OriginLongitude),原点纬度(OriginLatitude),
       标准纬度(StandardParallelOne)
   在城市GIS系统中均采用6度或3度分带的高斯-克吕格投影,因为一般城建坐标采用的是6度或3度分带的高斯-克吕格投影坐标。高斯-克吕格投影以6度或3度分带,每一个分带构成一个独立的平面直角坐标网,投影带中央经线投影后的直线为X轴(纵轴,纬度方向),赤道投影后为Y轴(横轴,经度方向),为了防止经度方向的坐标出现负值,规定每带的中央经线西移500公里,即东伪偏移值为500公里,由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,因此规定在横轴坐标前加上带号,如(4231898,21655933)其中21即为带号,同样所定义的东伪偏移值也需要加上带号,如21带的东伪偏移值为21500000米。
  假如你的工作区位于21带,即经度在120度至126度范围,该带的中央经度为123度,采用Pulkovo 1942基准面,那么定义6度分带的高斯-克吕格投影坐标系参数为:(8,1001,7,123,0,1,21500000,0)。
  那么当精度要求较高,实测数据为WGS1984坐标数据时,欲转换到北京54基准面的高斯-克吕格投影坐标,如何定义坐标系参数呢?你可选择WGS 1984(Mapinfo中代号104)作为基准面,当只有一个已知控制点时(见第2部分),根据平移参数调整东伪偏移、北纬偏移值实现WGS84到北京54的转换,如 8,104,7,123,0,1,21500200,-200),也可利用 AffineTransform坐标系变换对象,此时的转换系数(A、B、C、D、E、F)中A、B、D、E为0,只有X、Y方向的平移值C、F ;当有3个已知控制点时,可利用得到的转换系数(A、B、C、D、E、F)定义 AffineTransform坐标系变换对象,实现坐标系的转换,如:(8,104,7,123,0,1,21500000,0,map.AffineTransform),其中AffineTransform定义为AffineTransform.set(7,A、B、C、D、E、F)(7表示单位米);当然有足够多已知控制点时,直接求定7参数自定义基准面就行了。




台湾使用的是横墨卡脱投影,大陆使用的是高斯克吕格投影。

[ 本帖最后由 御天六龙 于 2006-8-7 03:28 PM 编辑 ]
御天六龙 (2006-8-07 11:37:58)
GIS中的坐标系定义与转换      
  作者:   青岛海洋地质研究所   戴勤奋      
        
  自“Mapinfo上的GIS系统开发”一文在计算机世界网上刊登后,有好几位读者向我询问坐标系定义与转换方面的问题,问题可归结为   (1)   地图在Mapinfo上显示得很好,但在MapX中却显示不出来或显示得不对;(2)   GPS定位得到的WGS84坐标怎么往北京54坐标地图上转。这些问题也是曾经困惑我的问题,在此我谈谈我个人的一些认识及经验,供各位读者参考,也希望相关方面的专业人士能给予纠正及补充。     
      
  1.   椭球体、基准面及地图投影     
      
  GIS中的坐标系定义是GIS系统的基础,正确定义GIS系统的坐标系非常重要。GIS中的坐标系定义由基准面和地图投影两组参数确定,而基准面的定义则由特定椭球体及其对应的转换参数确定,因此欲正确定义GIS系统坐标系,首先必须弄清地球椭球体(Ellipsoid)、大地基准面(Datum)及地图投影(Projection)三者的基本概念及它们之间的关系。     
      
  基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的基准面,我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国新的大地坐标系--西安80坐标系,目前大地测量基本上仍以北京54坐标系作为参照,北京54与西安80坐标之间的转换可查阅国家测绘局公布的对照表。   WGS1984基准面采用WGS84椭球体,它是一地心坐标系,即以地心作为椭球体中心,目前GPS测量数据多以WGS1984为基准。     
      
  上述3个椭球体参数如下:     
        
  椭球体                       Mapinfo中代号           年代         长半轴                     短半轴             1/扁率   
  Krassovsky                 3                                 1940             6378245                 6356863             298.3   
  IAG   75                       31                               1975             6378140                 6356755               298.25722101   
  WGS84                       28                               1984             6378137.000           6356752.314       298.257223563   
   
   
   
  椭球体与基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面,如前苏联的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的基准面定义肯定是不对的。因此,当系统精度要求较高时,一定要对所采用的参数进行检测、验证,确保坐标系定义的正确性。

3.   GIS中地图投影的定义     
      
  我国的基本比例尺地形图(1:5千,1:1万,1:2.5万,1:5万,1:10万,1:25万,1:50万,1:100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),又叫横轴墨卡托投影(Transverse   Mercator);小于50万的地形图采用正轴等角割园锥投影,又叫兰勃特投影(Lambert   Conformal   Conic);海上小于50万的地形图多用正轴等角园柱投影,又叫墨卡托投影(Mercator),我国的GIS系统中应该采用与我国基本比例尺地形图系列一致的地图投影系统。      
      
  在MapX中坐标系定义由基准面、投影两部分参数组成,方法如下:     
      
  CoordSys.Set(Type,   [Datum],   [Units],   [OriginLongitude],   [OriginLatitude],      
  [StandardParallelOne],   [StandardParallelTwo],   [Azimuth],   [ScaleFactor],      
  [FalseEasting],   [FalseNorthing],   [Range],   [Bounds],   [AffineTransform])      
      
  其中参数:Type表示投影类型,Type为1时地图坐标以经纬度表示,它是必选参数,它后面的参数都为可选参数;     
  Datum为大地基准面对象,如果采用非地球坐标(NonEarth)无需定义该参数;     
  Units为坐标单位,如Units为7表示以米为单位;     
  OriginLongitude、OriginLatitude分别为原点经度和纬度;     
  StandardParallelOne、StandardParallelTwo为第一、第二标准纬线;     
  Azimuth为方位角,斜轴投影需要定义该参数;     
  ScaleFactor为比例系数;     
  FalseEasting,   FalseNorthing为东伪偏移、北伪偏移值;     
  Range为地图可见纬度范围;     
  Bounds为地图坐标范围,是一矩形对象,非地球坐标(NonEarth)必须定义该参数;      
  AffineTransform为坐标系变换对象。     
      
  相应高斯-克吕格投影、兰勃特投影、墨卡托投影需要定义的坐标系参数序列如下:     
      
  高斯-克吕格:投影代号(Type),基准面(Datum),单位(Unit),     
  中央经度(OriginLongitude),原点纬度(OriginLatitude),     
  比例系数(ScaleFactor),     
  东伪偏移(FalseEasting),北纬偏移(FalseNorthing)     
      
  兰勃特:   投影代号(Type),基准面(Datum),单位(Unit),     
  中央经度(OriginLongitude),原点纬度(OriginLatitude),     
  标准纬度1(StandardParallelOne),标准纬度2(StandardParallelTwo),     
  东伪偏移(FalseEasting),北纬偏移(FalseNorthing)     
      
  墨卡托:   投影代号(Type),基准面(Datum),单位(Unit),     
  原点经度(OriginLongitude),原点纬度(OriginLatitude),     
  标准纬度(StandardParallelOne)     
      
  在城市GIS系统中均采用6度或3度分带的高斯-克吕格投影,因为一般城建坐标采用的是6度或3度分带的高斯-克吕格投影坐标。高斯-克吕格投影以6度或3度分带,每一个分带构成一个独立的平面直角坐标网,投影带中央经线投影后的直线为X轴(纵轴,纬度方向),赤道投影后为Y轴(横轴,经度方向),为了防止经度方向的坐标出现负值,规定每带的中央经线西移500公里,即东伪偏移值为500公里,由于高斯-克吕格投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,因此规定在横轴坐标前加上带号,如(4231898,21655933)其中21即为带号,同样所定义的东伪偏移值也需要加上带号,如21带的东伪偏移值为21500000米。     
      
  假如你的工作区位于21带,即经度在120度至126度范围,该带的中央经度为123度,采用Pulkovo   1942基准面,那么定义6度分带的高斯-克吕格投影坐标系参数为:(8,1001,7,123,0,1,21500000,0)。     
      
  那么当精度要求较高,实测数据为WGS1984坐标数据时,欲转换到北京54基准面的高斯-克吕格投影坐标,如何定义坐标系参数呢?你可选择WGS   1984(Mapinfo中代号104)作为基准面,当只有一个已知控制点时(见第2部分),根据平移参数调整东伪偏移、北纬偏移值实现WGS84到北京54的转换,如 8,104,7,123,0,1,21500200,-200),也可利用   AffineTransform坐标系变换对象,此时的转换系数(A、B、C、D、E、F)中A、B、D、E为0,只有X、Y方向的平移值C、F   ;当有3个已知控制点时,可利用得到的转换系数(A、B、C、D、E、F)定义   AffineTransform坐标系变换对象,实现坐标系的转换,如:(8,104,7,123,0,1,21500000,0,map.AffineTransform),其中AffineTransform定义为AffineTransform.set(7,A、B、C、D、E、F)(7表示单位米);当然有足够多已知控制点时,直接求定7参数自定义基准面就行了。
御天六龙 (2006-8-08 08:38:26)
常用地图投影算法- -
                          

转自 www.gissky.net

作者: 青岛海洋地质研究所 戴勤奋


最近几乎天天都有Email跟我要这样、那样的坐标系转换或投影转换公式,或问我编的投影程序公式是哪来的,有没有专门介绍投影公式的书等等,让我越来越觉得有必要就此方面写点东西,一来我自己总结一下,二来对那些我没有回Email的同行也有个交代,因为那些公式实在太难敲了。我在“海洋地质制图常用地图投影系列小程序”( http://www.gissky.net)中用的公式来自我原来的积累,同时参考了POSC( http://www.posc.org ,国际石油技术软件开放公司)的文献“Coordinate Conversions and Transformation including Formulas”,该文献由EPSG( http://www.epsg.org ,欧洲石油勘探组)编写,比较全面地介绍了各种地图投影与坐标系的转换方法及计算公式,而且最新更新到了2004年,是我目前看到的最全面、最新的相关文档了,只不过是英文的,我正在打算将它们翻成中文,到时与大家共享。
投影计算公式往往表达方式不止一种,有时很难分辨谁对谁错,我只把“墨卡托投影”、“高斯-克吕格投影”、“UTM投影”、“兰勃特等角投影” (1:100万地形图规范中称作正轴等角圆锥投影,GB/T 14512-93)的正反转换公式列出,因为我基本能保证这些公式的正确性。
“海洋地质制图常用地图投影系列小程序”( http://www.gissky.net)已升级,原下载者请注意下载更新版本。

1. 约定

本文中所列的转换公式都基于椭球体
a -- 椭球体长半轴
b -- 椭球体短半轴
f -- 扁率
e -- 第一偏心率
e’ -- 第二偏心率
N -- 卯酉圈曲率半径
R -- 子午圈曲率半径

B -- 纬度,L -- 经度,单位弧度(RAD)
-- 纵直角坐标, -- 横直角坐标,单位米(M)

2. 椭球体参数
我国常用的3个椭球体参数如下(源自“全球定位系统测量规范 GB/T 18314-2001”):
椭球体 长半轴 a(米) 短半轴b(米)
Krassovsky (北京54采用) 6378245 6356863.0188
IAG 75(西安80采用) 6378140 6356755.2882
WGS 84 6378137 6356752.3142


需要说明的是,在“海洋地质制图常用地图投影系列小程序”中,程序界面上的所谓“北京1954“西安1980”及“WGS 84”在实际计算中只涉及了相应的椭球体参数。

3. 墨卡托(Mercator)投影
3.1 墨卡托投影简介
墨卡托(Mercator)投影,是一种"等角正切圆柱投影”,荷兰地图学家墨卡托(Gerhardus Mercator 1512-1594)在1569年拟定, 假设地球被围在一中空的圆柱里,其标准纬线与圆柱相切接触,然后再假想地球中心有一盏灯,把球面上的图形投影到圆柱体上,再把圆柱体展开,这就是一幅选定标准纬线上的“墨卡托投影”绘制出的地图。
墨卡托投影没有角度变形,由每一点向各方向的长度比相等,它的经纬线都是平行直线,且相交成直角,经线间隔相等,纬线间隔从标准纬线向两极逐渐增大。墨卡托投影的地图上长度和面积变形明显,但标准纬线无变形,从标准纬线向两极变形逐渐增大,但因为它具有各个方向均等扩大的特性,保持了方向和相互位置关系的正确。
在地图上保持方向和角度的正确是墨卡托投影的优点,墨卡托投影地图常用作航海图和航空图,如果循着墨卡托投影图上两点间的直线航行,方向不变可以一直到达目的地,因此它对船舰在航行中定位、确定航向都具有有利条件,给航海者带来很大方便。
“海底地形图编绘规范”(GB/T 17834-1999,海军航保部起草)中规定1:25万及更小比例尺的海图采用墨卡托投影,其中基本比例尺海底地形图(1:5万,1:25万,1: 100万)采用统一基准纬线30°,非基本比例尺图以制图区域中纬为基准纬线。基准纬线取至整度或整分。
3.2 墨卡托投影坐标系
取零子午线或自定义原点经线(L0)与赤道交点的投影为原点,零子午线或自定义原点经线的投影为纵坐标X轴,赤道的投影为横坐标Y轴,构成墨卡托平面直角坐标系。
3.3 墨卡托投影正反解公式
墨卡托投影正解公式:(B,L)→(X,Y),标准纬度B0,原点纬度 0,原点经度L0

墨卡托投影反解公式:(X,Y) →(B,L),标准纬度B0,原点纬度 0,原点经度L0     

公式中EXP为自然对数底,纬度B通过迭代计算很快就收敛了。
4. 高斯-克吕格(Gauss-Kruger)投影和UTM(Universal Transverse Mercator)投影
4.1 高斯-克吕格投影与UTM投影异同
高斯-克吕格(Gauss-Kruger)投影与UTM投影(Universal Transverse Mercator,通用横轴墨卡托投影)都是横轴墨卡托投影的变种,目前一些国外的软件或国外进口仪器的配套软件往往不支持高斯-克吕格投影,但支持 UTM投影,因此常有把UTM投影当作高斯-克吕格投影的现象。从投影几何方式看,高斯-克吕格投影是“等角横切圆柱投影”,投影后中央经线保持长度不变,即比例系数为1;UTM投影是“等角横轴割圆柱投影”,圆柱割地球于南纬80度、北纬84度两条等高圈,投影后两条割线上没有变形,中央经线上长度比 0.9996。从计算结果看,两者主要差别在比例因子上,高斯-克吕格投影中央经线上的比例系数为1, UTM投影为0.9996,高斯-克吕格投影与UTM投影可近似采用 X[UTM]=0.9996 * X[高斯],Y[UTM]=0.9996 * Y[高斯],进行坐标转换(注意:如坐标纵轴西移了500000米,转换时必须将Y值减去500000乘上比例因子后再加500000)。从分带方式看,两者的分带起点不同,高斯-克吕格投影自0度子午线起每隔经差6度自西向东分带,第1带的中央经度为3°;UTM投影自西经180°起每隔经差6度自西向东分带,第1带的中央经度为-177°,因此高斯-克吕格投影的第1带是UTM的第31带。此外,两投影的东伪偏移都是500公里,高斯-克吕格投影北伪偏移为零,UTM北半球投影北伪偏移为零,南半球则为10000公里。
4.2 高斯-克吕格投影简介
高斯-克吕格(Gauss-Kruger)投影,是一种“等角横切圆柱投影”。德国数学家、物理学家、天文学家高斯(Carl Friedrich Gauss,1777一 1855)于十九世纪二十年代拟定,后经德国大地测量学家克吕格(Johannes Kruger,1857~1928)于 1912年对投影公式加以补充,故名。设想用一个圆柱横切于球面上投影带的中央经线,按照投影带中央经线投影为直线且长度不变和赤道投影为直线的条件,将中央经线两侧一定经差范围内的球面正形投影于圆柱面。然后将圆柱面沿过南北极的母线剪开展平,即获高斯一克吕格投影平面。
高斯一克吕格投影后,除中央经线和赤道为直线外,其他经线均为对称于中央经线的曲线。高斯-克吕格投影没有角度变形,在长度和面积上变形也很小,中央经线无变形,自中央经线向投影带边缘,变形逐渐增加,变形最大处在投影带内赤道的两端。由于其投影精度高,变形小,而且计算简便(各投影带坐标一致,只要算出一个带的数据,其他各带都能应用),因此在大比例尺地形图中应用,可以满足军事上各种需要,并能在图上进行精确的量测计算。
按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差,又要使带数不致过多以减少换带计算工作,据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。通常按经差6度或3度分为六度带或三度带。六度带自 0度子午线起每隔经差6度自西向东分带,带号依次编为第 1、2…60带。三度带是在六度带的基础上分成的,它的中央子午线与六度带的中央子午线和分带子午线重合,即自 1.5度子午线起每隔经差3度自西向东分带,带号依次编为三度带第 1、2…120带。我国的经度范围西起 73°东至135°,可分成六度带十一个,各带中央经线依次为75°、81°、87°、……、117°、123°、129°、135°,或三度带二十二个。
我国大于等于50万的大中比例尺地形图多采用六度带高斯-克吕格投影,三度带高斯-克吕格投影多用于大比例尺测图,如城建坐标多采用三度带的高斯-克吕格投影。
4.3 UTM投影简介
UTM投影全称为“通用横轴墨卡托投影”,是一种“等角横轴割圆柱投影”,椭圆柱割地球于南纬80度、北纬84度两条等高圈,投影后两条相割的经线上没有变形,而中央经线上长度比0.9996。UTM投影是为了全球战争需要创建的,美国于1948年完成这种通用投影系统的计算。与高斯-克吕格投影相似,该投影角度没有变形,中央经线为直线,且为投影的对称轴,中央经线的比例因子取0.9996是为了保证离中央经线左右约330km处有两条不失真的标准经线。
UTM投影分带方法与高斯-克吕格投影相似,是自西经180°起每隔经差6度自西向东分带,将地球划分为60个投影带。
我国的卫星影像资料常采用UTM投影。
4.4 高斯-克吕格投影与UTM投影坐标系
  高斯- 克吕格投影与UTM投影是按分带方法各自进行投影,故各带坐标成独立系统。以中央经线(L0)投影为纵轴X, 赤道投影为横轴Y,两轴交点即为各带的坐标原点。为了避免横坐标出现负值,高斯- 克吕格投影与UTM北半球投影中规定将坐标纵轴西移500公里当作起始轴,而UTM南半球投影除了将纵轴西移500公里外,横轴南移10000公里。由于高斯-克吕格投影与UTM投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,为了区别某一坐标系统属于哪一带,通常在横轴坐标前加上带号,如(4231898m,21655933m),其中21即为带号。   
4.5 高斯-克吕格投影与UTM投影正反解公式
高斯-克吕格投影和UTM投影公式从目前公开出版的教材、文献及网上我看到好几种版本,可归结为下列两组,我把原来教科书及国内文献上常见的一套公式列作高斯-克吕格投影公式,POSC(国际石油技术软件开放公司)及国外文献上见到的另一套公式列作UTM投影公式。常常能看到两套投影公式混用的文献资料,文中谈论的是UTM投影,但列出的公式却是国内教材上的高斯-克吕格投影公式,让我很困惑。为此,我设定比例因子都为1,用下列两组公式分别进行了同点的投影计算,计算结果在中高纬度时两套公式差异很小,小数后6位都是一致的;在低纬度时,投影结果差异拉大,横轴在小数第三位开始出现差异。假如精确到厘米级,上述试验说明两套公式混用应该没问题。不过,有可能会有其它极端的情况,毕竟是不同的投影公式。
高斯-克吕格投影正解公式:(B,L)→(X,Y),原点纬度 0,中央经度L0

上面公式中东纬偏移FE = 500000米 + 带号 * 1000000;
高斯-克吕格投影比例因子k0 = 1

UTM投影正解公式:(B,L)→(X,Y),原点纬度 0,中央经度L0

上面公式中东纬偏移 FE= 500000米 ;北纬偏移 FN北半球= 0,FN南半球= 10000000米;
UTM投影比例因子k0 = 0.9996,其它参数同高斯-克吕格投影正解公式

高斯-克吕格投影反解公式:(X,Y) →(B,L),原点纬度 0,中央经度L0     

UTM投影反解公式:(X,Y) →(B,L),原点纬度 0,中央经度L0
      
式中参数同高斯-克吕格投影反解公式

5. 兰勃特等角投影(Lambert Conformal Conic);
5.1 兰勃特等角投影简介
兰勃特等角投影,在双标准纬线下是一“等角正轴割圆锥投影”,由德国数学家兰勃特(J.H.Lambert)在1772年拟定。设想用一个正圆锥割于球面两标准纬线,应用等角条件将地球面投影到圆锥面上,然后沿一母线展开,即为兰勃特投影平面。兰勃特等角投影后纬线为同心圆弧,经线为同心圆半径。前面已经介绍的墨卡托(Mercator)投影是它的一个极端特例。
兰勃特投影采用双标准纬线相割,与采用单标准纬线相切比较,其投影变形小而均匀,兰勃托投影的变形分布规律是:a) 角度没有变形;b) 两条标准纬线上没有任何变形;c) 等变形线和纬线一致,即同一条纬线上的变形处处相等; d) 在同一经线上,两标准纬线外侧为正变形(长度比大于1),而两标准纬线之间为负变形(长度比小于1)。变形比较均匀,变形绝对值也比较小;e) 同一纬线上等经差的线段长度相等,两条纬线间的经纬线长度处处相等。
兰勃特投影常用于小比例尺地形图。“1:1000000地形图编绘规范及图式 GB/T 14515-93”中规定1:100万地形图采用正轴等角圆锥投影(兰勃特等角投影),并采用了国际地理学会规定的全球统一使用的国际百万分之一地图的分幅原则,按纬差4°从赤道向北、经差6°从-180°向东分幅,每个投影分幅单独计算坐标,每幅两条标准纬线,第一标准纬线为图幅南端纬度加30′的纬线,第二标准纬线为图幅北端纬度减30′的纬线。由于是纬差4°分带投影的,所以当沿着纬线方向拼接地图时,不论多少图幅,均不会产生裂隙;但是,当沿着经线方向拼接时,因拼接线分别处于上下不同的投影带,投影后的曲率不同,致使拼接时会产生裂隙。   
5.2 兰勃特等角投影坐标系
  以图幅的原点经线(一般是中央经线L0)作纵坐标X轴,原点经线与原点纬线(一般是最南端纬线)的交点作为原点,过此点的切线作为横坐标Y轴,构成兰勃特平面直角坐标系
5.3 兰勃特等角投影正反解公式
兰勃特等角投影正解公式:
(B,L)→(X,Y),原点纬度 B0,原点经度L0,第一标准纬线B1,第二标准纬线B2:
                        
兰勃特等角投影反解公式:
(X,Y) →(B,L),原点纬度 B0,原点经度L0,第一标准纬线B1,第二标准纬线B2:

式中参数同兰勃特等角投影正解公式
B通过迭代获取



一些公式图片没贴上来,请看这里
http://bornkilled.bokee.com/2084328.html
御天六龙 (2006-8-08 09:09:39)
一、只谈比较常用的几种:“墨卡托投影”、“高斯-克吕格投影”、“UTM投影”、“兰勃特等角投影”


1. 墨卡托(Mercator)投影

1.1  墨卡托投影简介

        墨卡托(Mercator)投影,是一种"等角正切圆柱投影”,荷兰地图学家墨卡托(Gerhardus Mercator 1512-1594)在1569年拟定, 假设地球被围在一中空的圆柱里,其标准纬线与圆柱相切接触,然后再假想地球中心有一盏灯,把球面上的图形投影到圆柱体上,再把圆柱体展开,这就是一幅选定标准纬线上的“墨卡托投影”绘制出的地图。
        墨 卡托投影没有角度变形,由每一点向各方向的长度比相等,它的经纬线都是平行直线,且相交成直角,经线间隔相等,纬线间隔从标准纬线向两极逐渐增大。墨卡托 投影的地图上长度和面积变形明显,但标准纬线无变形,从标准纬线向两极变形逐渐增大,但因为它具有各个方向均等扩大的特性,保持了方向和相互位置关系的正确。

        在地图上保持方向和角度的正确是墨卡托投影的优点,墨卡托投影地图常用作航海图和航空图,如果循着墨卡托投影图上两点间的直线航行,方向不变可以一直到达目的地,因此它对船舰在航行中定位、确定航向都具有有利条件,给航海者带来很大方便。

        “海底地形图编绘规范”(GB/T 17834-1999,海军航保部起草)中规定1:25万及更小比例尺的海图采用墨卡托投影,其中基本比例尺海底地形图(1:5万,1:25万,1:100万)采用统一基准纬线30°,非基本比例尺图以制图区域中纬为基准纬线。基准纬线取至整度或整分。

1.2  墨卡托投影坐标系

取零子午线或自定义原点经线(L0)与赤道交点的投影为原点,零子午线或自定义原点经线的投影为纵坐标X轴,赤道的投影为横坐标Y轴,构成墨卡托平面直角坐标系。

2. 高斯-克吕格(Gauss-Kruger)投影和UTM(Universal Transverse Mercator)投影

2.1  高斯-克吕格投影简介

高斯-克吕格(Gauss-Kruger)投影,是一种“等角横切圆柱投影”。德国数学家、物理学家、天文学家高斯(Carl Friedrich Gauss,1777一 1855)于十九世纪二十年代拟定,后经德国大地测量学家克吕格(Johannes Kruger,1857~1928)于 1912年对投影公式加以补充,故名。设想用一个圆柱横切于球面上投影带的中央经线,按照投影带中央经线投影为直线且长度不变和赤道投影为直线的条件,将中央经线两侧一定经差范围内的球面正形投影于圆柱面。然后将圆柱面沿过南北极的母线剪开展平,即获高斯一克吕格投影平面。

高斯一克吕格投影后,除中央经线和赤道为直线外,其他经线均为对称于中央经线的曲线。高斯-克吕格投影没有角度变形,在长度和面积上变形也很小,中央经线无变形,自中央经线向投影带边缘,变形逐渐增加,变形最大处在投影带内赤道的两端。由于其投影精度高,变形小,而且计算简便(各投影带坐标一致,只要算出一个带的数据,其他各带都能应用),因此在大比例尺地形图中应用,可以满足军事上各种需要,并能在图上进行精确的量测计算。

按一定经差将地球椭球面划分成若干投影带,这是高斯投影中限制长度变形的最有效方法。分带时既要控制长度变形使其不大于测图误差,又要使带数不致过多以减少换带计算工作,据此原则将地球椭球面沿子午线划分成经差相等的瓜瓣形地带,以便分带投影。通常按经差6度或3度分为六度带或三度带。六度带自0度子午线起每隔经差6度自西向东分带,带号依次编为第 1、2…60带。三度带是在六度带的基础上分成的,它的中央子午线与六度带的中央子午线和分带子午线重合,即自 1.5度子午线起每隔经差3度自西向东分带,带号依次编为三度带第 1、2…120带。我国的经度范围西起 73°东至135°,可分成六度带十一个,各带中央经线依次为75°、81°、87°、……、117°、123°、129°、135°,或三度带二十二个。

我国大于等于50万的大中比例尺地形图多采用六度带高斯-克吕格投影,三度带高斯-克吕格投影多用于大比例尺测图,如城建坐标多采用三度带的高斯-克吕格投影。

2.2  UTM投影简介

UTM投影全称为“通用横轴墨卡托投影”,是一种“等角横轴割圆柱投影”,椭圆柱割地球于南纬80度、北纬84度两条等高圈,投影后两条相割的经线上没有变形,而中央经线上长度比0.9996。UTM投影是为了全球战争需要创建的,美国于1948年完成这种通用投影系统的计算。与高斯-克吕格投影相似,该投影角度没有变形,中央经线为直线,且为投影的对称轴,中央经线的比例因子取0.9996是为了保证离中央经线左右约330km处有两条不失真的标准经线。

UTM投影分带方法与高斯-克吕格投影相似,是自西经180°起每隔经差6度自西向东分带,将地球划分为60个投影带。

我国的卫星影像资料常采用UTM投影。



2.3  高斯-克吕格投影与UTM投影异同

高斯-克吕格(Gauss-Kruger)投影与UTM投影(Universal Transverse Mercator,通用横轴墨卡托投影)都是横轴墨卡托投影的变种,目前一些国外的软件或国外进口仪器的配套软件往往不支持高斯-克吕格投影,但支持UTM投影,因此常有把UTM投影当作高斯-克吕格投影的现象。从投影几何方式看,高斯-克吕格投影是“等角横切圆柱投影”,投影后中央经线保持长度不变,即比例系数为1;UTM投影是“等角横轴割圆柱投影”,圆柱割地球于南纬80度、北纬84度两条等高圈,投影后两条割线上没有变形,中央经线上长度比0.9996。从计算结果看,两者主要差别在比例因子上,高斯-克吕格投影中央经线上的比例系数为1, UTM投影为0.9996,高斯-克吕格投影与UTM投影可近似采用 X[UTM]=0.9996 * X[高斯],Y[UTM]=0.9996 * Y[高斯],进行坐标转换(注意:如坐标纵轴西移了500000米,转换时必须将Y值减去500000乘上比例因子后再加500000)。从分带方式看,两者的分带起点不同,高斯-克吕格投影自0度子午线起每隔经差6度自西向东分带,第1带的中央经度为3°;UTM投影自西经180°起每隔经差6度自西向东分带,第1带的中央经度为-177°,因此高斯-克吕格投影的第1带是UTM的第31带。此外,两投影的东伪偏移都是500公里,高斯-克吕格投影北伪偏移为零,UTM北半球投影北伪偏移为零,南半球则为10000公里。

2.4  高斯-克吕格投影与UTM投影坐标系

     高斯- 克吕格投影与UTM投影是按分带方法各自进行投影,故各带坐标成独立系统。以中央经线(L0)投影为纵轴X,赤道投影为横轴Y,两轴交点即为各带的坐标原点。为了避免横坐标出现负值,高斯- 克吕格投影与UTM北半球投影中规定将坐标纵轴西移500公里当作起始轴,而UTM南半球投影除了将纵轴西移500公里外,横轴南移10000公里。由于高斯-克吕格投影与UTM投影每一个投影带的坐标都是对本带坐标原点的相对值,所以各带的坐标完全相同,为了区别某一坐标系统属于哪一带,通常在横轴坐标前加上带号,如(4231898m,21655933m),其中21即为带号。

二、分带方法

1.我国采用6度分带和3度分带:

  1∶2.5万及1∶5万的地形图采用6度分带投影,即经差为6度,从零度子午线开始,自西向东每个经差6度为一投影带,全球共分60个带,用1,2,3,4,5,……表示.即东经0~6度为第一带,其中央经线的经度为东经3度,东经6~12度为第二带,其中央经线的经度为9度。

  1∶1万的地形图采用3度分带,从东经1.5度的经线开始,每隔3度为一带,用1,2,3,……表示,全球共划分120个投影带,即东经1.5~ 4.5度为第1带,其中央经线的经度为东经3度,东经4.5~7.5度为第2带,其中央经线的经度为东经6度.我省位于东经113度-东经120度之间,跨第38、39、40共计3个带,其中东经115.5度以西为第38带,其中央经线为东经114度;东经115.5~118.5度为39带,其中央经线为东经117度;东经118.5度以东到山海关为40带,其中央经线为东经120度。

地形图上公里网横坐标前2位就是带号,例如:1∶5万地形图上的横坐标为20345486,其中20即为带号,345486为横坐标值。

2.当地中央经线经度的计算

六度带中央经线经度的计算:当地中央经线经度=6°×当地带号-3°,例如:地形图上的横坐标为20345,其所处的六度带的中央经线经度为:6°×20-3°=117°(适用于1∶2.5万和1∶5万地形图)。

三度带中央经线经度的计算:中央经线经度=3°×当地带号(适用于1∶1万地形图)。




Trackback: http://tb.blog.csdn.net/TrackBack.aspx?PostId=975558
御天六龙 (2006-8-08 09:12:42)
格林登投影
        格林登投影属任意投影。投影的条件是:首先使它角度和面积误差均不很大,使在赤道上的投影没有误差,并使地球投影在一正圆形内,在投影上经线和纬线都是圆 弧,赤道和中央经线为直线,且互相垂直,为大圆的直径。径线的间隔相等,纬线间隔自赤道向两极逐渐增大。赤道附近地区变形较小,越向高纬地区变形越大。因 此在绘制世界地图时,常采用矩形图廓,图内不表示两极。

等差分纬线多圆锥投影
       等差分纬线多圆锥投影是任意多圆锥投影的一种。是我国制图工作者根据我国领土的形状和位置,于1963年设计的。该投影在我国编制各种比例尺世界政区图及 其他类型世界地图中已得到较广泛的使用,并获得较好效果。该投影中纬线为对称于赤道的同轴圆圆弧,圆心位于中央经线上。中央经线为一直线,其它经线为对称 于中央经线的曲线,且离中央经线愈远,其经线间隔愈成比例地递减;极点表示为圆弧,其长度为赤道投影长度的二分之一,经纬网的图形有球形感。我国被配置在 地图中接近于中央的位置,使我国面积相对于同一条纬带上其它国家的面积不因面积变形而有所缩小,图形形状比较正确,图面图形完整,没有裂隙,也不出现重 复,保持太平洋完整,可显示我国与邻近国家的水陆联系。由于该投影的性质是接近等面积的任意投影,因此我国绝大部分地区面积变形小。在中学使用的世界地图 中的世界的气候类型和洋流图、世界政区图和自然带的分布图都是采用的这种投影图。

古德投影
        古德投影为美国科学家古德拟定,故名。是将伪圆柱投影的非制图区加以断裂,使制图区变形减少,编制成断裂地图的方法。它的作法是在地图上把几个主要制图区 的中央都定为中央经线,使地图分为几个部分,沿赤道连接在一起。不论是桑逊投影、毛尔威特投影以及其他伪圆柱投影都可以断裂。因为凡是纬线平行的投影,断 裂后都有可能相互对照。为了表示大陆的完整性则可在海洋部分断裂,如为了完整地表示海洋,则可在大陆部分断裂。还可采用桑逊投影和毛尔威特投影结合在一起 的断裂方法,在南北纬40°44′11.8°以内,采用桑逊投影以外采用毛尔威特投影。这样可减少变形,但缺点是在40°纬线附近有一折角,离中央经线越 远,折角越显著。

毛尔威特投影
       毛尔威特投影又称椭圆径线等面积伪圆柱投影,由德国人毛尔威特于1805年设计,故名,投影图上中央经线是直线,其他经线为对称于中央经线的椭圆弧。全图 呈椭圆状,在中央经线上从赤道向南向北纬线间距逐渐缩小,在赤道上从中央经线向东向西纬线间距相等。该图面积无变形,而长度和角度均有变形,中央经线和南 北纬40°的两交点为没有变形的点。该投影常用于编制小比例尺世界地图,在国外多用于地图集的封面。

彭纳投影
       彭纳投影即等积伪圆锥投影。为法国人彭纳所创。中央经线是直线,其他经线为对称于中央经线的曲线。纬线为同心圆弧。中央经线和标准纬线上没有变形,离开这两条线越远变形越大。图上所有纬线都保持长度不变,面积相等。彭纳投影常用作大洲图。

等角圆柱投影
       等角圆柱投影指保持角度、形状没有变形的圆柱投影。这是荷兰地图学家墨卡托于1569年创制的,又称墨卡托投影或等角正圆柱投影。该图上经纬线成互相直交 的平行直线,经线的间隔相等,纬线的间隔随纬度增高而加大。赤道处角度、形状没有误差,越向高纬度处误差越大。地面上的等方位角航线投影后为直线,故广泛 用于绘制航海图。但这种投影面积变形显著,在纬度60°地区经线和纬线比都扩大2倍,面积比例比实际扩大了4倍。到纬度80°附近,经线和纬线比例尺都扩 大将近6倍,面积扩大了33倍。所以在墨卡托投影上,纬度80°以上的地区就不绘出来了。中学使用的中国地图册中的时区图和世界地图册中的东南亚地图都是 采用这种投影绘制的。

等距圆锥投影
       等距圆锥投影指沿经线方向长度没有变形的圆锥投影。在这种图上,纬线间距相等,沿经线方向长度没有变形。除经线方向外其它方向的长度都有变形,面积和角度 也变形,但变形都不太大。这种投影适于编制各种教学用图和交通图。但在我国使用比较少。在中学使用的世界地图册中的苏联图是采用这种投影。

等积圆锥投影
       等积圆锥投影 指保持面积不变形的圆锥投影。正轴等面积割圆锥投影又称亚尔勃斯投影。这种投影在实践中应用较广。假设一个圆锥切割地球的某两条纬线为了保持等积条件,在 把地球经纬线投影到圆锥面时,两条标准纬线以内的经线局部比例尺要相应增大,而且在任一点上纬线缩小的比例和经线扩大的比例相等。两条标准纬线以外的纬线 比例尺大于主比例尺,经线比例尺小于主比例尺。在这种地图上纬线间距从中间向南向北逐渐缩小。这种投影因为没有面积变形,所以在图上保持正确的面积对比。 适于作经济地图或某些自然地图。我国行政区图、地形图、人口地图及社会经济等地图多采用这种投影,以25°N和47°N两条纬线为标准纬线相割。中学使用 的中国地图中的中国人口、民族、气候、水系、土地利用、农业、工业和矿产等分布图也都采用这种投影。

等角圆锥投影
        等角圆锥投影 指在地图上没有角度变形的圆锥投影。它是德国数学家兰勃脱所拟定,故又称兰勃脱正形圆锥投影,由于这种投影是一圆锥切割地球的两条标准纬线,又称双标准纬 线等角圆锥投影。很多中纬度国家和地区多采用这种投影来编制中、小比例尺地图。在图上,为了保持等角条件,必须使图上任一点的经线比与纬线比相等。圆锥面 展平后,经线为交于圆心的直线束,但经线之间的夹角小于纬线呈同心圆弧,纬线的间距从中间向南向北逐渐增大。同一纬线上的经线间隔相等。在绘制我国地图 时,多以25°N和47°N作为标准纬线,离开标准纬线愈远,面积变形愈大。我国的分省地图边多采用这种投影。中学使用的地图册中,中国地理的所有分区地 图,以及世界地图中的朝鲜、蒙古、日本、南亚、西亚、南欧、西欧、北欧、中欧、美国、墨西哥及西印度群岛等均用这种投影。

等积方位投影
       等积方位投影是使图上面积和相应的实际地面面积相等的方位投影。因地球面与投影面相切(或相割)的位置不同,分为正轴,横轴、斜轴投影。(1)等积正轴 (方位)投影中的经线表现为放射状直线,纬线表现为同心圆。从投影中心向外,纬线间隔不断缩小。这种投影主要适于绘制极地和南北半球图。如中学生使用的中 国地图册中的北半球和南半球图。(2)等积横轴(方位)投影又称赤道等积方位投影。在这种图上,通过投影中心的中央经线和赤道表现为直线,其他经纬线都表 现为曲线,在中央经线上从中心向南向北,纬线间隔逐渐缩小,在赤道上从地图中心向东向西,经线间隔逐渐缩小。我国所绘东西半球图,多用此投影,在中学生使 用的世界地图册中,东西半球图和非州图。(3)等积斜方位投影中央经线表现为直线,其他经纬线为曲线。在中央经线上从地图中心向上向下,纬线间隔逐渐缩 小。多用在地图集中做大洲图,各大洲面积便于对比。在中学使用的世界地图集中的陆半球和水半球。亚洲图、欧洲图、北美洲图、南美洲图、大洋洲及太平洋岛屿 等图均用此投影图(4)等距方位投影又称波斯托投影。沿一个主方向比例不变,在正投影中,经线不变,在横轴、斜轴投影中,沿垂直圈比例不变。经纬线形式和 等积方位投影相同,只是纬线间隔不同,当纬差相同时,在中央经线上纬线间隔距离相等。正轴投影主要用作极区地图,如我国出版的世界地图集中的北冰洋和南极 洲。

等距投影
        等距投影 是一种任意投影。沿某一特定方向之距离,投影之后保持不变,即沿该特定方向长度之比等于1。在实际应用中多把经线绘成直线,并保持沿经线方向距离相等,面积和角度有些变形,多用于绘制交通图。

等角方位投影
       等角方位投影 方位投影的一种。指保持角度正确的方位投影。因地球面与投影面相切或相割位置不同,分为正轴、横轴、斜轴投影。在正轴投影中,纬线是以极点为圆心的同心 圆,纬线间距从地图中心向外逐渐扩大。经线为由极点向外成放射状直线,经线间的夹角等于经度差。这种投影没有角度变形,但面积变形较大,到投影图的边缘, 面积变形为中心的四倍,在编制南北纬84°以上的地面1∶100万地图时,常采用等角正方位投影。

伪圆锥投影
伪圆锥投影是在圆锥投影的基础上,根据某些条件改变经线形状而成的。这类投影的纬线形状与圆锥投影类似,即纬线为同心圆弧,圆心位于中央经线上,但经线则不同,除中央经线为直线外,其余的经线均为对称于中央经线的曲线。

按投影的变形性质,伪圆锥投影和伪圆柱投影一样,没有等角投影,因为这种投影经纬线不直交。伪圆锥投影只有等积投影和任意投影。最常用的伪圆锥投影是等积伪圆锥投影。

等积伪圆锥投影又称彭纳投影,由法国水利工程师彭纳(Rigobert-Bonne)于1752年首先提出并应用于法国地形图而得名。彭纳投影的中央经线 为直线,其长度比等于1,即m0=1;纬线为同心圆弧,沿纬线长度比等于1,即n=1;图上面积与实际相应的面积相等,即P=1。在一条纬线上的经线间隔 相等,在中央经线上纬线间隔相等,中央经线与所有的纬线正交,中央纬线与所有的经线正交。

这个投影没有面积变形,中央经线和中央纬线是两条没有变形的线,离开这两条线愈远,变形愈大。

彭纳投影主要用于编制小比例尺的大洲图。例如,我国地图出版社出版的《世界地图图2-51彭纳投影及其最大角度变形集》中的亚洲政区图,英国《泰晤士世界地图集》中的澳大利亚与西南太平洋地图,都是采用的彭纳投影。

多圆锥投影
(一)多圆锥投影的概念

在切圆锥投影中,离开标准纬线愈远,变形愈大。如果制图区域包含纬差较大时,则在边缘纬线处将产生相当大的变形。因此,采用双标准纬线圆锥投影比 采用单标准纬线圆锥投影变形要小些。如果有更多的标准纬线,则变形会更小些,多圆锥投影就是由这样的设想建立起来的。假设有许多圆锥与地球面上的纬线相 切,将球面上的经纬线投影于这些圆锥面上,然后沿同一母线方向将圆锥剪开展成平面,如图2-48所示。由于圆锥顶点不是一个,所以纬线投影为同轴圆弧,其 圆心都在中央经线的延长线上,除中央经线为直线外,其余的经线投影为对称于中央经线的曲线。凡是经纬线形式符合上述特征的,均称为多圆锥投影。由于多圆锥 投影的经纬线系弯曲的曲线,具有良好的球形感,所以它常用于编制世界地图。

(二)普通多圆锥投影

普通多圆锥投影除了中央经线和每一条纬线的长度比等于1外,即m0=1,n=1其余经线长度比均大于1。这个投影在中央经线上纬线间隔相等,在每一条纬线 上经线间隔相等。普通多圆锥投影属于任意投影,中央经线是一条没有变形的线,离开中央经线愈远,变形愈大。这个投影适于作南北方向延伸地区的地图。美国海 岸测量局曾用此投影制作美国海岸附近地区的地图。

普通多圆锥投影的另一个用途是绘制地球仪用的图形。把整个地球按一定经差分为若干带,每带中央的经线都投影为直线,各带的投影图在赤道相接,将这样的图形贴于预制的球胎上,就成为一个地球仪。

(三)改良多圆锥投影

改良多圆锥投影是由普通多圆锥改良而成的。过去长时期国际上用它编绘百万分之一分幅地图,这是由1909年伦敦国际地理学会议决定的,故又名国际百万分之一地图投影。

国际百万分之一地图,在纬度0°—60°范围内,按纬差4°、经差6°分幅;在纬度60°—76°范围内,按纬差4°、经差12°分幅;在纬度76°— 88°范围内按纬差4°、经差24°分幅。每幅单独投影。每幅图的南北两条边纬线是同轴圆弧,其圆心位于中央经线的延长线上。将这两条纬线按经差1°等 分,过相应分点连成的直线即为各条经线。其他纬线是4等分各经线后,将相应分点联成的平滑曲线。

这个投影南北两条边纬线长度比等于1,其余纬线长度比均小于1,以中央纬线长度比为最小。在按经差6°的分幅中,距中央经线经差为±2°(在按经差12° 的分幅中,距中央经线经差为±4°,在按经差24°的分幅中,距中央经线经差为±8°)的经线长度比等于1,中间经线长度比小于1,边缘经线长度比大于 1。这种投影按变形性质来说属任意投影。由于每一幅图包括的范围不大,因而变形很小。在我国范围内长度变形不超过0.06%,面积变形不超过0.12%, 角度最大变形不超过5’。故总的来说,这种投影精度还是很高的。但因它不具有等角条件,故现已被等角圆锥投影所取代。

(四)等差分纬线多圆锥投影

这个投影是由我国地图出版社于1963年设计的一种不等分纬线的多圆锥投影。赤道和中央经线是互相垂直的直线,其他纬线为对称于赤道的同轴圆弧,其圆心均 在中央经线的延长线上;其他经线为对称于中央经线的曲线,各经线间的间隔,随离中央经线距离的增大而逐渐缩短,按等差递减。极点为圆弧,其长度为赤道的 1/2。

这种投影的变形性质属任意投影。我国绝大部分地区的面积变形在10%以内,面积比等于1的等变形线自东向西横贯我国中部;中央经线和纬度±44°交点处没 有角度变形,我国境内绝大部分地区的角度最大变形在10°以内,少数地区在13°左右。地图出版社用这一投影编制过数种比例尺的世界政区图和其他类型的世 界地图。

1976年地图出版社又拟定了另外一种不等分纬线多圆锥投影——正切差分纬线多圆锥投影,这个投影的经线间隔,由中央经线向东西两侧按与中央经线经差的正 切函数递减。正切函数随角度增加递增速度越来越快。因此,正切差分纬线多圆锥投影的经线间隔,在中央经线附近变化较小,在远离中央经线的地方,变化较大。 地图出版社1981年出版的1:1400万世界全图采用了这个投影。


Trackback: http://tb.blog.csdn.net/TrackBack.aspx?PostId=975788
0. (2006-8-13 17:21:21)
感谢分享!!!
sican (2006-8-13 22:36:26)
太专业了。。。:(

希望有短小精干、容易明白的。。呵呵~~

楼主辛苦了。。。

我来说两句 查看所有回复

<script language="javascript" src="http://www.100gps.com/include/js/ajax.js" type="text/javascript"></script> id="phpframe" name="phpframe" marginwidth="0" src="about:blank" frameborder="0" width="0" height="0"><script language="javascript" src="http://count22.51yes.com/click.aspx?id=224264587&logo=6" type="text/javascript"></script> 51YES网站统计系统 marginwidth="0" marginheight="0" src="http://count22.51yes.com/sa.aspx?id=224264587&refe=http%3A//www.google.com.sg/search%3Fhl%3Dzh-CN%26newwindow%3D1%26rls%3DGGLG%2CGGLG%3A2005-41%2CGGLG%3Azh-CN%26q%3Dutm%25E5%259D%2590%25E6%25A0%2587%26suggest%3D5%26sa%3DX%26oi%3Dcjkrefinements%26ct%3Dresult%26cd%3D6&location=http%3A//www.100gps.com/html/99/viewthread_tid_6399.html&color=16x&resolution=1024x768&returning=0&language=zh-cn&ua=Mozilla/4.0%20%28compatible%3B%20MSIE%206.0%3B%20Windows%20NT%205.0%3B%20.NET%20CLR%202.0.50727%3B%20.NET%20CLR%201.1.4322%29" frameborder="0" width="0" scrolling="no" height="0"> <script language="javascript" src="http://www.100gps.com/batch.html.php?uid=&itemid=&hash=cd7f22d3efb1b033dcab9b73de8c1734&p=.%2Fhtml%2F99%2Fviewthread_tid_6399.html&u=http%3A%2F%2Fwww.100gps.com%2Findex.php%3Faction%2Fviewthread%2Ftid%2F6399%2Fphp%2F1&t=viewthread" type="text/javascript"></script>  
  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值