24小时热门版块排行榜    

查看: 1698  |  回复: 7
【奖励】 本帖被评价4次,作者lk19920418增加金币 3.2 个

lk19920418

新虫 (小有名气)


[资源] ENVI和ERDAS中自定义坐标系的方法,很有用处,让我解决了困扰两周的问题

ENVI和ERDAS中自定义坐标系的方法
什么是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
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : 定义坐标系.docx
  • 2016-03-30 21:35:01, 15.26 K

» 收录本帖的淘帖专辑推荐

自然科学资源汇总

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kingsabre

木虫 (小有名气)


3楼2016-04-01 14:38:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guoyahui2010

金虫 (正式写手)


4楼2016-04-02 21:38:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
简单回复
bwbjlt2楼
2016-03-31 08:07   回复  
五星好评  顶一下,感谢分享!
absxwj5楼
2016-06-04 17:44   回复  
五星好评  顶一下,感谢分享!
dzqsyd6楼
2016-10-12 22:28   回复  
五星好评  顶一下,感谢分享!
fuyu80467楼
2016-11-04 23:36   回复  
五星好评  顶一下,感谢分享!
fuyu80468楼
2016-11-04 23:38   回复  
顶一下,感谢分享!
相关版块跳转 我要订阅楼主 lk19920418 的主题更新
☆ 无星级 ★ 一星级 ★★★ 三星级 ★★★★★ 五星级
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +11 xqy2009 2026-10-13 13/650 2026-10-15 21:32 by bear2007
[有机交流] 用NBS溴化时溶剂除了CCl4外还有什么 +3 18834002860 2026-10-12 6/300 2026-10-15 14:30 by 长沙摩尔信息科�
[基金申请] 申请科研项目纯靠实力能成功吗? 200+7 David1982 2026-10-13 11/550 2026-10-15 10:30 by zyqchem
[硕博家园] 极限拉扯,如何抉择 +11 Isak0320 2026-10-12 15/750 2026-10-14 21:55 by 猫鱼不是鱼
[论文投稿] 二审4个多月了还是Awaiting Reviewer Scores 10+3 纳斯达克ML 2026-10-11 4/200 2026-10-14 17:30 by babero
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 22:11 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:34 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
信息提示
请填处理意见