I'm writing a that needs to clip a multidimensional raster (Gridmet temperature

<a class="user-mention notranslate" data-hovercard-type="user" data-hovercard-url="/us

error running Clip raster function using Foreach loop about raster-functions HOT 2 OPEN

esri commented on May 28, 2024

error running Clip raster function using Foreach loop

from raster-functions.

Comments (2)

shenganzhan commented on May 28, 2024

@DavidAnderson-USFS
The clip tool internally converts a feature to a geometry JSON and use that as ClippingGeometry. For some reason we are investigating, the conversion is not done when used in Foreach. Please try this snippet as workaround:

import json

temp_max_file = r"D:\r3_analyst_work\Climate_Change_Modeling\climate_envelope\GridMetData_2020\tmmx_2020.nc"
max_temp_rast = arcpy.Raster(temp_max_file,True)
clip_frame= r"D:\r3_analyst_work\GIS_work\test_working_netcdf_climate_data\Default.gdb\small_area"

geo_list =[]
for row in arcpy.da.SearchCursor(clip_frame, ["OID@", "SHAPE@"]):
    geo_list.append(row[1])
geo_union = geo_list[0]
for i in range(1, len(geo_list)):
    geo_union = geo_union.union(geo_list[i])
geometry_inp = json.loads(geo_union.JSON)

max_temp_clip_raster = arcpy.sa.Foreach(max_temp_rast,"Clip",{"ClippingGeometry":geometry_inp, "ClippingType":1})

from raster-functions.

DavidAnderson-USFS commented on May 28, 2024

@shenganzhan I was able successfully implement this workaround. Thank you for rapid response in providing a solution.

@DavidAnderson-USFS The clip tool internally converts a feature to a geometry JSON and use that as ClippingGeometry. For some reason we are investigating, the conversion is not done when used in Foreach. Please try this snippet as workaround:

from raster-functions.

Recommend Projects