TM遥感影像解译技术手册
1、TM遥感影像的数据类型为DAT，我们拿到的一景影像中的原始数据有三个文件，一个是SCENE01文件夹，其中包括小程序制作教程视频TM影像的7个波段，band1到band7，还包括一个header.dat的头文件，共8个数据，另外的两个是LBL和SELF文件，可以用记事本打开，记录这景TM影像的一些基本信息，如时间、经纬度、像元的行列数等。
2、TM影像的打开
我们要打开这些TM影像有两个方法，一是用ERDAS打开，二是用ENVI打开，首先我们用ERDAS打开，因为TM原始数据为二进制DAT数据，首先将band1到band7TM的7个波段通过ERDAS的import模块将其转成ERDAS可以用的IMG格式文件，首先点击菜单条上的import然后在type(类型)中选择Generic Binary（二进制），在media里选择file，imput里选择我们要选择的TM影像波段，如band1，存出一个路径，点击OK按钮，我们转出的格式为ERDAS可以打开的IMG格式，只是在现在是单波影像，分为7个单独的文件，不利于我们以后的解译工作，需要将这7个波段组合成一个文件，点击菜单条上的interpreter模块，然后点击Utilities,在弹出的菜单条中选择layer stack,在弹出的窗口中选择需要合成的单波段文件，然后点击ad
d按钮，重复的将需要合成的7个单波段加入到对话框中，在output file中写入文件名和调整输出路径，输出数据类型为unsignde8 bit,波段组合选择union，然后选中c语言条件运算符实例Ignore Zero In Stats复选框（统计忽略0值），点击OK按钮，执行波段组合。这样我们就可以从ERDAS打开我们组合完成的IMG格式TM影像，包括7个波段。
现在利用ENVI打开TM影像的原始文件，打开TM影像的路径为，file→open external file→land sat→fast,51cto打开头文件header.dat文件，在弹出的对话框中选择RGB color,然后分别选择4、5、3三个波段，点击load band就可以打开TM影像，这种方式打开的TM影像是直接有经纬度的，只是由于投影和椭球体的原因使得经纬度与我们需要利用的有些偏差，需要下一步的将其校正，更改其投影和椭球体。
3、TM影像的校正
  TM影像的校正同样也有两种方法，一是利用ERDAS软件功能对其进行校正，再一个就是利用ENVI软件对其进行的校正，首先介绍利用ERDAS软件进行的校正，第一步：同时打开两个viewer,viewer1中打开需要校正的TM影像，viewer2中打开具有地理参考的校正过的影像，并将两个viewer平铺放置，（单击Session→Tile Viewers命令）。第二步：启动几
