topic Re: How to Query a Feature Layer with an Intersecting Point? in ArcGIS API for Python Questions

How to Query a Feature Layer with an Intersecting Point?

MikeBrown4 — Mon, 02 Apr 2018 15:41:55 GMT

I have a point, created with arcgis.geometry.Point, that I wish to use in a feature layer query to only return the records that intersect that point. I can do this via the REST API by setting the geometry parameters, but I am struggling to figure out how to do this with the Python API. I have read the reference guide and examples, but still not sure how to do this.

Any help would be much appreciated.

Thanks in advance.

Re: How to Query a Feature Layer with an Intersecting Point?

BalajiVeera — Mon, 02 Apr 2018 17:24:21 GMT

To do same in Python, create the PointGeometry first

point = arcpy.Point(x, y)

pointGeometry = arcpy.PointGeometry(point)

Apply the same PointGeometry in SelectLayerByLocation_management

arcpy.SelectLayerByLocation_management('query_feature_layer', 'CONTAINS', pointGeometry)

Re: How to Query a Feature Layer with an Intersecting Point?

MikeBrown4 — Mon, 02 Apr 2018 17:43:29 GMT

Thank you for the response.

I am trying to do this through the Python API, not with arcpy. Apologies if I did not make this clear in my question. I am wondering if there is a way to query the feature layer and include the point to only return back records that intersect. Something like:

query_result = feature_layer.query(where='1=1', some geometry parameter)

Using the REST API, I would construct a query like this to only get back the record that intersects the point: https://services6.arcgis.com/bxsVBCAOuftUsVxD/arcgis/rest/services/California_Cities/FeatureServer/0/query?where=1=1&geometryType=esriGeometryPoint&geometry=-104,35.6

Re: How to Query a Feature Layer with an Intersecting Point?

BalajiVeera — Mon, 02 Apr 2018 18:15:24 GMT

Sorry for confusion with arcpy instead Python API

As per Python API, we can use the parameter for geometry_filter

query(layer_defs_filter=None, geometry_filter=None, time_filter=None, return_geometry=True, return_ids_only=False, return_count_only=False, return_z=False, return_m=False, out_sr=None)

define geometry filter to use as parameter for feature_layer.query method

arcgis.geometry.filters.contains(point)

Re: How to Query a Feature Layer with an Intersecting Point?

JakeSkinner — Wed, 04 Aug 2021 17:17:05 GMT

Hi Mike,

Here's an example:

from arcgis.gis import GIS from arcgis.geometry import filters gis = GIS('https://www.arcgis.com', 'agol', 'gis123') ptFeatureLayer = gis.content.get('11ed3710e5a54d988d9e4f9d4068ca6f') polygonFeatureLayer = gis.content.get('48f9af87daa241c4b267c5931ad3b226') ptFeatures = ptFeatureLayer.layers[0] queryResult = ptFeatures.query(where='OBJECTID=10') ptFeat = queryResult.features[0].geometry polygonFeatures = polygonFeatureLayer.layers[0] queryResult = polygonFeatures.query(where='1=1', return_ids_only=True, geometry_filter=filters.intersects(ptFeat)) print(queryResult)‍‍‍‍‍‍‍‍‍‍‍‍‍‍

Re: How to Query a Feature Layer with an Intersecting Point?

MikeBrown4 — Mon, 02 Apr 2018 20:18:28 GMT

Thank you, this worked great. I used intersect to get the results I needed.

point = arcgis.geometry.Point({"x" : x, "y" : y, "spatialReference" : {"wkid" : 4326}})
gf = arcgis.geometry.filters.intersects(point)
test = cities_layer.query(where='1=1',geometry_filter=gf)