ERDAS西安80坐标系创建
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
ERDAS中设置西安80坐标系
1.打开ERDAS目录
一般情况下在C:\Leica Geosystems\Geospatial Imaging 9.2\etc中
找到spheroid.tab文件
选择记事本打开
在记事本中输入以下红色内容
"Xian_1980"{
75 6378140 6356755.2882
"D_xian_1980" 0 0 0 0 0 0 0 }
保存并关闭记事本
2.在ERDAS中配置投影参数
打开ERDAS
选择DataPrep
选择Reproject Images(重投影影像)
在Categories处选择自定义
在Custom选项卡中,首先选择投影类型
Gauss Kruger (高斯克鲁格投影)
Spheriod(椭球)中即可选择刚才添加的“Xian_1980”椭球
设置中央经线为东经99°,向东偏移33500000米
设置完成后保存
在Save as中输入Xian_1980_3_Degree_GK_Zone_33 In Category GCS_Xian_1980
双击回车
点击Cancel 退出自定义面板
便可以在Categories中选择刚刚创建的“GCS_Xian_1980”系列,以及33度带的投影了。
3.影像的重投影
选择需要重投影的源文件,并设置导出路径和文件名,点击左下角的OK即可对影像进行重投影
图中为将名为“mz-pinjie.img”的影像重投影到“GCS_Xian_1980”系列下的
“Xian_1980_3_Degree_GK_Zone_33”投影下。
附:(摘自网络)
RDAS中定义西安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