topic Re: Test for XY inside Polygon in Python Questions

Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:53:58 GMT

Im trying to do this and stuck

Take a JSON string with address and distance set those to variables
Geocode the Address to get an XY and set those to Variables
Buffer a polygon of the state using the distance variable created from the DataInput JSON
Define the State Boundary
Buffer the State Boundary and set into memory
Test ot see if the point is within that Buffered State Boundary
Return True if its within and False if Outside the Buffered State Boundary

ERROR:

Traceback (most recent call last):
File "E:\ArcGIS Server\ArcGIS\xxxxxxxx.py", line 87, in <module>
results = arcpy.SelectLayerByLocation_management(ptGeometry, 'INTERSECT', 'IN_MEMORY/BufferGeom')
File "C:\Program Files (x86)\ArcGIS\Desktop10.5\ArcPy\arcpy\management.py", line 7891, in SelectLayerByLocation
raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000368: Invalid input data.
Failed to execute (SelectLayerByLocation).

Can Anyone see anything obvious as to why I am tripping this error?

import arcpy
import requests
import json
import uuid

#dataInputParameter1 = arcpy.GetParameterAsText(0)
dataInputParameter1 = '{"Location":[{"address":"6188 Maddox Blvd, Chincoteague, VA, 23336","distance":"12" }]}'

        
# ADD ADDITIONAL VARIABLE FOR RETURN RESULT  
varResult = "true"
    
textDict = json.loads(dataInputParameter1)

langs = []
       
for SpecificLocation in textDict['Location']:
    varsearchAddress =  SpecificLocation["address"]
    varsearchDistance =  SpecificLocation["distance"]

    print ("address: " + varsearchAddress)
    print ("distance: " + varsearchDistance)
    print ""

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

    # GEOCODERS
    geoCodeUrl = r"https://gismaps.vita.virginia.gov/arcgis/rest/services/Geocoding/VGIN_Composite_Locator/GeocodeServer/findAddressCandidates"

    searchAddress = varsearchAddress
    searchDistance = varsearchDistance

    # CONVERT MILE TO METERS FOR SPATIAL BUFFER
    distanceMeters = int(varsearchDistance) * 1609.34

    def singleAdressGeocode(address, geoCodeUrl, outSR = "4326"):
    #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)


    arcpy.env.workspace = r"E:/ArcGISProjects/CountySelection/xxxxx.sde"

    # CREATE THE POINT FROM THE X AND Y
    point = arcpy.Point(coordinateX, coordinateY)
    ptGeometry = arcpy.PointGeometry(point)
    
    # DEFINE THE STATE FEATURE CLASS
    arcpy.MakeFeatureLayer_management(r"E:/ArcGISProjects/CountySelection/xxxxx.sde/Boundary", "state_lyr")

    # BUFFER THE STATE LAYER AND PUT IT IN MEMORY
    arcpy.Buffer_analysis('state_lyr', 'IN_MEMORY/BufferGeom', distanceMeters)

    # TEST TO SEE IF POINT IS INSIDE THE BUFFERED GEOMETRY
    results = arcpy.SelectLayerByLocation_management(ptGeometry, 'INTERSECT', 'IN_MEMORY/BufferGeom')

    #Test if there are any polygons points:  
    if int(arcpy.GetCount_management('IN_MEMORY/BufferGeom').getOutput(0))==0:  
        raise Exception("Selected points do not overlap any polygons.")  
        print "Selected points do not overlap any polygons."
        varResult = "true"‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

DanPatterson — Tue, 01 Sep 2020 22:14:19 GMT

Select Layer By Location (Data Management)—ArcGIS Pro | Documentation

you need to make a layer for your pntGeometry

Re: Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:54:01 GMT

OK that removed the error and an oversight by me.....thanks for the eyes....

For some reason I am getting true for any address....I think my code to test if its in or not is not working correctly.

If I put in an address for Wyoming it returns true....I set the variable to false at the beginning but for some reason I get a count from the results variable either way....

