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

合集下载

Erdas实习报告

Erdas实习报告

ERDAS遥感影像处理综合实习报告一:自定义坐标系(北京54、西安80、2000坐标系)操作步骤:1.1 添加椭球体修改文件为ellipse.txt,语法为<椭球体名称>,<长半轴>,<短半轴>。

这里的逗号为英文半角输入状态下的逗号,建议直接复制文件中已有的椭球体进行修改。

这里将下面三行加在ellipse.txt文件的末尾,保存关闭即可。

最终效果如下图所示。

• Krasovsky,6378245.0,6356863.0• IAG-75,6378140.0,6356755.3• CGCS2000,6378137.0,6356752.3注:ellipse.txt文件中已经有了克拉索夫斯基椭球,由于翻译原因,这里的英文名称是Krassovsky,为了让其他软件平台识别,这里新建一个Krasovsky椭球体。

1.2 添加基准面修改文件为datum.txt,语法为<基准面名称>,<椭球体名称>,<平移三参数>。

这里将下面三行添加在datum.txt文件末尾,保存关闭即可。

最终效果如下图所示。

• D_Beijing_1954, Krasovsky, -12, -113, -41• D_Xian_1980,IAG-75,0,0,0• D_China_2000,CGCS2000,0,0,01.3 定义坐标系操作步骤如下:(1)打开ENVI Classic,选择Map > Customize Map Projection工具;(2)在弹出的CustomizedMap ProjectionDefinition对话框内填写如图所示参数,其中Projection Name保持与ArcGIS中的名称一致;(3)选择Projection > Add New Project ion…,保存投影坐标系;(4)选择File > Save Projections…,在弹出对话框中点击OK,将新建坐标系保存在map_proj.txt文件内,以便下次启动ENVI后依然可以使用。

ENVI,ERDAS,PCI,ER_Mapper的特点与统计分类

ENVI,ERDAS,PCI,ER_Mapper的特点与统计分类

ENVI,ERDAS,PCI,ER_Mapper的特点与统计分类作者:李斯泽来源:《电脑知识与技术·学术交流》2008年第26期摘要:介绍了四种软件的主要功能并比较了在遥感图像处理中应用的特点。

PCI更适合于影像制图,ERDAS的数据融合效果最好,ENVI在针对像元处理的信息提取中功能最强大。

ER Mapper对于处理高分辨率影像效果较好。

关键词:遥感图像处理;统计分类;软件中图分类号:TP311文献标识码:A 文章编号:1009-3044(2008)26-1776-02The Features and Statistical Classification of ENVI, ERDAS, PCI, ER_MapperLI Si-ze(College of Engineering Remote Sensing Information, Wuhan University, Wuhan 430072, China)Abstract: The article introduces in detail the main function and the differences of four kinds of software which apply in remote sensing image. CPI is good at image mapping, ERDAS is skilled in data fusion, ENVI has powerful function in information acquisition when dealing with pixel; while for the software ER Mapper, the result of image manipulation of high resolution is better.Key words: GPS; application; positioning; direction of development1 四种软件1.1 ENVIENVI(The Environment for Visualizing Images)是一套功能齐全的遥感图像处理系统,是处理、分析并显示多光谱数据、高光谱数据和雷达数据的高级工具。

arcgis 创建自定义地理坐标变换 方法 参数值

arcgis 创建自定义地理坐标变换 方法 参数值

arcgis 创建自定义地理坐标变换方法参数值自定义地理坐标变换是在ArcGIS中进行空间数据分析和处理时不可或缺的功能。

通过创建自定义地理坐标变换,我们可以将数据从一个坐标系统转换到另一个坐标系统,实现不同坐标系统之间的数据转化和投影变换。

本文将介绍如何在ArcGIS中使用方法参数值创建自定义地理坐标变换,并探讨其方法和应用。

1. 什么是自定义地理坐标变换在ArcGIS中,自定义地理坐标变换是指用户通过指定一系列参数值和方法来定义自己的地理坐标转换方法,以实现特定的地理坐标转换需求。

这些参数值和方法可以是标准的地理坐标转换方法,也可以是用户自己定义的变换方法。

