draw a polygon with a python add-in and clip it

EstelleWuilliez1

Hello everybody, I really need some help for my project in ArcMap 10.2.

I want to create a tool with a python add-in, then click on it and draw a polygon on a map. Then I want to use this polygon to clip my map and to keep what I've selected with the polygon.

To do that, I've made this code in my Add-in :

Everything seems to work except one thing! The clip tool creates a empty output, I've got the error "Empty output generated", and I really don't understand why. My "Route_shapes.shp" and my "polygon.shp" have the same coordinate system.

Any suggestion?

Thank you!

XanderBakker

Interesting, but not what I was hoping for...

I suppose that you do have access to an Advanced license, right? I noticed that the "Feature To Polygon" tool requires that type of license. I do not have one available at the moment. I did rewrite the code to work as stand alone, so I can test it easier and to solve the "Feature To Polygon", in the code I create an empty featureclass (lines 17 - 25) and write the polygon to it (lines 28 - 29). The line geometry is created on lines 34 - 39.

import arcpy
arcpy.env.overwriteOutput = True

def testing(line_geometry):
    fc_routes = r"D:\Xander\GeoNet\ClipTool\ROUTE_SHAPE.shp"
    fc_graphics = r"D:\Xander\GeoNet\ClipTool\GRAPHICS.shp"
    fc_result = r"D:\Xander\GeoNet\ClipTool\POLYGON.shp"

    array = arcpy.Array()
    part = line_geometry.getPart(0)
    for pt in part:
        array.add(pt)
    array.add(line_geometry.firstPoint)

    polygon = arcpy.Polygon(array, fc_routes)

    # create empty shapefile
    geometry_type = "POLYGON"
    template = ""
    has_m = "DISABLED"
    has_z = "DISABLED"
    sr = arcpy.Describe(fc_routes).spatialReference
    ws_path, fc_out_name = os.path.split(fc_graphics)
    arcpy.CreateFeatureclass_management(ws_path, fc_out_name,
        geometry_type, template, has_m, has_z, sr)

    # FeatureToPolygon_management requires Advanced
    with arcpy.da.InsertCursor(fc_graphics, ("SHAPE@")) as curs:
        curs.insertRow((polygon,))

    arcpy.Clip_analysis(fc_routes, fc_graphics, fc_result)


# create a line geometry
array = arcpy.Array([arcpy.Point(467075, 211310),
                     arcpy.Point(467550, 214930),
                     arcpy.Point(470900, 212561)])
polyline = arcpy.Polyline(array)
line_geometry = arcpy.Polyline(array)


# call the clip function
testing(line_geometry)

This did work without problems:

Can you check your graphics shapefile to see if it has the desired content?

Kind regards, Xander