I have two shapefiles that I am a trying to compare to one another. The first set of points are various businesses. With assistance it was discovered that in order to get a subset category from that shapefile the following code should be used:
cat = 'ClothingStore' arcypy.SelectLayerByAttribute_management("lyr","NEW_SELECTIION",'"category"='+"'%s'"% cat)
After I have isolated the subset category from the shapefile I want to compare that subset of points to another set of points that are all within a single shapefile called Churches.shp. I want to find all the ClothingStores points that are within a 500 radius of Church points.
How could this be done through ArcPy?