zoukankan      html  css  js  c++  java
  • arcgis python xlstoshp

    import xlrd # must init xlrd 
    import arcpy
    
    # param
    arcpy.env.workspace = r"F:
    otepythonArcPy" #workspace
    excelPath = r"test.xlsx" # excel file path
    excelTableIndex = 0 # excel's table index
    outName = r"point.shp" # out file
    
    excel = xlrd.open_workbook('test.xlsx') # get excel
    table = excel.sheets()[0] # get table by sheets index
    nrows = table.nrows # number of table's row
    # get data
    pointGeometryList = [] # a list to hold the PointGeometry objects
    point = arcpy.Point() #create an empty Point object
    spRef = arcpy.SpatialReference('WGS84 ARC System Zone 18')
    
    for i in range(1,nrows): # get row once
        x = table.cell(i,0).value 
        y = table.cell(i,1).value
        point.X = float(x)
        point.Y = float(y)
        pointGeometry = arcpy.PointGeometry(point,spRef)
        pointGeometryList.append(pointGeometry)
    
    arcpy.CopyFeatures_management(pointGeometryList,outName) # save the shape file
    --------------------- 
    版权声明:本文为CSDN博主「PasserQi」的原创文章,遵循CC 4.0 by-sa版权协议,转载请附上原文出处链接及本声明。
    原文链接:https://blog.csdn.net/summer_dew/article/details/78116722

    =======================================修改===================================================
    import xlrd # must init xlrd
    import arcpy
    import os
    
    # param
    arcpy.env.workspace = r"F:kktemp" #workspace
    excelPath = r"sss.xlsx" # excel file path
    excelTableIndex = 0 # excel's table index
    outName = r"F:kktemppoint.shp" # out file
    
    fileName=r'F:kktempsss.xlsx'
    if os.path.exists(fileName):
        excel = xlrd.open_workbook(fileName) # get excel
        table = excel.sheets()[0] # get table by sheets index
        nrows = table.nrows # number of table's row
        # get data
        pointGeometryList = [] # a list to hold the PointGeometry objects
        point = arcpy.Point() #create an empty Point object
        #spRef = arcpy.SpatialReference('Xian_1980_3_Degree_GK_CM_111E')
    
        for i in range(1,nrows): # get row once
            x = table.cell(i,0).value
            y = table.cell(i,1).value
            point.X = float(x)
            point.Y = float(y)
            pointGeometry = arcpy.PointGeometry(point)
            pointGeometryList.append(pointGeometry)
    
        arcpy.CopyFeatures_management(pointGeometryList,outName) # save the shape file
  • 相关阅读:
    消失的 unittest.makeSuite()
    自定义错误信息在各个浏览器表现不同
    “创建Web应用程序”与“创建Web网站”的区别
    如何将Notepad++改造成Python开发工具
    64位服务器IIS不能识别32位framework版本。IIS没有Asp.net切换界面的解决办法。
    C++文件流,读入数据注意
    Chessboard(二分图完备匹配)
    courses二分图最大匹配
    Linux BASH Shell文件名匹配/输出重定向
    LinuxShell一些很少用到却很有用的指令
  • 原文地址:https://www.cnblogs.com/gisoracle/p/11354945.html
Copyright © 2011-2022 走看看