ENVI和ERDAS中自定义坐标系的方法.docx
- 文档编号:10537017
- 上传时间:2023-02-21
- 格式:DOCX
- 页数:19
- 大小:27.03KB
ENVI和ERDAS中自定义坐标系的方法.docx
《ENVI和ERDAS中自定义坐标系的方法.docx》由会员分享,可在线阅读,更多相关《ENVI和ERDAS中自定义坐标系的方法.docx(19页珍藏版)》请在冰豆网上搜索。
ENVI和ERDAS中自定义坐标系的方法
ENVI和ERDAS中自定义坐标系的方法
2008-01-2509:
20
ENVI和ERDAS中自定义坐标系的方法
什么是80西安坐标系?
1978年4月在西安召开全国天文大地网平差会议,确定重新定位,建立我国新的坐标系。
为此有了1980年国家大地坐标系。
1980年国家大地坐标系采用地球椭球基本参数为1975年国际大地测量与地球物理联合会第十六届大会推荐的数据。
该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇,位于西安市西北方向约60公里,故称1980年西安坐标系,又简称西安大地原点。
基准面采用青岛大港验潮站1952-1979年确定的黄海平均海水面(即1985国家高程基准)。
1980西安坐标系的投影椭球为IAG75,a=6378140,b=6356755.288投影类型为横轴墨卡托投影(transverse)。
ENVI中定义方法:
(1)选择Map>CustomizedMapProjection
(2)当出现CustomizedMapProjectionDefinition对话框时,键入一个新的投影名和参数。
具体参数设置如下:
西安80坐标系的定义:
ProjectionName:
比如:
xian80 --------投影名称
ProjectionType:
TransverseMercator --------投影类型
ProjectionEllipsoid:
UserDefineda=6378140,b=6356755.288--------基准椭球
FalseEasting:
500000 --------东伪偏移
FalseNorthing:
0.000000 ---------北伪偏移
Latitudeofprojectionorigin:
0.000000 ---------原点纬度
Longitudeofcentralmeridian:
117 ---------中央经度
Scale_Factor:
1.000000 ----------比例系数
北京54坐标系参数如下:
ProjectionName:
比如:
北京54 --------投影名称
ProjectionType:
TransverseMercator --------投影类型
ProjectionEllipsoid:
Krassovsky
FalseEasting:
500000 --------东伪偏移
FalseNorthing:
0.000000 ---------北伪偏移
Latitudeofprojectionorigin:
0.000000 ---------原点纬度
Longitudeofcentralmeridian:
117 ---------中央经度
Scale_Factor:
1.000000 ----------比例系数
北京本地独立坐标系参数如下:
ProjectionName:
比如:
beijinglocal --------投影名称
ProjectionType:
TransverseMercator --------投影类型
ProjectionEllipsoid:
Krassovsky
FalseEasting:
500000 --------东伪偏移
FalseNorthing:
300000 ---------北伪偏移
Latitudeofprojectionorigin:
395156.757 ---------原点纬度
Longitudeofcentralmeridian:
116210.9065---------中央经度
Scale_Factor:
1.000000 ----------比例系数
(3)一旦所有参数都被正确输入,选择Projection>AddNewProjection将投影添加到ENVI用的投影列表中。
当你关闭对话框时,可利用的投影将根据当前的ENVI设置发生改变。
系统将提问你是否将这一投影存到“map_proj.txt”文件中。
要存储新的或更改过的投影信息,选择File>SaveProjections。
在ENVI目录结构中的“map_proj.txt”文件将发生改变,包含新的投影。
这一文件也能用一些编辑器进行编辑达到上述改变定义的效果。
(4)选择Projection>LoadExistingProjection,出现MapProjectionSelection对话框,选择任意坐标系,可以查看参数设置,也可以更改参数设置。
(5)选择File>Cancel退出这一功能。
ERDAS中定义方法(以定义西安80坐标系为例)
ERDAS中添加椭球体和基准面:
ERDAS包含了一个能够自定义椭球体、基准面、投影方式的扩展库,通过这个扩展库,可以在ERDAS中添加任何可能存在的投影系统。
基于以上2点,就可以将IAG75椭球的参数添加到ERDAS中,并且应用这个椭球对栅格数据进行投影变换。
下面就来说明一下具体的添加过程。
1在ERDAS安装目录下的etc/spheroid.tab文件是用来记载椭球体和基准面参数的。
它是一个TXT文本文件,可以用文本编辑器对它进行修改,只要依照它的语法就可以任意添加自定义的椭球体和基准面参数。
基本语法为:
“椭球名称”{
“椭球序号”椭球体长半轴椭球体短半轴
“椭球名称”0000000
“基准面名称1”dx1dy1dz1rx1ry1rz1ds1
“基准面名称2”dx2dy2dz2rx2ry2rz1ds2
……….
}
其中:
“基准面名称”dxdydzrxrzds中,dx、dy、dz是x、y、z3个轴对于WGS84基准点的平移参数,单位为m。
rx、ry、rz是x、y、z、3个轴对于WGS84基准点的旋转参数,单位为rad。
Ds是对于WGS84基准点的比例因子。
在更多的情况下椭球的基准面是基于它本身的。
这时假定椭球的中心点是与没有经过任何平移或旋转的WGS84的基准面相重合,即这时椭球基准面的7个参数均为0,即这时椭球基准面的7个参数均为0。
我国在使用克拉索夫斯基椭球和IAG75椭球时就是用椭球体本身为基准。
在spheroid.tab文件末尾加入如下语句即可,假设spheroid.tab文件中最后一个椭球体序号为73(可以在文件最后一个椭球体中读出序号),则加入:
"IAG75"{
7463781406356755.2882
"xian80"0000000
}
经过以上的操作IAG75椭球就会出现在ERDAS的椭球选择列表中。
2.在Viewer中打开图像数据,Utility--->layerinfo,在projectioninfo栏中可以看到目前的数据投影信息还不完整。
点击edit菜单中的changemapmodel,在弹出窗口中将unite参数设为meters,projection,参数设为TansverseMercator。
接下来再点击edit菜单中的Add/Changeprojection,在弹出对话框中将原始投影参数添加进去。
Custom
ProjectionType:
TransverseMercator
SpheroidName:
IAG75
DatumName:
xian80
Scalefactoratcentralmeridian:
1.000000
Longitudeofcentralmeridian:
117:
00:
00.000000000000E
Latitudeoforiginofprojection:
0:
00:
00.000000000000N
Falseeasting:
39500000.00000000000meters
| 旧一篇:
中国主要城市经纬度表
一、只谈比较常用的几种:
“墨卡托投影”、“高斯-克吕格投影”、“UTM投影”、“兰勃特等角投影。
1.墨卡托(Mercator)投影
1.1墨卡托投影简介
墨卡托(Mercator)投影,是一种"等角正切圆柱投影”,荷兰地图学家墨卡托(GerhardusMercator1512-1594)在1569年拟定,假设地球被围在一中空的圆柱里,其标准纬线与圆柱相切接触,然后再假想地球中心有一盏灯,把球面上的图形投影到圆柱体上,再把圆柱体展开,这就是一幅选定标准纬线上的“墨卡托投影”绘制出的地图。
墨卡托投影没有角度变形,由每一点向各方向的长度比相等,它的经纬线都是平行直线,且相交成直角,经线间隔相等,纬线间隔从标准纬线向两极逐渐增大。
墨卡托投影的地图上长度和面积变形明显,但标准纬线无变形,从标准纬线向两极变形逐渐增大,但因为它具有各个方向均等扩大的特性,保持了方向和相互位置关系的正确。
在地图上保持方向和角度的正确是墨卡托投影的优点,墨卡托投影地图常用作航海图和航空图,如果循着墨卡托投影图上两点间的直线航行,方向不变可以一直到达目的地,因此它对船舰在航行中定位、确定航向都具有有利条件,给航海者带来很大方便。
“海底地形图编绘规范”(GB/T17834-1999,海军航保部起草)中规定1:
25万及更小比例尺的海图采用墨卡托投影,其中基本比例尺海底地形图(1:
5万,1:
25万,1:
100万)采用统一基准纬线30°,非基本比例尺图以制图区域中纬为基准纬线。
基准纬线取至整度或整分。
1.2墨卡托投影坐标系
取零子午线或自定义原点经线(L0)与赤道交点的投影为原点,零子午线或自定义原点经线的投影为纵坐标X轴,赤道的投影为横坐标Y轴,构成墨卡托平面直角坐标系。
2.高斯-克吕格(Gauss-Kruger)投影和UTM(UniversalTransverseMercator)投影
2.1高斯-克吕格投影简介
高斯-克吕格(Gauss-Kruger)投影,是一种“等角横切圆柱投影”。
德国数学家、物理学家、天文学家高斯(CarlFriedrichGauss,1777一1855)于十九世纪二十年代拟定,后经德国大地测量学家克吕格(JohannesKruger,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.2UTM投影简介
UTM投影全称为“通用横轴墨卡托投影”,是一种“等角横轴割圆柱投影”,椭圆柱割地球于南纬80度、北纬84度两条等高圈,投影后两条相割的经线上没有变形,而中央经线上长度比0.9996。
UTM投影是为了全球战争需要创建的,美国于1948年完成这种通用投影系统的计算。
与高斯-克吕格投影相似,该投影角度没有变形,中央经线为直线,且为投影的对称轴,中央经线的比例因子取0.9996是为了保证离中央经线左右约330km处有两条不失真的标准经线。
UTM投影分带方法与高斯-克吕格投影相似,是自西经180°起每隔经差6度自西向东分带,将地球划分为60个投影带。
我国的卫星影像资料常采用UTM投影。
2.3高斯-克吕格投影与UTM投影异同
高斯-克吕格(Gauss-Kruger)投影与UTM投影(UniversalTransverseMercator,通用横轴墨卡托投影)都是横轴墨卡托投影的变种,目前一些国外的软件或国外进口仪器的配套软件往往不支持高斯-克吕格投影,但支持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万地形图)。
GIS中坐标系定义是GIS系统的基础,GIS中的坐标系由基准面(Datum)和地图投影(Projection)两组参数确定。
地球椭球体:
地球是一个表面很复杂的球体,人们以假想的平均静止的海水面形成的“大地体”为参照,推求出近似的椭球体,理论和实践证明,该椭球体近似一个以地球短轴为轴的椭园而旋转的椭球面,这个椭球面可用数学公式表达,将自然表面上的点归化到这个椭球面上,就可以计算了。
常用的一些椭球及参数
海福特椭球(1910) 我国52年以前基准椭球
a=6378388mb=6356911.9461279mα=0.33670033670
克拉索夫斯基椭球(1940Krassovsky) 北京54坐标系基准椭球
a=6378245mb=6356863.018773mα=0.335232986921975年I.U.G.G推荐椭球(国际大地测量协会1975)西安80坐标系基准椭球 a=6378140mb=6356755.2881575mα=0.0033528131778
WGS-84椭球(GPS全球定位系统椭球、17届国际大地测量协会) WGS-84GPS基准椭球
a=6378137mb=6356752.3142451mα=0.00335281006247
Krasovsky_1940椭球及其相应参数
Alias:
Abbreviation:
Remarks:
AngularUnit:
Degree(0.017453292519943299)
PrimeMeridian(起始经度):
Greenwich(0.000000000000000000)
Datum(大地基准面):
D_Beijing_1954
Spheroid(参考椭球体):
Krasovsky_1940
SemimajorAxis:
6378245.000000000000000000
SemiminorAxis:
6356863.018773047300000000
InverseFlattening:
298.300000000000010000
地球椭球面上任一点的位置,可由该点的纬度(B)和精度(L)确定,即地面点的地理坐标值,由经线和纬线构成两组互相正交的曲线坐标网叫地理坐标网。
由经纬度构成的地理坐标系统又叫地理坐标系。
地理坐标分为天文地理坐标和大地地理坐标。
天文地理坐标是用天文测量方法确定的,大地地理坐标是用大地测量方法确定的。
我们在地球椭球面上所用的地理坐标系属于大地地理坐标系,简称大地坐标系
确定椭球的大小后,还要进行椭球定向,即把旋转椭球面套在地球的一个适当的位置,这一位置就是该地理坐标系的“坐标原点”,是全部大地坐标计算的起算点,俗称“大地原点”
基准面是利用特定椭球体对特定地区地球表面的逼近,因此每个国家或地区均有各自的基准面,我们通常称谓的北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。
我国参照前苏联从1953年起采用克拉索夫斯基(Krassovsky)椭球体建立了我国的北京54坐标系,1978年采用国际大地测量协会推荐的1975地球椭球体建立了我国新的大地坐标系--西安80坐标系,目前大地测量基本上仍以北京54坐标系作为参照,北京54与西安80坐标之间的转换可查阅国家测绘局公布的对照表。
椭球体与基准面之间的关系是一对多的关系,也就是基准面是在椭球体基础上建立的,但椭球体不能代表基准面,同样的椭球体能定义不同的基准面,如前苏联的Pulkovo1942、非洲索马里的Afgooye基准面都采用了Krassovsky椭球体,但它们的基准面显然是不同的。
GIS中地图投影的定义:
是为解决由不可展的椭球面描绘到平面上的矛盾,用几何透视方法或数学分析的方法,将地球上的点和线投影到可展的曲面(平面、园柱面或圆锥面)上,将此可展曲面展成平面,建立该平面上的点、线和地球椭球面上的点、线的对应关系。
我国的基本比例尺地形图(1:
5千,1:
1万,1:
2.5万,1:
5万,1:
10万,1:
25万,1:
50万,1:
100万)中,大于等于50万的均采用高斯-克吕格投影(Gauss-Kruger),又叫横轴墨卡托投影(TransverseMercator);小于50万的地形图采用正轴等角割园锥投影,又叫兰勃特投影(LambertConformalConic);海上小于50万的地形图多用正轴等角园柱投影,又叫墨卡托投影(Mercator),我国的GIS系统中应该采用与我国基本比例尺地形图系列一致的地图投影系统。
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
LinearUnit:
Meter(1.000000)
GeographicCoordinateSystem:
Name:
GCS_Beijing_1954
Alias:
Abbreviation:
Remarks:
AngularUnit:
Degree(0.017453292519943299)
PrimeMeridian:
Greenwich(0.000000000000000000)
Datum:
D_Beijing_1954
Spheroid:
Krasovsky_1940
SemimajorAxis:
6378245.000000000000000000
SemiminorAxis:
6356863.018773047300000000
InverseFlattening:
298.300000000000010000
高斯-克吕格直角坐标
高斯-克吕格投影是设想用一个椭圆柱横套在地球椭球的外面,并与设定的中央经线相切。
高斯-克吕格投影分带规定:
该投影是国家基本比例尺地形图的数学基础,为控制变形,采用分带投影的方法,在比例尺1:
2.5万-1:
50万图上采用6°分带,对比例尺为1:
1万及大于1:
1万的图采用3°分带。
6°分带法:
从格林威治零度经线起,每6°分为一个投影带,全球共分为60个投影
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ENVI ERDAS 自定义 坐标系 方法