Any ideas where this test for the result is wrong?

    # TEST TO SEE IF POINT IS INSIDE THE BUFFERED GEOMETRY
    results = arcpy.SelectLayerByLocation_management("point_lyr", 'INTERSECT', 'IN_MEMORY/BufferGeom')

    #Test if there are any polygons points:  
    if int(arcpy.GetCount_management(results).getOutput(0))==0:  
        #raise Exception("Selected points do not overlap any polygons.")  
        print "Selected points do not overlap any polygons."
        varResult = "true"‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

jaykapalczynski — Wed, 02 Sep 2020 10:56:56 GMT

Woke up at 3am and it dawned on my that I might be having issues with the projection of the Feature Class and the Point created from the Geocode...think once I fix that it might work....testing that this morning.

Re: Test for XY inside Polygon

jaykapalczynski — Wed, 02 Sep 2020 12:00:15 GMT

NAH that did not work....but at least the polygon state layer is WGS84 and the point being geocoded is output in WGS 84 so thats on the same page but I never return any features.....

Re: Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:54:04 GMT

dan patterson‌ OK I got it working but it seems the issue is around the point feature that I am creating from the Geocoded results...

Any thoughts? The data is in WGS84

If I point this to a point Feature Class as such this WORKS

arcpy.MakeFeatureLayer_management(r"E:/ArcGISProjects/CountySelection/xx.sde/Test", "testPoint")

results = arcpy.SelectLayerByLocation_management('testPoint', 'WITHIN', 'IN_MEMORY/BufferGeom') ‍‍‍‍‍‍‍‍‍

But I need to create a Point reference from the Geocoded XY this FAILS

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

# DEFINE THE POINT FEATURE CLASS
arcpy.MakeFeatureLayer_management(ptGeometry, 'point_lyr')

results = arcpy.SelectLayerByLocation_management('point_lyr', 'WITHIN', 'IN_MEMORY/BufferGeom') ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

RandyBurton — Sun, 12 Dec 2021 00:54:06 GMT

I assume BufferGeom is the polygon layer, so perhaps something like:

x = -75.369027 # longitude
y = 37.935822 # latitude

# WGS 1984: 4326
# GCS North American 1983: 4269 ## use the factory/authority code for BufferGeom layer
ptGeometry = arcpy.PointGeometry(arcpy.Point(x,y),arcpy.SpatialReference(4326)).projectAs(arcpy.SpatialReference(4269))

results = arcpy.SelectLayerByLocation_management('IN_MEMORY/BufferGeom', 'CONTAINS', ptGeometry, selection_type="NEW_SELECTION")‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:54:12 GMT

Think I might have it....Using the MakeFeatureLayer to take the ptGeometry and make a feature layer....

ptGeometry = arcpy.PointGeometry(arcpy.Point(coordinateX,coordinateY),arcpy.SpatialReference(4326)).projectAs(arcpy.SpatialReference(4269))

arcpy.MakeFeatureLayer_management(ptGeometry, "PointfromGeocode")‍‍‍

Also changed the Coord system to match that of my State Layer (webaux sphere)

ptGeometry = arcpy.PointGeometry(arcpy.Point(x,y),arcpy.SpatialReference(4326)).projectAs(arcpy.SpatialReference(3857))‍

For some reason not finding the XY coordinates keep telling me no matches...

    matchcount = int(arcpy.GetCount_management(results)[0]) 

    if matchcount == 0:
        print('no features matched spatial and attribute criteria')
    else:
        print('found something')
        print ""
        varResult = "true"‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:54:15 GMT

OK ... Puzzled here....I can get this to work if I reference the point that is in a Feature Class but

I cannot get it to work with the result from my Geocoded address and resulting XY

My original address that is geocoded is WGS84 I convert that to UTM Nad83 Zone 17N (26917)

The state State layer is stored in UTM Nad83 Zone17N (26917) Meters

The input from the parameters is Miles and I am converting that distance to Meters.

I know the distance thing is working because I can test that against a Point that is outside the state in a Feature Class and increase that distance and it eventually finds it...

Any thoughts as to why the Geocoded result is not being found...I am converting this to the Appropriate Spatial Reference.....ugggggg

