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,在d出窗口中将unite参数设为meters,projection,参数设为Tansverse Mercator。接下来再点击edit菜单中的Add/Change projection,在d出对话框中将原始投影参数添加进去。
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
sid文件是什么我不知道,一般地理编码都会在头文件里的。下面我说的方法不一定有用,但你可以参考下。
不知道你的文件可不可以在erdas里直接打开,如果不可以,先导入文件转换为img格式。
你说的坐标也不知道是什么坐标系下的坐标,还有你自己有没有GCP,如果只要地理坐标的话,在Google earth上找特征点,在erdas中进行影像的几何校正就可以了。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)