I have a rectangular raster image in ArcMap.
How can I produce a polygon which is just the extent of the raster image?
That is, I want to have a polygon layer with a single quadrilateral which is the edge of the image.
I've tried simply doing a Raster->Polygon conversion, but this tries to use the data in the image to generate the outline - all I want is the outline of the edge of the image.
I'm using ArcGIS Desktop 10.
If you have the "Spatial Analyst" then you can:
Another solution:
If your have 3D Analyst license then you can use the Raster Domain tool under 3D Analyst>Conversion>from raster (however it does not have the smooth function)
You can use the Raster Domain tool from 3D Analyst to create a polygon with the same extent as your raster data
I use this simple python script, it create a polygon featureclass with the extent of all the raster present in a folder. To use it, you can create a tool in ArcToolbox or you simply change the InFolder and Dest (Destination) variables.
import arcpy,os
InFolder = arcpy.GetParameterAsText(0)
Dest=arcpy.GetParameterAsText(1)
arcpy.env.workspace=InFolder
#The raster datasets in the input workspace
in_raster_datasets = arcpy.ListRasters()
arcpy.CreateFeatureclass_management(os.path.dirname(Dest),
os.path.basename(Dest),
"POLYGON")
arcpy.AddField_management(Dest,"RasterName", "String","","",250)
arcpy.AddField_management(Dest,"RasterPath", "String","","",250)
cursor = arcpy.InsertCursor(Dest)
point = arcpy.Point()
array = arcpy.Array()
corners = ["lowerLeft", "lowerRight", "upperRight", "upperLeft"]
for Ras in in_raster_datasets:
feat = cursor.newRow()
r = arcpy.Raster(Ras)
for corner in corners:
point.X = getattr(r.extent, "%s" % corner).X
point.Y = getattr(r.extent, "%s" % corner).Y
array.add(point)
array.add(array.getObject(0))
polygon = arcpy.Polygon(array)
feat.shape = polygon
feat.setValue("RasterName", Ras)
feat.setValue("RasterPath", InFolder + "\\" + Ras)
cursor.insertRow(feat)
array.removeAll()
del feat
del cursor
You can use the Build Footprint tool in the Data management Toolbox. You have to create a Mosaic Dataset out of your raster.
Create a new Mosaic in a File Geodatabase , add your raster.
Open your Mosaic in ArcMap and Extract the Footprint.
That's what I would do.
What about simply digitizing it? Click-click-click-click done.
Also, an easy way is to use Reclassify tool where you should click the button Classify and change the number of classes to 1. Click Ok. Now, the output raster can be easily converted using Raster to Polygon conversion Tool. ;)
I think that reclassifying and converting the raster to polygon can be way too time consuming. In my opinion, the easiest way to do it is:
If you have lots of rasters, this should be done in Python and can create hundreds of footprints in a matter of seconds.
In ArcCatalog you can create a raster catalog in a .gdb (Data Management Tools->Raster). Unmanaged is fine. Then right click and load (if the rasters are just sitting in a folder then "load from workspace". If you load the raster catalog into ArcMap you can right click-> data-> export footprint, as a shape file or feature class. The shape file will have a 'name' field populated with the name of the raster.
The Raster Catalog will automatically display the extents but I needed to send the extents to someone who wasn't using ArcGIS.
Use the "Raster to Polygon" tool in the conversion tools toolbox. Make sure to uncheck the simplify polygons box. This will create exactly what you want.
"Raster Domain" tool will create a polygon around an irregularly shaped raster. I believe it requires a license for the 3d analyst extension.
Using ArcCatalog is the best option. Create a Mosaic Dataset from available rasters. Define "NO VALUE" in the mosaic to avoid the background of raster from being included in the polygon. Build footprints using Radiometric method. Open the Mosaic dataset in ArcGIS and export Footprints to shp file. Great way to get the mosaic boundaries...
If you are trying to get the outline of a DEM raster then you will need to do the following
If you have Spatial Analyst and 3d Analyst then you can:
I read this post for advice on a similar need. I wanted to outline a raster that is very irregularly shaped, with a stair-step pattern around all the curves, which would have taken hours to digitize. All the other recommendations, mosaic, reclassify, raster calculator, raster domain, etc., all require higher level licenses that I don't have. I solved this problem by using the free program QGIS for parts of it. I loaded my raster into Q. I ran it through the tool Raster>Conversion>polygonize, which converted it into a polygon layer with many thousands of tiny polygons. I exported that layer as an ESRI shapefile, and brought that into Arcmap. Once in arcmap, I began an editing session, opened the attribute table, and selected all of the many thousands of polygons. The selection could have been more easily accomplished with the selector tool, but either way... From the Editor toolbar's main dropdown menu, I selected Union (different from the Union tool from the toolbox). It created the outline I wanted, already highlighted. From the attribute table, I clicked on switch selection, and then right clicked the selected polygons and deleted the selection, so all I had left in the layer was the outline. Then saved edits and stopped editing.
You can do it easily with "arcpy" module in Python: The important thing is that you have to read the spatial georeference system of the input raster, and create a polygon feature with the same spatial georeference system. You read the corners of the raster and add points to a polygon shape. Then you add the shape to your feature class. This is the code for doing it.
import arcpy
import sys
import os
fp_inp = "C:\\raster\\raster.tif"
fp_out_dir, fp_out_name = "C:\\raster", "polygon.shp"
fp_out = os.path.join(fp_out_dir, fp_out_name)
if arcpy.Exists(fp_inp):
desc = arcpy.Describe(fp_inp)
spatialName, spatRef = "Unkown", None
if hasattr(desc, 'spatialReference'):
spatRef = desc.spatialReference
spatialName = desc.spatialReference.name
if spatialName != "" and spatialName != "Unkown":
arcpy.env.workspace = fp_out_dir
try:
arcpy.CreateFeatureclass_management(fp_out_dir, fp_out_name, 'POLYGON', spatial_reference=spatRef)
polygonSuc = True
except:
print("ERROR (126401236004)")
print(arcpy.GetMessages())
fieldsOk = False
if polygonSuc:
try:
arcpy.AddField_management(fp_out, "src", "TEXT", "", "", 254, "src",
"NULLABLE")
fieldsOk = True
except:
msg01 = "\n" + arcpy.GetMessages()
print(msg01)
corners = ["lowerLeft", "lowerRight", "upperRight", "upperLeft"]
ras = arcpy.Raster(fp_inp)
points_array = arcpy.Array()
for corner in corners:
pt = arcpy.Point(getattr(ras.extent, "%s" % corner).X, getattr(ras.extent, "%s" % corner).Y)
points_array.add(pt)
polygon_shape = arcpy.Polygon(points_array)
if fieldsOk:
with arcpy.da.InsertCursor(fp_out,
["SHAPE@", "src"]) as cursor:
cursor.insertRow([polygon_shape, "your raster name goes here"])
print("Inserted polygon into shapefile: " + fp_out)
points_array.removeAll()
array = None
else:
print("Input raster does not exist.")