何校正模块，在Viewer #1中单击Raster→Geometric Correction,打开对话框选择Polynomial，单击OK。弹出Polynomial Model Properties对话框，在Polynomial Order选择2，单击Close,弹出GCP对话框→Existing viewer→OK，单击参考图象→OK。第三步：采集地面控制点，可以看见GCP工具对话框（GCP Tool），利用配准图象校正TM原始影像，选取具有代表性地物作为控制点（GCP），如山顶点，公路交叉口，城市道路等，选取6个控制点后计算机将根据前6个点来自动选取第7个点，这个时候一般是不准确的，误差很大，需要手动继续调整，一景TM影像一般选取20个GCP点就可以做到校正，但是这些点必须选取均匀，尤其在我们利用的城市周围及农作物耕作土地上，尽量选取多的控制点，这样就不会引起图象与实际地区的几何形变，利于我们下一步的工作。第四步：选取控制点后计算转换模型，单击Geo Correction Tools对话框方块图标，开始计算，弹出Resample对话框，写入输出图象文件名，定义输出像元大小（Output Cell Sizes）X/Y值都为30，选中Ignore Zero in Stats复选框，单击OK。这样TM影像的校正完成。
下面介绍利用ENVI校正TM影像过程，File→open External File→landsat→hdf打开TM原始影像中的头文件header.dat,打开后存成TIFF格式文件,我们将以TIFF格式文件为校正对象,打开TIFF格式文件后,单击MAP→Convert Map Projectionsstruct函数使用→选择需要校正的TIFF格式文件,
单击OK,弹出Convert Map Projection Paramerers对话框,在Solect output map projection列表中选择我们需要的投影类型albers,确定输出文件名并选择输出到File.选择输出文件路径,点击OK.校正时间较长,一般TM影像需要1小时左右,校正后的文件需要再次存成TIFF格式,然后通过ERDAS软件Input模块转换成ERDAS可以利用的IMG格式.进行下一步工作.一般情况下,利用ENVI软件校正后的影像与真实情况有一定的偏移,需要我们利用ERDAS软件对其进行平移,而且由于ENVI和ERDAS软件的不通用性,使得TM影像在TIFF转成IMG格式时,影像的彩及亮度都有变化,使得TM影像看起来不太清晰,但是不影响栅格文件的灰度值,不影响工作.
4、TM遥感影像的分类
TM遥感影像的分类我们只利用ERDAS软件,因为ERDAS软件的分类比较准确并且操作简单,第一步,启动非监督分类,Classification→Unsupervised Classification打开分监督分类对话框,第二步,进行非监督分类设置及进行非监督分类,确定输入文件名(input raster file),确定输出文件名(output file),取消输出文类模版文件(output signature set),确定聚类参数(initialize from statistics),确定分类数,一般我们需要的TM影像分类都为20类(number of clas
ses),定义最大循环次数(maximum iterations)为24,设置循环收敛阈值(convergence threshold)为0.95,单击OK.这样分类开始,一般只需要十几分钟就可以完成分类.第三步,进行分类后处理,分类后处理是为了消除计算机的盲目性,包括两步,聚类统计和去处分析,我们先进行聚类统计,工具条interpreter→GIS analysis→Clump命令,确定输入文件,定义输出文件名,确定聚类统计邻域大小(Connect Neighbors)为8,单击OK.聚类统计需要较长时间,特别当邻域为paramecium8时,但是聚类为8时效果较好,一般由于TM影像也是比较小,所以推荐用8.聚类统计完成后应用去除分析,工具条interpreter→GIS analysis→Eliminate命令,打开对话框,确定输入文件,确定输出文件名,坐标类型为Map,确定最小图斑大小(Minimum)为16pixels,输出数据类型为unsigned 4 Bit,单击OK.这样就完成了分类图象的后处理.
5、TM影像的分割,一般情况下都要求对TM影像进行分割,一般都是按照旗县来分割TM影像,我们在对TM影像分割推荐使用ERDAS的Mask来对其进行掩模,当然ENVI软件也可以对其进行分割,但是分割后的图像为HDR格式文件,还需要进行转换成上交时要求的IMG格式,在转换的过程中会出现一些图象显示上的错误.所以我们推荐使用ERDAS软件进行下一步工作.首先我们把各个旗县的SHP格式文件转换成COVERAGE格式,因为ERDAS软件只支持矢量的coverageproperties文件用什么打开格式文件,由于不知道的原因,在转换成coverage格式后,矢量文件缺失椭
球文件,所有转换完成后需要在属性表里手动加入托球体文件.开始进行分割,工具条interpreter→utilities→mask,弹出对话框,确定输入文件,确定掩模文件,确定输出文件,我们不对其进行设置,默认情况下就可以对其进行掩模了.
6、面积的提取
按照旗县切割完的图象利用光标查询就可以了解分类完成后的像元数目,因为TM影像的空间分辨率为30,通过像元数就可以乘以30m2来计算我们需要的农作物面积,如果分类为20类,那么第8或者第9类就为我们需要的作物,如果分类为10类,那么一般情况下是第2类为农作物.通过光标查询就可以看到我们需要的像元数目,我们一般是通过地面调查的样地来进行分类,通过ARCVIEW软件来进行观察,看分类后哪种类型为我们需要的农作物面积.之后把我们需要的通过EXCEL表格来进行统计后就可以知道我们需要的面积.这样,我们的工作就完成了.

版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。