Como dividir pontos por polígono usando o ArcGIS Desktop?

Estou procurando uma maneira de dividir pontos por uma classe de recurso de polígono do condado. A imagem anexada mostra duas classes de recurso: uma composta por pontos e a outra uma classe de recurso de polígono do município. A ferramenta de divisão seria ideal se permitisse pontos. Eu preferiria uma saída limpa sem campos adicionais adicionados aos atributos de ponto. Além disso, gostaria de ter o nome do condado da classe de recurso de polígono que define as novas classes de recurso de ponto (por exemplo, Kiowa, Clark, Comanche). Agradeço quaisquer soluções e conselhos.

Neste exemplo, o produto final deve ser de três classes de feições pontuais denominadas "Kiowa", "Clark" e "Comanche", produzidas a partir de uma classe de feições pontuais maior. Uma solução automatizada seria ideal, pois tenho muitos FCs de pontos mesclados em dezenas de municípios para trabalhar (abrangendo aproximadamente o estado de KS).

Depois de usar Junção espacial ou Intersecção para obter o atributo Nome do condado em cada ponto, tente usar a ferramenta Camada dividida por atributos de Dan Patterson, disponível na Galeria de ferramentas de script e modelo de geoprocessamento .

Como alternativa, você pode usar o ModelBuilder para automatizar isso usando uma abordagem diferente envolvendo Selecionar camada por local (clique na miniatura para obter uma imagem completa).

Retirado desse segmento nos fóruns da ESRI: Seleção em lote por local e exportação de shapefiles

Muito mais fácil de usar o comando intersect.
Desative os campos indesejados.
1. Cruze para uma nova saída.
2. Junte a saída com o oid original.
3. Calcule um novo campo com o nome cnty.

Confira o ettools para obter algumas versões aprimoradas de junção espacial e divida por local.

Pelo que li, eu executaria uma sobreposição (ArcToolbox Analysis Tools - Overlay - e Identity, Intersect ou Spatial Join) de seus pontos com um conjunto de dados de polígonos do condado. No ArcMap, você poderia limitar os campos exibidos na camada do condado ao campo com o nome do condado, e somente esse campo seria adicionado ao resultado. Se mesmo esse campo adicionado for muito grande, calcule seus valores no campo correspondente na tabela de pontos e exclua o campo adicionado.

Usando Python, você pode fazer um loop em "Select Layer By Location" para pegar os pontos que se sobrepõem a cada país e salvá-lo como um Shapefile. Seria algo como isto:

import arcgisscripting
# starts geoprocessing
gp = arcgisscripting.create()
gp.OverWriteOutput = 1

# Variable iniciation
points = u"Path for your point shape"
counties = u"Path for counties shape"
outDir = u"path for output directory"

#Load points as a layer
gp.MakeFeatureLayer(points, "points")

# Go county by county
rows = gp.searchcursor(counties)
row = rows.next()

# loop County by county overlap
while row.countyName: #adapt to your countie table of attributes
    # Make a layer from the feature class
    gp.MakeFeatureLayer(counties, "counties lyr", "[countieName]='" + row.countieName + "'")
    # Select all points that intersect the current countie polygon
    gp.SelectLayerByLocation(points, "intersect", "counties_lyr", 0, "new_selection")
    outSHP = outDir + row.counties + u".shp"
    gp.CopyFeatures_management(points, outSHP)
    print(outSHP)
    row = rows.next()
 # End of loop

 del rows, row, gp

 #END

Eu não testei, você provavelmente precisa adaptá-lo ao seu caso.

Edit: Alterado Select_analysis por CopyFeatures_management seguindo @ blah238 conselhos.

Agora existe um tutorial em vídeo do eLearning (não gratuito, mas relativamente barato) chamado Usando os Cursores do ArcPy para Dividir por Atributo e Local que inclui mostrar como:

Escreva um script Python para dividir uma classe de recurso usando uma classe de recurso de polígono em uma classe de recurso por recursos encontrados em cada um de seus locais de polígono

Eu acho que aborda os requisitos desta pergunta.

Isenção de responsabilidade : sou o autor e o apresentador deste tutorial em vídeo do eLearning