topic Re: Expand extent for purposes of clipping in Python Questions

Expand extent for purposes of clipping

RebeccaStrauch__GISP — Fri, 10 Dec 2021 21:21:14 GMT

I'm wanting to clip a hydro layer (with a reselect) based on the extent of my study area extent....but expanded a bit (% to be determined). I have code that works, but it's a bit clunky. I grab the extent, create a poly, buffer, grab-new-extent, create the poly and then do my clip/query. It's the getting the extent->poly->buffer->extent->poly that seems a bit overkill. I'm wondering if I am missing an obvious arcpy option for expanding the extent for the clip. Trying to keep it to arcpy and not get deep into the geometry of this one...but all suggestions will be considered.

Thanks.

This works:

import arcpy
from arcpy import env
import os
arcpy.env.overwriteOutput = "True"

inStudy =  r'C:\Prep.gdb\study_boundary'
srDesc = arcpy.Describe(inStudy).spatialReference
extent = arcpy.Describe(inStudy).extent
arcpy.env.outputCoordinateSystem = srDesc.PCSCode

# ---> Ugly way to do this, but it works and fairly fast.  Ok for now.
# get extent, a buff distance (1% of X), and create a temp polygon
extBuffDist = ((int(abs(extent.lowerLeft.X - extent.lowerRight.X))) * .01)

# Array to hold points for the bounding box for initial extent
origExtentPts = arcpy.Array()
origExtentPts.add(extent.lowerLeft)
origExtentPts.add(extent.lowerRight)
origExtentPts.add(extent.upperRight)
origExtentPts.add(extent.upperLeft)
origExtentPts.add(extent.lowerLeft)  # ensures polygon closes
# Create and buffer the polygon object
polygonTmp1 = arcpy.Polygon(origExtentPts)
arcpy.Buffer_analysis(polygonTmp1, "polyBuff", extBuffDist, "OUTSIDE_ONLY")
origExtentPts.removeAll()

#get extent of buffered poly
extent = arcpy.Describe("polyBuff").extent
# Array to hold points for the buffer bounding box
newExtentPts = arcpy.Array()
newExtentPts.add(extent.lowerLeft)
newExtentPts.add(extent.lowerRight)
newExtentPts.add(extent.upperRight)
newExtentPts.add(extent.upperLeft)
newExtentPts.add(extent.lowerLeft) # ensures polygon closes
polygonTmp2 = arcpy.Polygon(newExtentPts)

# save to disk
#extentPoly = "extentPoly" 
arcpy.CopyFeatures_management(polygonTmp2, "extentPoly")
del polygonTmp1, polygonTmp2

hydroNHDSource = r"\\myServer\NHDH_AK.gdb\Hydrography\NHDFlowline"
hydroNHDQuery = "(GNIS_Name >= '1' AND ( FType = 566 OR FType = 460 OR FType = 558))"

arcpy.Clip_analysis(hydroNHDSource, "extentPoly", "hydroExtent")   
arcpy.MakeFeatureLayer_management("hydroExtent", "HydroReselect", hydroNHDQuery)
arcpy.Buffer_analysis("HydroReselect", "hydroBuff", 500, "FULL", "ROUND", "ALL" )
hydroBuffStudy = arcpy.Clip_analysis("hydroBuff", inStudy)

Re: Expand extent for purposes of clipping

DarrenWiens2 — Fri, 10 Dec 2021 21:21:17 GMT

You can make use of the buffer method of the polygon object to make a new buffered polygon geometry:

polygonTmp1 = arcpy.Polygon(origExtentPts)
buffPoly = polygonTmp1.buffer(extBuffDist)
newExtent = buffPoly.extent # this returns an 'extent' object. May have to convert to 'polygon'

Re: Expand extent for purposes of clipping

DanPatterson_Retired — Wed, 09 Dec 2015 20:09:03 GMT

Darren, just grab the points from the extent to create the polygon, unfortunately, one would assume that they were a variant of a list but you have to go through the process Rebecca of getting them by name to create the polygon. This hasn't even changed in Pro Extent—ArcPy Classes | ArcGIS for Desktop but at least your buffer cuts out a step.

Re: Expand extent for purposes of clipping

RebeccaStrauch__GISP — Fri, 10 Dec 2021 21:21:19 GMT

Darren's answer was helpful, but as you mention Dan, I think I would still have to go thru the process of grabbing the LL, LR, etc corners into an array to make the poly. ?? I guess it's the extent->poly that I was hoping to simplify, i.e. without the array. Maybe not possible?