This Does NOT work

    ## CONVERT MILE TO METERS FOR SPATIAL BUFFER
    # Value received from input parameter in JSON string
    distanceMeters = int(varsearchDistance) * 1609.34
   arcpy.MakeFeatureLayer_management(r"E:/ArcGISProjects/Selection/xx.sde/Virginia_UTM_Meters", "state_lyr")

    x = -77.434769 # longitude
    y = 37.541290 # latitude

    # convert from WGS84 to UTM zone 17N
    ptGeometry = arcpy.PointGeometry(arcpy.Point(x,y),arcpy.SpatialReference(4326)).projectAs(arcpy.SpatialReference(26917))

    arcpy.MakeFeatureLayer_management(ptGeometry, "PointfromGeocode")

    distance = int(distanceMeters)
       
    results = arcpy.SelectLayerByLocation_management('PointfromGeocode', 'WITHIN_A_DISTANCE', 'state_lyr', distance, selection_type='NEW_SELECTION')  ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

This works - only real difference is that I am pointing to a Feature Class point instead of usign the XY to create the point geometry

    ## CONVERT MILE TO METERS FOR SPATIAL BUFFER
    # Value received from input parameter in JSON string
    distanceMeters = int(varsearchDistance) * 1609.34
   arcpy.MakeFeatureLayer_management(r"E:/ArcGISProjects/Selection/xx.sde/Virginia_UTM_Meters", "state_lyr")

arcpy.MakeFeatureLayer_management(r"E:/ArcGISProjects/Selection/xx.sde/TestPoint", "testPoint")

    distance = int(distanceMeters)
       
    results = arcpy.SelectLayerByLocation_management('testPoint', 'WITHIN_A_DISTANCE', 'state_lyr', distance, selection_type='NEW_SELECTION')  ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

HOW CAN I Test this new Layer? I tried the below but I dont get anything returned.

    ptGeometry = arcpy.PointGeometry(arcpy.Point(x,y),arcpy.SpatialReference(4326)).projectAs(arcpy.SpatialReference(26917))
    arcpy.MakeFeatureLayer_management(ptGeometry, "PointfromGeocode")

    print "start"

    fields = ['SHAPE@XY']

    with arcpy.da.SearchCursor('PointfromGeocode', fields) as cursor:
        for row in cursor:
            print(u'{0}'.format(row[0]))

    print "end of first"
‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

RandyBurton — Sun, 12 Dec 2021 00:54:18 GMT

Regarding testing the new layer:

I was unable to turn ptGeometry into a layer as you are attempting to do in line 2 of your last code segment. However, I was able to use ptGeometry as the select_features parameter in Select Layer By Location. To test the values being used for x and y, you can print the JSON containing the x and y values:

# x,y for 6188 Maddox Blvd, Chincoteague, VA, 23336
x = -77.434769 # longitude
y = 37.541290 # latitude

ptGeometry = arcpy.PointGeometry(arcpy.Point(x,y),arcpy.SpatialReference(4326)).projectAs(arcpy.SpatialReference(26917))

print(ptGeometry.JSON)
''' Output:
{"x":815010.20711794402,"y":4160898.4941948522,"spatialReference":{"wkid":26917,"latestWkid":26917}}
''' ‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

I am a bit confused with your state layer, buffer and distance. Does the layer have only one feature, the state of Virginia? Or does it contain other states or Virginia counties. Is the buffer being used to see if the address is close to the border of Virginia, possibly in a nearby state?

Re: Test for XY inside Polygon

DanPatterson — Sun, 12 Dec 2021 00:54:23 GMT

Why not make things easier and get the geometry objects for both. Then you don't have to make featurelayers etc etc

As an example, a point geometry being tested against a polygon

You can project geometry if needed. To find out more, just do a 'dir' on the objects

