ENVI实验报告

1、一、实验目的ENVI 是一套功能齐全的遥感图像处理系统，是处理、分析并显示多光谱数据、高光谱数据和雷达数据的高级工具。此次实习主要是学习一些关于ENVI 的基本操作，如：图像预处理，影像分析，图像增强，几何校正，监督分类以及专题制图等步骤。二、实验数据LE71440292000268SGS00.tar.gzELEVATION_SOURCE = "GLS2000"PROCESSING_SOFTWARE = "LPGS_9.1"EPHEMERIS_TYPE = "DEFINITIVE"SPACECRAFT_ID = "Lands

2、at7"SENSOR_ID = "ETM+"SENSOR_MODE = "SAM"ACQUISITION_DATE = 2000-09-24WRS_PATH = 144BAND_COMBINATION = "123456678"PRODUCT_UL_CORNER_LAT = 45.5786828PRODUCT_UL_CORNER_LON = 84.0750064PRODUCT_UR_CORNER_LAT = 45.6157964PRODUCT_UR_CORNER_LON = 87.2821725PRODUCT_LL_CORN

3、ER_LAT = 43.5718357PRODUCT_LL_CORNER_LON = 84.1739972PRODUCT_LR_CORNER_LAT = 43.6064525PRODUCT_LR_CORNER_LON = 87.2726073PRODUCT_UL_CORNER_MAPX = 271800.000PRODUCT_UL_CORNER_MAPY = 5051400.000PRODUCT_UR_CORNER_MAPX = 522000.000PRODUCT_UR_CORNER_MAPY = 5051400.000PRODUCT_LL_CORNER_MAPX = 271800.000PR

4、ODUCT_LL_CORNER_MAPY = 4828200.000PRODUCT_LR_CORNER_MAPX = 522000.000PRODUCT_LR_CORNER_MAPY = 4828200.000三、实验内容1.遥感数据下载2.数据格式转化及多波段的合成3.图像边框的裁剪4.图像的镶嵌6.图像几何校正7.裁剪8.监督分类9.专题出图四、实验步骤1.遥感数据下载输入地名及日期，搜索符合自己通途的数据，查看并下载2.数据格式转化及多波段的合成由于下载图像的格式tif 格式，需要转成img 格式。并将单波段融合成为多波段，这一步在erdas中完成。tif 转 img单波段融合成多波段3

5、.图像边框的裁剪通过主图像窗口 -tools-region of interest-roi tool ，然后自己在图像上定义出一块感兴趣区 （实际中难于确定严格准确的裁剪界限）。4.图像的镶嵌Map mosaicking georeferenced启动镶嵌模板6.图像几何校正设置参数，参数在原始数据中的txt 中可以查看启动校正模块，并利用google地球进行校正7.裁剪分别打开需要裁剪的图像和所在区域的shapefile 文件将矢量文件转成ROI利用 ROI 裁剪图像选择主菜单 basic tools subset data via ROIS，选择裁剪图像在 spatial subset v

6、ia ROI parameters中，选择由矢量生成的ROI，在“ Mask pixels outside of ROI ？”项中选择 yes，Mask BackgroundValue 为 08.监督分类打开裁剪后的图像， 以 RGB ：543 显示，在 display 中选择 Overlay Region of interest样本选择完成后，在Option Compute ROI Separability ，选择分类图像，查看分离性Input File:裁剪ROI Name: (Jeffries-Matusita, Transformed Divergence)水体 Blue 36084

7、points:农田 Green 18122 points: (2.00000000 2.00000000)沙漠 Sienna 200333 points: (1.99973178 2.00000000)未利用土地Orchid 57805 points: (1.99999494 2.00000000)山体 Purple 355228 points: (1.98296692 2.00000000)农田 Green 18122 points:水体 Blue 36084 points: (2.00000000 2.00000000)沙漠 Sienna 200333 points: (1.9995482

8、0 1.99999924)未利用土地Orchid 57805 points: (2.00000000 2.00000000)山体 Purple 355228 points: (1.98744710 1.99928798)沙漠 Sienna 200333 points:水体 Blue 36084 points: (1.99973178 2.00000000)农田 Green 18122 points: (1.99954820 1.99999924)未利用土地Orchid 57805 points: (1.98700177 2.00000000)山体 Purple 355228 points: (

9、1.84635561 1.96524858)未利用土地Orchid 57805 points:水体 Blue 36084 points: (1.99999494 2.00000000)农田 Green 18122 points: (2.00000000 2.00000000)沙漠 Sienna 200333 points: (1.98700177 2.00000000)山体 Purple 355228 points: (1.99987319 2.00000000)山体 Purple 355228 points:水体 Blue 36084 points: (1.98296692 2.000000

10、00)农田 Green 18122 points: (1.98744710 1.99928798)沙漠 Sienna 200333 points: (1.84635561 1.96524858)未利用土地Orchid 57805 points: (1.99987319 2.00000000)Pair Separation (least to most);沙漠 Sienna 200333 points and山体 Purple 355228 points - 1.84635561水体 Blue 36084 points and山体 Purple 355228 points - 1.9829669

