I am attempting to access ArcGIS Map Services and query a dataset to only return records that intersect the boundary of a custom geometry I created. I have queried and downloaded data before using this script, but never using geometry. Do i simply add my geometry_json to the Where clause?
#My geometry (a simply polygon)geometry_json = {"geometryType" : "esriGeometryPolygon","spatialReference" : {"wkid" : 4326,"latestWkid" : 4326 },"features" : [{"geometry":{"curveRings":[ [[-123.77335021899995,49.353870065000081],{"a":[[-123.77335021899995,49.353870065000081],[-123.88831213030559,49.338864996255367],0,1]}]]}}]}The script to access and download the data
import arcpy # Import the arcpy module, which provides access to ArcGIS functionality.import urllib.request # Import the urllib.request module to work with URLs.import json # Import the json module for JSON data handling.# Setuparcpy.env.overwriteOutput = True baseURL = r"https://gisp.dfo-mpo.gc.ca/arcgis/rest/services/FGP/MSDI_Dynamic_Current_Layer/MapServer/0" # Base URLfields = "*" # Return all fieldsoutdata = r"out/path" #output location# Get record extract limiturlstring = baseURL +"?f=json" # Construct the URL to request JSON metadata from the service.j = urllib.request.urlopen(urlstring) # Open the URL and get the response.js = json.load(j) # Load the JSON response into a Python dictionary.maxrc = int(js["maxRecordCount"]) # Extract the maximum record count from the JSON response.print ("Record extract limit: %s" % maxrc) # Print the maximum record count.# Get object ids of featureswhere = "1=1" # Define a where clause to retrieve all records.urlstring = baseURL +"/query?where={}&returnIdsOnly=true&f=json".format(where) # Construct the URL to request object IDs.j = urllib.request.urlopen(urlstring) # Open the URL and get the response.js = json.load(j) # Load the JSON response into a Python dictionary.idfield = js["objectIdFieldName"] # Extract the name of the object ID field.idlist = js["objectIds"] # Extract the list of object IDs.idlist.sort() # Sort the list of object IDs.numrec = len(idlist) # Get the number of records.print ("Number of target records: %s" % numrec) # Print the number of records.# Gather featuresprint ("Gathering records...") # Print a message indicating feature gathering.fs = dict() # Create an empty dictionary to store features.for i in range(0, numrec, maxrc): # Loop over the object IDs in chunks based on the maximum record count. torec = i + (maxrc - 1) # Calculate the end index of the chunk. if torec > numrec: # Adjust the end index if it exceeds the number of records. torec = numrec - 1 fromid = idlist[i] # Get the starting object ID of the chunk. toid = idlist[torec] # Get the ending object ID of the chunk. where = "{} >= {} and {} <= {}".format(idfield, fromid, idfield, toid) # Define a where clause for the chunk. print (" {}".format(where)) # Print the where clause. urlstring = baseURL +"/query?where={}&returnGeometry=true&outFields={}&f=json".format(where,fields) # Construct the URL to request features. fs[i] = arcpy.FeatureSet() # Create a new empty FeatureSet object. fs[i].load(urlstring) # Load features from the URL into the FeatureSet.# Save featuresprint ("Saving features...") # Print a message indicating feature saving.fslist = [] # Create an empty list to store FeatureSets.for key,value in fs.items(): # Iterate over the dictionary of FeatureSets. fslist.append(value) # Append each FeatureSet to the list.arcpy.Merge_management(fslist, outdata) # Merge all FeatureSets into a single feature class at the specified output location.print ("Done!") # Print a message indicating completion.