topic Re: arcrest geometryfilter? in Python Questions

arcrest geometryfilter?

RichardWebb2 — Sat, 11 Dec 2021 16:31:06 GMT

Is there a better way of getting that geometry into the second query? Not happy with line 28

from __future__ import print_function
from arcrest.ags import FeatureLayer
import arcpy
import json

if __name__ == "__main__":

    arcpy.env.overwriteOutput = True

    lgaUrl = "http://services.thelist.tas.gov.au/arcgis/rest/services/Public/CadastreAndAdministrative/MapServer/4"

    flLGA = FeatureLayer(
        url=lgaUrl)

    lgaWhereClause = "NAME = 'Launceston'"
    outSR = '28355'
##    outFC="C:\Scripts\data.gdb\LGA"
##    # get lga geometry
##    lgaGeom = flLGA.query(where=lgaWhereClause,out_fields='shape',returnGeometry=True,returnFeatureClass=True,out_fc=outFC,outSR=outSR)
    lgaGeom = flLGA.query(where=lgaWhereClause,out_fields='shape',returnGeometry=True,outSR=outSR)
    geomText = lgaGeom.toJSON
    geomRings = json.loads(geomText)

    for key, value in geomRings.iteritems():
        if key == 'features':
            features = str(value)

    features = features.replace("[{u'geometry': ","").replace(", u'attributes': {}}]","").replace("{u'rings':","{'rings':")

    # change to Crown Licenses
    crownUrl = "http://services.thelist.tas.gov.au/arcgis/rest/services/Public/CadastreAndAdministrative/MapServer/3"

    flCrown = FeatureLayer(
    url=crownUrl)

    crownWhereClause = "1=1"
    outFC="C:\Scripts\data.gdb\CrownLicenses"
    geometryType = "esriGeometryPolygon"

    # get crownLicenses
    crownLicenses =

38.

flCrown.query(where=crownWhereClause,geometry=features,out_fields='objectid,pid,shape',returnGeometry=Tr

e,returnFeatureClass=True,out_fc=outFC,geometryType=geometryType,inSR=outSR,outSR=outSR)

Re: arcrest geometryfilter?

JoshuaBixby — Tue, 31 May 2016 14:31:48 GMT

Have you inquired about this at the ArcREST GitHub site? If so, what was the response? If not, it might be worth asking there as well while you wait for feedback here.

Re: arcrest geometryfilter?

AndrewChapkowski — Sat, 11 Dec 2021 16:31:09 GMT

Try this:

import arcrest
if __name__ == "__main__":
  fl_url = "https://site.com/arcgis/rest/services/reference/adminboundaries/FeatureServer/0"
  token_url = "https://site.com/arcgis/admin/generateToken"
  org_url="https://site.com/arcgis"
  username = "admin"
  pw = "some password"
  where = "ISO_CC = 'US'"
  sh = arcrest.AGSTokenSecurityHandler(username=username,
  password=pw,
  org_url=org_url,
  token_url=token_url)
  fl = arcrest.ags.FeatureLayer(url=fl_url,securityHandler=sh)
  fs = fl.query(where=where)
  print ("Found %s features with the where clause: %s" % (len(fs), where) )
  if len(fs.features) > 0:
     feature = fs.features[0]
     geom = feature.geometry
     geomFilter = arcrest.filters.GeometryFilter(geomObject=geom)
     result = fl.query(where=where, geomtryFilter=geomFilter, as_json=True) # set as_json=False for the featureset object
     print ("Found %s features with the where clause: %s and geometry filter" % (len(result['features']), where) )