ARCGIS实战（2），由遥感数据和道路、河流数据，生成生态廊道指标（指标为连续的栅格影像）。

本文总览：
生态廊道可以简单理解为：河流、道路两侧的植被（本文仅考虑道路，河流处理方式与道路类似）
现有数据：landsat8遥感影像、道路线数据
目标数据：代表生态廊道好坏的连续的栅格数据
难点一：道路为极长且相互交叉的地物，一个道路要素的起端和末端的生态廊道情况差别很大，多个道路相互交叉，几何关系混乱。
难点二：道路两侧的生态廊道未能覆盖整个区域，如何是目标数据为连续的栅格数据。
技术路线图：

具体操作：
1、植被提取
（1）landsat8遥感影像（可用高分影像替代）由地理空间数据云下载。
（2）如果一副遥感影像不需要进行纠正处理。
（3）arcgis打开landsat8红波段（band4）和近红外波段（band5）

（4）打开栅格计算器输入公式：(NIR-R)/(NIR+R)。特别说明，如下图，需要在公式前方加1.0*，保证生成数据为浮点型。

（5）本研究研究区为淄博张店城区，在此进行一步掩膜操作。


（6）在遥感影像上寻找纯植被的NDVI值NDVIveg和纯裸地的NDVI值NDVIsoil，将大于NDVIveg的栅格赋值为NDVIveg，小于NDVIsoil赋值为NDVIsoil，可通过栅格计算器con条件函数实现。下图为将小于NDVIsoil赋值为NDVIsoil的实现步骤。（大于NDVIveg的栅格赋值为NDVIveg，操作类似）

（7）使用栅格计算器，对上一步结果进行归一化处理，具体公式：f=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)。

（8）完成植被提取。

（9）由于（8）结果的空间分辨率为30m，为与道路两侧10m面要素匹配，将（8）植被栅格数据空间分辨率变为3m。
2、道路数据处理
（1）22年高德路网下载，可通过python爬取，也可以通过某宝获取。

（2）根据fclass属性或者oneway属性确定缓冲半径，通过属性选择工具和缓冲工具对不同级别的道路线进行缓冲，生成道路面数据。（道路面数据也可以通过矢量化或者其他途径获取）


（3）将道路面再此进行缓冲处理，缓冲半径可以根据自己需要所定。本文探究对象为道路两侧10m内的生态廊道情况，故缓冲半径为10m。


（4）新建渔网线，将渔网线合并处理，利用渔网线对道路面缓冲结果进行裁剪。并使用高级工具下的拆分多部件要素工具对结果进行打散处理。

（5）打开属性表，添加浮点型字段，使用计算几何下的计算面积工具，计算廊道的面积。

3、道路与植被交互处理
（1）使用arctoolbox，空间分析工具，区域分析，以表格显示分区统计。将植被覆盖栅格与廊道矢量面关联，关联规则为每一个廊道面的植被覆盖均值。最后生成表格。

（2）通过属性连接功能，对表格数据与廊道数据进行连接，最后如图所示，每个廊道均有平均植被覆盖率属性。

（3）打开廊道属性表，使用字段计算器将廊道的面积和廊道植被均值相乘，生成总植被覆盖度字段。

（3）新建与2（4）相同大小的渔网面，并以渔网面为目标要素，廊道面为连接要素，廊道面的面积字段、总植被覆盖度字段的合并规则为总和。空间连接后，在属性表中添加廊道绿地率字段，使用字段计算器为其赋值，公式为：总植被覆盖度字段/廊道面的面积字段。得到基于网格的廊道绿地率数据。

（4）通过面转点、插值分析得到反映廊道植被情况的连续的栅格影像。