arcrest geometryfilter?

3399
2
Jump to solution
05-30-2016 07:09 PM
RichardWebb2
New Contributor II

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)

0 Kudos
1 Solution

Accepted Solutions
AndrewChapkowski
Esri Regular Contributor

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) )

View solution in original post

2 Replies
JoshuaBixby
MVP Esteemed Contributor

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.

0 Kudos
AndrewChapkowski
Esri Regular Contributor

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) )