赫子淏-人造覆盖提取操作流程资料

1、土地分类居民地(无参考数据)提取软件操作流程： (1) 生成坡度（SLOPE）数据 a.在已拼接好的DEM数据中裁切出待分类影像所对应的DEM。为了减少分类过程中的计算量，一般要对影像上不参与运算的区域进行掩膜。具体步骤是：在ENVI里打开一景待分类的TM影像，点击主菜单Basic Tools-Masking-Build Mask，选择要处理的显示窗口后，打开Mask Definition窗口。点击Options-Import data range.，打开选择TM影像点击OK后，出现Input for Data Range Mask窗口，参数设置如下图：点击OK后回到Mask Definiti

2、on窗口，点击Choose按钮选择输出mask文件的路径和名称，完成mask文件的制作。Memory1为制作好的掩膜文件，选中Memory1右键，Open File in ArcMap，在ArcMap中选择Conversion Tools-From Raster-Raster to Polygon，转成shp。（注意 在栅格转shp时，Simplify polygons（optional）前面的勾要取消）开始编辑生成的shp文件，删除外部边框，引入拼接好的DEM，Date Management Tools-Raster-Rester Processing-Clip，运算后得到对应TM影像区域的

3、DEM。b.把DEM数据转化为SLOPE坡度数据。 由于拼接好的DEM数据是经纬度投影，像元单位是度。而计算坡度的时候使用的单位是米，因此在计算坡度数据之前，要把DEM的经纬度投影转换为对应TM影像的UTM n48(n49)投影。打开ENVI软件，并打开拼接好的DEM文件和待处理的TM影像。点击主菜单Map-Convert Map Projection，选择要转换投影的拼接好的DEM数据，可以打开Convert Map Projection Parameters窗口。点击Options-Match existing file.，选择待处理的TM影像，点击OK后回到Convert Map Pro

4、jection Parameters窗口。在Resampling中设置参数为Bilinear，点击Enter Output Filename边上的Choose按钮，选择输出转换成UTM投影后的DEM文件存放的路径和文件名，点击OK等待完成。（输出的DEM文件区域位置大小是与待处理的TM影像对应的）下面生成SLOPE数据。打开转换为UTM投影的DEM数据，点击ENVI主菜单Topographic-Topographic Modeling菜单，选择DEM文件后，打开Topo Model Parameters窗口，参数设置如下图：点击Choose按钮，选择输出slope坡度数据的路径和文件名，点击O

5、K，运算后得到对应TM影像区域的SLOPE数据。(2) 对原始影像进行去水处理。 利用ENVI中的辅助工具-矢量转栅格，把原始TM影像的水体（shp文件格式）转为栅格文件，点击ENVI主菜单中的Basic Tools-Band Math-在Enter an expression：中输入公式 b1 eq 1 and b2 eq 0-点击Add to list-选中Previous Band Math Expressions：下的b1 eq 1 and b2 eq 0-OK b1选择步骤(1)a中得到的原始TM对应的掩膜文件，b2选择原始TM对应的水体栅格文件，点击Choose按钮，选择输出路径和

6、文件名，点击OK，得到在原始TM掩膜文件中去水的掩膜文件。点击主菜单Basic Tools-Masking-Apply Mask，打开Apply Mask Input File窗口选择原始TM影像，然后点击Select Mask Band按钮选择掩膜波段，选择刚生成mask波段，最后点击OK，并保存。(3) 样本选择。 人造覆盖指以人工结构覆盖为主的地表，由沥青、混凝土、沙石、砖瓦、玻璃以及其他人工建材覆盖的地表。人造覆盖以居民地为主，同时也包含道路、机场、厂矿用地、露天采矿场、工矿堆积地、建筑用地和城市空地等。影像特征： 6,4,3波段呈现粉红色，平原地区居民地呈团状，山区居民地的分布受地形

7、限制，多呈带状、丁字状、Y字状、弧状等。点击 ENVI主菜单->Open Image File，打开去水后的原始TM点击以上窗口中的“RGB Color”,然后再点击窗口左下角“Load Band”,弹出下面窗口，点击ENVI主菜单- Basic Tools-Region Of Interest-ROI Tool，弹出以下窗口，这个窗口用来选取样本。以上窗口中，在ROI Name这一列中，列出需要选取样本的几个类型，如在本次居民地分类中，需要选取两类样本，居民地和非居民地，因此，在选取样本时需要列两类，artificial和other,窗口如下，在Color这一列中对每种类型的样本颜色进

8、行选择，在“window”后面的选项中，选取“zoom”窗口在这个窗口中选择样本，开始选择样本时，需要鼠标选中要选择样本的这一个类型，然后在“zoom”窗口中选择样本，选择结果如下所示，红色的是水体的居民地，绿色是非居民地在样本选取完成后，需要保存样本，如下图，在“#1ROI Tool->Save ROIs”,在弹出的右边的窗口中点击“Select All Items”，这样选择完成需要保存的样本类型，在下面点击“choose”设置保存路径保存样本。(4) 进行SVM监督分类 在ENVI主菜单中选择Classification-Supervised-Support Vector Mach

9、ine选择去水的原始影像-OK点击Select All Items，Output Rule Images选择No，点击Choose选择输出路径和命名。点击OK。(5) TM影像的分割 利用ecogniton软件，设定不同的尺度，进行分割。得到不同分割尺度的shape文件和栅格影像（TIF格式）。(6) 使用分割文件与坡度数据辅助分类。点击分类工具-农田分类-专家知识，如下图：点击后可以打开专家知识窗口，如下所示：点击SVM分类结果Choose按钮，选择步骤(4)中得到的SVM分类结果；点击分割影像文件choose按钮，选择步骤(5)得到的分割影像（tif格式）；点击坡度影像文件choose按钮

10、，选择步骤(1)b的坡度文件；Choose按钮，选择输出路径和命名，程序运行完毕后，可以得到最终居民地提取结果。(7) 分类结果的后处理。由于步骤(6)得到的分类结果是ENVI标准文件形式的栅格格式，为了便于下一步的人工编辑后处理，需要转换为矢量格式。而后在Arcgis里面进行人工修改编辑，完成之后保存。右键，Open File in ArcMap，在ArcMap中选择Conversion Tools-From Raster-Raster to Polygon，转成shp格式。（注意 在栅格转shp时，Simplify polygons（optional）前面的勾要取消）而后在Arcgis里面

11、进行人工修改编辑，完成之后保存。(8) 把人工编辑好的为矢量（shape）分类结果再转换回栅格形式 利用ENVI辅助工具-矢量转栅格得到最终结果。土地分类居民地(有参考数据)提取软件操作流程： (1)在Arcmap中导入和原始数据对应的水体矢量，开始编辑，删除shp文件中的全部内容。目的是得到与原始TM影像投影等信息都相同的空shp文件。 (2)对原始数据进行制作掩膜，导入Arcmap中转为矢量(shp)格式，开始编辑，去除外边框，得到原始TM的边界文件。 (3)引入参考数据。 (4)根据位置进行选择：Selection-Select By Location 在the following layer(s):窗口下选择参考数据，