topic Re: How to Get Polygon of Raster Extent in ArcGIS Pro in ArcGIS Pro Questions

How to Get Polygon of Raster Extent in ArcGIS Pro

dgray — Thu, 09 Sep 2021 20:33:20 GMT

Hi,

I would like to know a way, or several methods if they exist, of generating a rectangular polygon representing the footprint of a raster using ArcGIS Pro. I don't want to go through the extra time to convert a raster to polygon as is, because this is a potentially very time-consuming step (especially if the raster is an image), whereas I only want a single-valued output.

If this method does not require higher licensing levels or extensions, that would be even better, but if not that's ok too.

If there's a custom Python script that will do that job, or a snippet of Python I could insert into a script, please let me know.

Thanks!

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

DanPatterson — Wed, 08 Sep 2021 01:47:12 GMT

If you can get the extent of the raster, you can create a polygon then add that to a new featureclass

Extent—ArcGIS Pro | Documentation

If you prefer a simple raster representing the extent, then

Create Constant Raster (Spatial Analyst)—ArcGIS Pro | Documentation

using a value of 1 and the cell size and extent of your current raster. A constant raster... even if temporary... will convert to a polygon featureclass quite quickly since it only has one value

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

Elijah — Wed, 08 Sep 2021 05:20:48 GMT

I think Raster Domain 3D tool can do it. It outputs a polygon/line geometry of the input raster.

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

DrVSSKiran — Wed, 08 Sep 2021 18:42:59 GMT

Mosaic datset also generqte the raster foot print based on raster extent wothout any add on extension licenses.

Here is the link of arcmap, but similar you can try for pro. I have done in pro 2.0 and later versions

https://desktop.arcgis.com/en/arcmap/10.3/manage-data/raster-and-images/mosaic-dataset-footprints.htm

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

dgray — Wed, 08 Sep 2021 21:35:55 GMT

Thanks for your suggestions so far. I will test the mosaic method, as I'm looking specifically for a way to not worry about licensing levels and extensions.

So far, I've created a working Python script (derived from another script meant to composite satellite imagery), which generates extent shapefiles of a satellite imagery raster based on min and max coordinates. The biggest problem is that it assumes all rasters are oriented perfectly upright. So it did what I asked it to do, but it's not actually the real extent, as my source rasters are angled.

Is there a way to get the corner coordinates of my raster's extent, not the mins and maxes? Please take a look at my script (remove the .txt from the file) and provide your thoughts.

Thanks!

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

DanPatterson — Wed, 08 Sep 2021 22:10:35 GMT

xmin = extent.XMin xmax = extent.XMax ymin = extent.YMin ymax = extent.YMax lowerLeft = arcpy.Point(xmin, ymin) lowerRight = arcpy.Point(xmax, ymin) # etcetera

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

dgray — Wed, 08 Sep 2021 22:33:11 GMT

Yeah, but if you have an angled raster (e.g. diamond shape), the lower left corner will not have a coincidence of xmin and ymin, etc. You'll have a ymin, but the x value will be somewhere between the min and max extents of x. Same with y...the upper right corner will have the ymax value, but the x value isn't necessarily at its max.

What I'm really wanting to do is create a polygon that shows the exact footprint of the raster, regardless of its orientation. I guess I'm saying the word extents, but it's the footprint I'm after, if I'm using terminology correctly.

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

DanPatterson — Wed, 08 Sep 2021 23:00:42 GMT

That can't be derived from the extent obviously, you will have to convert the nodata collar to a polygon, which you could get by classifying the nodata collar to 1 and everything else to 0. Your probably multipart polygon would then have to be queried for each part's extent

A simpler solution would be to return the intersection points of the nodata and data collars

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

dgray — Wed, 08 Sep 2021 23:09:37 GMT

Hmm, ok, can this be done without extra licensing levels or extensions? How do you get the data and nodata collars?

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

dgray — Thu, 09 Sep 2021 20:04:18 GMT

So I ended up creating a script using the Mosaic dataset footprint generation. It creates extra intermediate steps but it isn't too bad--I just set a scratch workspace and delete the content later.

However, I'm still interested in learning more about the data and nodata collar method.

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

ClaireMcHenry2 — Thu, 19 Oct 2023 21:19:59 GMT

Did anything ever come from this? I am attempting to batch export the extent of multiple historic maps (all raster). I suppose I am looking for the footprint as well and not the extent itself.

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

JohnGaiot — Fri, 16 Feb 2024 21:27:11 GMT

It seems like a basic function but in my experience, it's not. There is one function that comes very close and it's a quick one-step solution. It's called "Raster Domain".. However, the Downsides: it requires 3D Analyst, and more importantly, the final polygon goes through the centroid of each cell at the data extent of the raster, so it's not an exact 'footprint' as it falls a half cell short on all sides of features represented. But if you want something quick and approximate, this would be the ideal solution.

For something more exact, the simplest method would be to use Raster Calculator (which requires a Spatial or Image Analyst license). Set the "NoData" Environment variable to something other than Zero. If you have a floating point or non-discrete raster like a DEM, you will first need to convert your raster to integer. Second, multiply the raster by Zero, so you get a single value for the entire raster where there were valid values. Third, convert the Zero raster to a polygon. Now you have your extents based on features in a raster.

Re: How to Get Polygon of Raster Extent in ArcGIS Pro

Town_ofCBS_GIS — Fri, 05 Apr 2024 16:34:01 GMT

Hello, I was also trying to do this so I could label all my rasters. To do this I had to create a polygon layer. Make sure the project you are working on is saved and closed, then open the python notebook in a different project. Also, I only had one 'Map' - when I had multiple it would not work. Below is my python code:

import arcpy

# Set the workspace environment
arcpy.env.workspace = r"C:\path\to\your\geodatabase.gdb" # Update the path to your geodatabase

# Input map project file (.aprx)
map_project = r"C:\path\to\your\map_project.aprx" # Update the path to your map project file (.aprx)

# Reference the map project
aprx = arcpy.mp.ArcGISProject(map_project)

# Get all the maps in the project
maps = aprx.listMaps()

# Create a polygon feature class to store layer names as polygons
polygon_fc_name = "LayerNamesPolygon"
polygon_fc_path = arcpy.management.CreateFeatureclass(arcpy.env.workspace, polygon_fc_name, "POLYGON").getOutput(0)

# Add a field to store layer names in the polygon feature class
arcpy.management.AddField(polygon_fc_path, "Layer_Name", "TEXT")

# Iterate through each map in the project
for map_obj in maps:
# Iterate through each layer in the map object
for lyr in map_obj.listLayers():
if lyr.isRasterLayer:
# Get the raster layer name
raster_layer_name = lyr.name

# Get the extent of the raster layer
desc = arcpy.Describe(lyr)
extent = desc.extent

# Create a polygon from the extent
array = arcpy.Array([extent.lowerLeft, extent.lowerRight, extent.upperRight, extent.upperLeft])
polygon = arcpy.Polygon(array)

# Insert the polygon with layer name into the feature class
with arcpy.da.InsertCursor(polygon_fc_path, ["SHAPE@", "Layer_Name"]) as cursor:
cursor.insertRow([polygon, raster_layer_name])

# Save changes to the map project
aprx.save()

# Delete the ArcGISProject object to release locks
del aprx