Feature Layer query for multiple coordinates

Question

Hi,

I have a FeatureLayer that I'm trying to query for a bunch of coordinates. Is there any way to query the FeatureLayer once (from the WWW to my local machine) and then query it once per coordinate, instead of doing an HTTP request per coordinate?

Thanks,

Gilles

sigccil · Answer

Did you ever get follow up on this issue?

OSINT_ENGINEER · Answer

Whenever we try to use Multipoint directly against the query operation of a hosted feature service endpoint, it seems the feature service always returns the first matching feature, not all matching features. Would you be so kind and validate that maybe our Multipoint geometry is malformed?

AnttiKajanus · Answer

I totally get that you want to do all in one query to avoid the ridiculous amount of HTTP calls when doing these kind of queries one by one. I think that what you are looking here is to use a correct spatial relationship. Feature Service supports following spatial relationships on the query * esriSpatialRelIntersects —Part of an input feature is contained in the query geometry * esriSpatialRelContains —Part or all of an input feature is contained within the query geometry * esriSpatialRelCrosses —An input feature crosses the query geometry * esriSpatialRelEnvelopeIntersects —The envelope of an input feature intersects the envelope of the query geometry * esriSpatialRelOverlaps —An input feature overlaps the query geometry * esriSpatialRelTouches —An input feature touches the border of the query geometry * esriSpatialRelWithin—An input feature is completely enclosed by the query geometry So instead of using contains you should use intersects instead. That should/could/maybe do the trick. from arcgis.gis import GIS from arcgis.features import FeatureLayerCollection from arcgis.geometry import * feature_service_url = 'https://gis.wri.org/server/rest/services/LandMark/comm_ind_Documented/MapServer' feature_service_root = FeatureLayerCollection(feature_service_url) feature_layer = feature_service_root.layers[1] feature_layer # multipoint_json = {"points":[ [-109.271351,36.357824], [127.484020,-14.681117] ], "spatialReference":{"wkid":4326 }}; multipoint = MultiPoint(multipoint_json) print(multipoint.is_valid()) # True query_result = feature_layer.query( geometry_filter = filters.intersects(multipoint), return_geometry='false', out_fields = 'NAME') print("Found {} results.".format(len(query_result.features))) for feature in query_result.features: print(feature) # Found 3 results. # {"attributes": {"OBJECTID": 30989, "Name": "Navajo Nation Reservation"}} # {"attributes": {"OBJECTID": 106183, "Name": "FORREST RIVER"}} # {"attributes": {"OBJECTID": 107290, "Name": "Balanggarra (Combined)"}}‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍‍