# -*-coding:utf-8-*-
import arcpy
import pyodbc
import os
# 指定工作空间
arcpy.env.workspace = r"E:\shp"
# 2000坐标系
fc_2000 = "YDBP_2000.shp"
cursor_2000 = arcpy.da.InsertCursor(fc_2000, ["FLOWSN", "XMMC", "KCID", "DKID", "COOR_TYPE", "[email protected]"])
# 80坐标系
fc_80 = "YDBP_80.shp"
cursor_80 = arcpy.da.InsertCursor(fc_80, ["FLOWSN", "XMMC", "KCID", "DKID", "COOR_TYPE", "[email protected]"])
# 未知坐标系
fc_unknown = "YDBP_unknown.shp"
cursor_unknown = arcpy.da.InsertCursor(fc_unknown, ["FLOWSN", "XMMC", "KCID", "DKID", "COOR_TYPE", "[email protected]"])
# 根据需要设置变量：流程号、项目名称、勘测定界ID、地块ID、坐标系类型
FLOWSN = ""
XMMC = ""
KCID = ""
DKID = ""
COOR_TYPE = ""
count = 0
# 遍历每个mdb文件
path = r"E:\LYMDB"
mdblist = os.listdir(path)
try:
for mdbfile in mdblist:
# 获取项目名称
XMMC = mdbfile[:-4]
print XMMC.decode("gb2312") + " is beginning..."
# 获取mdb文件的完整路径
mdbfullpath = os.path.join(path, mdbfile)
# 根据本机的驱动程序编写连接字符串
# mdbsource = r‘DRIVER={Microsoft Access Driver (*.mdb,*.accdb)};DBQ=‘ + mdbfullpath + ‘;‘
mdbsource = r‘DRIVER={Microsoft Access Driver (*.mdb)};DBQ=‘ + mdbfullpath.decode("gb2312") + ‘;‘
with pyodbc.connect(mdbsource) as dbconn:
cur = dbconn.cursor()
# 根据勘测定界技术说明，获得勘测ID与坐标系的对应关系
xysql = "select * from DELIMITATION_TAC order by SB_SB_ID;"
cur.execute(xysql)
xyalldata = cur.fetchall()
xynum = len(xyalldata)
xydict = {}
if xynum > 0:
for xyid in range(xynum):
xydict[xyalldata[xyid][26]] = xyalldata[0][23]
# 根据勘测定界表，获取FLOWSN，勘测ID与地块ID的对应关系
kcidsql = "select * from PLOT order by PL_ID;"
cur.execute(kcidsql)
kcidalldata = cur.fetchall()
kcidnum = len(kcidalldata)
kcdict = {}
if kcidnum > 0:
FLOWSN = kcidalldata[0][14]
for id in range(kcidnum):
kcdict[kcidalldata[id][13]] = kcidalldata[id][15]
# 根据坐标点顺序升序排序,并获取界址点成果表
sql = "select * from PNT_COORD order by PL_ID,PNT_SERIAL;"
cur.execute(sql)
alldata = cur.fetchall()
# 计算成所有坐标行数，包括表头
rownum = len(alldata)
# 存放坐标点的列表
polygonarray = arcpy.Array()
# 根据需要，设置字段
PNT_SERIAL = ""
PL_ID = ""
# 遍历读取坐标表中的坐标
for i in range(1, rownum, 1):
row = alldata[i]
X_COORD = row[1]
Y_COORD = row[2]
PNT_SERIAL = row[7]
PL_ID = row[9]
# 对DKID第一次赋值
if DKID == "":
DKID = row[9]
# 生成坐标点
pnt = arcpy.Point()
pnt.ID = int(PNT_SERIAL)
pnt.X = float(Y_COORD)
pnt.Y = float(X_COORD)
# 如果地块ID发生变化，说明包含多个地块
if DKID == PL_ID:
polygonarray.add(pnt)
else:
# 生成面要素
poly = arcpy.Polygon(polygonarray)
# 根据地块ID找到对应的勘测ID
KCID = kcdict[DKID]
# 根据勘测ID获取坐标系类型
COOR_TYPE = xydict[KCID]
# 追加到要素类中
if COOR_TYPE == u"80国家大地坐标系":
cursor_80.insertRow([FLOWSN, XMMC, KCID, DKID, COOR_TYPE, poly])
elif COOR_TYPE == u"2000国家大地坐标系":
cursor_2000.insertRow([FLOWSN, XMMC, KCID, DKID, COOR_TYPE, poly])
else:
cursor_unknown.insertRow([FLOWSN, XMMC, KCID, DKID, COOR_TYPE, poly])
# 将新的地块ID赋值给fid
DKID = row[9]
# 清除掉之前的坐标点，并将当前坐标点添加到polygonarray中
polygonarray.removeAll()
polygonarray.add(pnt)
if polygonarray.count > 0:
# 生成面要素
poly = arcpy.Polygon(polygonarray)
# 根据地块ID找到对应的勘测ID
KCID = kcdict[DKID]
# 根据勘测ID获取坐标系类型
COOR_TYPE = xydict[KCID]
# 根据不同的坐标系，追加到不同的图层中
if COOR_TYPE == u"80国家大地坐标系":
cursor_80.insertRow([FLOWSN, XMMC, KCID, DKID, COOR_TYPE, poly])
elif COOR_TYPE == u"2000国家大地坐标系":
cursor_2000.insertRow([FLOWSN, XMMC, KCID, DKID, COOR_TYPE, poly])
else:
cursor_unknown.insertRow([FLOWSN, XMMC, KCID, DKID, COOR_TYPE, poly])
# 计数，便于统计
count = count + 1
print mdbfile.decode("gb2312")
# 关闭mdb连接
dbconn.close()
# 将变量初始化
FLOWSN = ""
XMMC = ""
KCID = ""
DKID = ""
COOR_TYPE = ""
print "一共上图{0}个...".format(count)
except Exception as e:
print e.args
# 删除图层游标，释放资源
if cursor_2000:
del cursor_2000
if cursor_80:
del cursor_80
if cursor_unknown:
del cursor_unknown
截图如下：
原文地址：https://www.cnblogs.com/apromise/p/11595860.html