DEM与土地利用叠加分析

       最近要做DEM与土地利用数据叠加分析，类似下图，最终得出区域尺度景观特征，看起来应该会很简单很容易完成的样子，但操作过程出现了很多问题，现将操作过程总结，希望能帮到更多的同学。

一、准备数据

       DEM、土地利用数据、研究范围

二、操作

1、获取数据

       先从公网上获取DEM和土地利用数据，分别查看数据源的像元大小，两个数据尽量保持像元一致，本次数据不一致，所以需要将数据进行处理，保持像元一致，可在数据裁剪（提取掩膜）的时候设置好像元大小。

2、按掩膜提取，裁剪保留研究区域范围

        对DEM与土地利用数据进行裁剪，只保留与研究范围一致大小即可，像元可以在裁剪时设置。

       （1）点击Space Analyst Tools（空间分析）——Extraction（提取）——Extract by Mask（掩膜提取）工具，在弹出的掩膜提取面板位置4输入需要裁剪的DEM，位置5输入用于裁剪的矢量范围数据，在位置6选择输出路径，点位置7的环境，设置Output Coordinate System（输出坐标系）、Processing Extent（处理范围）和Raster Analysls（像元设置），详细参数设置如下图：

        1）坐标系设置：输出坐标系可以根据需要，选跟自己需要输出坐标系相同的图层，这里输出跟研究范围一致的坐标系即可，整个分析过程尽量保持坐标系一致，如下图所示：

环境——输出坐标系设置

        2）输出范围设置：可以在位置1设置输出范围，本次输出范围与研究范围一致大小即可，所以这里选择“与研究范围一样”，位置2选择输出的栅格像元跟裁剪前保持一致，所以选择dem栅格，如下图所示：

环境——处理范围设置

        3）像元设置：栅格分析参数是设置像元大小，如果要保持裁剪前后保持一样，就直接选择裁剪栅格，像元大小就自动在下面的框内显示，如果不想要这个像元，可直接在像元框内输入像元大小，如下图所示：

环境——输出像元设置

        （2）提取掩膜完成后，结果如下图所示，输出与研究范围一致的栅格，右击栅格，在属性——源，可以查看像元大小，确保数据是否准确。

3.对栅格数据进行重分类

        （1）对DEM按高程范围重分类，便于将高程范围值划分界定，选择Space Analyst Tools（空间分析）——Reclass（重分类）——Reclassify（重分类）工具，在位置4输入重分类栅格数据，在位置5下方也就是框起来的范围输入重分类范围，位置6输出路径，点确定即可，如下图：

        （2）重分类完成后得到一个重分类好的DEM数据如下图：

4.栅格转面

        （1）将DEM栅格数据转面，修改属性表，将分类信息整理至属性表。单击工具箱Conversion Tools（转换工具）——From Raster（来自栅格）——Raster to Polygon（栅格转面），在弹出的栅格转面框中，位置4输入需要转面的栅格，位置5输入分类字段，位置6填写输出路径，最后点OK即可执行；

        （2）修改属性值：打开转出来的面图层属性表，点击文件——添加字段，用于存放高程分类信息，分别将1,2,3,4,5,6,7对应数据筛选出来，然后在新建字段中输入对应分类；

添加字段

        （3）点击属性表的文件——按属性筛选，在弹出框中，按步骤2,3,4,5完成一类数据的筛选，右击新字段，选择字段计算器，在位置7输入分类信息，点ok即可完成一类数据信息的标识，其他类按以上步骤给值，直至完成7类数据属性信息修改；

注意：本次操作中DEM、土地利用数据均为栅格数据，两个数据均需执行文档中操作2、3、4步骤，本次只详细写了DEM的详细操作过程，土地利用省略了，当执行到步骤4的时候应该再重新执行一遍2、3、4步骤对土地利用数据进行操作。

5.求DEM和土地利用交集

        （1）单击工具箱Analysis Tools（分析工具）——Overlay（覆盖）——Intersect（相交）如下图在位置1选择要求相交的两个矢量数据文件，然后点确定即可。

        （2）编辑属性表：右击打开新生成图层属性表，在属性表中新建type字段用于存放两者交集分类信息，右击type字段，在弹出的计算器中将DEM分类信息字段2和土地利用分类字段3并一起，组成两个图层的分类信息；

        （3）符号设置：右击图层，在符合系统页面位置3选择用于展示的分类字段，然后点击位置4添加所有分类，最后点击确定即可完成分类渲染，操作如下符号设置图，可以点击每个分类前的符号重新设置颜色；

        （4）符号完成后，结果如下图所示。