11、2沙漠 Sienna 200333 points and未利用土地Orchid 57805 points - 1.98700177农田 Green 18122 points and山体 Purple 355228 points - 1.98744710农田 Green 18122 points and沙漠 Sienna 200333 points - 1.99954820水体 Blue 36084 points and沙漠 Sienna 200333 points - 1.99973178未利用土地Orchid 57805 points and山体 Purple 355228 points -

12、 1.99987319水体水体农田Blue 36084 points andBlue 36084 points andGreen 18122 points and未利用土地Orchid 57805 points - 1.99999494农田 Green 18122 points - 2.00000000未利用土地Orchid 57805 points - 2.00000000分类统计结果分类精度评价混淆矩阵Confusion Matrix: H:123456jianduOverall Accuracy = (646026/667344) 96.8055%Kappa Coefficient =

13、0.9486Ground Truth (Pixels)Class水体农田沙漠未利用土地山体Unclassified00000水体 Blue 335922000286农田 Green751790748019沙漠 Sienna346119581353615580未利用土地 O402166572690山体 Purple4915423060339115Total360841812220033357805355000Ground Truth (Pixels)ClassTotalUnclassified0水体 Blue 336208农田 Green18049沙漠 Sienna212024未利用土地 O59

14、439山体 Purple341624Total667344Ground Truth (Percent)Class水体农田沙漠未利用土地山体Unclassified0.000.000.000.000.00水体 Blue 399.550.000.000.000.08农田 Green0.2198.810.020.000.01沙漠 Sienna0.090.3497.740.934.39未利用土地 O0.010.001.0899.070.00山体 Purple0.140.851.150.0095.53Total100.00100.00100.00100.00100.00Ground Truth (Per

15、cent)ClassTotalUnclassified0.00水体 Blue 35.43农田 Green2.70沙漠 Sienna31.77未利用土地 O8.91山体 Purple51.19Total100.00Class CommissionOmissionCommissionOmission(Percent) (Percent)(Pixels)(Pixels)水体 Blue 30.790.45286/36208162/36084农田 Green0.791.19142/18049215/18122沙漠 Sienna7.652.2616211/2120244520/200333未利用土地 O3

16、.650.932170/59439536/57805山体 Purple0.734.472509/34162415885/355000Class Prod. Acc. User Acc.Prod. Acc.UserAcc.(Percent) (Percent)(Pixels)(Pixels)水体 Blue 399.5599.2135922/3608435922/36208农田 Green98.8199.2117907/1812217907/18049沙漠 Sienna97.7492.35195813/200333195813/212024未利用土地 O99.0796.3557269/578055

17、7269/59439山体 Purple95.5399.27339115/355000339115/3416249.专题出图1、虚拟边框设置1）在主显示窗口菜单栏中选择File Preferences，打开Display Parameters对话框，设置虚拟边框的边界值和颜色。2）点击 OK完成虚拟边框的设置。2、公里网设置ENVI支持同时显示像素公里网、地图坐标公里网以及地理坐标（纬度 /经度）网。添加或修改地图影像公里网：1）主显示窗口菜单栏中选择Overlay Grid Lines。打开 Grid Line Parameters对话框中，显示默认的公里网设置。2）设置公里网属性参数。选择

18、Options Edit Map Grid Attributes 或Edit GeographicAttributes获 Edit Pixel Attributes 修改所选公里网的属性3）点击 OK，完成参数设置。4）在Grid Line Parameters对话框中，点击 Apply 将新的公里网应用到影像图中。3、注记要素操作1）在主显示窗口菜单栏中选择Overlay Annotation，打开Annotation对话框。2）从 Annotation对话框菜单栏的 Object下拉菜单中选择所需的注记要素。3）选择 Image单选按钮指定注记放置的窗口。4）在主显示窗口鼠标左键点击注记要素放置的位置，点击鼠标右键锁定注记的位置。5）编辑注记要素。a、移动注记要素菜单中选择 Object Selection/Edit，用鼠标左键拖画出一个矩形框选择待移动的注记要素。 通过点击小圆柄并拖放可以重新设置注记要素新位置。b、修改注记要素的属性菜单中选择 Object Selection/Edit，用鼠标左键拖画出一个矩形