nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Get smallest extent from a list of rasters - different approach

Andreas

Original thread located here: Get smallest extent from a list of rasters

I had issues with that method, if some extents do not overlap, arcgis extent is "empty" while the above minof is (always?) still a rectangle!
In the docs it is said thats minof is the intersection of inputs. But that seems not to be a Intersect_analyst-equivalent. I tried with three extent objects, two do not overlap and one overlaps all other.

Message was edited by: Xander Bakker (branched to new thread)

extents.png

Find more posts tagged with

Accepted answers

XanderBakker

If the calculation was done by collecting all the minx, miny, maxx and maxy and use the resulting mins and maxs as points to construct the extent, the result will be as you showed in your image. This is however, as you concluded yourself, not the correct result.

Let's take a look at the code below that shows what happens:

import arcpy

def main():
    ext1 = arcpy.Extent(1, 0, 9, 5)
    ext2 = arcpy.Extent(0, 1, 10, 10)
    ext3 = arcpy.Extent(0, 6, 6, 12)

    print "getIntersectingExtent", getIntersectingExtent([ext1, ext2, ext3])
    print "getMinsMaxsExtentsAsExtent", getMinsMaxsExtentsAsExtent([ext1, ext2, ext3])

def getIntersectingExtent(lst):
    ext = lst[0]
    for e in lst[1:]:
        if ext.overlaps(e) or ext.contains(e) or ext.within(e) or ext.equals(e):
            ext = getOverlapExents(ext, e)
        else:
            return None
            break
    return ext

def getOverlapExents(ext1, ext2):
    return getMinsMaxsExtentsAsExtent([ext1, ext2])

def getMinsMaxsExtentsAsExtent(lst):
    xmin = max([e.XMin for e in lst])
    xmax = min([e.XMax for e in lst])
    ymin = max([e.YMin for e in lst])
    ymax = min([e.YMax for e in lst])
    return arcpy.Extent(xmin, ymin, xmax, ymax)

if __name__ == '__main__':
    main()

The result of the code is:

getIntersectingExtent None
getMinsMaxsExtentsAsExtent 1 5 6 6 NaN NaN NaN NaN

Both use the same procedure where the min and max of the extents is determined. However, the first result will test is there is overlap and then determine the minimums and maximums. The second will simply determine the minimums and maximums of the extents and when provided to the arcpy.Extent() method this will correct for minimums being larger than maximums.

All comments

DanPatterson_Retired

not sure

yellow blue intersect produces max-X and max-Y for MIN-OF

purplish blue intersect produces min-X and min-Y for MIN-OF

IF min-of is your left-over bit

Now I am wondering if the order affects the results. overlay blue last...

order as yellow, purple then blue

compare to

blue, then yellow, purple

Andreas

Thank you for your contribution,

I got an extent "between" the non-overlapping rasters, but it should be none, i think.

The Problem araised with a tool that shoud calculate a hexagonal fishnet, where the extent is provided by the user, with an extent parameter in a script tool. Now my question was how to calculate the extent, if the parameter value is MINOF and was glad to see your solution. If i use a tool, like Copy_management and an MINOF in environments extent, the tool creates no output, because the two of three extents do not overlab. But this isn't possible to get with the approach above?

I attached a new sketch to show what i got.

Andreas

extents2.png

XanderBakker

True, and that's what the script from the original post did. That's why I would like to hear what the desired output should be, to see if something can be created to match that...

DanPatterson_Retired

From environment settings Output Extent (Environment setting)—Help | ArcGIS for Desktop

With intersection of inputs...

Intersection of Inputs—The extent where all input features or rasters overlap (intersect one another). Note that it is possible that none of the features overlap and that a null extent (zero width and height) may result. In such cases, no features will be processed.

XanderBakker

Could you elaborate on what the desired output should be based on the three rasters you mentioned?

Join the Conversation

Signed in members can post, follow updates, and more. New here? Register a free account.

Quick Links

Community Help
Find useful guides, FAQs, and documents to help you navigate and make the most of Esri Community.