dir(pnt)
['JSON', 'WKB', 'WKT', ... snip...
, '_fromGeoJson', '_go', '_passthrough', '_repr_svg_', 'angleAndDistanceTo', 'area',
'boundary', 'buffer', 'centroid', 'clip', 'contains', 'convexHull', 'crosses', 'cut', 
'densify', 'difference', 'disjoint', 'distanceTo', 'equals', 'extent', 'firstPoint', 
'generalize', 'getArea', 'getGeohash', 'getLength', 'getPart', 'hasCurves', 
'hullRectangle', 'intersect', 'isMultipart', 'labelPoint', 'lastPoint', 'length', 
'length3D', 'measureOnLine', 'overlaps', 'partCount', 'pointCount', 
'pointFromAngleAndDistance', 'positionAlongLine', 'projectAs', 'queryPointAndDistance',
 'segmentAlongLine', 'snapToLine', 'spatialReference', 'symmetricDifference', 
'toCoordString', 'touches', 'trueCentroid', 'type', 'union', 'within']


# ---- ditto
dir(poly)
['JSON', 'WKB', 'WKT', '__add__',... snip ..., 'within']

# ---- pnt.within(poly)  # ---- point in polygon checks
‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

RandyBurton — Sun, 12 Dec 2021 00:54:26 GMT

The ptGeometry item should have a spatial reference; it should not be null (Line 7 in the first code block in your previous message). I think the first parameter in SelectLayerByLocation should be a feature layer, and ptGeometry (which is only geometry for a point) is not a feature. I have been experimenting with some code in the Python window of Desktop 10.5 that may help you get the geocoding address coordinates into point geometry and then into a feature. The test code did create a point feature and it appears to be correctly located using the Streets basemap (and Google Maps). Here's my test code:

# create the point geometry

address: 1520 Split Oak Ln, Henrico, VA, 23229
distance: 25

>>> geocodeResult
{u'candidates': [{u'attributes': {}, u'score': 100, u'location': {u'y': 4167142.8740720972, u'x': 805193.5919402852}, u'address': u'1520 SPLIT OAK LN, HENRICO, VA, 23229'}], u'spatialReference': {u'wkid': 26917, u'latestWkid': 26917}}
>>> addressResult = geocodeResult["candidates"][0]["address"]
>>> addressResult
u'1520 SPLIT OAK LN, HENRICO, VA, 23229'
>>> ptGeometry = arcpy.PointGeometry(arcpy.Point(
...     geocodeResult["candidates"][0]["location"]["x"],
...     geocodeResult["candidates"][0]["location"]["y"]),
...     arcpy.SpatialReference(geocodeResult["spatialReference"]["wkid"]))
...     
>>> ptGeometry.JSON
u'{"x":805193.59194028517,"y":4167142.8740720972,"spatialReference":{"wkid":26917,"latestWkid":26917}}'
>>>

# create an in_memory feature class to hold the point

>>> fc = arcpy.CreateFeatureclass_management("in_memory","points","POINT",
...                                          spatial_reference=arcpy.SpatialReference(geocodeResult["spatialReference"]["wkid"]))
... 
>>> fc
<Result 'in_memory\\points'>
>>> cursor = arcpy.da.InsertCursor(fc,["SHAPE@"])
>>> cursor.insertRow([ptGeometry])
1L
>>> del cursor

# delete the in_memory feature when done

>>> arcpy.Delete_management('in_memory/points')
<Result 'true'>‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

RandyBurton — Sun, 12 Dec 2021 00:54:29 GMT

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.

Re: Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:54:32 GMT

OK I am here again....still having issues with the first parameter ....Im trying to create the Feature and put it into memory as you suggested....but getting the same error where the parameters are incorrect.

I am not sure about this line in your example

fcPoint<Result 'in_memory\\points'>

ERROR:

Traceback (most recent call last):
File "E:\ArcGIS Server\ArcGIS Published MXDs\Projects\NWTL\Python\NWTL_In_Virginia.py", line 132, in <module>
results = arcpy.SelectLayerByLocation_management('IN_MEMORY\\points', 'WITHIN_A_DISTANCE', 'state_lyr', distance, selection_type='NEW_SELECTION')
File "C:\Program Files (x86)\ArcGIS\Desktop10.5\ArcPy\arcpy\management.py", line 7891, in SelectLayerByLocation
raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000368: Invalid input data.
Failed to execute (SelectLayerByLocation).

