[续]基于高斯两步移动搜寻法（2SFCA）的城市绿地可达性分析[格网]

原文链接

[续]基于高斯两步移动搜寻法（2SFCA）的城市绿地可达性分析[格网]https://mp.weixin.qq.com/s/BkjZVqCAORrKH5bsyOSTbg

之前我用立方学社公开课的数据做了一次2SFCA，传送门：基于高斯两步移动搜寻法（2SFCA）的城市绿地可达性分析http://mp.weixin.qq.com/s?__biz=MzkzNjA5NjA3OA==&mid=2247487203&idx=1&sn=4c7598e608e172616df15251100d92f0&chksm=c2a2b6c2f5d53fd49ba9d09b2f3aa8e88952019065985d6b7633e6e52bf02a6f873fda71858b&scene=21#wechat_redirect。不过当时做的时候小区数据是没有人口的，也有人问过这个问题，里面用的小区人口是根据每个小区的户数估算出来的，所以准确度就低了一些，修人口学了不能白学是吧[跟练]基于七普修正Worldpop人口栅格数据(附2020年worldpop100m人口栅格)http://mp.weixin.qq.com/s?__biz=MzkzNjA5NjA3OA==&mid=2247491255&idx=1&sn=493ff37dc588bfe8e8aae4516d571453&chksm=c2a2a696f5d52f8042cb7a29229429db27442a35db33b67a44959f077fe0dc7d2cb641a67fd9&scene=21#wechat_redirect...这次嘛就用worldpop的人口栅格就好啦。

我在想着用渔网格代替小区，就是下面这篇论文啊

使用[创建渔网]工具创建渔网格，格网单元大小为200m×200m；输出范围和输出坐标系在环境里设置一下（在所有数据都经过投影变换的前提下）（渔网格命名为finsh200，以下都用此名代替）

然后用创建的渔网格通过[以表格显示分区统计]统计每个网格内的人口，得到sum_pop

接着根据FID字段将sum_pop中汇总的人口数量连接至finsh200中。连接进去后外围的那些渔网格肯定是没有值的，需要把这些无值的处理掉，右键finsh200-图层属性-定义查询，选择人口＞0的。

绿地计算一下面积，公式什么的也不再说了，可以看之前那篇基于高斯两步移动搜寻法（2SFCA）的城市绿地可达性分析或者自己翻翻论文。这样需求端和供给端都处理好了接下来进行操作。

第一步，绿地至小区（就是计算这个公式）

[邻域分析]-[生成近邻表]，距离阈值还是用上一次取的2.5km。

得出的表中，IN_FID就是绿地，NEAR_FID就是格网，NEAR_DIST是距离，这个生成近邻表生成的是直线距离，所以有些不足，最好还是结合路网进行分析。然后通过IN_FID连接绿地面积，通过NEAR_FID连接人口。将人口和绿地面积都放在一个表里。

计算那个G，高斯函数，GIS里我把公式贴在这，根据自己的数修改距离字段和阈值即可。

(math.exp((-1/2)*(!NEAR_DIST!/2500)*(!NEAR_DIST!/2500))-math.exp(-1/2))/(1-math.exp(-1/2))

 G有了，接着计算G*Pk（人口）

字段名写错了，是GPk，叫PkDkBk都无所谓，代表的是人口，大家注意下就行

接着是按照公式，∑汇总GPk，以绿地（也就是IN_FID为汇总字段汇总GPk）

将sum_GPk这个表连接回绿地表，再用绿地面积/sum_GPk就得到了Rj

 至此，第一步完成。

第二步，小区至绿地（计算这个公式）

还是上一步操作[生成近邻表]，距离阈值仍然是2500m，只是将输入要素和邻近要素换一下，IN_FID就是格网，NEAR_FID就是绿地。

第一步计算出来了Rj，根据NEAR_FID（绿地）与sum_GPk这个表中的IN_FID进行连接，将Rj连接过来到这个表中。

根据公式，还得需要计算这个G，跟第一步一样。完了接着计算G*Rj就行了，不再演示

接着是按照公式，∑汇总G*Rj，以网格（也就是IN_FID为汇总字段汇总GRj）,得到Ai。

至此，第二步完成。

最后

最后便是将Ai根据IN_FID连接回finsh200，在符号系统中进行可视化

操作思路大概就是这么个操作，结果大概就这样子

之前的练习数据可以公众号【一点规划】回复【两步】，人口栅格数据可以回复【2020】