通过自定义地理坐标变换,我们可以实现高精度、高效率的数据转换和投影变换,满足实际应用中各种不同坐标系统的需求。

2. ArcGIS中的自定义地理坐标变换方法和参数值在ArcGIS中,创建自定义地理坐标变换需要指定一些方法和参数值,来定义地理坐标的转换过程。

一般包括以下几个步骤:- 第一步,选择合适的坐标系统。

根据实际需求,选择待转换数据的坐标系统和目标坐标系统。

在选择坐标系统时,需要考虑数据的地理范围、精度要求和应用环境等因素。

- 第二步,选择合适的转换方法。

根据待转换数据的空间特征和应用需求,选择相应的坐标转换方法。

常用的转换方法包括七参数转换、仿射转换和多项式转换等。

- 第三步,确定转换参数值。

根据所选转换方法,确定相应的参数值。

这些参数值可以根据实际数据和转换需求进行设置。

对于七参数转换,需要确定七个参数的数值。

- 第四步,验证和优化转换。

创建自定义地理坐标转换后,需进行验证和优化。

可以使用样本数据进行验证,检查转换效果和精度。

如果发现转换效果不理想,可以适当调整参数值和方法,进行优化。

3. 自定义地理坐标变换的应用与价值自定义地理坐标变换在实际应用中具有广泛的应用和价值。

它可以满足不同坐标系统之间数据转换的需求。

当我们需要将一个国家的坐标数据转换到另一个国家的坐标系统时,可以通过创建自定义地理坐标变换来实现。

在ENVI中添加beijing54坐标系

在ENVI中添加beijing54坐标系

如何在ENVI里自定义高斯投影(北京54标)与转换投影baicai0114一.添加北京54坐标通过修改ENVI相应的系统文档,可以将北京54坐标系添加进ENVI软件系统,具体方法和步骤如下:1、找到ENVI安装目录下相应的坐标系存储文档,如“C:\RSI\IDL63\products\envi43\map_proj\datum.txt”,并打开。

2、将以下描述语句添加在文档最后一行:“Beijing-54, Krassovsky, -12, -113, -41”,其中“Beijing-54”是新坐标的名称;“Krassovsky”是它所使用的椭球体名称;“-12, -113, -41”是新坐标系与WGS-84坐标系的坐标差。

3、将修改后的文档保存,关闭。

4、重新启动ENVI。

新的坐标系添加成功。

二、自定义高斯-克吕格(Kauss Kruger)投影具体方法和步骤如下:1.选择Map> Customize Map Projection,或在任何地图投影选择对话框中,点击“New”按钮。

2、当出现Customized Map Projection Definition 对话框时,假设为定义高斯投影的13带,13带的中央经度为75度,具体设置参数如下(上图所示):Projection name :GK Zone 13(Beijing-54)Projection Type:Transverse MercatorProjection Datum:Beijing-54False easting:500000False northing:0Latitude of projection origin: 0Longitude of central meridian: 75Scale factor:13.选择Projection > Add New Projection将投影添加到ENVI 所用的投影列表中。

ENVI地图调整

ENVI地图调整

5.1 地形图的调整5.1.1 ENVI中定义坐标系(1)地形图坐标系本地图采用西安80坐标系,相关参数如下:(2)修改ENVI中的坐标定义文件ENVI中的坐标定义文件存放在 \ITT\IDL708\products\envi46\map_proj文件夹下。

添加椭球体:把“IAG-75,6378140.0,6356755.3”(西安80椭球体)加入ellipse.txt末端。

添加基准面:把“Xi’an-80,IAG-75,0,0,0”加入datum.txt末端。

添加坐标系:打开ENVI4.6。

选择Map->Customize Map Projections, 根据地图填入参数。

从地图上可看出,本地图位于39度带,采用3度分带投影,中央经线为东经117度。

填入参数如图1:图1 添加坐标系选择Projection -> Add New Projection,将投影添加到ENVI所有的投影中。

再选择File->Save Projections,保存添加好的投影。

(3)为地图添加Map Info打开地图3.tif。

选择Map->Registration->Select GCPs:Image To Map。

选择之前创建好的投影 Xi’an-80 3 Degree GK CM117E,并把X Pixel Size和Y Pixel Size改为1米。

