如何使用Python和ArcGIS绘制重庆板块地图？

本文共计1040个文字，预计阅读时间需要5分钟。

使用Python和ArcGIS绘制重庆市区块地图：数据来源：[此处插入数据来源]完整代码：pythonimport arcpy

设置工作环境arcpy.env.workspace=C:/GIS/重庆市区块地图.gdb

输入数据input_feature_class=C:/GIS/重庆市区块.shp

输出数据output_feature_class=C:/GIS/重庆市区块地图.shp

创建新的要素类arcpy.CreateFeatureclass_management(arcpy.env.workspace, output_feature_class, POLYGON)

添加字段arcpy.AddField_management(output_feature_class, 区块名称, TEXT)

添加图形要素with arcpy.da.InsertCursor(output_feature_class, [SHAPE@, 区块名称]) as cursor: with arcpy.da.SearchCursor(input_feature_class, [SHAPE@, 区块名称]) as search_cursor: for row in search_cursor: cursor.insertRow([row[0], row[1]])

代码解读1. 导入arcpy模块2. 设置工作环境3. 输入数据4. 输出数据5. 创建新的要素类6. 添加字段7. 添加图形要素8. 代码解读利用Python的arcpy模块绘制重庆市区块地图，完整代码如上所示。

python for ArcGIS 绘制重庆市板块地图

• ​​原数据​​
• ​​完整代码​​
• ​​代码解读​​

利用python的arcpy模块绘制出重庆市板块地图如下

原数据

完整代码

# -*- coding: utf-8 -*-
"""
Project_name:drawing plate for Hangzhou
@author: 帅帅de三叔
Created on Mon Dec 09 15:01:15 2019
"""
import sys
arcpy_path = [r'D:\Python27\ArcGIS10.6\Lib\site-packages',
r'D:\Program Files (x86)\ArcGIS\Desktop10.6\arcpy',
r'D:\Program Files (x86)\ArcGIS\Desktop10.6\bin',
r'D:\Program Files (x86)\ArcGIS\Desktop10.6\ArcToolbox\Scripts']
sys.path.extend(arcpy_path)
stdi,stdo,stde=sys.stdin,sys.stdout,sys.stderr
reload(sys) #通过import引用进来时,setdefaultencoding函数在被系统调用后被删除了，所以必须reload一次
sys.stdin,sys.stdout,sys.stderr=stdi,stdo,stde
sys.setdefaultencoding('utf-8')

import arcpy #导入地理处理模块
from arcpy import env #导入环境类
env.workspace=r'D:\python for ArcGIS\绘制重庆板块' #设置环境变量
env.overWriteOutput=True #覆盖
import pandas as pd #导入数据分析模块
plate_data=pd.read_excel(u"重庆板块.xlsx") #读取数据
rows,cols=plate_data.shape #确定数据尺寸
plate_name=plate_data['板块名称'].unique() #板块名称
lng=plate_data[u'经度'] #板块边界经度
lat = plate_data[u'纬度'] #板块边界纬度
spatial_ref=arcpy.SpatialReference('China Geodetic Coordinate System 2000') #设置空间参考系

polygonPoints=arcpy.Array() #用来存放多边形的折点
polygonGeometryList=[] #用来存放多边形几何对象
for name in plate_name: #按板块名称循环
lng_data=plate_data[plate_data[u'板块名称'] == name][u'经度']
lat_data=plate_data[plate_data[u'板块名称'] == name][u'纬度']
for lng, lat in zip(lng_data, lat_data):
point=arcpy.Point()
point.id=name.index; point.X=float(lng); point.Y=float(lat) #构造折点对象
polygonPoints.add(point) #点成线
polygon=arcpy.Polygon(polygonPoints,spatial_ref,"","") #线成面
polygonGeometryList.append(polygon) #多边形组
polygonPoints.removeAll() #移除折点
arcpy.CopyFeatures_management(polygonGeometryList, r"D:\python for ArcGIS\绘制重庆板块\plate_Chongqing.shp","POLYGON")

代码解读

整个实现过程包括3步

读数，即读取板块边界经纬度数据

拆点，即采用split()函数将经纬度数据分割，构成折点

连线，即将折点连起来形成封闭多边形，即板块

若还有不明白的，可以来“三行科创”微信公众号交流群。

1，​​python for ArcGIS 绘制上海市板块地图​​

2，​​python for ArcGIS 绘制上海市环线地图​​

3，​​python for ArcGIS 绘制北京市板块地图​​

4，​​python for ArcGIS 绘制广州市板块地图​​

5，​​python for ArcGIS 绘制深圳市板块地图​​

6，​​python for ArcGIS 绘制成都市板块地图​​

7，​​python for ArcGIS 绘制武汉市板块地图​​

8，​​python for ArcGIS 绘制杭州市板块地图​​

9，​​python for ArcGIS 绘制重庆市板块地图​​

10，​​python for ArcGIS绘制合肥市板块地图​​

11，​​python for ArcGIS 绘制南京市板块地图​​

12，​​python for ArcGIS绘制宁波市板块地图​​

13，​​python for ArcGIS 绘制沈阳市板块地图​​

14，​​python for ArcGIS 绘制苏州市板块地图​​

15，​​python for ArcGIS 绘制天津市板块地图​​

16，​​python for ArcGIS 绘制西安市板块地图​​

17，​​python for ArcGIS 绘制郑州市板块地图​​