《遥感数字图像处理》01次实习报告

（1）以多波段组合方式将GeoTIFF格式的白银市TM原始数据转换为ENVI Standard 格式：

利用Basic Tools/Layer Stacking弹出对话框然后Import File，弹出对话框，导入GeoTIFF格式的TM原始数据，选择波段1、2、3、4、5和7，

点击OK，利用Choose选择输出路径及文件名，同时可以利用Reorder Files对输入的文件根据自己的需要进行调换顺序，点击OK输出ENVI Standard格式的数据。

（2）查询并记录影像文件的基本信息、投影信息，以及各个波段直方图信息，然后编辑头文件：

利用Basic Tools/Resize Data弹出对话框里面选择要查看的影像，左

边会出现其基本信息，如图所示：也有投影信息，既可以用来看单波段的也可以看合成后整个影像的信息。在对话框下，合成影像的名字上右击，选择Quick Statistics弹出对话框，在此对话框中点击Select

Plo下拉菜单，选择单波段或者多波段的直方图，相应的对话框中会出现直方图（在结果与分析中记录），还可以右击选择edit修改横、纵坐标的单位。

同样的在合成影像的名字上右击，选择Edit Head，弹出对话框

然后点击Edit Attributes/Band Name弹出对话框，选中波段输入修改

后名字，点击OK即可进行波段名字的修改。点击Edit Attributes/Wavelengths弹出进行相应的波长的修改。

（3）在View视窗中，利用影像缩小、放大、漫游工具识别影像中的土地利用/土地覆盖类型:

可以结合当地的google earth上高分辨率的遥感影像，进行识别，利用Viewer视窗下Tools/SPEAR/Google Earth/Jump to Location可以在google earth上显示View主视窗中相应选中地物对应的位置。

（4）利用Viewer视窗打开影像，分别选取4、3、2和7、4、2波段组合进行假彩色合成，观察实习内容中所要求地物的色调变化：

利用File/Open Image File，选择第1步合成的ENVI Standard 格式的数据，弹出对话框，在其中选择RGB Color，将R、G、B分别设为4、3、2波段，点击Load

Band，在Viewer#1中出现了4、3、2波段组合的假彩色图像，再在此窗口中，点击Display/New Display,弹出Viewer#2，选择RGB Color,将R、G、B分别设为7、4、2波段，点击Load Band，在Viewer#2中出现了7、4、2波段组合的假彩色图像，在Viewer窗口中右击选择Link

Displays，弹出对话框，点击OK，可以把两个窗口中同一位置进行连接起来，

即其中一个窗口放大、缩小、漫游到某个位置，另外一个也跟着漫游到其相对应的位置。这样可以进行地物色调变化的对比。

（5）提取6种地物在不同波段的数值（Digital Number，DN），做光谱剖面图：

在Viewer视窗中Tools/Profile/Z Profile(Spectrum)弹出对话框,在其

Options下拉菜单中勾选Plot Key,对话框中出现了Viewer视窗中选中的目标地物的X，Y坐标，然后勾选Collect Spectra，鼠标箭头变为十字箭头，在目标地物中取九个点（本来图上就有一个，总共是十个点），然后在选择File/Save Plot As/ASCII弹出对话框

，点击Select All Items，利用Choose选择输出路径和文件名，点击

OK，将其保存为.txt格式。选六种地物，重复以上操作，提取不同波段的数值（Digital Number，DN）。将.txt格式的文件用excel打开，然后用插入函数中的average函数求出每种地物的平均DN值，然后做出光谱剖面（光谱图如结果与分析中所示）。

（6）使用Excel制作6种地物的样本特征光谱统计表：

在Excel中分别使用插入函数中的AVERAGE、VAR、STDEV、MAX和MIN函数求出各地物样本DN值在各个波段的平均值、方差、标准差、最大值和最小值。然后，在07版Excel 的“Microsoft Office 按钮”，单击“Excel 选项”。“加载项”，然后在“管理”框中，选择“Excel 加载项”，单击“转到”弹出“加载宏”，在弹出来的对话框中选择“分析工具库”，并点击确定。然后从“工具”中找到“数据分析”，从“数据分析”对话框中选择“协方差”，并导入某种地物需求协方差的数据区域并选择“逐行”进行，最后选择数据输出区域并确定，则可得该地物的协方差矩阵。同理，在从“数据分析”对话框中选择“相关系数”，进行相应操作，可求得相关系数矩阵。（在结果与分析中附有个地物的样本特征光谱统计表）（7）制作散点图：

在Excel中，打开6种地物的样本DN数据（5步骤产生的），选择band2和band4做散

点图。从插入图表中选择散点图，然后点击“选择数据”，弹出对话框，进行加载系列和

数据：如图所示：按“添加”后弹出

进行系类名称和X、Y轴数据的选择，X为Band2，Y为Band4，如此重复完成六个系类的添加，点击确定，即可生成散点图。

（8）从可见光/近红外的6个波段影像中选取512×512子区并保存：

利用Basic Tools/Resize Data,弹出对话框，选择第1步生成的ENVI Standard格式的数据，点击Spatial Subset,弹出对话框，然后选择Image，

弹出如下对话框：，选择自己感兴趣的地方，并在samples 和lines后输入512，连续点击OK，最后在对话框中，利用choose选择输出路径和名字，点击OK即可生成裁剪的图。

结果与分析

（1）影像的基本信息：

查看ENVI Standard格式数据（即经过转化后的数据）的信息：宽、高、层数和储存形式: 8031 x 7151 x 6 [BSQ]，大小：344,578,086 bytes，文件类型: ENVI Standard，投影: UTM, Zone 48 North，分辨率: 30 Meters，基准面（datum）: WGS-84。

直方图信息：