15.如何在ENVI里快速添加高斯投影(北京54坐标)yang1026
- 格式:doc
- 大小:34.00 KB
- 文档页数:2
ENVI实验2:影像地理坐标定位和配准主要介绍在ENVI中对影像进行地理校正,添加地理坐标,以及如何使用ENVI进行影像到影像的配准和影像到地图的校正。
一、ENVI中带地理坐标的影像ENVI对带地理坐标的影像提供了全面的支持,它能够对许多预定义的地图投影进行处理,这些地图投影可以采用UTM或State Plane投影方式。
此外,ENVI的用户自定义地图投影功能能够创建自定义的地图投影,它允许使用6种基本投影类型,超过35种的不同的椭球体以及100多种的基准数据集(Datum)来满足大多数地图投影的需要。
ENVI的地图投影参数存储在ASCII码文本文件Map_proj.txt中,该文本文件能够被ENVI地图投影工具修改,或者直接被用户编辑。
这个文件中的信息会被影像相应的头文件所使用,而且ENVI允许使用已知的地图投影坐标来简单地指定相Magic Pixel(地图坐标系统的起始点)。
然后,选择的ENVI函数就能够使用该信息,在带地理坐标的数据空间中进行操作处理。
ENVI的影像配准和几何纠正工具允许用户将基于像素的影像定位到地理坐标上,然后对它们进行几何纠正,使其匹配基准影像的几何信息。
使用全分辨率(主影像窗口)和缩放窗口来选择地面控制点(GCPs),进行影像到影像和影像到地图的配准。
基准影像和未校正影像的控制点坐标都会显示出来,同时由指定的校正算法所得的误差也会显示出来。
地面控制点预测功能能够使对地面控制点的选取简单化。
将使用重采样、缩放比例和平移(PST),以及多项式函数(多项式系数可以从1到n)或者Delaunay三角网的方法,来对影像进行校正。
所支持的重采样方法包括最近邻法(nearest-neighbor)、双线性内插法(bilinear interpolation)和三次卷积法(Cubic convolution)。
使用ENVI的多重动态链接功能对基准影像和校正后的影像进行比较,可以快速评估配准的精度。
如何在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,全称为Environment,是一种遥感图像处理软件。
以下是ENVI软
件的基本操作手册:
1. 打开影像:在ENVI中,可以通过File菜单的Open选项或者直接点击工具栏上的Open按钮来打开影像。
在弹出的对话框中选择要打开的影像文件,然后点击OK即可。
2. 显示影像:打开影像后,可以通过点击工具栏上的Display按钮来显示影像。
在弹出的对话框中可以选择显示方式、颜色方案、透明度等参数,然后点击OK即可。
3. 调整影像大小:可以通过工具栏上的Zoom和Pan按钮来调整影像的大
小和位置。
Zoom按钮可以放大或缩小影像,Pan按钮可以平移影像。
4. 创建ROI(感兴趣区域):在ENVI中,可以通过工具栏上的ROI按钮来创建感兴趣区域。
在创建ROI时,可以选择不同的形状、大小和位置,并
且可以在ROI上添加标签和注解。
5. 提取光谱信息:在ENVI中,可以通过Spectral Analysis工具来提取光
谱信息。
可以选择不同的光谱分析方法,如光谱曲线、光谱角、光谱匹配等,并可以自定义波段和阈值等参数。
6. 生成图像地图:在ENVI中,可以通过Map工具来生成图像地图。
可以
选择不同的地图投影和坐标系,并可以添加图层、标注、符号等元素。
7. 导出数据:在ENVI中,可以将处理后的数据导出为多种格式,如TIFF、JPEG、BMP等。
在导出数据时,可以选择导出的范围、数据类型、分辨率等参数,并可以设置输出文件的格式和质量。
以上是ENVI软件的基本操作手册,希望对您有所帮助。
ENVI和ERDAS中自定义坐标系的方法什么是80西安坐标系?1978年4月在西安召开全国天文大地网平差会议,确定重新定位,建立我国新的坐标系。
为此有了1980年国家大地坐标系。
1980年国家大地坐标系采用地球椭球基本参数为1975年国际大地测量与地球物理联合会第十六届大会推荐的数据。
该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇,位于西安市西北方向约60公里,故称1980年西安坐标系,又简称西安大地原点。
基准面采用青岛大港验潮站1952-1979年确定的黄海平均海水面(即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下影像处理详细步骤:1、坐标转换1)定义投影:提交的原始影像为经纬度坐标,应用ArcGIS中的投影定义工具定义影像坐标系,如xian80或北京54;2)投影转换:使用ArcGIS高级工具箱中的Project Raster工具对影像进行投影转换,得到系统需要的平面坐标成果。
如图1中所示,重采样选用BILINEAR法。
图12、影像裁切影像经坐标转换后,原矩形影像将会变成其它不规则图形,需要进行影像裁剪再得到矩形影像,裁切方法如下:ENVI软件下裁切:1)打开影像:Image File,打开要裁切的影像;2)裁切:Basic Tools—Resize Data (Spatial/Spectral),打开Resize Data Input File对话框,如图2所示。
选中要裁切的影像,点击Spatial Subset按钮,弹出Select Spatial Subset对话框,如图3所示。
然后点击Image按钮,选择按影像裁切方式,弹出Subset by Image对话框,见图4,用红框画定要保留的影像范围。
接下来点三次OK分别确定三个对话框的设置,在最后弹出的Resize Data Parameters对话框中设置重采集方式“Bilinear”,点击Choose按钮确定输出影像的路径和名称后点OK确定即运行影像载切。
图2 对话框Resize Data Input File图3 对话框Select Spatial Subset图4 Subset by Image对话框图5 Resize Data Parameters对话框3、格式转换应用ENVI裁切输出的影像为ENVI的默认格式,需将其转换成ArcGIS能应用的*.img格式。
方法:ENVI下—ERDAS IMAGING将影像另存为*.img格式即可。
注;格式转换后再在ArcGIS下重新定义一下投影。
与地图投影有关的软件使用说明1.建立自定义的地图投影菜单位置:Map——Customized Map ProjectionENVI 软件中包含了一些标准的地图投影(记录在文件"map_proj.txt"中,该文件一般在map_proj子目录下,例如Universal Transverse Mercator (UTM和高斯克吕格投影。
通过这一菜单命令,用户不仅利用已知的投影,而且可以建立自定义的投影。
一般情况下,需要输入的投影参数包括:a投影类型(Transverse Mercator, Lamberts Conformal Conic, Lambert Azimuthal EqualArea, Oblique Mercator (A and B, Stereographic, Albers Equal Area, and Polyconic;b椭球体及其参数(Ellipsoid, Datum,主要是长半轴和短半轴。
c经度和纬度方向上的偏置(False Easting and Northingd起始纬度e中央经线f放缩比例因子等注意:这个对话框根据用户选择的投影类型不同,具有不同的内容。
如要了解地图投影更加详尽的信息,请访问以下INTERNET站点/depts/grg/gcraft/notes/datum/datum.html/jbanta/根据已有投影的修改是一个比较常用的方法。
首先装载标准投影,对话框的菜单命令Projection >Load existing project,该投影的参数自动出现在编辑框中,用户可以根据自己的需要修改这些参数。
包括输入“投影名称”等。
需要注意的是,不同地图投影所需要的参数是不同的。
ENVI 软件会根据投影的类别,自动显示和隐藏参数的编辑框。
例如横轴墨卡托投影的参数有中央经线、基准经线等;用户自定义椭球体时,需要输入定义椭球体的长半轴和短半轴。
ENVI中如何建立北京54或西安80的地图投影ENVI中如何建立北京54或西安80的地图投影?首先要资料收集,比如西安80的椭球体参数等。
在map里自定义投影,可以选择系统定义的投影参数加以修改。
如系统中GK投影,用的是普尔科沃1942坐标系,克拉索夫斯基椭球体,这就是北京1954坐标系。
不过它的平移参数为18500000,加了带号。
我要与mapgis 的图框相套合,所以我改成了500000米。
西安中的椭球体是IUGG1975推荐椭球体,可以查它的半径。
ENVI 里只能输入小数点后1位(还是2位,我记不清了),其它的设置好像差不多一样.平移,中央经度,投影带号等.RSI\IDL60\products\envi40\map_proj目前的遥感处理软件中并没有我国常用的北京54坐标系,在ENVI中涉及了2个文件--datum.txt和map_proj.txt。
datum.txt文件中写入基准面信息,map_proj.txt文件中写入各投影带的具体参数。
datum.txt:Gauss-Kruger (Beijing1954 3-degree zone), Krassovsky, 0, 0, 0Gauss-Kruger (Beijing1954 6-degree zone), Krassovsky, 0, 0, 0map-proj.txt:3, 6378245.0, 6356863.0, 0.000000, 63.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 11 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 69.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 12 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 75.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 13 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 81.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone14 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 87.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 15 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 93.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 16 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 99.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 17 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 105.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 18 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 111.000000, 500000.0, 0.0, 1.000000,3, 6378245.0, 6356863.0, 0.000000, 117.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 20 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 123.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 21 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 129.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 22 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 135.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 6-degree zone), GK Zone 23 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 75.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 25 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 78.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 26 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 81.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 27 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 84.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 28 (Beijing1954) 3,6378245.0, 6356863.0, 0.000000, 87.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 29 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 90.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 30 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 93.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 31 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 96.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 32 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 99.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 33 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 102.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 34 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 105.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 35 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 108.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 36 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 111.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 37 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 114.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 38 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 117.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 39 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 120.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 40 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 123.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone41 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 126.000000, 500000.0, 0.0, 1.000000,3, 6378245.0, 6356863.0, 0.000000, 129.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 43 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 132.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 44 (Beijing1954) 3, 6378245.0, 6356863.0, 0.000000, 135.000000, 500000.0, 0.0, 1.000000, Gauss-Kruger (Beijing1954 3-degree zone), GK Zone 45 (Beijing1954)。
如何在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 Zone
13 (Beijing-54)
3, 6378245.0, 6356863.0, 0.000000, 81.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone
14 (Beijing-54)
3, 6378245.0, 6356863.0, 0.000000, 87.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone
15 (Beijing-54)
3, 6378245.0, 6356863.0, 0.000000, 93.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone
16 (Beijing-54)
3, 6378245.0, 6356863.0, 0.000000, 99.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone
17 (Beijing-54)
3, 6378245.0, 6356863.0, 0.000000, 105.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone
18 (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 Zone
20 (Beijing-54)
3, 6378245.0, 6356863.0, 0.000000, 123.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK6 Zone
21 (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 Zone
23 (Beijing-54)
他们是高斯-克吕格(Kauss Kruger)投影的13-23带相关参数,分别是:序号,椭球体长短半轴,向北偏移量,纬度,中央经度,向东偏移,坐标名称,投影名称。
3、将修改后的文档保存,关闭。
4、重新启动ENVI。
新的坐标系添加成功。
三、添加高斯-克吕格(Kauss Kruger)投影(3度带)
3度带的添加方法和6度带的添加方法基本一样,只用改变每一个带的中央经线即可。
如25带中央经线是75度。
添加内容如下:
3, 6378245.0, 6356863.0, 0.000000, 75.000000, 500000.0, 0.0, 1.000000, Beijing-54, GK3 Zone
25 (Beijing-54)
其他带依此类推。
best wishes!
--------------
仰满荣 (Miss Yang)。