Latest Contributions by JohannesLindner

It sounds like you only used the buffers as input fot the tool. If you use the buffers and the parcels as input, the tool will return a feature class that has each intersection between buffers and parcels. If Buffer_1 intersects the whole Parcel_1 and Buffer_2 intersects only a corner of Parcel_1, there will be two entries in the result feature class: whole geometry of Parcel_1 and the attributes of Parcel_1 & Buffer_1 only the corner of Parcel_1 with the attributes of Parcel_1 & Buffer_2 You can then use the buffer attributes (eg ObjectID or another primary key) to filter out the intersections belonging to a certain buffer. With a Python script, it would look like this: parcel_layer = "Parcels" buffer_layer = "Buffers" output_class = "IntersectParcelsBuffers" # intersect result = arcpy.analysis.Intersect([parcel_layer, buffer_layer], output_class) # create a layer of this feature class for each buffer buffer_id_field = "FID_Buffers" # ID field in the result table buffer_ids = {row[0] for row in arcpy.da.SearchCursor(result, [buffer_id_field])} for buffer_id in sorted(buffer_ids): arcpy.management.MakeFeatureLayer(result, f"Buffer {buffer_id}", f"{buffer_id_field} = {buffer_id}")     If you really want to use Clip, you can select each buffer and run it separately. This will take much longer and save a result feature class for each buffer (as opposed to the one result table for Intersect). parcel_layer = "Parcels" buffer_layer = "Buffer" output_gdb = "" # default gdb, change if you want buffer_id_field = "OBJECTID" buffer_ids = {row[0] for row in arcpy.da.SearchCursor(buffer_layer, [buffer_id_field])} import os for buffer_id in buffer_ids: arcpy.management.SelectLayerByAttribute(buffer_layer, "NEW_SELECTION", f"{buffer_id_field} = {buffer_id}") arcpy.analysis.Clip(parcel_layer, buffer_layer, os.path.join(output_gdb, f"Clip_Buffer_{buffer_id}"))   ... View more

You don't save your snow_values anywhere, so snow_values is just the date part of the last raster in location_a, which is why the comparison in line 28 returns False. You have location_a and location_b the wrong way around (or you mislabeled your screenshots). Your variable names are really misleading. Here's what you want to do (untested): # build a dictionary {date: raster} for the snow rasters # I'm going by your screenshots: snow rasters in location_b snow_raster_dict = dict() arcpy.env.workspace = location_b snow_raster_list = arcpy.ListRasters("", "TIF") for snow_name in snow_raster_list: snow_raster = Raster(snow_name) #print(type(snow_raster)) snow_name_parts = snow_name.split("_") snow_date = snow_name_parts[2] snow_raster_dict[snow_date] = snow_raster #print(snow_raster_dict) # get the names of airtemp rasters # again, going by your screenshots -> location_a arcpy.env.workspace = location_a airtemp_raster_list = arcpy.ListRasters("", "TIF") #print(airtemp_raster_list) for airtemp_name in airtemp_raster_list: airtemp_raster = Raster(airtemp_name) airtemp_name_parts = airtemp_name.split("_") airtemp_date = airtemp_name_parts[0] # try to get the corresponding snow raster try: snow_raster = snow_raster_dict[airtemp_date] # do the calculation output_raster = Con(snow_raster==1, 0, airtemp_raster) output = os.path.join(location,airtemp_name.split("_")[0] +'_airtemp'+ '.tif') output_raster.save(output) #print(f"created new airtemp raster for {airtemp_date}") # catch KeyErrors (date not found in snow_raster_dict) except KeyError: print(f"no snow raster found for {airtemp_date}")   ... View more

from pathlib import Path def get_folder(dataset): parent = Path(dataset).parent if "." in parent.name: # .gdb, .sde return str(parent.parent) return str(parent) ... View more

Your variable is named width_cat, but you're using width in the When() function. ... View more

This is absolutely possible. Instead of specifying the field for the rule and returning a value, you leave the field empty and return a dictionary. This dictionary has to have a defined structure, which you can find here: https://pro.arcgis.com/en/pro-app/latest/help/data/geodatabases/overview/attribute-rule-dictionary-keywords.htm   A simple example:   // Calculation Attribute Rule // field: empty // triggers: insert return { edit: [ {className: "OtherTable", adds: [ {attributes: {TextField: "Value", IntegerField: 5}} ] } ] }     EDIT: And you can do lots of other things, too. A few examples are listed in this blog post: Advanced Attribute Rules - Editing features on another class with attribute rules (esri.com)   ... View more

Your error is in the if statement. Most programming languages use "=" as assignment (var x = 3) and "==" as check for equal (x == 3 --> true). if($feature["Data_record_complete"] == "YES") { return Now() }   While the default return value of Date() is the current date and time, I would use Now() here. It returns the same value but when reading the code it's clearer what you want to do, as you can construct arbitrary Dates with Date(). Also, depending on your field, you might want to use Today(), which only returns the date without the time.   ... View more