ArcGIS教程：1997-2013年北京用地许可热度图（一）——热度图生成

图1   2013年北京用地许可热度图

    某天，笔者看到BCL发布的北京用地许可动态热度图，遂想能否在ArcGIS中实现，正好将热度图的生成方法分享给大家。

    下载了数据会发现，仅仅是张xlsx表，但是里面有x、y坐标信息，这样就能把每个用地许可证的建设位置定位到地图上。

1.生成点

    右键Create Feature Class—From XY table

图2

    经纬度，笔者选的阿里云坐标。坐标系统GCS_WGS_1984。

图3

    这样，我们得到了好多点，如下。

图4

2.每年提取一份数据

    上面得到的点涵盖了1997-2013年所有的用地许可，我们希望每年单独生成一份数据，将来可以做逐年的动态展示。

    属性表里有两列“证书编号”、“核发日期”都包含年份信息，我们根据其中一列就能提取年份。

图5

    添加一个text型字段“year”，取“证书编号”前4位为其赋值。借助mid函数，year= mid( [证书编号] ,1,4)，如下。

图6

图7

    有年份列，我们就能逐年提取用地许可了。用select工具即可，表达式"year" = '2013'。

3.生成热度图

    对提取出来的2013年用地许可点，使用核密度工具Kernel Density生成热度图。

    核密度工具根据点或线要素计算每单位面积的量值，并将各个点或线拟合为光滑锥状表面。

    输出象元大小和搜索半径可以根据出图的平滑效果调整。

图8

    出来的效果可能如下，默认的符号化方式不好看，但是能看到点密集程度的雏形。

图9

    那我们来换种拉伸色阶：

图10

    最后可以从ArcGIS Onlie上添加北京底图，就能展示出如图1的2013年用地许可热度图了。