ENVI和ERDAS中自定义坐标系的方法.docx

ENVI和ERDAS中自定义坐标系的方法 （转）ENVI和ERDAS中自定义坐标系的方法 什么是80西安坐标系？1978年4月在西安召开全国天文大地网平差会议，确定重新定位，建立我国新的 坐标系。为此有了1980年国家大地坐标系。1980年国家大地坐标系采用地球椭球基本参数为1975年国际大地测量与地球物理联合会第十六届大会推荐的 数据。该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇，位于西安市西北方向约60公里，故称1980年西安坐标系，又简称西安大地原点。基准面采用 青岛大港验潮站19521979年确定的黄海平均海水面（即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: 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 -投影名称Projection 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 Projection 将投影添加到 ENVI 用的投影列表中。当你关闭对话框时，可利用的投影将根据当前的 ENVI 设置发生改变。系统将提问你是否将这一投影存到 “map_proj.txt” 文件中。要存储新的或更改过的投影信息，选择 File Save Projections 。在 ENVI 目录结构中的 “map_proj.txt” 文件将发生改变，包含新的投影。这一文件也能用一些编辑器进行编辑达到上述改变定义的效果。（4）选择ProjectionLoad 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.2882xian 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 meridian:1.000000Longitude of central meridian:117:00:00.000000000000 E