单击OK。

在地图左上角选择一点,并输入这点经纬度,如图2。

图2 选点并输入该点坐标然后单击Add Point。

再重复上述选点方法选则2个点。

所选点坐标如图3。

图3 第二个点坐标图4 第三个点坐标选择Options->Warp File。

在弹出的窗口中选择3.tif,单击OK。

保存成新的文件,命名为3-1。

(4)调整地图地理坐标打开3-1.tif。

a.选择Overlay->Grid Line。

因为本图经度跨度3分45秒,维度跨度2分30秒,所以设置格网间隔1分15秒,使格网与地图四角经纬度数相一致。

ENVI基本影像处理流程操作ppt课件

ENVI基本影像处理流程操作ppt课件
7
1.1ENVI简介——立体像对高程提取扩展模块— DEM Extraction
• 快速从ALOS PRISM, ASTER, CARTOSAT-1, FORMOSAT-2,
GeoEye-1, IKONOS, KOMPSAT-2, OrbView-3, QuickBird, WorldView-1, SPOT 1-5等以及航空影像立体像对中提取 DEM。
3
1.快速认识ENVI
•1.1 ENVI简介 •1.2 安装目录结构 •1.3 栅格文件系统和储存 •1.4 数据输入 •1.5 数据显示 •1.6 常见系统设置
4
1.1ENVI简介——ENVI/IDL体系结构
扩展模块 主模块 开发语言
Atmospheric Correction
大气校正模 块
Feature Extraction
5
1 55 2 3-5 3
1.1-3.7 26 2-3
2每5天
其他卫星
传感器 Landsat1~7
发射时间 72~99
SPOT4
1999
中巴资源卫星-01/02 1999
13
1.3栅格文件系统和储存
• ENVI栅格文件格式:ENVI使用的是通用栅格数据格式,包含一
个简单的二进制文件( a simple flat binary )和一个相关 的ASCII(文本)的头文件。
–ENVI头文件包含用于读取图像数据文件的信息,它通常创建于一个数据文件第一次被 ENVI读取
时。单独的ENVI头文本文件提供关于图像尺寸、嵌入的头文件(若存在)、数据格式及其它相 关信息。所需信息通过交互式输入,或自动地用“文件吸取”创建,并且以后可以编辑修改。 您可以在ENVI之外使用一个文本编辑器生成一个ENVI头文件

15.如何在ENVI里快速添加高斯投影(北京54坐标)yang1026

15.如何在ENVI里快速添加高斯投影(北京54坐标)yang1026

如何在ENVI里快速添加高斯投影(北京54坐标)一.添加北京54坐标通过修改ENVI相应的系统文档,可以将北京54坐标系添加进ENVI软件系统,具体方法和步骤如下:1、找到ENVI安装目录下相应的坐标系存储文档,如“C:\RSI\IDL63\products\envi43\map_proj\datum.txt”,并打开。

2、将以下描述语句添加在文档最后一行:“Beijing-54, Krassovsky, -12,-113, -41”,其中“Beijing-54”是新坐标的名称;“Krassovsky”是它所使用的椭球体名称;“-12, -113, -41”是新坐标系与WGS-84坐标系的坐标差。

3、将修改后的文档保存,关闭。

二、添加高斯-克吕格(Kauss Kruger)投影(6度带)具体方法和步骤如下:1、找到ENVI安装目录下相应的坐标系存储文档,如“C:\RSI\IDL63\products\envi43\map_proj\map_proj.txt”,并打开。

