请选择 进入手机版 | 继续访问电脑版

遥感集市论坛

 
登录  立即注册
查看: 1444|回复: 0

ENVI和ERDAS中自定义坐标系的方法添加80坐标

[复制链接]

33

主题

51

帖子

151

积分

版主

Rank: 7Rank: 7Rank: 7

积分
151
发表于 2016-2-29 15:16:57 | 显示全部楼层 |阅读模式
ENVI和ERDAS中自定义坐标系的方法添加80坐标


什么是80西安坐标系?
1978年4月在西安召开全国天文大地网平差会议,确定重新定位,建立我国新的坐标系。为此有了1980年国家大地坐标系。1980年国家大地坐标系采用地球椭球基本参数为1975年国际大地测量与地球物理联合会第十六届大会推荐的数据。该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇,位于西安市西北方向约60公里,故称1980年西安坐标系,又简称西安大地原点。基准面采用青岛大港验潮站1952-1979年确定的黄海平均海水面(即1985国家高程基准)。1980西安坐标系的投影椭球为IAG75, a=6378140,b=6356755.288 投影类型为横轴墨卡托投影(transverse)。


ENVI中定义方法:
(1) 选择 Map> Customized Map Projection
(2) 当出现 Customized Map Projection Definition 对话框时,键入一个新的投影名和参数。具体参数设置如下


西安80坐标系的定义:
Projection Name: 比如:xian80 --------投影名称
Projection Type: Transverse Mercator --------投影类型
Projection Ellipsoid: User Defined a=6378140,b=6356755.288 --------基准椭球
False Easting: 500000 --------东伪偏移
False Northing: 0.000000 ---------北伪偏移
Latitude of projection origin: 0.000000 ---------原点纬度
Longitude of central meridian: 117 ---------中央经度
Scale_Factor: 1.000000 ----------比例系数


北京54坐标系参数如下:
Projection Name: 比如:北京54 --------投影名称
Projection Type: Transverse Mercator --------投影类型
Projection Ellipsoid: Krassovsky
False Easting: 500000 --------东伪偏移
False Northing: 0.000000 ---------北伪偏移
Latitude of projection origin: 0.000000 ---------原点纬度
Longitude of central meridian: 117 ---------中央经度
Scale_Factor: 1.000000 ----------比例系数


北京本地独立坐标系参数如下:
Projection Name: 比如:beijinglocal --------投影名称
Projection Type: Transverse Mercator --------投影类型
Projection Ellipsoid: Krassovsky
False Easting: 500000 --------东伪偏移
False Northing: 300000 ---------北伪偏移
Latitude of projection origin: 39 51 56.757 ---------原点纬度
Longitude of central meridian: 116 21 0.9065 ---------中央经度
Scale_Factor: 1.000000 ----------比例系数


(3) 一旦所有参数都被正确输入,选择 Projection > Add New Projection 将投影添加到 ENVI 用的投影列表中。
当你关闭对话框时,可利用的投影将根据当前的 ENVI 设置发生改变。系统将提问你是否将这一投影存到 “map_proj.txt” 文件中。
要存储新的或更改过的投影信息,选择 File > Save Projections 。
在 ENVI 目录结构中的 “map_proj.txt” 文件将发生改变,包含新的投影。这一文件也能用一些编辑器进行编辑达到上述改变定义的效果。


(4)选择Projection>Load Existing Projection,出现Map Projection Selection对话框,选择任意坐标系,可以查看参数设置,也可以更改参数设置。


(5) 选择 File > Cancel 退出这一功能。


ERDAS中定义方法(以定义西安80坐标系为例)
ERDAS中添加椭球体和基准面:
ERDAS包含了一个能够自定义椭球体、基准面、投影方式的扩展库,通过这个扩展库,可以在ERDAS中添加任何可能存在的投影系统。基于以上2点,就可以将IAG75椭球的参数添加到ERDAS中,并且应用这个椭球对栅格数据进行投影变换。下面就来说明一下具体的添加过程。


1 在ERDAS安装目录下的etc/spheroid.tab文件是用来记载椭球体和基准面参数的。它是一个TXT文本文件,可以用文本编辑器对它进行修改,只要依照它的语法就可以任意添加自定义的椭球体和基准面参数。
基本语法为:
“椭球名称”{
“椭球序号”椭球体长半轴 椭球体短半轴
“椭球名称” 0 0 0 0 0 0 0
“基准面名称1” dx1 dy1 dz1 rx1 rz1 ds1
“基准面名称2” dx2 dy2 dz2 rx2 rz1 ds2
……….
}
其中:“基准面名称” dx dy dz rx rz ds中,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(可以在文件最后一


个椭球体中读出序号),则加入:
"IAG 75" {
74 6378140 6356755.2882
"xian 80" 0 0 0 0 0 0 0
}




经过以上的操作IAG75 椭球就会出现在ERDAS的椭球选择列表中。
2. 在Viewer中打开图像数据,Utility--->layer info,在projection info 栏中可以看到目前的数据投影信息还不完整。点击edit菜单中的change map model,在弹出窗口中将unite参数设为meters,projection,参数设为Tansverse Mercator。接下来再点击edit菜单中的Add/Change projection,在弹出对话框中将原始投影参数添加进去
Custom
Projection Type:Transverse Mercator
Spheroid Name:IAG 75
Datum Name:xian80
Scale factor at central meridian:1.000000
Longitude of central meridian:117:00:00.000000000000 E
Latitude of origin of projection:0:00:00.000000000000 N
False easting:39500000.00000000000 meters
False northing:0.00000000000000000 meters

回复

使用道具 举报

发表回复

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

其它相关

2 人签到
将意见发给我们
  • 通过email将您的想法和建议发给我们

    编辑部:bbs@rscloudmart.com

联系我们
  • 论坛官方QQ群:128651389

    论坛客服QQ:2739549450

遥感集市服务热线
  • 400-674-7498

QQ|小黑屋|遥感集市论坛 ( 粤ICP备12003192号-2

GMT+8, 2018-10-16 00:41 , Processed in 0.114980 second(s), 33 queries, Gzip On.

Powered by Discuz! X3.1 Copyright
© 2001-2013 Comsenz Inc.    All Rights Reserved.

快速回复 返回顶部 返回列表