孔庙地址来自于cnrds，公司注册地址来自国泰安。半径内孔庙数量由这两个地址算出。

数据年份跨度:2000_2020（国泰安现在还没有更新21年的公司经纬度，所以21年的为空）

文件带有stata代码，可以自己算。

r_5代表5km内孔庙数量，r_10代表10㎞内，以此类推。