前情提要：主页的另一篇代码计算SVF的结果有较大纰漏，所以学了新的方法。

SAGA下载链接SAGA GIS - Advanced Geodata Analysis Software

建筑数据依旧来自：3D-GloBFP: the first global three-dimensional building footprint dataset

ESSD - 3D-GloBFP：全球首个三维建筑覆盖区数据集

一、arcgis前置处理

先把建筑数据shapefile转化成栅格，栅格值就是height字段，我这边设置的分辨率为3m，大家分辨率可以任意选。

因为要计算SVF，所以平地部分也需要处理。打开栅格计算器， Con( IsNull("bfp_raster") ,0,"bfp_raster")，这样得到的最终栅格影像平地部分的值就是0。

二、SAGA

下载SAGA，【file-grid】将建筑栅格导入。

打开工具箱，选择SVF工具

这个工具要求输入的参数不多，只必须输入一个栅格影像，具体参数可以看说明文档。

我基本默认设置。然后点击okay就可以开始运行了。差不多一个半小时跑出结果。

结果可以右键save保存到本地，也可以在软件内可视化。其实蛮简单的。

三、arcgis后置处理

把前面的bfp栅格用栅格计算器Con(IsNull("Building"),0,1)，变成一张二值化影像Building_mask。然后再SetNull("Building_mask" == 1, "SVF")，把svf中属于房顶的部分挖空。

最后得到的结果就是研究区内的SVF了。