当前很多应用场景都需要对城市进行网格化，比如销售分析、人群分析、店铺选址等等。在开展业务的时候，在遇到如下需求，
- 根据城市的行政区划边界，生成100*100米的网格
- 对生成的网格进行编码，编码规则，从上到下，从做到右
- 获取网格中心点坐标以及四个定点的坐标。
- arcgis
- qgis
备注：单独使用arcgis就可以，只是本人中途换了mac电脑，用不了arcgis，所以中途换了qgis，小伙伴如果不想用qgis的话可以对照步骤去找arcgis中的工具进行操作
- 厦门市行政边界.shp
使用acrgis中的渔网工具
qgis也可以生成渔网，但是qgis中定义网格大小只能用经纬度，无法用米做单位
要求（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
在QGIS中完成（在arcgis中同样可以完成，本人更换了苹果使用不了arcgis，改成QGIS）
基于“网格划分与编码”的结果对网格进行编码后，生成网格是中心点和四个边界点
3.1、获取网格中心点
3.2、获取中心点的坐标
3.3、修改中心点的字段名称
3.4、将中心点坐标给原来的面数据
3.5、删除多余字段
3.6、提取网格四个顶点的坐标
3.7、删除多余字段