批量生成采矿_挖空矿区_xls字符串

# -*-coding:utf-8-*-
import arcpy
import xlrd
import re
 
# 采矿权坐标格式举例
# 1,4,1,4145796.42,37500324.70,2,4145822.42,37500429.70,3,4145717.42,37500429.70,4,4145713.42,37500349.70,1167,1090,,1,
#2,4,1,3950410.72,19466190.96,2,3950212.72,19466470.96,3,3950144.72,19466415.96,4,3950344.72,19466140.96,750,550,,1,4,1,3949944.71,19466065.96,2,3949804.71,19466265.96,3,3949254.71,19465995.96,4,3949304.71,19465715.95,750,550,,1,
 
# 读取xls文件
xlsfile = xlrd.open_workbook(r"E:\ck.xls", "r")
mysheet1 = xlsfile.sheet_by_name("ck")
rownum = mysheet1.nrows
colnum = mysheet1.ncols
 
# 生成的shp图形存放目录
arcpy.env.workspace = r"E:\shp"
 
# 2000坐标系
fc_2000_19 = "ck2000_19.shp"
fc_2000_37 = "ck2000_37.shp"
fc_2000_38 = "ck2000_38.shp"
 
# 80坐标系
fc_80_19 = "ck80_19.shp"
fc_80_37 = "ck80_37.shp"
fc_80_38 = "ck80_38.shp"
 
# 54坐标系
fc_54_19 = "ck54_19.shp"
fc_54_37 = "ck54_37.shp"
fc_54_38 = "ck54_38.shp"
 
# 未知坐标系
fc_unknown_19 = "unknown_19.shp"
fc_unknown_37 = "unknown_37.shp"
fc_unknown_38 = "unknown_38.shp"
 
# 检查工作空间下是否存在2000坐标系的要素类
isexist_fc_2000_19 = arcpy.Exists(fc_2000_19)
isexist_fc_2000_37 = arcpy.Exists(fc_2000_37)
isexist_fc_2000_38 = arcpy.Exists(fc_2000_38)
 
# 检查工作空间下是否存在80坐标系的要素类
isexist_fc_80_19 = arcpy.Exists(fc_80_19)
isexist_fc_80_37 = arcpy.Exists(fc_80_37)
isexist_fc_80_38 = arcpy.Exists(fc_80_38)
 
# 检查工作空间下是否存在54坐标系的要素类
isexist_fc_54_19 = arcpy.Exists(fc_54_19)
isexist_fc_54_37 = arcpy.Exists(fc_54_37)
isexist_fc_54_38 = arcpy.Exists(fc_54_38)
 
# 检查工作空间下是否存在未知坐标系的要素类
isexist_fc_unknown_19 = arcpy.Exists(fc_unknown_19)
isexist_fc_unknown_37 = arcpy.Exists(fc_unknown_37)
isexist_fc_unknown_38 = arcpy.Exists(fc_unknown_38)
 
# 记录没有匹配的记录
str_notmatch = ""
 
if not (
        isexist_fc_2000_19 and isexist_fc_2000_37 and isexist_fc_2000_38 and isexist_fc_80_19 and isexist_fc_80_37 and isexist_fc_80_38 and isexist_fc_54_19 and isexist_fc_54_37 and isexist_fc_54_38 and isexist_fc_unknown_19 and isexist_fc_unknown_37 and isexist_fc_unknown_38):
    print "要素类不存在!"
    exit()
 
# 创建2000坐标系各要素类的插入游标,许可证号xkzh,名称name,区域坐标coord
cursor_2000_19 = arcpy.da.InsertCursor(fc_2000_19, ["xkzh", "name", "coord", "SHAPE@"])
cursor_2000_37 = arcpy.da.InsertCursor(fc_2000_37, ["xkzh", "name", "coord", "SHAPE@"])
cursor_2000_38 = arcpy.da.InsertCursor(fc_2000_38, ["xkzh", "name", "coord", "SHAPE@"])
 
# 创建80坐标系各要素类的插入游标,许可证号xkzh,名称name,区域坐标coord
cursor_80_19 = arcpy.da.InsertCursor(fc_80_19, ["xkzh", "name", "coord", "SHAPE@"])
cursor_80_37 = arcpy.da.InsertCursor(fc_80_37, ["xkzh", "name", "coord", "SHAPE@"])
cursor_80_38 = arcpy.da.InsertCursor(fc_80_38, ["xkzh", "name", "coord", "SHAPE@"])
 
# 创建54坐标系各要素类的插入游标,许可证号xkzh,名称name,区域坐标coord
cursor_54_19 = arcpy.da.InsertCursor(fc_54_19, ["xkzh", "name", "coord", "SHAPE@"])
cursor_54_37 = arcpy.da.InsertCursor(fc_54_37, ["xkzh", "name", "coord", "SHAPE@"])
cursor_54_38 = arcpy.da.InsertCursor(fc_54_38, ["xkzh", "name", "coord", "SHAPE@"])
 