2、将以下描述语句添加在文档最后一行:3, 6378245.0, 6356863.0, 0.000000, 75.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone13 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 81.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone14 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 87.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone15 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 93.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone16 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 99.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone17 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 105.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone18 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 111.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone 19(Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 117.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone20 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 123.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone21 (Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 129.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone 22(Beijing-54)3, 6378245.0, 6356863.0, 0.000000, 135.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone23 (Beijing-54)他们是高斯-克吕格(Kauss Kruger)投影的13-23带相关参数,分别是:序号,椭球体长短半轴,向北偏移量,纬度,中央经度,向东偏移,坐标名称,投影名称。

ENVI基本操作之自定义坐标

ENVI基本操作之自定义坐标

IAG-75,6378140.0,6356755.3
影像预处理
分别是:基准面名称、 椭球名称、和三个参数。 都用半角豆号隔开。
影像预处理
这就是加进去 的椭球参数
影像预处理
• 第二步:定义基准面。 • 将基准面参数拷贝到datum.txt 文件夹中。 • 同样是半角豆号隔开。
附:北京54和西安80的基本面参数。 Beijing-54,Krasovsky,-12,-113,-41 Xi′an-80,IAG-75,0,0,0
影像预处理
• 第一步:定义椭球体。打开安装目录下
影像预处理
• 点开ellipse.txt,有三个参数:
椭球体名称
短半轴
长半轴
影像预处理
• 要增加椭球体和参数可直接在这个文件夹中
加上。这里要注意,名称和参数都是半角豆 号。
附:自定义投影参数。可选用 Krasovsky,6378245.0,6356863.0
选择已有的北 京54坐标基准 面
影像预处理
坐标偏移,这保证X 为正值,将中央经线 向东偏移500km 南偏为0 中央纬线为0
比例系数,0.9996 是个经验值。也可输 入1。
中央经线根据自己 的图幅位置而定
注意输入值后要回车以表示确定。
影像预处理
上述完成 后,将定 义的坐标 参数确定 下来。
影像预处理
ENVI基本操作之
自定义坐标
影像预处理
自定义坐标系
• 坐标系原理:地理坐标系是以经纬度为单位
的地球坐标系统,地理坐标系中有2个重要 部分,即地球椭球体(Spheroid)和大地基 准面(datum)。
大地基准面指目前参考椭球与WGS84参考椭球间 的相对位置关系(3个平移,3个旋转,1个缩 放),可以用其中3个、4个或者7个参数来描述 它们之间的关系,每个椭球体都对应一个或多 个大地基准面。
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

ENVI和ERDAS中自定义坐标系的方法1 地理投影的基本原理常用到的地图坐标系有2种,即地理坐标系和投影坐标系。

地理坐标系是以经纬度为单位的地球坐标系统,地理坐标系中有2个重要部分,即地球椭球体(spheroid)和大地基准面(datum)。

由于地球表面的不规则性,它不能用数学公式来表达,也就无法实施运算,所以必须找一个形状和大小都很接近地球的椭球体来代替地球,这个椭球体被称为地球椭球体,我国常用的椭球体如下大地基准面指目前参考椭球与WGS84参考椭球间的相对位置关系(3个平移,3个旋转,1个缩放),可以用其中3个、4个或者7个参数来描述它们之间的关系,每个椭球体都对应一个或多个大地基准面。

投影坐标系是利用一定的数学法则把地球表面上的经纬线网表示到平面上,属于平面坐标系。

数学法则指的是投影类型,目前我国普遍采用的是高斯——克吕格投影,在英美国家称为横轴墨卡托投影(Transverse Mercator)。

高斯克吕格投影的中央经线和赤道为互相垂直,分带标准分为3度带和6度带。

美国编制世界各地军用地图和地球资源卫星像片所采用的全球横轴墨卡托投影(UTM)是横轴墨卡托投影的一种变型。

高斯克吕格投影的中央经线长度比等于1,UTM投影规定中央经线长度比为0.9996。

我国规定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年北京坐标系。

为了适应大地测量的发展,我国于1978年采用国际大地测量协会推荐的IAG-75地球椭球体建立了我国新的大地坐标系,并在1986年宣布在陕西省泾阳县设立了新的大地坐标原点,由此计算出来的各大地控制点坐标,称为1980年大地坐标系。

对于地理坐标,只需要确定两个参数,即椭球体和大地基准面。

对于投影坐标,投影类型为Gauss Kruger(Transverse Mercator),除了确定椭球体和大地基准面外,还需要确定中央经线。

大地基准面的确定关键是确定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*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 --------基准椭球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)选择Projection>Load 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.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 meridian:1.000000Longitude of central meridian:117:00:00.000000000000 ELatitude of origin of projection:0:00:00.000000000000 NFalse easting:39500000.00000000000 metersFalse northing:0.00000000000000000 meters。

相关文档
最新文档