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

ENVI基本影像处理流程操基本影像处理流程操作作邓书斌E-mail:dengsbesrichina-ESRI中国（北京）有限公司 遥感事业部主要内容主要内容 1.快速认识ENVI 2.影像一般预处理 3.影像信息基本提取方法 4.制图与三维可视化数据的输入输出图像显示与分析图像预处理（几何纠正、融合、镶嵌等）影像信息提取（人工解译、自动分类、特征提取、动态检测等）专题制图/三维可视化分析（集成GIS现有数据)成果报告（GIS分析/共享）遥感图像处理的一般流程遥感图像处理的一般流程1.快速认识ENVI1.1 ENVI简介1.2 安装目录结构1.3 栅格文件系统和储存1.4 数据输入1.5 数据显示1.6 常见系统设置Feature Extraction空间特征提取模块DEMExtraction立体像对高程提取模块Atmospheric Correction大气校正模块1.1ENVI简介ENVI/IDL体系结构ENVIIDL 扩展模块 主模块 开发语言IDL Advanced 数学与统计扩展 工具包IDL DataMiner 数据库连接工 具包Orthorectification正射校正模块NITFNITF数据支持模块1.1ENVI简介ENVI从图像中获得您所需的信息Image-InformationEXTEND IDL,C+,Java,.NetPanchromatic 全色 Multispectral 多光谱 Hyperspectral 高光谱TRadarLidarhermal TerrainGPSVector雷达激光雷达热量数据地 形 数 据 位 置 数 据 矢量数据制图工具矢量工具地形工具高光谱工具信息提取工具预处理工具雷达工具读取显示工具三维可视分析制图输出GeodatabaseGIS分析1.1ENVI简介简介大气校正扩展模块大气校正扩展模块(Atmospheric Correction)可 以 有 效 地 去 除 水 蒸 气可 以 有 效 地 去 除 水 蒸 气,气溶胶散射，漫反射的邻域效气溶胶散射，漫反射的邻域效应。应。获 得 地 物 反 射 率 和 辐 射 率、地 表 温 度 等 真 实 物 理 模 型 参 数。获 得 地 物 反 射 率 和 辐 射 率、地 表 温 度 等 真 实 物 理 模 型 参 数。使影像变得使影像变得“清晰清晰”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。全面支持RPC模型参数，尽可能少的控制点以达到有效的精度 使用DEM编辑工具对提取的DEM做局部编辑 交互量测特征地物的高度和收集3D信息并导出为3D Shapefile文件格式提供面向对象方法、易于使用的向导操作流提供面向对象方法、易于使用的向导操作流程从高分辨程从高分辨 率全色和多光谱数据中提取地物率全色和多光谱数据中提取地物信息。包括：信息。包括：交 通 工 具交 通 工 具飞机飞机,坦克坦克,汽车汽车,船只船只 建 筑 物建 筑 物建筑物轮廓建筑物轮廓,屋顶屋顶 基 础 设 施基 础 设 施道路道路,桥梁桥梁,机场机场,海港码头海港码头 自 然 要 素自 然 要 素河流河流,湖泊湖泊,森林森林,田地田地 云和雾云和雾1.1ENVI简介简介空间特征提取扩展模块空间特征提取扩展模块Feature Extraction（FX）1.1ENVI简介简介正射校正扩展模块正射校正扩展模块Orthorectification 由 瑞 典 的由 瑞 典 的Spacemetric公司开发公司开发 采用的正射校正方法具有可靠和高精度的特采用的正射校正方法具有可靠和高精度的特点，并且该方点，并且该方 法被行业所认可。法被行业所认可。支 持 大 区 域 范 围 内 的 多 幅 影 像、多 传 感 器 的 一 次 正 射 校 正支 持 大 区 域 范 围 内 的 多 幅 影 像、多 传 感 器 的 一 次 正 射 校 正 具 有 镶 嵌 结 果 的 功 能，并 提 供 接 边 线 和 颜 色 平 衡 辅 助 工 具具 有 镶 嵌 结 果 的 功 能，并 提 供 接 边 线 和 颜 色 平 衡 辅 助 工 具 采 用 流 程 化 的 向 导 式 操 作 方 式 和 工 程 化 管 理。采 用 流 程 化 的 向 导 式 操 作 方 式 和 工 程 化 管 理。自 定 义 传 感 器 模 型自 定 义 传 感 器 模 型 提供接口函数，便于扩展功能。提供接口函数，便于扩展功能。1.1ENVI简介简介NITF数据支持扩展模块 NITF 读 写、显 示 标 准读 写、显 示 标 准NITF格式文件格式文件 JPEG2000编 码 压 缩编 码 压 缩NITF格式文件格式文件 支 持支 持NITF2.0、NITF2.1和和NSIF1.0之间的转换之间的转换 读 写 从 商 业 卫 星、读 写 从 商 业 卫 星、NCDRD和和 第二图像格式（第二图像格式（NSIF）中获得中获得政府标准数据政府标准数据 广 泛 支 持广 泛 支 持NSDE的分类或未分类的的分类或未分类的TREs，也包括，也包括自定义的自定义的TREs1.1ENVI简介简介ENVI/IDL的特点的特点1.先进、可靠的影像信息提取工具全套影像信息智 能化提取工具,全面提升影像的价值。2.专业的光谱分析高光谱分析一直处于世界领先地 位。3.随心所欲扩展新功能底层的IDL语言可以帮助用户 轻松地添加、扩展ENVI的功能，甚至开发定制自己的 专业遥感平台。4.流程化图像处理工具ENVI将众多主流的图像处理 过程集成到流程化（Workflow）图像处理工具中，进 一步提高了图像处理的效率。5.与ArcGIS的整合从2007年开始，与ESRI公司的全 面合作,为遥感和GIS的一体化集成提供了一个典型的 解决方案。1.2安装目录结构安装目录结构 Bin：相应的ENVI运行目录。Data：数据目录，保存一矢量文件夹（一些矢量数据）、两个TM5栅格 数据、两个DEM数据和一个高光谱数据。Filt_func：ENVI常规传感器的光谱库文件。例如：aster、modis、spot、tm等。Help：ENVI的帮助文档。Lib：IDL生成的可编译的程序，用于二次开发。Map_proj：影像的投影信息，文本格式，客户可以进行定制。Menu：ENVI菜单文件，可以进行中、英文菜单互换。Save：应用IDL可视化语言编译好的、可执行的ENVI程序。Save_add：客户自主开发的、可执行程序，比如各种补丁程序。Spec_lib：波谱库，不同地区可以有不同的波谱库，用户可以自定义。1.3栅格文件系统和储存栅格文件系统和储存 ENVI栅格文件格式：ENVI使用的是通用栅格数据格式，包含一 个简单的二进制文件（a simple flat binary）和一个相关 的ASCII（文本）的头文件。ENVI头文件包含用于读取图像数据文件的信息，它通常创建于一个数据文件第一次被 ENVI读取 时。单独的ENVI头文本文件提供关于图像尺寸、嵌入的头文件（若存在）、数据格式及其它相 关信息。所需信息通过交互式输入，或自动地用“文件吸取”创建，并且以后可以编辑修改。您可以在ENVI之外使用一个文本编辑器生成一个ENVI头文件通用栅格数据都会存储为二进制的字节流，通常它将以BSQ（按波段顺序）、BIP（波段按像元交叉）或者BIL（波段按行交叉）的方式进行存储。储存窗口菜单界面FileSave File As，将影像按照需要的格式进行存储，保存的为原始数据，没有拉伸。主影像窗口FileSave Images As，将影像按照需要的格式进行存储，存储的影像是显示的影像样式。FileSave Zoom As，将Zoom窗口显示的影像按照需要的格式进行存储。其他窗口下的文件存储例如：MapMosaicking的镶嵌窗口下：Apply；Save Template等；Classification等功能下：Output Result to等。1.4 数据显示数据显示 波段列表每次打开的文件都显示在Available Bands List中，列表中可以完成当前 在ENVI中打开的或存储在内存中的文件的信息，还可以进行包括：打开新 文件、关闭文件、将内存数据项保存到磁盘，以及编辑ENVI头文件等操作。三视窗显示当你打开一个图像文件时，会在一个ENVI的三视窗图像显示中，其中包括 主图像窗口，缩放窗口和滚动窗口（应用于大的图像），如图1.5所示，目 前大部分的ENVI图像处理操作都在这个窗口中完成。ENVI ZOOM显示将图层管理、图像显示、鼠标信息等集中在一个窗体中，目前只有部分 ENVI图像处理操作在这个窗口中完成，如面向对象的特征提取、Pan sharping、异常检测等，在新的软件版本中会有更多的功能集成在此窗体 中完成。1.4 数据显示1.5 数据输入数据输入一般数据的打开一般数据的打开AVHRRHDF SeaWiFSMrSIDBMPJPEGNLAPSERMapper,PCI (.pix)JPEG 2000PDSERDAS 7.x(.lan)Landsat7Fast(.fst)RADARSATERDAS IMAGINE8.x(.img)Landsat 7 HDFSRFGeoTIFFMAS-50TIFFHDFMRLC(.dda)使用Open Image File打开 ENVI 图像文件或其它已知格式的 二进制图像文件。ENVI 自动地识别和读取下列类型的文件：1.5 数据输入数据输入特定数据的打开特定数据的打开 对于特定的已知文件类型，利用内部或外部对于特定的已知文件类型，利用内部或外部的头文件信息的头文件信息 通常会更加方便。使用通常会更加方便。使用 Open External File 选项，选项，ENVI 能够读取一些标准能够读取一些标准文件类型的若干格式，包括精选的遥感文件类型的若干格式，包括精选的遥感 格式、格式、军事格式、数字高程模型格式、图像处理软件格军事格式、数字高程模型格式、图像处理软件格式式 及通用图像格式。及通用图像格式。ENVI 从内部头文件读取从内部头文件读取必要的参数，因必要的参数，因 此不必在此不必在Header Information对话框中输入任何信息。对话框中输入任何信息。1.6常见系统设置常见系统设置用户自定义文件用户自定义文件 这里可以选择自定义的图形颜色文件、颜色表文件、ENVI的菜单文件（ENVI Menu File、Display Menu、Shortcut File）、地图投影文件 等，需要重启ENVI。1.6常见系统设置常见系统设置默认文件目录默认文件目录 设置一些ENVI默认打开的文件夹，如数据目录、临时文件目录、输出 文件目录、ENVI补丁文件、光谱库文件、备用头文件目录等，需要重 启ENVI。1.6常见系统设置常见系统设置显示设置显示设置 可以设置三窗口中各个分窗口的显示大小，窗口显示式样等。其中可 以设置数据显示拉伸方式（Display Default Stretch），默认为2%线 性拉伸。1.6常见系统设置常见系统设置综合设置综合设置 这个选项设置的是一些杂项，值得设置的为制图单位（Page Unit），默认为英寸（Inches），可设置为厘米（Centimeters）；还有缓冲大 小（cache size），可以设置为物理内存的50-75%左右，文件碎片大 小设置为cache size的1/10。2.影像预处理2.1 一般预处理流程介绍一般预处理流程介绍2.2 预处理中基础知识预处理中基础知识2.3 自定义坐标系自定义坐标系2.4 ENVI中的几何校正中的几何校正2.5 ENVI中的图像融合、镶嵌、裁剪中的图像融合、镶嵌、裁剪校正模型选择影像参考源控制点选取误差检查不 符 合图像融合图像镶嵌图像裁剪配准影像其它影像同名点选取（人工/自动）影像配准影像重采样校正影像2.1数据预处理一般流程数据预处理一般流程2.2预处理中基础知识预处理中基础知识常见商业高分辨率卫星常见商业高分辨率卫星传感器发射时间国家多光谱波段空间分辨率(米)重返周期IKONOS1999美国红、绿、蓝、近红外全色：1多光谱：41.5-2.9SPOT52001法国红、绿、近红外、中远红外全色：5 或2.5（超模式多光谱:1026Quick Bird（快鸟2001美国红、绿、蓝、近红外全色：0.61多光谱：2.441-3.5FORMOSAT II2004中国台湾红、绿、蓝、近红外全色：2多光谱：41EROS-B2006以色列/全色：0.7（立体）55CartoSAT-1（P5）2005印度/全色：2.5（立体）ALOS2005日本红、绿、蓝、近红外全色：2.5（立体）多光谱：102北京一号小卫星2005中国红、绿、近红外全色：4多光谱：323-5KOMPSAT-22006韩国红、绿、蓝、近红外全色：1 多光谱：43WorldView-1/22008美国红、绿、蓝、近红外 红边、海岸、黄、近红外2全色：0.5 多光谱：2.41.1-3.7资源应用卫星-2B星2008中国/全色：2.37多光谱：19.526GeoEye-12008美国红、绿、蓝、近红外全色：0.41（0.5）多光谱：1.652-3RapidEye2008德国蓝、绿、红、红边、近红外5每天其他卫其他卫星星传感器发射时间国家多光谱波段空间分辨率(米)Landsat177299美国蓝、绿、红、近红外、短波红外、热红外15、30、60、80、120SPOT41999法国绿、红、近红外、中远红外全色：10 多光谱:20中巴资源卫星-01/021999中国蓝、绿、红、近红外多光谱：19.5Resourcesat（P6）2003印度绿、红、近红外、短波红外多光谱 24 米全色 5.8米ALOS2005日本微波、立体像对、多光谱2.5米立体像对、10米多光谱、3米RadarTerraSAR-X2007德国微波1 m Radar、3m、5mCOSMO-SkyMed2007意大利微波3米、15米RADARSAT II2008加拿大微波3m 超细化模式1m景观光线模式NOAA气象卫星/美国红、近红外、中红外和两个热红外1.1km风云系列卫星/中国可见光4个，近红外2个，中远红外2个，热红外2个。1.1kmMODIS/美国36个波段250m、500m和1000m减小卫星A、B星2008中国多光谱近中红外(4波段)、高光谱(111波段)多光谱：30米 高光谱：100米Hyperion/EO-12000美国0.42.5 m共有220波段30米2.2 预处理中基础知识数据源的选择图像选择经济成本专题目的专题地域环境专题图比例尺空间分辨率时间分辨率波谱分辨率2.2 预处理中基础知识预处理中基础知识影像格式影像格式 传 感 器 文 件 格 式传 感 器 文 件 格 式不同的卫星传感器研发或运行机构一般会给所分发的卫星不同的卫星传感器研发或运行机构一般会给所分发的卫星数据设计数据设计 一种分发格式，如一种分发格式，如Landsat系列的系列的Fast格式、格式、EOS系列卫星的系列卫星的HDF格格 式等。式等。商 业 软 件 文 件 格 式商 业 软 件 文 件 格 式商业化的图像处理软件都会开发出软件本身的图像格式，商业化的图像处理软件都会开发出软件本身的图像格式，如如ENVI的的 Hdr&img格式，格式，Erdas的的IMG格式，格式，PCI的的pix格格式等。式等。通 用 图 像 文 件 格 式通 用 图 像 文 件 格 式很多图像格式成为国际通用，被大多数软件所支持。如很多图像格式成为国际通用，被大多数软件所支持。如TIFF、JPEG2000、BMP等。等。2.2预处理中基础知识预处理中基础知识引起图像畸变因素引起图像畸变因素 系 统 误 差系 统 误 差有规律的、可预测的。比如扫描畸变有规律的、可预测的。比如扫描畸变 非 系 统 误 差非 系 统 误 差无规律的无规律的如传感器平台的高度、经纬度、速度和姿态的不稳，地球如传感器平台的高度、经纬度、速度和姿态的不稳，地球曲率及空曲率及空 气折射，地形影响等气折射，地形影响等2.3预处理中基础知识预处理中基础知识几何校正中的几个概念几何校正中的几个概念 几何校正：纠正系统和非系统因素引起的几何畸变。图像配准（Registration）：同一区域里一幅图像（基准 图像）对另一幅图像校准，以使两幅图像中的同名像素配 准。图像纠正（Rectification）：借助一组控制点，对一幅图像进行地理坐标的校正。又叫地理参照（Geo-referencing）图像地理编码（Geo-coding）：特殊的图像纠正方式，把图像矫正到一种统一标准的坐标系。图像正射校正（Ortho-rectification）：借助于地形高程 模型（DEM），对图像中每个像元进行地形的校正，使图像 符合正射投影的要求。2.4预处理中基础知识预处理中基础知识实际中的概念实际中的概念 几 何 粗 校 正几 何 粗 校 正校正系统误差，地面站完成校正系统误差，地面站完成 几 何 精 校 正几 何 精 校 正包括图像纠正、地理编码和部分包括图像纠正、地理编码和部分图像配准图像配准 图 像 配 准图 像 配 准 正射校正正射校正2.2预处理中基础知识预处理中基础知识卫星影像的校正卫星影像的校正 根 据 卫 星 轨 道 参 数，包 括 位 置、姿 态、轨 道 及 扫 描 特 征，根 据 卫 星 轨 道 参 数，包 括 位 置、姿 态、轨 道 及 扫 描 特 征，校 正 影 像校 正 影 像(有有时加入时加入DEM)。地 面 控 制 点 校 正地 面 控 制 点 校 正+校正模型校正模型 轨 道 参 数轨 道 参 数+地面控制点地面控制点+DEM2.2预处理中基础知识预处理中基础知识多项式模型多项式模型 多 项 式 模 型多 项 式 模 型x=a0+a1x+a2Y+a3x2+a4xy+a5y2+y=b0+b1x+b2Y+b3x2+b4xy+b5y2+最 少 控 制 点 个 数最 少 控 制 点 个 数N=（n+1）*（n+2）/2 误 差 计 算误 差 计 算RMSEerror=sqrt（x-x）2+（y-y）2）2.2预处理中基础知识预处理中基础知识控制点获得途径控制点获得途径 基础数据基础测绘数据数字线画图（DLG）数字栅格图（DRG）影像数据正射影像（DOM）实地测量2.2预处理中基础知识预处理中基础知识控制点质量控制控制点质量控制 图 像 选 点 原 则图 像 选 点 原 则选取图像上易分辨且较精细的特征点：道路交叉点，河流弯选取图像上易分辨且较精细的特征点：道路交叉点，河流弯曲或分曲或分叉处，海岸线弯曲处，飞机场，城廓边缘等叉处，海岸线弯曲处，飞机场，城廓边缘等特征变化大的地区需要多选特征变化大的地区需要多选图像边缘部分一定要选取控制点图像边缘部分一定要选取控制点尽可能满幅均匀选取尽可能满幅均匀选取 数 量 原 则数 量 原 则在图像边缘处，在地面特征变化大的地区，需要增加控制点在图像边缘处，在地面特征变化大的地区，需要增加控制点保证一定数量的控制点，不是控制点越多越好。如一景保证一定数量的控制点，不是控制点越多越好。如一景TM的的控制点控制点数量在数量在30-50左右。左右。2.2预处理中基础知识预处理中基础知识重采样方法（插值算法）重采样方法（插值算法）最近邻法取与所计算点（x,y）周围相邻的 4个点，比较它们与被计算点的距 离，哪个点距离最近，就取哪个 亮度值作为（x,y）点的亮度值简单易用，计算量小，图像的亮 度具有不连续性，精度差2.2预处理中基础知识预处理中基础知识重采样方法（插值算法）重采样方法（插值算法）双线性内插法取（x,y）点周围的4邻点，在y方向内插二次，再在x方向 内插一次，得到（x,y)点的亮 度值f（x,y)双线性内插法比最近邻发虽 然计算量有所增加，但精度明 显提高，特别是对亮度不连续 现象或线状特征的块状化现象 有明显的改善。内插法会对图像起到平滑作 用，从而使对比度明显的分界 线变得模糊。x y双线内插算法原理示意图原始图像2.2预处理中基础知识重采样方法（插值算法）三次卷积内插法进一步提高内插精度的一 种方法，通过增加邻点来 获得最佳插值函数取与计算点周围相邻的16 个点，先在某一方向内插，再根据计算结果在另一个 方向上内插，得到一个连 续内插函数计算量大，精度高，细节 表现更为清楚，对控制点 要求较高 x y123452.2预处理中基础知识预处理中基础知识图像融合图像融合 图 像 融 合图 像 融 合将低分辨率的多光谱影像与高分辨率的单波段影像重采样生将低分辨率的多光谱影像与高分辨率的单波段影像重采样生成成一成成一副高分辨率多光谱影像遥感的图像处理技术。副高分辨率多光谱影像遥感的图像处理技术。关 键 技 术关 键 技 术两个影像配准在亚像元范围内两个影像配准在亚像元范围内融合方法选择融合方法选择运算速度和交换缓冲空间运算速度和交换缓冲空间2.2预处理中基础知识图像融合方法彩色合成数学运算彩色技术空间配准高分辨率多光/谱图像变换HIS变换加 与 乘差值 比值主 成 分 分 析滤波分析小波分析HSV变换2.2预处理中基础知识预处理中基础知识图像镶嵌图像镶嵌 图 像 镶 嵌图 像 镶 嵌指在一定地数学基础控制下，把多景相邻遥感影像拼接成一指在一定地数学基础控制下，把多景相邻遥感影像拼接成一个大范个大范围的影像图的过程。围的影像图的过程。关 键 技 术关 键 技 术颜色的平衡颜色的平衡接边处理接边处理运算速度和交换缓冲空间运算速度和交换缓冲空间2.2预处理中基础知识预处理中基础知识图像裁剪图像裁剪 图 像 裁 减图 像 裁 减图像裁剪的目的是将研究之外的区域去除，常用的是按照行图像裁剪的目的是将研究之外的区域去除，常用的是按照行政区划政区划边界或自然区划边界进行图像的裁剪。边界或自然区划边界进行图像的裁剪。关 键 技 术关 键 技 术裁剪区的确定裁剪区的确定无数据区处理无数据区处理2.3 自定义坐标系自定义坐标系坐标系原理坐标系原理 地理坐标系是以经纬度为单位的地球坐标系统，地理坐标 系中有2个重要部分，即地球椭球体（spheroid）和大地基 准面（datum）。大地基准面指目前参考椭球与WGS84参考椭球间的相对位置关系（3 个平移，3个旋转，1个缩放），可以用其中3个、4个或者7个参数 来描述它们之间的关系，每个椭球体都对应一个或多个大地基准面。投影坐标系是利用一定的数学法则把地球表面上的经纬线 网表示到平面上，属于平面坐标系。数学法则指的是投影 类型，目前我国普遍采用的是高斯克吕格投影，在英 美国家称为横轴墨卡托投影（Transverse Mercator）。2.3 自定义坐标系自定义坐标系北京北京54与西安与西安80坐标系坐标系 都是投影直角坐标系 北京54坐标系、西安80坐标系实际上指的是我国的两个大地基准面。坐标名称投影类型椭球体基准面北京54Gauss Kruger（Transverse Mercator）Krasovsky北京54西安80Gauss Kruger（Transverse Mercator）IAG75西安80椭球体名称年代长半轴（米）短半轴（米）扁率WGS8419846378137.06356752.31：298.257克拉索夫斯基（Krasovsky19406378245.06356863.01：298.3IAG-7519756378140.06356755.31：298.2572.3 自定义坐标系自定义坐标系 ENVI中自定义中自定义坐标系坐标系 ENVI中的坐标定义文件存放在HOME ITTIDLxxproductsenvixx map_proj 文件夹下，三个 文件记录了坐标信息：ellipse.txtdatum.txtmap_proj.txt椭球体参数文件 基准面参数文件 坐标系参数文件 在ENVI中自定义坐标系分三步：定义椭球体、基准面和定 义坐标参数2.3 自定义坐标系自定义坐标系定义椭球体定义椭球体 语法为,。这里将 “Krasovsky，6378245.0，6356863.0”和“IAG-75，6378140.0，6356755.3”加入ellipse.txt末端。注：ellipse.txt文件中已经有了克拉索夫斯基椭球，由于 翻译原因，这里的英文名称是Krassovsky，为了让其他软 件平台识别，这里新建一个Krasovsky椭球体。2.3 自定义坐标系自定义坐标系添加基准面与添加基准面与定义坐标系定义坐标系 语法为,。这里 将“Beijing-54,Krasovsky,-12,-113,-41”和 “Xian-80,IAG-75,0,0,0”加入datum.txt 末端。注：有的时候为了与其他软件平台兼容，基准面的名称直 接写成所用的椭球体名称。在ENVI任何用到投影坐标的功能模块中都可以新建坐标系 (在任何地图投影选择对话框中，点击“New”按钮。)，或 者直接选择主菜单-Map-Customize Map Projection2.3 自定义坐标系自定义坐标系坐标系定义练坐标系定义练习习 利 用 自 定 义 坐 标 系 将 一 幅 北 京利 用 自 定 义 坐 标 系 将 一 幅 北 京54坐标系转化为西安坐标系转化为西安80坐坐标标 系。系。试 验 的 栅 格 数 据 情 况 为：一 幅 北 京 坐 标 系 的 栅 格 数 据，投试 验 的 栅 格 数 据 情 况 为：一 幅 北 京 坐 标 系 的 栅 格 数 据，投影参数如下：影参数如下：投影类型：投影类型：Transverse Mercator 椭球：椭球：Krasovsky 基准面：基准面：Krasovsky（自定义）（自定义）中央经线：中央经线：117东向偏移：东向偏移：500000m2.4 ENVI中的几何校正中的几何校正传感器参数校传感器参数校正正 传 感 器传 感 器(带有地理定位文件带有地理定位文件)SPOT1-4SeaWiFSASTERAVHRRENVISATMODISRADARSAT 自 定 义 地 理 定 位 文 件自 定 义 地 理 定 位 文 件GLTIGM2.4 ENVI中的几何校正中的几何校正传感器参数校正传感器参数校正练习练习 数 据 源数 据 源Modis传感器的传感器的2级数据（级数据（“1-Modis”文件夹内文件夹内)EOS-HDF格式储存格式储存 处 理 过 程处 理 过 程利用自带地理定位文件进行几何校正利用自带地理定位文件进行几何校正 输 出输 出几何校正结果几何校正结果2.4 ENVI中的几何校正中的几何校正几何精校正流几何精校正流程程显示图像文件采集地面控制点计算误差选择几何模型检验校正结果开始结束重采样输出误差 太 大2.5 ENVI中的几何校正中的几何校正几何精校正练几何精校正练习习 数 据 源数 据 源已经做过几何校正的已经做过几何校正的SPOT4全色全色10米分辨率影像（米分辨率影像（“2-几何几何精校正精校正”文件夹）文件夹）待校正的待校正的Landsat5 TM 30米分辨率影像（米分辨率影像（“2-几何精校正几何精校正”文件文件夹）夹）处 理 过 程处 理 过 程用用SPOT4影像作为基准影像，选择控制点来校正影像作为基准影像，选择控制点来校正TM影像。影像。输 出输 出校正结果校正结果2.5 ENVI中的影像融合、镶嵌、裁剪中的影像融合、镶嵌、裁剪融合融合 ENVI中的融合方法HSV变换、主成分分析、Brovey变换，CN SpectralSharpening，高保真的Gram-Schmidt，Pansharping，小波融合（补丁）自定义 有地理坐标和无地理坐标都可以融合 操作方式主模块流程化操作ENVI ZOOM中 数据源数据源已经做过几何校正的已经做过几何校正的SPOT4全色全色10米分辨率影像（米分辨率影像（“3-影像影像融合融合”文件夹）文件夹）待校正的待校正的Landsat5 TM 30米分辨率影像（米分辨率影像（“3-影像融合影像融合”文件文件夹）夹）处 理 过 程处 理 过 程用主成分分析、用主成分分析、HSV颜色变换等方法融合两个影像颜色变换等方法融合两个影像 输 出输 出融合结果融合结果2.5 ENVI中的影像融合、镶嵌、裁剪中的影像融合、镶嵌、裁剪镶嵌镶嵌 基 于 像 素 镶 嵌 和 基 于 地 理 坐 标 镶 嵌基 于 像 素 镶 嵌 和 基 于 地 理 坐 标 镶 嵌 自 动 颜 色 平 衡，边 缘 直 方 图 匹 配自 动 颜 色 平 衡，边 缘 直 方 图 匹 配 接 边 线、接 边 羽 化接 边 线、接 边 羽 化 虚 拟 镶 嵌虚 拟 镶 嵌 运 算 速 度 快运 算 速 度 快 占用非常少的虚拟内存空间占用非常少的虚拟内存空间2.5 ENVI中的影像融合、镶嵌、裁剪中的影像融合、镶嵌、裁剪镶嵌练习镶嵌练习 数 据 源数 据 源两幅已经校正好的两幅已经校正好的SPOT4 10米全色影像（米全色影像（“4-影像镶嵌影像镶嵌”文文件夹中）件夹中）处 理 过 程处 理 过 程用注记工具勾画两影像接边线，用羽化和颜色校正等使两幅用注记工具勾画两影像接边线，用羽化和颜色校正等使两幅影像镶影像镶嵌在一起。嵌在一起。输 出输 出镶嵌结果镶嵌结果2.5 ENVI中的影像融合、镶嵌、裁剪中的影像融合、镶嵌、裁剪裁剪裁剪 空 间 裁 减空 间 裁 减基于感兴趣区（基于感兴趣区（ROI）的裁减）的裁减基于矢量基于矢量/栅格数据文件的裁剪栅格数据文件的裁剪自定义裁剪自定义裁剪 波谱裁剪波谱裁剪2.5 ENVI中的影像融合、镶嵌、裁剪中的影像融合、镶嵌、裁剪裁剪练习裁剪练习 数 据 源数 据 源一幅一幅TM影像、影像区域的影像、影像区域的Shapefile矢量文件（矢量文件（“5-影像裁剪影像裁剪”文文件夹中）件夹中）处 理 过 程处 理 过 程用用ROI工具在工具在TM影像上绘制不规则的多边形感兴趣区域，影像上绘制不规则的多边形感兴趣区域，后利用这后利用这 个感兴趣区域裁剪个感兴趣区域裁剪TM影像影像利用利用Shapefile矢量文件裁剪矢量文件裁剪TM影像影像 输 出输 出裁剪结果裁剪结果3.影像信息基本提取方法3.1 影像信息提取技术概述3.2 影像增强处理3.3 监督分类3.4 非监督分类3.1 影像信息提取技术概述影像信息提取技术概述 遥感影像通过亮度值或像元值的高低差异（反映地物的光 谱信息）及空间变化（反映地物的空间信息）来表示不同 地物的差异，这是区分不同影像地物的物理基础。遥感影像分类就是利用计算机通过对遥感影像中各类地物 的光谱信息和空间信息进行分析，选择特征，将图像中每 个像元按照某种规则或算法划分为不同的类别，然后获得 遥感影像中与实际地物的对应信息，从而实现遥感影像的 分类。3.1 影像信息提取技术概述影像信息提取技术概述遥感分类技术的发展遥感分类技术的发展 可分为四个阶段人工解译基于光谱 计算机自 动分类基于专家 知识的决 策树分类面向对象 特征自动 提取四种方法并存3.2图像增强处理 图像变换主成分（PCA）独立主成分最小噪声分离(MNF)颜色空间变换（HSV,HLS）穗帽变换波段运算 图像拉伸去相关、饱和度、彩色直方图（匹配、拉伸）滤波卷积形态学、纹理、自适应、自定义频率域局部增强智能数字化工具智能数字化工具 智 能 数 字 化 工 具智 能 数 字 化 工 具提高数字化的效率提高数字化的效率增强增强ENVI矢量的功能矢量的功能计算提取的线状地物的长度计算提取的线状地物的长度过头去除以及两线相交过头去除以及两线相交ENVI中基于光谱分类方法中基于光谱分类方法非监督分类ISODATAK-Means监督分类基于传统统计分析分类器 平行六面体 最小距离 马氏距离 最大似然基于人工智能分类器 神经网络基于模式识别分类器 支持向量机 模糊分类类别定义/特征判别影像分类分类器选择样本选择分类后处理结果验证平行六面体 最小距离 马氏距离 最大似然 波谱角 二进制编码 光谱信息散度 神经网络 支持向量机分类 模糊分类3.3 监督分类监督分类3.2 监督分类监督分类练习练习数据源以Landsat TM为数据源（“6-监督与非监督分类”文件夹内）。处理过程选择样本,后选择一种分类器进行分类。分类后处理 类后处理Majority/Minority 分析ClumpSieve 精度分析 生成随机样本 混淆矩阵结果分类结果3.4 非监督分类非监督分类 非监督分类：也称为聚类分析或点群分类。在多光谱图像中搜寻、定义其自然相似光谱集群的过程。分类器选择ISODATAK-means其他影像分类类别定义/类别合并码分类后处理 结果验证影像分析3.4 非监督分类非监督分类练习练习数据源以Landsat TM为数据源（“6-监督与非监督分类”文件夹内）处理过程分类器选择ISODATA或者K-mean对TM进行分类。分类后处理 类别定义 类后处理Majority/Minority 分析、Clump、Sieve 重新组合类别 精度分析 生成随机样本 混淆矩阵结果分类结果4.制图与三维可视化4.1 ENVI的快速制图4.2 三维可视化4.1 ENVI的快速制图的快速制图QuickMap随意选择边界颜色随意选择边界颜色设置所有格网线属性设置所有格网线属性颜色、样式、字体等颜色、样式、字体等整饰和必要的标注作为图像显示整饰和必要的标注作为图像显示自动标注自动标注地图综合地图综合三维场景构建与三维场景构建与3D 曲面飞行曲面飞行 快 速 利 用 影 像 及快 速 利 用 影 像 及DEM构建三维场景构建三维场景 叠 其 他 数 据，如 矢 量叠 其 他 数 据，如 矢 量 用 户 自 定 义 背 景用 户 自 定 义 背 景 可 对可 对DEM和影像进行重采样，提高浏览速度和影像进行重采样，提高浏览速度 放 大 时 逐 步 增 强 细 节放 大 时 逐 步 增 强 细 节 设置飞行路线，可将飞行录制成视频设置飞行路线，可将飞行录制成视频ENVI高级预处理主要内容主要内容 1、GLT几何校正几何校正 2、单 景 高 分 辨 率 影 像 的 正 射 校 正、单 景 高 分 辨 率 影 像 的 正 射 校 正 3、影像自动配准、影像自动配准1、GLT几何校正几何校正 GLT几何校正法利用输入的几何文件生成一个地理位置查找 表文件（geographic lookup table，GLT），从该文件中 可以了解到某个初始像元在最终输出结果中实际的地理位 置。地理位置查找表文件是一个二维图像文件，文件中所包含 两个波段：地理校正影像的行和列，文件对应的灰度值表 示原始影像每个像素对应的地理位置坐标信息，用有符号 整型储存，它的符号说明输出像元是对应于真实的输入像 元，还是由邻近像元生成的填实像元（infill pixel）。符号为正时说明使用了真实的像元位置值；符号为负时说 明使用了邻近像元的位置值，值为0说明周围7个象元内没 有邻近像元位置值。GLT校正风云三号气象卫星校正风云三号气象卫星安装安装补丁补丁 先安装ENVI的HDF5读取补丁，将文件解压放 在ITTIDL70productsenvi45save_add下，打 开ITTIDL70productsenvi45menuenvi.men文件，在0 File1 Open Image File open envi fileenvi_menu_event下添加以下菜单1 Open HDF5 File not used open_hdf5_event启动ENVI,在菜单File下新增一个菜单Open HDF5 File。GLT校正风云三号气象卫星校正风云三号