topic Re: Rectangle for Clip_management tool in arcpy in Geoprocessing Questions

Rectangle for Clip_management tool in arcpy

RichNauman — Mon, 02 Aug 2010 15:57:15 GMT

Hello,
I have been struggling to get this part of a script to work. If I hard code the four coordinate pairs: like this:

rectangle = "-121.677618762735 42.7497214280392 -119.35205017646 45.3659860875993"

it works but we have to remember to change this every time we change study areas. When I run the script pasted below I get the following error:

Traceback (most recent call last):
File "C:\Python26\ArcGIS10.0\Lib\site-packages\pythonwin\pywin\framework\scriptutils.py", line 325, in RunScript
    exec codeObject in __main__.__dict__
File "C:\Documents and Settings\Administrator\Desktop\clip_test.py", line 13, in <module>
    arcpy.Clip_management (raster, rectangle, outraster)
File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\management.py", line 7692, in Clip
    raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000840: The value is not a Envelope.
Failed to execute (Clip).

Here is my code:

import arcpy
arcpy.env.workspace = r"C:\Temp"

raster = "some_raster"
outraster = r"C:\output_raster"
extent = r"C:\GIS_Data\test_extent2.shp"

desc = arcpy.Describe(extent)
rectangle = desc.extent

print rectangle

# Clip_management (in_raster, rectangle, out_raster, {in_template_dataset}, {nodata_value}, {clipping_geometry})
arcpy.Clip_management (raster, rectangle, outraster)

print "Success!"

What am I missing?
Thanks!

Re: Rectangle for Clip_management tool in arcpy

ChrisSnyder — Mon, 02 Aug 2010 17:00:44 GMT

Try this:

import arcpy
arcpy.env.workspace = r"C:\Temp"
raster = "some_raster"
outraster = r"C:\output_raster"
extentFC = r"C:\GIS_Data\test_extent2.shp"
dsc = arcpy.describe(extent)
arcpy.Clip_management (raster, dsc.extent, outraster)
###or if you want to mask it using extentFC as the outline:
##arcpy.Clip_management (raster, dsc.extent, outraster, extentFC, "", "CLIPPINGGEOMETRY")

Re: Rectangle for Clip_management tool in arcpy

RichNauman — Mon, 02 Aug 2010 17:20:54 GMT

Thanks but I still get the "The value is not a Envelope." error:

Traceback (most recent call last):
File "C:\Python26\ArcGIS10.0\Lib\site-packages\pythonwin\pywin\framework\scriptutils.py", line 325, in RunScript
    exec codeObject in __main__.__dict__
File "C:\Documents and Settings\Administrator\Desktop\clip_test.py", line 15, in <module>
    arcpy.Clip_management (raster, dsc.extent, outraster)
File "C:\Program Files\ArcGIS\Desktop10.0\arcpy\arcpy\management.py", line 7692, in Clip
    raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000840: The value is not a Envelope.
Failed to execute (Clip).

Re: Rectangle for Clip_management tool in arcpy

ChrisSnyder — Mon, 02 Aug 2010 17:25:03 GMT

Hmm... Maybe a bug in v10?

Here's an excerpt of some v9.3 code that works for me:

#Process: Creates a blank GDB
fgdbName = "oesf_score_gnn_raster"
fgdbPath = root + "\\" + fgdbName + ".gdb"
gp.CreateFileGDB_management(root, fgdbName); showGpMessage()

#Process: Extract the OESF
hcpUnitFC = r"\\Snarf\am\div_lm\ds\gis\tools\sde_connections_read\ropa_gis_layer_user.sde\ROPA.HCPUNIT"
oesfFC = fgdbPath + "\\oesf"
gp.Select_analysis(hcpUnitFC, oesfFC, "HCPUNIT_NM = 'OESF'"); showGpMessage()

#Process: Buffer the oesf
oesfBufferUnits = "FEET"
oesfBufferDistance = "35000" #NOTE: This MUST be the same buffer distance as in 'nso_make_path_distance_grids_v93.py'
oesfBufferFC = fgdbPath + "\\oesf_buff_" + str(oesfBufferDistance) + oesfBufferUnits.lower()
gp.Buffer_analysis(oesfFC, oesfBufferFC, str(oesfBufferDistance) + " " + oesfBufferUnits.upper(), "FULL", "ROUND", "NONE", ""); showGpMessage()

#Process: Clip the GNN data to our area of interest
gnnReadGrd = r"\\snarf\am\div_lm\ds\gis\tools\sde_connections_read\raster_gis_layer_user.sde\SHARED_LM.GNN_20100303"
gp.extent = oesfBufferFC
gp.extent = snapExtentToRaster(gp.extent, gnnReadGrd)
gp.mask = oesfBufferFC
gp.cellsize = gnnReadGrd
gnnClipGrd = fgdbPath + "\\gnn_clip"
gp.Clip_management(gnnReadGrd, gp.extent, gnnClipGrd, oesfBufferFC, "", "CLIPPINGGEOMETRY")

