by
by
by
by
by
RandyBurton
Another approach...
The Spatial Join (Analysis) will compute distances in certain cases. As a test, a Virginia layer (containing only one feature - an outline of the state in UTM Zone 17N) was opened in Desktop 10.5. I geocoded the two addresses previously used with some others that lie in and outside the state. I pasted code into the Python window, ran the spatial join and had the following results:
The DistanceTo field was generated by the spatial join. A distance of 0 indicates the point is inside the state. If it lies outside, a distance in meters is calculated. This distance can be compared to the buffer in the Distance field. If the DistanceTo is less than the Distance buffer, the point is within the buffer. A select layer by attribute can use this in a where clause.
Here's the test code, which may need modification depending on your settings or if run outside ArcMap. In_memory was used in creating the points and spatialJoin layers, so they should be deleted at some point.
import arcpy, requests, json
# ['searchAddress', searchDistance (in miles)],
addressList = [['1520 Split Oak Ln, Henrico, VA, 23229', 25],
['6188 Maddox Blvd, Chincoteague, VA, 23336', 12],
['318 Market St, Peterstown, WV 24963', 1],
['109 Federal St, Rich Creek, VA 24147',1],
['439 Rockford St, Mt Airy, NC 27030',2] ]
KMile = 1609.34
def singleAdressGeocode(address):
geoCodeUrl = r"https://gismaps.vita.virginia.gov/arcgis/rest/services/Geocoding/VGIN_Composite_Locator/GeocodeServer/findAddressCandidates"
outSR = "26917"
r = requests.get(geoCodeUrl, params={'SingleLine': address, 'outSR' : outSR, 'maxLocations': 1, 'f': 'pjson'}) # requests will url encode params
data = r.json()
if data["candidates"]:
return data
else:
return None # indicates error
def insertResult(address, distance, result, cursor):
geoAddress = geocodeResult["candidates"][0]["address"]
score = geocodeResult["candidates"][0]["score"]
pointX = geocodeResult["candidates"][0]["location"]["x"]
pointY = geocodeResult["candidates"][0]["location"]["y"]
ptGeometry = arcpy.PointGeometry(arcpy.Point(
geocodeResult["candidates"][0]["location"]["x"],
geocodeResult["candidates"][0]["location"]["y"]),
arcpy.SpatialReference(geocodeResult["spatialReference"]["wkid"]))
cursor.insertRow([ptGeometry, address, geoAddress, score, pointX, pointY, distance])
# create in_memory feature to save address points
fc = arcpy.CreateFeatureclass_management("in_memory","points","POINT", spatial_reference=arcpy.SpatialReference(26917))
arcpy.AddField_management(fc,'Address','TEXT',"#","#",100,'Address',"NULLABLE","NON_REQUIRED")
arcpy.AddField_management(fc,'geoAddress','TEXT',"#","#",100,'Geocoded Address',"NULLABLE","NON_REQUIRED")
arcpy.AddField_management(fc,'geoScore','Double',"#","#","#",'Geocode Score',"NULLABLE","NON_REQUIRED")
arcpy.AddField_management(fc,'POINT_X','Double',"#","#","#",'Point X',"NULLABLE","NON_REQUIRED")
arcpy.AddField_management(fc,'POINT_Y','Double',"#","#","#",'Point Y',"NULLABLE","NON_REQUIRED")
arcpy.AddField_management(fc,'Distance','Double',"#","#","#",'Distance',"NULLABLE","NON_REQUIRED")
# start an insert cursor
cursor = arcpy.da.InsertCursor(fc,["SHAPE@", 'Address', 'geoAddress', 'geoScore', 'POINT_X', 'POINT_Y', 'Distance'])
# loop through addresses and insert in points feature
for address, distance in addressList:
geocodeResult = singleAdressGeocode(address)
if geocodeResult:
insertResult(address, int(distance)*KMile, geocodeResult, cursor)
else:
print('ERROR: Address ({}) did not geocode'.format(address))
del cursor # done with insert cursor
# a spatial join can be used to compute distances
# Virginia layer contains only 1 feature, an outline of VA using NAD_1983_UTM_Zone_17N
# field_mapping is optional, only fields in points feature were used
arcpy.SpatialJoin_analysis(target_features="points",
join_features="Virginia",
out_feature_class="in_memory/spatialJoin",
join_operation="JOIN_ONE_TO_ONE",
join_type="KEEP_ALL",
field_mapping='''
Address "Address" true true false 100 Text 0 0 ,First,#,points,Address,-1,-1;
geoAddress "Geocoded Address" true true false 100 Text 0 0 ,First,#,points,geoAddress,-1,-1;
geoScore "Geocode Score" true true false 0 Double 0 0 ,First,#,points,geoScore,-1,-1;
POINT_X "Point X" true true false 0 Double 0 0 ,First,#,points,POINT_X,-1,-1;
POINT_Y "Point Y" true true false 0 Double 0 0 ,First,#,points,POINT_Y,-1,-1;
Distance "Distance" true true false 0 Double 0 0 ,First,#,points,Distance,-1,-1''',
match_option="CLOSEST_GEODESIC",
distance_field_name="DistanceTo")
# a search cursor can print a results table
fields = ['Address', 'geoAddress', 'geoScore', 'POINT_X', 'POINT_Y', 'Distance', 'DistanceTo']
print('\t'.join(fields)+'\tLocation')
with arcpy.da.SearchCursor("spatialJoin", fields) as cursor:
for a, g, s, x, y, d, dt in cursor:
if dt == 0:
inout = 'Inside Virginia'
elif dt < d:
inout = 'Within Buffer'
else:
inout = 'Outside Buffer'
print('{}\t{}\t{}\t{}\t{}\t{}\t{}\t{}'.format(a, g, s, x, y, d, dt, inout))
# delete the in_memory feature when done
# arcpy.Delete_management('in_memory/points')
# arcpy.Delete_management('in_memory/spatialJoin')
''' # address points can be selected where the distance is greater than the buffer
arcpy.SelectLayerByAttribute_management(in_layer_or_view="spatialJoin",
selection_type="NEW_SELECTION",
where_clause="DistanceTo > Distance")
'''Hope this helps.
