Selecione e copie recursos no ArcMap usando a ferramenta de suplemento Python

Estou tentando adicionar uma ferramenta a um toobar de suplemento para selecionar recursos de uma classe de recurso existente e copiá-los para outra classe de recurso. Preciso de duas ferramentas, uma para um ponto e outra para um retângulo. As ferramentas de suplemento lidam apenas com retângulos, então estou tentando gerar a coordenada do ponto fora do evento onMouseDown e simplesmente usar o feedback X e Y que recebo. O problema é que não consigo passar essas coordenadas para algo útil para extrair os recursos. Eu tentei usar extensões de ambiente e, em seguida, executar um CopyFeatures simples, mas isso não está funcionando, pois copia todo o fc e parece não respeitar as configurações de ambiente definidas.

Posso passar coordenadas para a ferramenta Selecionar por local de alguma forma ou existe outra maneira de passar coordenadas para algo a ser extraído nessa extensão?

Isto é o que eu tenho agora:

import arcpy
import pythonaddins
arcpy.overWriteOutput = True

# Replace this with the SDE layer once it's available.
global ELAtemplate, fc
ELAtemplate = r'C:\Data\nsw_map_units.shp'
fc = ""

class DefineUnitsbyPoint(object):
    """Implementation for DefineUnitsbyPoint.tool (Tool)"""
    def __init__(self):
        self.enabled = True
        self.shape = "Rectangle" # Use onMouseDown to get initial extent of the rectangle.
    def onMouseDownMap(self, x, y, button, shift):
        # fc can be assigned from a combo box selection in another class
        global ELAtemplate, fc
        if fc == "":
            pythonaddins.MessageBox('Choose a layer to select from.', 'Choose a Layer', 0)
        else:
            mxd = arcpy.mapping.MapDocument("CURRENT")
            pointGeom = arcpy.PointGeometry(arcpy.Point(x,y), mxd.activeDataFrame.spatialReference)
            arcpy.SelectLayerByLocation_management(ELAtemplate, "INTERSECT", pointGeom, "", "ADD_TO_SELECTION")
            arcpy.CopyFeatures(ELAtemplate, fc)

Você provavelmente desejará usar, em onMouseDownMapvez de onMouseDownretornar a localização nas coordenadas do mapa, não nas coordenadas da janela.

Além disso, certifique-se de passar um SpatialReferenceobjeto válido para o PointGeometryconstrutor, caso contrário, provavelmente não funcionará. No exemplo abaixo, uso a referência espacial do quadro de dados ativo.

Por fim, convém especificar um search_distanceno seu, SelectLayerByLocationpara que os recursos de ponto e linha possam ser selecionados sem clicar exatamente neles. No ArcObjects, você normalmente usaria a tolerância de seleção do ArcMap em pixels e expandiria o envelope do seu ponto por esse valor nas coordenadas do mapa. Não consegui encontrar uma maneira de acessar a configuração de tolerância de seleção do ArcMap no arcpy, mas se você quiser usar o padrão de 3 pixels (ou passar por conta própria), poderá passar a saída da função nesta resposta como um search_distance( em polegadas) para SelectLayerByLocation.

def onMouseDownMap(self, x, y, button, shift):
    mxd = arcpy.mapping.MapDocument("CURRENT")
    pointGeom = arcpy.PointGeometry(arcpy.Point(x, y), mxd.activeDataFrame.spatialReference)
    searchdistance = getSearchDistanceInches(mxd.activeDataFrame.scale)
    lyr = arcpy.mapping.ListLayers(mxd)[0] # assumes you want to select features from 1st layer in TOC
    arcpy.SelectLayerByLocation_management(lyr, "INTERSECT", pointGeom, "%d INCHES" % searchdistance)
    arcpy.RefreshActiveView()

Suponho que você esteja usando as técnicas descritas em Resposta a É possível obter as coordenadas de um ponto clicado pelo usuário em um MXD atual com o ArcGIS 10.0 ArcPy? hoje cedo ou algo semelhante.

Depois de criar sua geometria de ponto ou polígono (retângulo) usando algo como o código abaixo, você poderá usá-lo como recursos selecionados para Selecionar camada por local.

pt=arcpy.PointGeometry(arcpy.Point(x,y))

A partir daí, os Recursos de cópia devem obter os recursos que você procura.