topic Re: Address, buffer, select in ArcGIS API for Python Questions

Address, buffer, select

jaykapalczynski — Mon, 23 Mar 2020 16:04:05 GMT

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?

Re: Address, buffer, select

JakeSkinner — Mon, 23 Mar 2020 16:21:54 GMT

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.

Re: Address, buffer, select

jaykapalczynski — Mon, 23 Mar 2020 16:28:30 GMT

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?

Re: Address, buffer, select

jaykapalczynski — Mon, 23 Mar 2020 16:40:38 GMT

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

Re: Address, buffer, select

jaykapalczynski — Mon, 23 Mar 2020 19:18:12 GMT

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

Re: Address, buffer, select

JakeSkinner — Sun, 12 Dec 2021 08:39:43 GMT

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")‍‍‍‍

Re: Address, buffer, select

JakeSkinner — Tue, 24 Mar 2020 09:33:11 GMT

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

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

Re: Address, buffer, select

jaykapalczynski — Sun, 12 Dec 2021 08:39:48 GMT

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))‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Address, buffer, select

jaykapalczynski — Tue, 24 Mar 2020 17:56:35 GMT

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?

Re: Address, buffer, select

jaykapalczynski — Tue, 24 Mar 2020 19:54:24 GMT

I believe it was the projection on my Counties....passing DD and data was in UTM....

think I got it...

will be posting all code shortly

Re: Address, buffer, select

jaykapalczynski — Sun, 12 Dec 2021 08:39:51 GMT

Well here is the final result.....and I will try and explain what is going on.....

Right now its set up to grab a few parameters that are going to be passed to it once it is published as a Rest Endpoint service....I have those 3 parameters hardcoded for testing purposes (address, distance, and unique ID)

# PARAMETERS PASSED FROM CALL TO THE SERVICE
searchAddress = "86 Rockland Street West Davenport, NY 13860"
searchDistance = 20
searchid = "123"‍‍‍‍

The script takes the address and geocodes it to get XY

The XY are then used to create a temp Buffer at the user specified distance

this buffer is then used to select counties that INTERSECT it

this returned data is then packaged up to be send back to the requester...

Any questions please feel free to ask.....

THANK YOU ALL FOR YOUR HELP

OH one more thing....I was tripping up on the spatial reference....my Counties were in UTM and the XY being returned from the Geocode was Decimal Degrees....it could never find anything in the spatial select. Knowing that my data was in UTM I made sure that I set the output of the Geocode result to UTM as well

def singleAdressGeocode(address, geoCodeUrl, outSR = "26917"):‍

Last thing my data was also in meters so when the Miles were sent as parameters from the requester I had to convert those to meters to feed the buffer parameters

searchDistance = 20
searchid = "123"

# CONVERT MILE TO METERS FOR SPATIAL BUFFER
distanceMeters = searchDistance * 1609.34‍‍‍‍‍

FINAL CODE

import arcpy
import requests

# PARAMETERS PASSED FROM CALL TO THE SERVICE
searchAddress = "86 Rockland Street West Davenport, NY 13860"
searchDistance = 20
searchid = "123"

# GEOCODERS
# esri geocoding service
geoCodeUrl = "https://geocode.arcgis.com/arcgis/rest/services/World/GeocodeServer/findAddressCandidates"

# CONVERT MILE TO METERS FOR SPATIAL BUFFER
distanceMeters = searchDistance * 1609.34

# SEARCH PARAMETERS
print " - SEARCH PARAMETERS -"
print "search address: " + searchAddress
print "passed unique ID : " + searchid
print "search miles: " + str(searchDistance)
print "search meters: " + str(distanceMeters)
print ""
print "-------------------"

# LOCAL VARIABLES
coordinates = ""
countList = []   
countyField = "FIPS2"  
countList2 = []
finalList = []
txt_list = ""
txt_list2 = ""

def singleAdressGeocode(address, geoCodeUrl, outSR = "26917"):
  ##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
geocodeResult = singleAdressGeocode(searchAddress, geoCodeUrl)
addressResult = geocodeResult["candidates"][0]["address"]
coordinateX = geocodeResult["candidates"][0]["location"]["x"]
coordinateY = geocodeResult["candidates"][0]["location"]["y"]
strcoordinateX = str(coordinateX)
strcoordinateY = str(coordinateY)


print "gecode result: " + str(geocodeResult)
print "-------------------"
print "geocoded address: " + addressResult
print "-------------------"
print "X coordinate: " + str(coordinateX)
print "-------------------"
print "Y coordinate: " + str(coordinateY)
print "-------------------"

# CREATE THE POINT FROM THE X AND Y
point = arcpy.Point(coordinateX, coordinateY)
ptGeometry = arcpy.PointGeometry(point)

# FEED THE BUFFER ANALYSIS THE POINT GEOMETRY THE DISTANCE PARAMETERS
arcpy.Buffer_analysis(ptGeometry, 'IN_MEMORY/BufferGeom', distanceMeters)

arcpy.env.workspace = "C:\Work\CountySelection\Data\SelectCounties.gdb"
arcpy.MakeFeatureLayer_management('CountiesDissolved', 'counties_lyr')

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

# GET THE FIPS CODES
rows = arcpy.SearchCursor(results,"","","FIPS2")
for row in rows:  
 neighborVal = str(row.FIPS2)  
 #print(neighborVal)  
 countList.append(neighborVal)  
 txt_list = ','.join(countList)
print "FIPS codes: " + str(txt_list)
print "-------------------"

# GET THE JURIS NAME
rows2 = arcpy.SearchCursor(results,"","","FIRST_JURI")
for row in rows2:  
 neighborVal2 = str(row.FIRST_JURI)  
 #print(neighborVal)  
 countList2.append(neighborVal2)  
 txt_list2 = ','.join(countList2)
print "County Name: " + str(txt_list2)
print "-------------------"

# GET BOTH THE FIPS AND JURIS NAME
fields = ['FIPS2', 'FIRST_JURI']
with arcpy.da.SearchCursor(results, fields) as cursor:
    for row in cursor:
        finalList.append('[{0}-{1}]'.format(row[0], row[1]))
        #print('{0}, {1}'.format(row[0], row[1]))
print "[" + "[" + strcoordinateX + "]," + "[" + strcoordinateY + "]," + "[" + str(searchDistance) + "]," + "[" + searchid + "]," + "[" + txt_list + "]," + "[" + txt_list2 + "], [" + str(finalList) + "] ]"
print "-------------------"


# 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 were located within {1} miles of your location'.format(
                                                  matchcount, searchDistance))
print "-------------------"
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Address, buffer, select

jaykapalczynski — Sun, 12 Dec 2021 08:39:46 GMT

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))



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