zoukankan      html  css  js  c++  java
  • ArcGIS将shp按照属性字段进行分割为多个polygon矢量

     将一个shp中所有的要素都导出为shp

    流域.shp属性表:

    import arcpy
    import os
    from arcpy import env
    
    shp ="D:/01RiverPro/01DATA/流域.shp" #要处理的shp文件路径
    out_path="D:/01RiverPro/01DATA/子流域"#输出文件夹的路径
    
    with arcpy.da.SearchCursor(shp, ["SHAPE@",'ID']) as cursor:
    #SHAPE@指代单个要素,ID是一个字段,该字段也是我们想要作为每个polygon命名的值,也可以改为其他的字段
        for row in cursor:
            out_name=str(row[1])+'.shp'#输出文件名 确保均为字符串类型,注意命名是唯一的,这样就能将该shp数据中的所有要素都导出 
            arcpy.FeatureClassToFeatureClass_conversion (row[0],out_path,out_name)

     注意:

    row[0]代表的就是polygon要素

    row[1]表示的是ID字段,最好确保该字段的唯一性

    修改的时候,不需要该row[0] row[1] 只需要改‘ID’为其他即可

    结果:

    地图展示:

    另外一种方式:

     参考:https://www.cnblogs.com/yhpan/p/13556106.html

    import re
    import arcpy
    from arcpy import env
    
    def validateTitle(title):
        rstr = r"[/\:*?"<>|\,. ]"  
        new_title = re.sub(rstr, "_", title)  
        return new_title
    
    
    env.workspace = "D:/01RiverPro/01DATA/01Headwater/hill_watershed3/test"
    
    in_features = "D:/01RiverPro/01DATA/01Headwater/hill_watershed3/test/三级河网.shp"
    
    rows = arcpy.SearchCursor(in_features, fields="GRID_CODE")
    
    #遍历每一行 对符合要求的polygon进行分割
    for row in rows:
        selected_field_value = row.getValue("GRID_CODE")
        where_clause = '"GRID_CODE" =3'
        #outfilename = validateTitle(selected_field_value)
        out_feature_class = "D:/01RiverPro/01DATA/01Headwater/hill_watershed3/test/main.shp"
        arcpy.Select_analysis(in_features, out_feature_class, where_clause)
        print(selected_field_value+' has done!!!')
    
    #将 GRID_CODE字段作为名字 如果是“name” 则需改为%s
    for row in rows:
        selected_field_value = row.getValue("GRID_CODE")
        where_clause = '"GRID_CODE" = %x'%(selected_field_value)
        #outfilename = validateTitle(selected_field_value)
        out_feature_class = "D:/01RiverPro/01DATA/01Headwater/hill_watershed3/test/%x.shp"%(selected_field_value)
        arcpy.Select_analysis(in_features, out_feature_class, where_clause)
        #print(selected_field_value+' has done!!!')

    参考:

    https://blog.csdn.net/qq_36196621/article/details/90171166

  • 相关阅读:
    关于springMVC转换json出现的异常
    jQuery实现,动态自动定位弹窗。JS分页,Ajax请求
    servlet为什么要配置web.xml
    Jmeter系列(4)- Jmeter 脚本录制
    后缀数组模板
    NOIP2016 玩脱记
    TERSUS无代码开发(笔记21)-流程执行顺序思考(转载)
    ===>===>===>特色思TERSUS常用功能整理
    TERSUS无代码开发(笔记20)-本地开发测试mysql数据库连接
    TERSUS无代码开发(笔记19)-mysql-connector-java-5.-bin.jar下载方法
  • 原文地址:https://www.cnblogs.com/icydengyw/p/13556201.html
Copyright © 2011-2022 走看看