# 创建未知坐标系各要素类的插入游标,许可证号xkzh,名称name,区域坐标coord
cursor_unknown_19 = arcpy.da.InsertCursor(fc_unknown_19, ["xkzh", "name", "coord", "SHAPE@"])
cursor_unknown_37 = arcpy.da.InsertCursor(fc_unknown_37, ["xkzh", "name", "coord", "SHAPE@"])
cursor_unknown_38 = arcpy.da.InsertCursor(fc_unknown_38, ["xkzh", "name", "coord", "SHAPE@"])
 
try:
    # 遍历xls文件中所有行
    for i in range(rownum):
        row = mysheet1.row_values(i)
 
        # 读取基本信息:许可证号、名称、坐标系类型、区域坐标,这些基本信息主要根据xls文件来定
        xkzh = row[0]
        name = row[1]
        xytype = str(row[2])
        coord = row[3]
        # print (name + ",is beginging...")
 
        # 调试代码用
        if i == -1:
            print "!!!"
 
        # 外圈坐标,即外圈矿山
        _xypolylist = []
        # 挖空坐标,即挖空矿山
        _wkpolylist = []
        # 临时坐标
        _tempxylist = arcpy.Array()
        # 坐标点的编号
        bh = 0
        # 投影坐标的带号
        dh = 0
 
        # 获取区块数num_qk,字符串第1个逗号前的数字表示区块数
        index_qk = coord.find(",")
        num_qk = int(coord[:index_qk])
 
        # 除过区块数,剩余的字符串
        coordlist = coord[index_qk + 1:]
 
        # 根据正则表达式,获取各区块结束的字符串标识       
        pattern = re.compile(
            r"([-]{0,1}[0-9]{1,4}[\.]{0,1}[0-9]{0,3}[,][-]{0,1}[0-9]{1,4}[\.]{0,1}[0-9]{0,3}[,][.]{0,}[\d]{0,}[\u4E00-\u9FFF]{0,}[,][-]{0,1}[0-1],)",
            re.U)
        reglist = re.findall(pattern, coordlist)       
 
        # 如果没有匹配出坐标结束标识,则做出记录
        num_reglist = len(reglist)
        if num_reglist == 0:
            # print name + " is not matching..."
            print xkzh
            str_notmatch = r"'" + xkzh + "'," + str_notmatch
            continue
 
        # 如果匹配出的坐标结束标识的数量与区块数不一致,则做出记录
        if num_reglist != num_qk:
            # print name + " qks is not equal ..."
            print xkzh
            str_notmatch = r"'" + xkzh + "'," + str_notmatch
            continue
 
        for j in range(num_qk):
 
            # 获得区块结束的字符串标识
            str_reg = str(reglist[j])
 
            # 如果匹配出的标识符不是4个,则说明匹配错了
            num_dot = str_reg.count(",")
            if num_dot != 4:
                break
 
            # 获得索引值index_regend
            index_regend = coordlist.find(str_reg)
            qk_coordlist = coordlist[:index_regend]
 
            # 获取坐标点数num_coord,字符串中第1个逗号前的数值表示坐标点数。(已经去掉了区块数)
            index_coord = qk_coordlist.find(",")
            num_coord = int(qk_coordlist[:index_coord])
 
            # 获取是否挖空的标识
            wk_info = str_reg.split(",")
            wk_mark = wk_info[3]
 
            # 获得坐标点
            qk_coordlist = qk_coordlist[index_coord + 1:]
            xylist = qk_coordlist.split(",")
 
            # 编号、投影带号、横坐标、纵坐标
            bh = 0
            dh = 0
            _x = 0
            _y = 0
 
            for k in range(num_coord * 3):
 
                # 坐标点标识,纵坐标(7位整数)、横坐标(8位整数)
                if k % 3 == 1:
 
                    # 检查纵坐标是否满足7位
                    y = xylist[k].format("0000000.###")
                    yindex_dot = y.find(".")
                    if yindex_dot == -1:
                        if len(y) != 7:
                            print xkzh + " ylength is not enough..."
                            break
                    else:
                        if len(y[:yindex_dot]) != 7:
                            print xkzh + " ylength is not enough..."
                            break
                    _y = float(y)
 
                elif k % 3 == 2:
 
                    # 检查横坐标是否满足8位
                    x = xylist[k].format("00000000.###")
                    xindex_dot = x.find(".")
                    if xindex_dot == -1:
                        if len(x) != 8:
                            print xkzh + " xlength is not enough..."
                            break
                    else:
                        if len(x[:xindex_dot]) != 8:
                            print xkzh + " xlength is not enough..."
                            break
                    _x = float(x)
 
                    # 获取投影带号
                    dh = str(x[:2])
 
                    # 读取坐标值
                    pnt = arcpy.Point()
                    # 依次为坐标点编号、纵坐标、横坐标
                    bh = bh + 1
                    pnt.ID = bh
                    pnt.X = _x
                    pnt.Y = _y
                    _tempxylist.append(pnt)
                else:
                    continue
 
            # 如果坐标点数为0,则说明该坐标有问题(不满足前7后8),继续一个矿权
            if len(_tempxylist) == 0:
                break
 
            # 根据挖空标识不同,将外圈、内圈坐标存放在不同的列表中
            _xytemppolygon = arcpy.Polygon(_tempxylist)
            if wk_mark == "1" or wk_mark == "0":
                _xypolylist.append(_xytemppolygon)
            elif wk_mark == "-1":
                _wkpolylist.append(_xytemppolygon)
            else:
                print name + "is error..."
 
            # 清楚临时列表
            _tempxylist.removeAll()
 
            # 剩余的坐标串
            len_reg = len(str_reg)
            coordlist = coordlist[index_regend + len_reg:]
 
        # 获取外圈的数量,如果等于0则继续下一个循环
        xypolynum = len(_xypolylist)
        if xypolynum == 0:
            continue
 
        wkpolynum = len(_wkpolylist)
        # 如果外圈矿山只有1个,挖空矿山1个或者多个,则执行裁剪,即交集取反
        if xypolynum == 1 and wkpolynum >= 1:
            poly = _xypolylist[0]
            for p in range(wkpolynum):
                poly = poly.symmetricDifference(_wkpolylist[p])
 
            # 2000坐标系
            if xytype == "3":
                if dh == "37":
                    cursor_2000_37.insertRow([xkzh, name, coord, poly])
                elif dh == "38":
                    cursor_2000_38.insertRow([xkzh, name, coord, poly])
                elif dh == "19":
                    cursor_2000_19.insertRow([xkzh, name, coord, poly])
                else:
                    print xkzh + " dh is error..."
 
            # 80坐标
            elif xytype == "2":
                if dh == "37":
                    cursor_80_37.insertRow([xkzh, name, coord, poly])
                elif dh == "38":
                    cursor_80_38.insertRow([xkzh, name, coord, poly])
                elif dh == "19":
                    cursor_80_19.insertRow([xkzh, name, coord, poly])
                else:
                    print xkzh + " dh is error..."
 
            # 54坐标系
            elif xytype == "1":
                if dh == "37":
                    cursor_54_37.insertRow([xkzh, name, coord, poly])
                elif dh == "38":
                    cursor_54_38.insertRow([xkzh, name, coord, poly])
                elif dh == "19":
                    cursor_54_19.insertRow([xkzh, name, coord, poly])
                else:
                    print xkzh + " dh is error..."
 
            # 坐标系未知
            else:
                if dh == "37":
                    cursor_unknown_37.insertRow([xkzh, name, coord, poly])
                elif dh == "38":
                    cursor_unknown_38.insertRow([xkzh, name, coord, poly])
                elif dh == "19":
                    cursor_unknown_19.insertRow([xkzh, name, coord, poly])
                else:
                    print xkzh + " dh is error..."
            # print name + " is finished!"
            continue
 
        # 对于多个外圈矿山,1个或者多个挖空矿山,无法判断对哪个外圈矿山挖空
        if xypolynum > 1 and wkpolynum >= 1:
            print (name + ",无法判断挖空矿山!")
            continue
 
        # 遍历形成外圈地块,一个矿山无论是否存在多个区块,只能是一种坐标系
        for q in range(xypolynum):
 
            # 2000坐标系
            if xytype == "3":
                if dh == "37":
                    cursor_2000_37.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "38":
                    cursor_2000_38.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "19":
                    cursor_2000_19.insertRow([xkzh, name, coord, _xypolylist[q]])
                else:
                    print xkzh + " dh is error..."
 
            # 80坐标
            elif xytype == "2":
                if dh == "37":
                    cursor_80_37.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "38":
                    cursor_80_38.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "19":
                    cursor_80_19.insertRow([xkzh, name, coord, _xypolylist[q]])
                else:
                    print xkzh + " dh is error..."
 
            # 54坐标系
            elif xytype == "1":
                if dh == "37":
                    cursor_54_37.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "38":
                    cursor_54_38.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "19":
                    cursor_54_19.insertRow([xkzh, name, coord, _xypolylist[q]])
                else:
                    print xkzh + " dh is error..."
 
            # 坐标系未知
            else:
                if dh == "37":
                    cursor_unknown_37.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "38":
                    cursor_unknown_38.insertRow([xkzh, name, coord, _xypolylist[q]])
                elif dh == "19":
                    cursor_unknown_19.insertRow([xkzh, name, coord, _xypolylist[q]])
                else:
                    print xkzh + " dh is error..."
 
        # print "It is finished!"
except Exception as err:
    # print (err.args[0]).decode("gbk")
    print err.message
else:
    print "全部完成!"
    print str_notmatch
    del cursor_2000_37
    del cursor_2000_38
    del cursor_2000_19
    del cursor_80_37
    del cursor_80_38
    del cursor_80_19
    del cursor_54_37
    del cursor_54_38
    del cursor_54_19
    del cursor_unknown_37
    del cursor_unknown_38
    del cursor_unknown_19
 

转载于:https://www.cnblogs.com/apromise/p/11496404.html

你可能感兴趣的:(批量生成采矿_挖空矿区_xls字符串)