As a side note, the code Darren provided did allow me to remove two lines...

arcpy.Buffer_analysis(polygonTmp1, "polyBuff", extBuffDist, "OUTSIDE_ONLY")
#...
extent = arcpy.Describe("polyBuff").extent

Re: Expand extent for purposes of clipping

DanPatterson_Retired — Fri, 10 Dec 2021 21:21:22 GMT

Yes you still have to get the values since they aren't variants of lists, but so you can't slice them but have to get them by name.

for row in arcpy.da.SearchCursor(in_FC, ["SHAPE@"]):
    extent = row[0].extent
    L,B,R,T = [extent.XMin, extent.YMin, extent.XMax, extent.YMax]
    frmt = "arcpy extent object {}\nScraped values for ... L {} B {} R {} T {}"
    print(frmt.format(extent, L, B, R, T))

but then at least you add/subtract from L,B,R,T to produce new extent objects and formulate a new point if you wanted to. or just use that after buffering the r[0] object above

Re: Expand extent for purposes of clipping

DarrenWiens2 — Wed, 09 Dec 2015 21:04:01 GMT

I suppose since you're expanding the extent by some number in all directions, you could just add/subtract that number to/from your original extent, rather than running a buffer process at all (I'm not sure if that's what Dan was getting at. If so, carry on).

Re: Expand extent for purposes of clipping

RebeccaStrauch__GISP — Fri, 10 Dec 2021 21:21:25 GMT

I took bits from both you and Dan and was able to get it to something that works (and give the same results, as I needed). I didn't need a cursor as Dan pointed out (overkill for this project) but reassigning the extent.lowerLeft etc so shorter variables and then creating the array with a for loop worked for me. I wanted to use a buffer since that was the easiest/cleanest to move the extent out equal, in all directions.

I could use the output of the temp polygon object created from the extent of the buffer in my clip, so all works fine. Final snippet of code..

import arcpy
from arcpy import env
import os
arcpy.env.overwriteOutput = "True"

inStudy =  r'C:\Prep.gdb\study_boundary' # r"c:\temp\USA.gdb\States"
srDesc = arcpy.Describe(inStudy).spatialReference
arcpy.env.outputCoordinateSystem = srDesc.PCSCode

# get extent, a buff distance (1% of X), and create a temp polygon
extent = arcpy.Describe(inStudy).extent
ll, lr, ur, ul = [extent.lowerLeft, extent.lowerRight, extent.upperRight, extent.upperLeft]
extBuffDist = ((int(abs(ll.X - lr.X))) * .01)
origExtentPts = arcpy.Array()
# Array to hold points for the bounding box for initial extent
for coords in [ll, lr, ur, ul, ll]:
    origExtentPts.add(coords)

polygonTmp1 = arcpy.Polygon(origExtentPts)
# buffer the temporary poly by 1% of width of extent as caluculated above
buffPoly = polygonTmp1.buffer(extBuffDist)
newExtent = buffPoly.extent

# extent points for polygonTmp2
ll2, lr2, ur2, ul2 = [newExtent.lowerLeft, newExtent.lowerRight, newExtent.upperRight, newExtent.upperLeft]
newExtentPts = arcpy.Array()
for coords in [ll2, lr2, ur2, ul2, ll2]:
    newExtentPts.add(coords)

polygonTmp2 = arcpy.Polygon(newExtentPts)

hydroNHDSource = r"\\myServer\NHDH_AK.gdb\Hydrography\NHDFlowline"
hydroNHDQuery = "(GNIS_Name >= '1' AND ( FType = 566 OR FType = 460 OR FType = 558))"

arcpy.Clip_analysis(hydroNHDSource, polygonTmp2, "hydroExtent") 
arcpy.MakeFeatureLayer_management("hydroExtent", "HydroReselect", hydroNHDQuery)
arcpy.Buffer_analysis("HydroReselect", "hydroBuff", 500, "FULL", "ROUND", "ALL" )
hydroBuffStudy = arcpy.Clip_analysis("hydroBuff", inStudy)

At least the array is easier to read. So, I think Dan get's the correct answer for pointing out that I couldn't really get around the array for the coordinates.

For anyone that cares, this is a part of a much larger project, but I needed to buffer a subset of the rivers within a study area, but since the NHD hydro in our area is still a bit fuzzy, and many of our study areas use the topo maps to create the boundaries (many times based on hydro), I needed to buffer the rivers before I clipped.