# USE THE GEOCODED LOCATION
    ptGeometry = arcpy.PointGeometry(arcpy.Point(\
        geocodeResult["candidates"][0]["location"]["x"],\
        geocodeResult["candidates"][0]["location"]["y"]),\
        arcpy.SpatialReference(geocodeResult["spatialReference"]["wkid"]))

    print(ptGeometry.JSON)
    print ""

    out_path = "IN_MEMORY"
    out_name = "points"
    geometry_type = "POINT"
    template = ""
    has_m = ""
    has_z = ""

    fcPoint = arcpy.CreateFeatureclass_management(out_path,out_name,geometry_type,template,has_m,has_z,\
        spatial_reference=arcpy.SpatialReference(geocodeResult["spatialReference"]["wkid"]))

    #fcPoint<Result 'in_memory\\points'>
    cursor = arcpy.da.InsertCursor(fcPoint,["SHAPE@"])
    cursor.insertRow([ptGeometry])
    del cursor
            
## TEST TO SEE IF POINT IS INSIDE THE BUFFERED GEOMETRY
    results = arcpy.SelectLayerByLocation_management('IN_MEMORY/points', 'WITHIN_A_DISTANCE', 'state_lyr', distance, selection_type='NEW_SELECTION')    

# If features matched criteria, write them to a new feature class
    matchcount = int(arcpy.GetCount_management(results)[0]) 

    if matchcount == 0:
        print('no features matched spatial and attribute criteria')
    else:
        print('found something')
        print ""
        varResult = "true"

# DELETE THE FEATURE WHEN DONE
    arcpy.Delete_management('IN_MEMORY/points')‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

jaykapalczynski — Tue, 08 Sep 2020 13:28:22 GMT

Think I might have it.....but not sure still testing...took the fcPoint we created and used MakeFeatureLayer and used that

arcpy.MakeFeatureLayer_management(fcPoint, "PointfromGeocode")

results = arcpy.SelectLayerByLocation_management('PointfromGeocode', 'WITHIN', 'state_lyr')

Re: Test for XY inside Polygon

jaykapalczynski — Tue, 08 Sep 2020 19:04:40 GMT

Dan

I am interested in this approach...is there a better example? Trying to get a better understanding of this.

Re: Test for XY inside Polygon

BlakeTerhune — Wed, 09 Sep 2020 15:31:21 GMT

I think he's referring to the within() method for Geometry objects.

https://community.esri.com/thread/197968-re-how-to-use-arcpy-geometry-measureonline

Re: Test for XY inside Polygon

DanPatterson — Wed, 09 Sep 2020 16:21:19 GMT

As Blake suggests, since you have a mix of geometry origins (json, arcpy), and your query is essentially simple, it makes less sense now to be creating a bunch of featureclasses and/or featurelayers. Just get the geometry at its root (that is, arcpy geometry objects) and do the check. You basically need a boolean... if True, then do something, if False, move on

Re: Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:54:35 GMT

OK I think I am understanding now....so I would just read in my geometry....

Do I have this correct....

read in my Polygons coordinates instead of the squares XYs
create the array from list of coords
create my point with my XY coordinates
Test and point.within(polygon) will return true if within and false if outside

# READ IN THE COORDINATES OF MY POLY INSTEAD OF THE SQUARE BELOW
    coords = [[0, 0], [0, 2], [2, 2], [2, 0], [0, 0]]

# CREATE ARRAY OF THOSE COORDS
    polygon = arcpy.Polygon(arcpy.Array([arcpy.Point(*a) for a in coords]), sr)
    print "polygon area:", polygon.area

# READ IN MY POINT GEOMETRY
    point = arcpy.PointGeometry(arcpy.Point(X Coord, Y Coord), sr)

# TEST IF INSIDE WILL RETURN TRUE IF YES AND FALSE IF NO?
    print "point within polygon:", point.within(polygon)‍‍‍‍‍‍‍‍‍‍‍‍

Re: Test for XY inside Polygon

jaykapalczynski — Sun, 12 Dec 2021 00:54:38 GMT

Last questions.....Can you do?

distance = 25 + " Mile"
point.with-distance-of(polygon, distance)‍‍‍‍