《遥感图像处理讲课教案》由会员分享，可在线阅读，更多相关《遥感图像处理讲课教案（122页珍藏版）》请在金锄头文库上搜索。
1、遥感图像处理,一、遥感图像基本操作 二、遥感图像的校正 三、遥感图像信息增强 四、遥感图像融合处理 五、遥感图像计算机分类及应用,一、遥感图像基本操作,（一）图像的输入与输出 （二）图像的裁剪 （三）图像的重采样,（一）图像的输入与输出,（一）图像的输入与输出,（一）图像的输入与输出,（一）图像的输入与输出,在“Available Bands List”菜单中选择“Available Files List”还可以可以显示出遥感图像的基本信息，具体如图所示,（一）图像的输入与输出,点击“MAP Info ”也可以看到部分图像信息,（一）图像的输入与输出,若要保存图像，需要在图像所在窗口中选择Fi
2、leSave Image AsImage File，弹出Output Display to Image File对话框：,对于单波段图像，选择8-bit Color，而多波段彩色合成图像则选择24-bit Color。 图像的保存方式有两种：一是直接保存为文件；二是选择“Memory”，记忆在“Available Bands List”菜单中。,（二）图像的裁剪,对遥感图像进行裁剪是图像处理工作的基础，本例针对图像can_tmr.img对进行裁剪。 1首先打开图像can_tmr.img，选择主菜单中的“Basic Tools” “Resize Data”：,（二）图像的裁剪,弹出Resize
3、Data Input File对话框 ，选择你要裁剪的图像之后，点击Spatial Subset，弹出Select Spatial Subset对话框，就可以进行裁剪了。,（二）图像的裁剪,（二）图像的裁剪,（二）图像的裁剪-image方式,选择Image方式对图像进行裁剪，它是一种根据图像裁剪的方法，拖拽红色框的左上角来选择裁剪范围，图像中的红框大小可以通过输入数值改变其大小,（二）图像的裁剪-file方式,还可以选择File方式进行裁剪，它是根据另一个已知文件来选择裁剪范围的方法，点击File之后，弹出Subset By File Input Filed对话框，选择已经存在的图像之后进行裁
4、剪。,（二）图像的裁剪-map方式,还可以根据地图坐标裁剪。在Select Spatial Subset对话框中，选择map方式，弹出Spatial Subset by Map Coordinates对话框，分别输入左上角和右下角的地理经纬度坐标，即可完成对图像的裁剪操作，保存即可,（三）图像的重采样,常用的方法有最邻近法（Nearest Neighbor）、双线性内插法（Bilinear）、以及立方卷积法（Cubic Convolution）,（三）图像的重采样,图像裁剪设置完毕后，点击OK，就可以进行重采样设置。,（三）图像的重采样,在Resize Data Parameters对话框的相
5、应文本中，输入所需要重采样的样本数和行数；或在“xfac”和“yfac”中输入你要放大或者缩小的倍数。,（三）图像的重采样,二、遥感图像的校正,（一）辐射校正 （二）几何校正,（一）辐射校正,辐射校正包括成像时期照明条件的差异、大气影响、系统噪声、仪器误差等的改正和补偿。 1、大气校正 2、系统噪声消除,1、大气校正,对遥感图像进行大气校正指的是去掉大气散射作用造成的路径辐射附加在地物中的部分。,式中，LAi为卫星观测值；Ri为地面反射率； a和b为回归系数；,1、大气校正,（1）公式法 与卫星扫描同步进行野外波谱测试，将地面测量结果与卫星影像对应像元亮度值进行回归分析，回归方程为：,1、大气
6、校正,（2）回归分析法 作法：在不受大气影响的波段（如TM5）和待校正的某一波段图像中，选择由最亮至最暗的一系列目标，将每一目标的两个待比较的波段灰度值提取出来进行回归分析。 例如：,a1,b1计算如下：,T1、T5表示TM1与TM5波段灰度值，,为TM1波段校正后的灰度值。,辐射回归分析图：,（3）直方图校正法 通过灰度直方图对比找出校正量。,1、大气校正,打开一幅遥感影像的第二、第四波段，由于第四波段属于近红外波段，波长较长，受大气影响较小，可以利用波长较长的第四波段来对波长较短的第二波段进行校正。,1、大气校正,首先把两幅图像动态连接，在较暗部比较像元亮度值的大小，估计亮度值差值的一个准
7、确值。,1、大气校正,按右图设置动态连接的选项。,1、大气校正,如果二波段相同位置像元的亮度值比四波段大，这个差距是由大气的辐射影响造成的，可以利用波段运算的方法对第二波段的亮度值减去他们的差值，然后保存就可以消除第二波段中大气辐射的影响。,1、大气校正,还可以用直方图对比的方法进行大气校正。其原理如下图所示。,1、大气校正,如上面的方法一样，将zoom图像放置于比较暗的地方，对第二波段和第四波段的图像进行连接，打开第二、四波段的直方图，即可找出两个波段最暗部分的亮度值差，然后可以按照上面的方法来消除大气辐射的影响。,1、大气校正,用同样的方法，显示出第四波段的直方图，可以比较两个波段最小亮度
8、值的差别，即可完成大气校正。,1、大气校正,此外，ENVI 中还有特定的功能对各种数据进行大气校正工作， 选择菜单Basic ToolsPreprocessingCalibration Utilities，这个工具可以对各种不同的数据，进行自动的大气校正工作。,2、系统噪声消除,图像中的主要的噪声是周期性噪声、条带与坏线现象及尖峰噪声等。 （1）周期性噪声一般重叠在原图像上，成为周期性的干涉图形，与尖峰噪声一道可以通过滤波的方式进行处理。 （2）遥感图像中通常会出现与扫描方向平行的条带，还有一些与辐射信号无关的条带噪声，一般称为坏线，它的出现和分布具有随机性质。,（2）去除条带及坏线,去除条带
9、以及坏线处理，既可以在原空间域（即原图像）中进行，也可以通过傅立叶变换进行滤波处理来实现。这里讲述的是如何在原图像中进行去除坏线。,（2）去除条带及坏线,打开图像，可以看到很多条带噪声。,（2）去除条带及坏线,选择Basic ToolsPreprocessingGeneral PurposeReplace Bad Lines。出现Bad Lines Input File 对话框，选择一个输入文件，也可以通过Spectral Subset 来选择单一波段进行去除。,（2）去除条带及坏线,（2）去除条带及坏线,在“Half Wildth to Average”文本框中，键入要参与计算平均值的邻近行
10、数。在要替代的行周围数值应是对称的。例如，值为“2”意味着每边各有两行将参与平均值计算。,（2）去除条带及坏线,对于图像坏线的准确行数，可以在打开的图像的上点击右键，选择Pixel Locator，然后在zoom图像窗口中，准确确定坏线的行数。,（二）几何校正,遥感图像几何校正包括光学校正和数字纠正两种方法。 数字纠正是通过计算机对图像每个像元逐个地解析纠正处理完成的，其包括两方面：一是像元坐标变换；二是像元灰度值重新计算（重采样）。,（二）几何校正,1、图像对图像的校正 2、图像对地图的校正,1、图像对图像的校正,1、在实习数据路径下打开SPOT数据文件bldr_sp.img，并加载到显示窗
11、口，它已具有地理信息，展开波段列表中的Map Info选项，如下图所示。,1、图像对图像的校正,1、图像对图像的校正,由于刚才打开的SPOT图像具有正确的坐标信息，可以用它来校正另一幅同一地区的不具有地理信息的TM图像。,1、图像对图像的校正,一旦两幅图像都已经显示，在ENVI主菜单中选择MapRegistrationSelect GCPs:Image to Image,1、图像对图像的校正,，出现Image to Image Registration对话框时，在“Base Image:”下面点击选择基图像（参照图像）的显示（Display #1，SPOT图像）。在 “Warp Image:”
12、 下方点击选择被纠正的图像的显示(Display #2，TM图像)。,1、图像对图像的校正,单击OK，出现 Ground Control Points Selection 对话框。在该对话框里添加地面控制点（GCP）。,1、图像对图像的校正,1、图像对图像的校正,选择好控制点后，如左图所示，在Ground Control Points Selection 对话框中选择OptionsWarp file，,1、图像对图像的校正,选择bldr_tm.img文件，点击OK,1、图像对图像的校正,出现Registration Parameters对话框，纠正方法选择RST，重采样方法采用最邻近采样：,1
13、、图像对图像的校正,右图为校正后的TM影像。,2、图像对地图的校正,图像对地图配准要求至少由一幅图像是打开的。图像中地面控制点(GCPs)由缩放窗口中的光标选择。也可以选择亚像元坐标。相应的地图坐标被手工输入或从矢量窗口输入。一旦选择了足够的点定义一个纠正多项式，在纠正图像中的GCP位置就能预测。,2、图像对地图的校正,1、在实习数据路径下打开TM图像bldr_tm.img，图像默认加载到显示窗口中（图像打开时的默认显示可以在头文件里进行设置）,2、图像对地图的校正,选择Map Registration Select GCPs: Image to Map，出现Image to Map Regi
14、stration对话框,2、图像对地图的校正,在该对话框中，在 “Select Registration Projection” 列表中，点击类型选择要输出的投影。这里选择UTM投影，在标有 “Zone” 的文本框里输入区域号13，像元分辨率接受默认的30m，点击OK,2、图像对地图的校正,出现 Ground Control Points Selection 对话框，在该对话框中选择地面控制点，地面控制点的选择类似于图像-图像配准,2、图像对地图的校正,在纠正图像的缩放窗口中定位要添加的控制点，它的相应的像素坐标出现在“Image X” 和 “Image Y”框里（可能为亚像元）。然后在分别在
15、标有“E”和“N”的文本框里，键入该点对应的地图坐标（可从地图上读取）。,2、图像对地图的校正,一旦选择完图像中需要的像元，且地图坐标已经输入，点击 Ground Control Points Selection 对话框中的 “Add Point” 按钮，将点添加到 GCPs 的列表里。,2、图像对地图的校正,用同样的方法添加另外的GCPs。 按前述方法执行纠正。,三、遥感图像信息增强,遥感图像信息增强包括空间信息增强和波谱信息增强，这里主要介绍最为基础的空间信息增强。,三、遥感图像信息增强,（一）反差扩展 （二）彩色合成,（一）反差扩展,反差扩展是一种通过拉伸或扩展图像的亮度数据分布，使其占
16、满整个动态范围（0-255），以达到扩大地物之间亮度差异，分出更多的亮度范围等级的一种处理技术。反差扩展只能在单波段的图像上进行。,（一）反差扩展,从主菜单的File中，打开TM 1-7波段数据，在窗口中打开任意波段的数据。系统在打开图像时默认进行2%的线性扩展。 所谓2%的线性扩展指的是分别去除2%的小值和大值，这样拉伸显示异常值基本都能去掉。,（一）反差扩展,从上图可以看出，由于进行了2%的线性拉伸，使得图像层次感下降，可以选择主图像菜单中的enhance scroll linear菜单对图像进行反差扩展。,（一）反差扩展,要想改变打开文件时，系统默认的2%线性扩展，可以从主菜单File中，打开PreferencesDisplay Default，将%Linear中的2.0改为0.0，选择OK后。关闭对话框，保存生成的新文件再打开，就形成新的从0开始的线性变换。,（一）反差扩展,主图像菜单中选择EnhanceInteractive Stretching 可以进行交互式的反差扩展。Strech_Type中可以选择各种扩展
