• 博客园logo
  • 会员
  • 众包
  • 新闻
  • 博问
  • 闪存
  • 赞助商
  • HarmonyOS
  • Chat2DB
    • 搜索
      所有博客
    • 搜索
      当前博客
  • 写随笔 我的博客 短消息 简洁模式
    用户头像
    我的博客 我的园子 账号设置 会员中心 简洁模式 ... 退出登录
    注册 登录

gisoracle

  • 博客园
  • 联系
  • 订阅
  • 管理

公告

View Post

arcgis python 把一个图层数据合并在一起

def getLayerGeomtry(inFeature):
   
    num=getCount(inFeature)
    
    mytemp=inFeature
    if num>1:#融合在一起
        mytemp="in_memory/gisoracle"
        arcpy.Dissolve_management(inFeature, mytemp,"", "", "MULTI_PART")
        
    desc = arcpy.Describe(mytemp)
    shapefieldname = desc.ShapeFieldName
    
    rows = arcpy.SearchCursor(mytemp)
   
    for row in rows:
        pgeometry=row.getValue(shapefieldname)
      
    if row:
        del row
    if rows:
        del rows
        
    return pgeometry
#到Pro中使用
def getLayerGeomtryPro(inFeature):
    desc = arcpy.Describe(inFeature)
    shapefieldname = desc.ShapeFieldName
    rows = arcpy.SearchCursor(inFeature)
    i=0
    for row in rows:
        if i<1:
            pgeometry=row.getValue(shapefieldname)
        else:
            pgeometry=pgeometry.union(row.getValue(shapefieldname))
            
        i=i+1
    if row:
        del row
    if rows:
        del rows
        
    return pgeometry

 

pgeometry=getLayerGeomtry(inFeature)
yl_tmp = arcpy.CreateUniqueName("yl_tmp") #临时数据
arcpy.Select_analysis(pgeometry, yl_tmp, '')

posted on 2019-06-04 11:33  gisai  阅读(716)  评论(0)    收藏  举报

刷新页面返回顶部
 
博客园  ©  2004-2025
浙公网安备 33010602011771号 浙ICP备2021040463号-3