Address, buffer, select

[jaykapalczynski]

I am looking for the best way to eventually publish a GP tool that will accept an address, geocode the address, buffer the location X amount of distance and then select all the features it intersects.  Resulting in a list of intersected features.

Any examples?  

[JakeSkinner]

Hey Jay,

Take a look at the Public Notification widget:

Public Notification widget—Web AppBuilder for ArcGIS | Documentation 

You have the option to write the result to a CSV.

[jaykapalczynski]

could I stand this up and then call it from .Net C#, passing it an address and buffer distance and get a JSON list back?

[jaykapalczynski]

this is one cog in a larger wheel here and need to do this grammatically, which is why I am asking if I can call it, pass parameters, and get JSON back....this is the workflow that I need

[jaykapalczynski]

OK its really easy to Geocode an address I send to a python script...which yeilds an XY coordinate.

I am now looking at the Buffer analysis and it appears that I have to input a Feature Class.....Is there a way that I can pass it XY coordinates instead.  I am not finding a reason to create a FC....

I want to Take an XY and create a Temp Buffer and select features from another FC without creating a Point Feature and a Buffer Feature.....Can I do this with temp features?  Or do I have to create a point in a FC when its goecoded, then create a buffer FC from that point FC, then select the features from another FC and get a list of unique IDs.

Only to have to delete the point FC and Buffer FC everything I run the python script

[JakeSkinner]

You could create a point geometry from the XY coordinate, and then create an IN_MEMORY output for the buffer:

point = arcpy.Point(35.723, -75.706)
ptGeometry = arcpy.PointGeometry(point)

arcpy.Buffer_analysis(ptGeometry, "IN_MEMORY/BufferGeom", "100 Feet")‍‍‍‍

[jaykapalczynski]

This is what I have right now.....and it is working...ALTHOUGH I am using my own buffer FC....I still have not figure out how to use the in memory buffer:

arcpy.Buffer_analysis(ptGeometry, "IN_MEMORY/BufferGeom", "50 miles")

I am returning the correct number of results with my fake buffer FC so thats good....I just need to figure out how to implement the in memory buffer that I am trying to create.

I try and add the testBuffer instead of the direct call to my buffer FC....to this but it errors out....

testBuffer = arcpy.Buffer_analysis(ptGeometry, "IN_MEMORY/BufferGeom", "50 miles")

Successful in using FC buffer  

finalResults = arcpy.SelectLayerByLocation_management('counties_lyr', 'INTERSECT', 'buffer_1')

Trying to Use IN memory buffer

finalResults = arcpy.SelectLayerByLocation_management('counties_lyr', 'INTERSECT', 'testBuffer')

import arcpy
import requests
geoCodeUrl = "https://geocode.arcgis.com/arcgis/rest/services/World/GeocodeServer/findAddressCandidates"

searchAddress = "1520 split oak ln, Henrico, VA 23229"
searchDistance = "50"
coordinates = ""

def singleAdressGeocode(address, geoCodeUrl, outSR = "4326"):
  #clean up the address for url encoding
  address = address.replace(" ", "+")
  address = address.replace(",", "%3B")

  #send address to geocode service
  lookup = requests.get(geoCodeUrl + "?SingleLine=" + address + "&outSR=" + outSR + "&maxLocations=1&f=pjson")
  data = lookup.json()
  
  if data["candidates"]:
    #coords = data["candidates"][0]["location"]
    #return coords
    return data
  else:
    #no results
    return "Address not geocoded: " + address


# DEFINE VARIABLES FOR THE ADDRESS AND X AND Y VALUES
#print singleAdressGeocode("1520 split oak ln, Henrico, VA 23229", geoCodeUrl)
geocodeResult = singleAdressGeocode(searchAddress, geoCodeUrl)
addressResult = geocodeResult["candidates"][0]["address"]
coordinateX = geocodeResult["candidates"][0]["location"]["x"]
coordinateY = geocodeResult["candidates"][0]["location"]["y"]

print geocodeResult
print addressResult
print coordinateX
print coordinateY

# CREATE THE POINT FROM THE X AND Y
#point = arcpy.Point(35.723, -75.706)
point = arcpy.Point(coordinateY,coordinateX)
ptGeometry = arcpy.PointGeometry(point)
searchdistancevar = searchDistance + " miles"

print point
print searchdistancevar

# FEED THE BUFFER ANALYSIS THE POINT GEOMETRY THE DISTANCE PARAMETERS
testBuffer = arcpy.Buffer_analysis(ptGeometry, "IN_MEMORY/BufferGeom", "50 miles")

arcpy.env.workspace = "C:/Users/SelectCounties.gdb"
arcpy.MakeFeatureLayer_management('Counties', 'counties_lyr') 
finalResults = arcpy.SelectLayerByLocation_management('counties_lyr', 'INTERSECT', 'buffer_1')

# If features matched criteria write them to a new feature class
matchcount = int(arcpy.GetCount_management('counties_lyr')[0]) 
if matchcount == 0:
    print('no features matched spatial and attribute criteria')
else:
    print('{0} counties that matched criteria written to {0}'.format(
                                                  matchcount))



‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

[JakeSkinner]

Here is an example on how to call the in memory buffer:

finalResults = arcpy.SelectLayerByLocation_management('counties_lyr', 'INTERSECT', 'IN_MEMORY/BufferGeom')‍

[jaykapalczynski]

Thanks....although something is not right....below is what I am doing....

If I use the FC buffer (comment out below I get 23 features returned)

If I use the IN memory Buffer I get 0 and the 'no features matched spatial and attribute criteria'

Can you see anything I am doing wrong?

point = arcpy.Point(37.723, -77.706)
ptGeometry = arcpy.PointGeometry(point)

arcpy.Buffer_analysis(ptGeometry, "IN_MEMORY/BufferGeom", "50 miles")

arcpy.env.workspace = "C:/Users/Data/SelectCounties.gdb"

arcpy.MakeFeatureLayer_management('CountiesJurisdictions', 'counties_lyr')

# USING FC BUFFER
#results = arcpy.SelectLayerByLocation_management('counties_lyr', 'INTERSECT', 'buffer')

# USING IN MEMORY BUFFER
results = arcpy.SelectLayerByLocation_management('counties_lyr', 'INTERSECT', 'IN_MEMORY/BufferGeom')

# GET THE COUNT OF THE RETURNED VALUES    
# If features matched criteria write them to a new feature class
matchcount = int(arcpy.GetCount_management('counties_lyr')[0]) 
if matchcount == 0:
    print('no features matched spatial and attribute criteria')
else:
    print('{0} counties that matched criteria written to {0}'.format(
                                                  matchcount))‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

[jaykapalczynski]

I've tried the below.....got rid of the " " and replaced with ' ' which I believe is the proper syntax....BUT still nothing...

arcpy.Buffer_analysis(ptGeometry, 'IN_MEMORY/BufferGeom', '50 miles')

arcpy.Buffer_analysis(ptGeometry, 'IN_MEMORY/BufferGeom', '50')

Maybe a spatial reference issue?  I think the County Layer its trying to select from is in UTM....going to project it to decimal degrees and see if that works?