Re: Rectangle for Clip_management tool in arcpy

DavidWynne — Sat, 11 Dec 2021 17:34:05 GMT

Hi Richard,
You won't be able to use a Extent object directly as input to Clip, but since you can use geometry objects, Polygon, etc. So, really it just a matter of converting that extent into a polygon.

Working from your original code, something like below should do it.

import arcpy
arcpy.env.workspace = r"C:\Temp"

raster = "some_raster"
outraster = r"C:\output_raster"
desc = arcpy.Describe(r"C:\GIS_Data\test_extent2.shp")
extent = desc.extent

rectangle = arcpy.Polygon(arcpy.Array([extent.lowerLeft, extent.upperLeft, extent.upperRight, extent.lowerRight, extent.lowerLeft]),
                          desc.spatialReference)

print rectangle

# Clip_management (in_raster, rectangle, out_raster, {in_template_dataset}, {nodata_value}, {clipping_geometry})
arcpy.Clip_management(raster, rectangle, outraster)

print "Success!"

-Dave

Re: Rectangle for Clip_management tool in arcpy

ChrisSnyder — Tue, 03 Aug 2010 17:05:26 GMT

Looks like I will be rewriting a lot of v9.3 code to make it arcpy compatible...

Re: Rectangle for Clip_management tool in arcpy

DavidWynne — Sat, 11 Dec 2021 17:34:08 GMT

Looks like I will be rewriting a lot of v9.3 code to make it arcpy compatible...

Sorry, I was overthinking that example way too much. It could also be done by converting the extent to a string like so:

desc = arcpy.Describe(r"C:\GIS_Data\test_extent2.shp")
extent = desc.extent

# Clip_management (in_raster, rectangle, out_raster, {in_template_dataset}, {nodata_value}, {clipping_geometry})
arcpy.Clip_management(raster, str(extent), outraster)

Hopefully, that is more agreeable
-Dave

Re: Rectangle for Clip_management tool in arcpy

ChrisSnyder — Tue, 03 Aug 2010 17:14:55 GMT

I like that better!

Re: Rectangle for Clip_management tool in arcpy

RichNauman — Tue, 03 Aug 2010 20:19:15 GMT

Thanks for the help but I am still struggling.

This script works:

import arcpy
raster = r"Y:\GIS_Data\raster1"
outraster = r"Y:\GIS_Data\delete_me"

desc = arcpy.Describe(r"Y:\GIS_Data\test_extent2.shp")
frame = str(desc.extent)
frame = str(frame [:-16])
frame = str('"' + frame + '"')
print frame

arcpy.Clip_management (raster, "-121.677618762735 42.7497214280392 -119.35205017646 45.3659860875993", outraster)
print "Done!"

But this one does not:

import arcpy
raster = r"Y:\GIS_Data\raster1"
outraster = r"Y:\GIS_Data\delete_me"

desc = arcpy.Describe(r"Y:\GIS_Data\test_extent2.shp")
frame = str(desc.extent)
frame = str(frame [:-16])
frame = str('"' + frame + '"')
print frame

arcpy.Clip_management (raster,frame, outraster)
print "Done!"

I copied and posted the output from the "print frame" statement into the clip statement to get it to work.

Any ideas?

Re: Rectangle for Clip_management tool in arcpy

DavidWynne — Sat, 11 Dec 2021 17:34:11 GMT

desc = arcpy.Describe(r"Y:\GIS_Data\test_extent2.shp")
frame = str(desc.extent)
frame = str(frame [:-16])
frame = str('"' + frame + '"')

Don't bother with the extra quotes in this last line, I'm pretty sure this is the cause of your error.

Instead, just do this (leaving out your last frame = line)...

frame = str(desc.extent)
frame = frame[:-16]

or, if your intention is to strip out all the NaN's at the end, you could also do this...

frame = str(desc.extent).strip(" NaN")

Re: Rectangle for Clip_management tool in arcpy

RichNauman — Sat, 11 Dec 2021 17:34:13 GMT

Thanks for the help!!!

I'm not sure where I went off track but the following works:

import arcpy
raster = r"c:\GIS_Data\some_raster"
outraster = r"c:\GIS_Data\delete_me"

desc = arcpy.Describe(r"c:\GIS_Data\test_extent2.shp")
frame = str(desc.extent)

arcpy.Clip_management (raster, frame, outraster)
print "Done!"

Re: Rectangle for Clip_management tool in arcpy

JenniferCorcoran — Sat, 19 Mar 2011 21:48:29 GMT

Have you tried applying this code to several rasters from different resolution source imagery? Do you get the same extents in the long run? I am having major trouble clipping my rasters to specific extents and having them come out the EXACT same. I keep getting some of them the same as I am trying to use and some with extra decimals (ex. Raster Extent Interested: 524002.75 5167611.25 542557.75 5188001.25; Raster Extent Result: 524002.74999979 5167611.25012111 542557.74999979 5188001.25012111).

Thoughts?