Arcgis、QGIS网格划分、编码与坐标获取

目的

当前很多应用场景都需要对城市进行网格化，比如销售分析、人群分析、店铺选址等等。在开展业务的时候，在遇到如下需求，

1. 根据城市的行政区划边界，生成100*100米的网格
2. 对生成的网格进行编码，编码规则，从上到下，从做到右
3. 获取网格中心点坐标以及四个定点的坐标。

工具

1. arcgis
2. qgis

备注：单独使用arcgis就可以，只是本人中途换了mac电脑，用不了arcgis，所以中途换了qgis，小伙伴如果不想用qgis的话可以对照步骤去找arcgis中的工具进行操作

数据

1. 厦门市行政边界.shp

步骤

1、生成渔网

使用acrgis中的渔网工具

qgis也可以生成渔网，但是qgis中定义网格大小只能用经纬度，无法用米做单位

2、给网格进行编码

要求（x,y）编码，其中从左到右，从上到下的顺序

2.1 获取每个网格左上角的坐标

基本公式

!shape.extent.Xmin!：计算要素最小外包矩形的左上角X坐标;

!shape.extent.XMax!：计算要素最小外包矩形的右下角X坐标;

!shape.extent.Ymax!：计算要素最小外包矩形的左上角Y坐标;

!shape.extent.YMin!：计算要素最小外包矩形的右下角Y坐标;

!shape.area!：计算要素面积;

!shape.length!：计算要素周长。

用Pathon语句分别计算出这2个字段，计算语法为：!shape.extent.Xmin!和!shape.extent.Ymax!

 2.2、 使用排序工具

第一次排序

针对排序结果进行编码

#顺序编码
rec=0
def numberxp():
global rec
pStart=1
pInterval=1
if(rec==0):
rec=pStart
else:
rec=rec+pInterval
return rec

第二次排序

进行编码

rec=0
def numberxp():
global rec
pStart=1
pInterval=1
if(rec==0):
rec=pStart
else:
rec=rec+pInterval
return rec

修改编码

对于x来说，第一列为1，第二列为2，第三列为3。

对于y来说，第一行为1，第二行为2，第三行为3

def get_hang(id):
if (id%539==0):
rec= id // 539-1
else:
rec= id // 539
return rec

3、网格坐标点获取

在QGIS中完成（在arcgis中同样可以完成，本人更换了苹果使用不了arcgis，改成QGIS）

基于“网格划分与编码”的结果对网格进行编码后，生成网格是中心点和四个边界点

3.1、获取网格中心点

3.2、获取中心点的坐标

 

3.3、修改中心点的字段名称

3.4、将中心点坐标给原来的面数据

3.5、删除多余字段

3.6、提取网格四个顶点的坐标

3.7、删除多余字段

3.8、提取四个顶点的坐标