erdas中坐标的定义

1、ENVI和ERDAS中自定义坐标系的方法1 地理投影的根本原理常用到的地图坐标系有2种，即地理坐标系和投影坐标系。地理坐标系是以经纬度为单位的地球坐标系统，地理坐标系中有2个重要局部，即地球椭球体spheroid和大地基准面datum。由于地球外表的不规那么性，它不能用数学公式来表达，也就无法实施运算，所以必须找一个形状和大小都很接近地球的椭球体来代替地球，这个椭球体被称为地球椭球体，我国常用的椭球体如下大地基准面指目前参考椭球与WGS84参考椭球间的相对位置关系3个平移，3个旋转，1个缩放，可以用其中3个、4个或者7个参数来描述它们之间的关系，每个椭球体都对应一个或多个大地基准面。投影坐标系

2、是利用一定的数学法那么把地球外表上的经纬线网表示到平面上，属于平面坐标系。数学法那么指的是投影类型，目前我国普遍采用的是高斯克吕格投影，在英美国家称为横轴墨卡托投影Transverse Mercator。高斯克吕格投影的中央经线和赤道为互相垂直，分带标准分为3度带和6度带。美国编制世界各地军用地图和地球资源卫星像片所采用的全球横轴墨卡托投影UTM是横轴墨卡托投影的一种变型。高斯克吕格投影的中央经线长度比等于1，UTM投影规定中央经线长度比为0.9996。    我国规定1：1万、1：2.5万、1：5万、1：10万、1：25万、1：50万比例尺地形图，均采用高斯克吕格投影。1

3、：2.5至1：50万比例尺地形图采用经差6度分带，1：1和1：2.5万比例尺地形图采用经差3度分带。大地坐标，在地面上建立一系列相连接的三角形，量取一段精确的距离作为起算边，在这个边的两端点，采用天文观测的方法确定其点位经度、纬度和方位角，用精密测角仪器测定各三角形的角值，根据起算边的边长和点位，就可以推算出其他各点的坐标。这样推算出的坐标，称为大地坐标。   我国1954年在北京设立了大地坐标原点，由此计算出来的各大地控制点的坐标，称为1954年北京坐标系。为了适应大地测量的开展，我国于1978年采用国际大地测量协会推荐的IAG-75地球椭球体建立了我国新的大地坐标系，并

4、在1986年宣布在陕西省泾阳县设立了新的大地坐标原点，由此计算出来的各大地控制点坐标，称为1980年大地坐标系。对于地理坐标，只需要确定两个参数，即椭球体和大地基准面。对于投影坐标，投影类型为Gauss KrugerTransverse Mercator，除了确定椭球体和大地基准面外，还需要确定中央经线。 大地基准面确实定关键是确定7个参数或者其中几个参数，北京54基准面可以用三个平移参数来确定，即“-12,-113,-41,0,0,0,0，很多软件近似为Krasovsky(0,0,0,0,0,0,0)基准面；西安80的7参数比拟特殊，各个区域不一样。一般有两个途径：一是直接从测绘部门获取；二

5、是根据三个以上具有西安80坐标系与其他坐标系的同名点坐标值，利用软件来推算，有一些绿色软件具有这个功能，如Coord MG。中央经线获取可有以下两种方法，第一种根据带号计算，6度带用6*N-3，3度带用3*N；第二种方法是根据经度从图1中查找。ENVI中定义方法：1 选择 Map> Customized Map Projection2 当出现 Customized Map Projection Definition 对话框时，键入一个新的投影名和参数。具体参数设置如下：西安80坐标系的定义：Projection Name: 比方：xian80 -投影名称Projection Type:

6、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 -投影名

7、称Projection Type: Transverse Mercator -投影类型Projection Ellipsoid: KrassovskyFalse Easting: 500000 -东伪偏移False Northing: 0.000000 -北伪偏移Latitude of projection origin: 0.000000 -原点纬度Longitude of central meridian: 117 -中央经度Scale_Factor: 1.000000 -比例系数北京本地独立坐标系参数如下：Projection Name: 比方：beijinglocal -投影名称Pro

8、jection Type: Transverse Mercator -投影类型Projection Ellipsoid: KrassovskyFalse 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 Proje

9、ction 将投影添加到 ENVI 用的投影列表中。当你关闭对话框时，可利用的投影将根据当前的 ENVI 设置发生改变。系统将提问你是否将这一投影存到 “map_proj.txt 文件中。要存储新的或更改正的投影信息，选择 File > Save Projections 。在 ENVI 目录结构中的 “map_proj.txt 文件将发生改变，包含新的投影。这一文件也能用一些编辑器进行编辑到达上述改变定义的效果。4选择Projection>Load Existing Projection，出现Map Projection Selection对话框，选择任意坐标系，可以查看参数设置，

10、也可以更改参数设置。5 选择 File > Cancel 退出这一功能。ERDAS中定义方法以定义西安80坐标系为例ERDAS中添加椭球体和基准面：ERDAS包含了一个能够自定义椭球体、基准面、投影方式的扩展库，通过这个扩展库，可以在ERDAS中添加任何可能存在的投影系统。基于以上2点，就可以将IAG75椭球的参数添加到ERDAS中，并且应用这个椭球对栅格数据进行投影变换。下面就来说明一下具体的添加过程。1 在ERDAS安装目录下的etc/spheroid.tab文件是用来记载椭球体和基准面参数的。它是一个TXT文本文件，可以用文本编辑器对它进行修改，只要依照它的语法就可以任意添加自定义

11、的椭球体和基准面参数。根本语法为：“椭球名称“椭球序号椭球体长半轴 椭球体短半轴“椭球名称 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的

12、基准面相重合，即这时椭球基准面的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，在弹出对话框中将原始投影参数添加进去。CustomProjection Type:Transverse MercatorSpheroid Name:IAG 75Datum Name:xian80Scale factor at central meri