怎样通过居民地面图层判断地块类型

l 问题背景

项目中需要根据一个点的位置坐标判断该点所在的地块类型是城市、农村或城乡结合部。

现只能提供一个居民地面矢量数据图层，且图层仅包含图形信息（即：位置、长度、面积），没有可用的属性信息（如：地上层数、地下层数、结构材质、用途）。

居民地面图层数据（黄色）在ArcMap中叠加谷歌影像地图后显示如下：

 

属性信息如下：

 

l 解决方案

1. 版本一：

u 操作步骤

查询缓冲半径（如：200米）范围（圆形）内有无房屋，无房屋则为农村，有房屋则根据房屋面积再次判断，即：若有面积大于指定数值（如：500平方米）的房屋，则为城市，若面积均小于该数值，则为城乡结合部。

u 识别效果

不理想。如缓冲半径过小，则即便是城市地区，位于房屋之间的部分及公园、道路等无房屋的地方，均不能正确识别。即便缓冲半径合适，由于城乡结合地区也有不少大于500平方米的房屋，故不能正确区分城市和城乡结合部。

2. 版本二：

u 操作步骤

查询固定的缓冲半径（如：200米）范围（圆形）内有无房屋，无房屋则为农村，有房屋则根据房屋个数、面积再次判断，即：分别统计面积大于、小于指定数值（如：500平方米）的房屋个数，再求二者的比值。如该比值大于某个数（如：0.5），则说明面积较大的房屋相对较多，为城市，否则为城乡结合部。

缓冲半径必须大于指定值，不允许设置过小的缓冲半径。

u 识别效果

稍好些。但对有只有非常少量的面积较小的区域，本应为农村，却被识别成了城乡结合部，即对农村地区的识别正确率不够高。

3. 版本三：

u 操作步骤

查询固定的缓冲半径（如：200米）范围（矩形）内有无房屋，无房屋则为农村，有房屋则根据房屋个数、面积再次判断，即：统计该范围内的所有房屋的面积之和，并取其占矩形范围面积的比例。

对结果比例设置两个间断值（如：10%、50%）。若比例小于较小间断值，则为农村；高于较大间断值，则为城市；介于二者之间，则为城乡结合部。

u 识别效果

较好，但依然存在三个影响识别结果精度的地方。一是两个间断值如何设置比较合理，二是缓冲半径如何设置比较合理，三是接近于间断值的情况如何处理。

4. 版本四：

u 操作步骤

根据给定点的位置和缓冲半径，生成九宫格，对每一个宫格（矩形范围），分别判断其地块类型。判断过程如下：

1. 查询固定的缓冲半径（如：200米）范围（矩形）内有无房屋，无房屋则为农村，有房屋则根据房屋个数、面积再次判断，即：统计该范围内的所有房屋的面积之和，并取其占矩形范围面积的比例。

2. 对结果比例设置两个间断值（如：10%、50%）。若比例小于较小间断值，则为农村；高于较大间断值，则为城市；介于二者之间，则为城乡结合部。

3. 求取这9个宫格的类型后，分别统计三种类型的个数，由类型较多者决定该点所在的地块类型。若出现两种或三种类型相同的情况，则根据该九宫格的矩形范围求取房屋面积占比，按其落入的间断值的区间判断地块类型。

此外，先根据对全图范围按不同半径值运行的结果（房屋占比统计），人为设定间断值的默认值及最佳缓冲半径值。如：设置初始半径值为50米，将全图划分为若干个100m*100m的矩形鱼网，统计各鱼网范围内的房屋面积占比；半径值每次递增50米，重新划分鱼网，统计房屋面积占比；半径值递增至1000米停止运算。统计结果结合人眼对地物的判别，确定最佳半径值和间断值。

u 识别效果

比较理想。半径值的选取非常关键，初测结果显示，750米较为理想。

l 正确的识别结果示例

如下（中心点位置为传入的坐标参数）：

   

左：城市      

  中：农村             

右：城乡结合部