《2022年ERDAS中如何自定义坐标系 .pdf》由会员分享，可在线阅读，更多相关《2022年ERDAS中如何自定义坐标系 .pdf（6页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、ERDAS 中如何自定义坐标系ENVI和 ERDAS 中自定义坐标系的方法1 地理投影的基本原理常用到的地图坐标系有2 种，即地理坐标系和投影坐标系。地理坐标系是以经纬度为单位的地球坐标系统，地理坐标系中有2 个重要部分，即地球椭球体（ spheroid）和大地基准面（ datum） 。由于地球表面的不规则性， 它不能用数学公式来表达， 也就无法实施运算， 所以必须找一个形状和大小都很接近地球的椭球体来代替地球，这个椭球体被称为地球椭球体，我国常用的椭球体如下大地基准面指目前参考椭球与WGS84参考椭球间的相对位置关系（3 个平移，3个旋转， 1 个缩放） ， 可以用其中 3 个、 4 个或者
2、 7 个参数来描述它们之间的关系，每个椭球体都对应一个或多个大地基准面。投影坐标系是利用一定的数学法则把地球表面上的经纬线网表示到平面上，属于平面坐标系。数学法则指的是投影类型，目前我国普遍采用的是高斯克吕格投影，在英美国家称为横轴墨卡托投影（Transverse Mercator ） 。高斯克吕格投影的中央经线和赤道为互相垂直，分带标准分为3 度带和 6 度带。美国编制世界各地军用地图和地球资源卫星像片所采用的全球横轴墨卡托投影（UTM）是横轴墨卡托投影的一种变型。高斯克吕格投影的中央经线长度比等于1，UTM 投影规定中央经线长度比为0.9996。我国规定 1：1 万、1：2.5 万、1：5
3、 万、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 年大地坐标系。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 6 页 - - - - - - - - - 对于地理坐标，只需要确定两个参数，即椭球体和大地基准面。对于投影坐标，投影类型为Gauss Kruger（Transverse Mercator） ，除了确定椭球体
5、和大地基准面外，还需要确定中央经线。大地基准面的确定关键是确定7 个参数（或者其中几个参数） ，北京 54 基准面可以用三个平移参数来确 定 ， 即 “ -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
6、*N；第二种方法是根据经度从图1 中查找。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 - 基准椭球Fa
7、lse Easting: 500000 - 东伪偏移False Northing: 0.000000 - 北伪偏移Latitude of projection origin: 0.000000 -原点纬度Longitude of central meridian: 117 -中央经度Scale_Factor: 1.000000 - 比例系数名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 6 页 - - - - - - - - - 北京 54 坐标系参数如下：Project
8、ion Name: 比如：北京 54 -投影名称Projection Type: Transverse Mercator - 投影类型Projection Ellipsoid: Krassovsky False Easting: 500000 - 东伪偏移False Northing: 0.000000 - 北伪偏移Latitude of projection origin: 0.000000 -原点纬度Longitude of central meridian: 117 -中央经度Scale_Factor: 1.000000 - 比例系数北京本地独立坐标系参数如下：Projection Na
9、me: 比如： beijinglocal -投影名称Projection Type: Transverse Mercator - 投影类型Projection Ellipsoid: Krassovsky False 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） 一旦所有参数都被
10、正确输入，选择Projection > Add New Projection 将投影添加到ENVI 用的投影列表中。当你关闭对话框时， 可利用的投影将根据当前的ENVI 设置发生改变。 系统将提问你是否将这一投影存到“map_proj.txt” 文件中。要存储新的或更改过的投影信息，选择File > Save Projections 。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 6 页 - - - - - - - - - 在 ENVI 目录结构中的“map
11、_proj.txt” 文件将发生改变，包含新的投影。这一文件也能用一些编辑器进行编辑达到上述改变定义的效果。（4）选择 Projection>Load Existing Projection ，出现 Map Projection Selection对话框，选择任意坐标系，可以查看参数设置，也可以更改参数设置。（5） 选择 File > Cancel 退出这一功能。ERDAS 中定义方法（以定义西安80坐标系为例）ERDAS 中添加椭球体和基准面：ERDAS 包含了一个能够自定义椭球体、基准面、投影方式的扩展库，通过这个扩展库，可以在ERDAS中添加任何可能存在的投影系统。基于以上2
12、 点，就可以将 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 ,. 其中：
13、 “基准面名称”dx dy dz rx rz ds中，dx、dy、dz 是 x、y、z3个轴对于 WGS84基准点的平移参数，单位为m。rx、ry、rz 是 x、y、z、3 个轴对于 WGS84基准名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 6 页 - - - - - - - - - 点的旋转参数，单位为rad。Ds是对于 WGS84基准点的比例因子。在更多的情况下椭球的基准面是基于它本身的。这时假定椭球的中心点是与没有经过任何平移或旋转的WGS84的基准面相重合，即
14、这时椭球基准面的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->
15、;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 c
16、entral 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 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 6 页 - - - - - - - - - 名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 6 页 - - - - - - - - -