Temporary layers..

and the final result (study area and buffers)

Thanks Darren and Dan for you thoughts on this. It helped a lot.

Re: Expand extent for purposes of clipping

DarrenWiens2 — Wed, 09 Dec 2015 21:50:55 GMT

I see you've got this figured out, but for this, why fight with the extents at all? Just buffer the study area, clip those streams in buffered area, buffer streams, clip back to study area. Perhaps it takes more time to buffer the study area than to clip a bunch of extra streams. Maybe it's six one way, half dozen the other, but I'd guess dealing with the fewest streams possible would be a good idea.

Re: Expand extent for purposes of clipping

RebeccaStrauch__GISP — Wed, 09 Dec 2015 22:01:54 GMT

Ahh, yes good thoughts, and in many cases that would be better, but I'm using the same clip boundary for a bunch of other layers I'm clipping in the same script (topo maps, DEM, etc) that I want to be that same extent+, so this was just a piece of what I needed. I just simplified the code to keep that part out, but worth noting for others.

Re: Expand extent for purposes of clipping

DanPatterson_Retired — Fri, 10 Dec 2021 21:21:27 GMT

Rebecca Strauch, GISP alternate ways to get extent parameters to simplify code

>>> import numpy as np
>>> ext = extent.__str__()
>>> L, B, R, T, *ZZMM = np.fromstring(ext,dtype=float,sep= " ")
>>> l, b, r, t, *zzmm = [float(i) for i in ext.split(" ")]

Where lines 2 and 4 are the key lines are for pure python implementations and 2 and 3 are if numpy is used.... of course the letters represent their obvious values and *zzmm variants are the extra z, m min, max values if they are needed.

Re: Expand extent for purposes of clipping

JoshuaBixby — Sat, 12 Dec 2015 19:14:27 GMT

I am not sure if the OP stated whether this is for ArcGIS for Desktop or ArcGIS Pro. Since extended iterable unpacking (PEP 3132) wasn't introduced until Python 3.0, neither lines 3 nor 4 will work with ArcGIS for Desktop since it is stuck with Python 2.x.

Re: Expand extent for purposes of clipping

RebeccaStrauch__GISP — Sat, 12 Dec 2015 19:22:42 GMT

Joshua, it is initially in a toolbox (so could be for either), but will be part of a python addin which aren't support by Pro. I'm trying to keep it simple, and in arcpy as much possible. Improved speed of numpy isn't necessary in this case, since this is a pretty minor part....the clipping is the main part.

Re: Expand extent for purposes of clipping

JoshuaBixby — Fri, 10 Dec 2021 21:21:30 GMT

Not sure what version of ArcGIS you are using, but starting at ArcGIS 10.3 the "Extent object now supports JSON and polygon properties."

>>> ext = arcpy.Describe(inStudy).extent
>>> new_ext_poly = ext.polygon.buffer(int((ext.XMax - ext.XMin)*0.01)).extent.polygon
>>>

If you don't mind just a little rounding of the extent corners, you could drop .extent.polygon from the end of line 2.

Re: Expand extent for purposes of clipping

DanPatterson_Retired — Sun, 13 Dec 2015 00:17:02 GMT

ooh yea.... you would have to expand the *zzmm term out in full... such a nice enhancement when you only need the first few and can't remember the rest...but it does save all that XMin...Xmax explicit calling stuff though. This would be a an area that would be nice to see arcpy go...reflecting more generic objects than proprietary

Re: Expand extent for purposes of clipping

JoshuaBixby — Sun, 13 Dec 2015 16:59:01 GMT

Extended iterable unpacking is a fantastic enhancement, I really wish it could have made its way backwards to 2.7. Of course, Esri could help the situation by moving ArcGIS for Desktop to Python 3.x, but I am becoming suspicious that will never happen.

I am sure I have lamented about this before, but I think Esri really missed an opportunity with ArcGIS Pro to take ArcPy to the next level. In some sense, that was done out of necessity with the ArcPy Mapping module (arcpy.mapping => arcpy.mp) because the interface is so different between Desktop and Pro, but the changes to all of the other modules and base package were minimal. As has been brought up in GeoNet before, it would be great to see the ArcPy Geometry classes expanded and made much more Pythonic, wouldn't have the change to ArcGIS Pro been a great break in the product line to introduce such new functionality? Oh well, c'est la vie in Esri-land.