erdas中坐标的定义

1、ENVI和ERDAS中自定义坐标系的方法1地理投影的根本原理常用到的地图坐标系有2种，即地理坐标系和投影坐标系。地理坐标系是以经纬度为单位的地球坐标系统，地理坐标系中有2个重要部分，即地球椭球体spheroid和大地基准面datum。由于地球外表的不规那么性，它不能用数学公式来表达，也就无法实施运算，所以必须找一个形状和大小都很接近地球的椭球体来代替地球，这个椭球体被称为地球椭球体，我国常用的椭球体如下(JK)(jT&inO3U 2K1510停第車舁。|.|J(140 0<15(755 Jh 29ft 257大地基准面指目前参考椭球与 WGS8参考椭球间的相对位置关系3个平移，

2、3个旋转，1个缩放，可以用其中3个、4个或者7个参数来描述它们之间的 关系，每个椭球体都对应一个或多个大地基准面。投影坐标系是利用一定的数学法那么把地球外表上的经纬线网表示到平面上， 属于平面坐标系。数学法那么指的是投影类型，目前我国普遍采用的是高斯一一克吕格投影，在英美国家称为横轴墨卡托投影Transverse Mercator 。高斯克吕格投影的中央经线和赤道为互相垂直，分带标准分为3度带和6度带。美国编制世界各地军用地图和地球资源卫星像片所采用的全球横轴墨卡托投影UTM是横轴墨卡托投影的一种变型。高斯克 吕格投影的中央经线长度比等于1，UTM投影规定中央经线长度比为 0.9996。我国规

3、定1： 1万、1： 2.5万、1： 5万、1： 10万、1： 25万、1： 50万比 例尺地形图，均采用高斯克吕格投影。1： 2.5至1： 50万比例尺地形图采用经 差6度分带，1： 1和1： 2.5万比例尺地形图采用经差3度分带。大地坐标，在地面上建立一系列相连接的三角形， 量取一段精确的距离作为 起算边，在这个边的两端点，采用天文观测的方法确定其点位经度、纬度和 方位角，用精密测角仪器测定各三角形的角值，根据起算边的边长和点位，就 可以推算出其他各点的坐标。这样推算出的坐标，称为大地坐标。我国1954年在北京设立了大地坐标原点，由此计算出来的各大地控制点的 坐标，称为1954年北京坐标系。

4、为了适应大地测量的开展，我国于1978年米用 国际大地测量协会推荐的IAG-75地球椭球体建立了我国新的大地坐标系，并在1986年宣布在陕西省泾阳县设立了新的大地坐标原点，由此计算出来的各大地 控制点坐标，称为1980年大地坐标系。北肃 XGuiFKiMCtf CTiwonK 冷Kfahvr朗r北栩电IMAM皆俅对于地理坐标，只需要确定两个参数，即椭球体和大地基准面。对于投影坐标，投影类型为 Gauss KrugerTransverse Mercator，除了 确定椭球体和大地基准面外，还需要确定中央经线。大地基准面确实定关键是确定7个参数或者其中几个参数，北京 54基准面可以用三个平移参数来

5、确定，即“-12,-113,- 41,0,0,0,0 ，很多软件近似为Krasovsky(0,0,0,0,0,0,0)基准面；西安80的7参数比拟特殊，各个区域不一样。一般有两个途径：一是直接从测绘部门获取；二是根据三个以上具有西 安80坐标系与其他坐标系的同名点坐标值，利用软件来推算，有一些绿色软件 具有这个功能，如Coord MG中央经线获取可有以下两种方法，第一种根据带号计算，6度带用6*N-3，3度带用3*N；第二种方法是根据经度从图1中查找。ENVI中定义方法：1选择 Map> Customized Map Projection2当出现 Customized MapProjec

6、tion Definition对话框时，键入一个新的投影名和参数。具体参数设置如下：西安80坐标系的定义：Projection Name:比方：xian80 投影名称Projectio n Type: Tran sverse Mercator投影类型Projection Ellipsoid: User Defi ned a=6378140,b=6356755.288准椭球False East ing: 500000东伪偏移False Northing: 0.000000北伪偏移Latitude of projection origin: 0.000000原点纬度Lon gitude of ce

7、n tral meridia n: 117中央经度Scale Factor: 1.000000比例系数北京54坐标系参数如下:Projection Name: 比方：北京 54投影名称Projectio n Type: Tran sverse Mercator投影类型Projectio n Ellipsoid: KrassovskyFalse East ing: 500000东伪偏移False Northing: 0.000000北伪偏移Latitude of projection origin: 0.000000原点纬度Lon gitude of cen tral meridia n: 11

8、7中央经度Scale Factor: 1.000000比例系数北京本地独立坐标系参数如下:Projection Name:比方：beijinglocal投影名称Projectio n Type: Tran sverse Mercator投影类型Projectio n Ellipsoid: KrassovskyFalse East ing: 500000东伪偏移False Northi ng: 300000北伪偏移Latitude of project ion orig in: 39 51 56.757原点纬度Lon gitude of central meridia n: 116 21 0.9

9、065中央经度Scale Factor: 1.000000比例系数> Add NewProject ion3一旦所有参数都被正确输入，选择Projection投影添加到ENVI用的投影列表中。当你关闭对话框时，可利用的投影将根据当前的ENVI设置发生改变。系统将提 问你是否将这一投影存到“map_proj.txt 文件中。要存储新的或更改正的投影信息，选择File > Save Projectio ns 。在ENVI目录结构中的“map_proj.txt 文件将发生改变，包含新的投影。这一文件也能用一些编辑器进行编辑到达上述改变定义的效果。4选择 Projection>Loa

10、d Existing Projection ，出现 Map ProjectionSelection对话框，选择任意坐标系，可以查看参数设置，也可以更改参数设置。5选择File > Cancel退出这一功能。ERDA中定义方法以定义西安 80坐标系为例ERDA中添加椭球体和基准面：ERDA包含了一个能够自定义椭球体、基准面、投影方式的扩展库，通过这个扩 展库，可以在ERDA中添加任何可能存在的投影系统。基于以上2点，就可以将IAG75椭球的参数添加到ERDA中,并且应用这个椭球对栅格数据进行投影变换。 下面就来说明一下具体的添加过程。1在ERDA安装目录下的etc/spheroid.tab

11、文件是用来记载椭球体和基准面参数的。它是一个TXT文本文件，可以用文本编辑器对它进行修改，只要依照它的语法就可以任意添加自定义的椭球体和基准面参数。根本语法为：“椭球名称 “椭球序号椭球体长半轴椭球体短半轴“椭球名称 0 0 0 0 0 0 0“基准面名称 T 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个轴对于 WGS8基准点的平移参数，单位为mb rx、ry、rz是x、y、z、3个轴对于 WGS84 基准点的旋转参数，单位为rad。Ds

12、是对于WGS8基准点的比例因子。在更多的情况下椭球的基准面是基于它本身的。这时假定椭球的中心点是与没有经过任何平移或旋转的WGS8的基准面相重合，即这时椭球基准面的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

13、经过以上的操作IAG75椭球就会出现在ERDAS勺椭球选择列表中。2.在 Viewer 中翻开图像数据，Utility->layer info，在 projection info 栏 中可以看到目前的数据投影信息还不完整。点击edit菜单中的cha nge mapmodel，在弹出窗口中将 unite 参数设为 meters ,projection ，参数设为 Tansverse Mercator。接下来再点击 edit菜单中的 Add/Change projection ，在弹出对话 框中将原始投影参数添加进去。CustomProjectio n Type:Tra nsverse MercatorSpheroid Name:IAG 75Datum Name:xia n80Scale factor at central merid