# Import the required ArcGIS API for Python modules import arcgis from arcgis.gis import GIS from arcgis.geoanalytics import manage_data
# Connect to your ArcGIS Enterprise portal and check that GeoAnalytics is supported portal = GIS("https://myportal.domain.com/portal", "gis_publisher", "my_password", verify_cert=False) if not portal.geoanalytics.is_supported():
    print("Quitting, GeoAnalytics is not supported")    exit(1)   
# Find the big data file share dataset you're interested in using for analysis search_result = portal.content.search("", "Big Data File Share")
# Look through search results for a big data file share with the matching name bd_file = next(x for x in search_result if x.title == "bigDataFileShares_LandUse")
# Look through the big data file share for grazing land grazing_land = next(x for x in bd_file.layers if x.properties.name == "GrazingLand")
# Find a feature layer named "Watersheds" in your ArcGIS Enterprise portal and apply a filter watersheds = portal.content.search("Watersheds", "Feature Layer") watersheds_layer = layer_result[0].layers[0] watersheds_layer.filter = "region = 'Missouri'"
# Set the tool environment settings arcgis.env.verbose = True arcgis.env.defaultAggregations = True
# Run the Overlay Layers tool overlay_result = manage_data.overlay_data(input_layer = grazing_land,                                           overlay_layer = watersheds_layer, 
                                          overlay_type = "Intersect",                                           output_name = "Watershed_intersections")
# Visualize the tool results if you are running Python in a Jupyter Notebook processed_map = portal.map('Missouri', 6) processed_map.add_layer(overlay_result)