¿Obtener los límites de una imagen rasterizada como polígono en ArcGIS Desktop?

Question

Tengo una imagen rasterizada rectangular en ArcGIS. Mi pregunta es: ¿cómo puedo producir un polígono que es sólo la extensión de la imagen de trama? Es decir, quiero tener una capa de polígono con un solo cuadrilátero que es el borde de la imagen.

He intentado simplemente hacer una conversión Raster->Polígono, pero esto intenta utilizar los datos de la imagen para generar el contorno - todo lo que quiero es el contorno del borde de la imagen.

¿Alguna idea?

Estoy utilizando ArcGIS Desktop 10.

Answer 1

Si tiene Spatial Analyst, puede hacerlo:

1. multiplicar el raster por 0 para crear un raster de valor constante utilizando la Calculadora de Raster.
2. convierta el resultado del paso 1 en polígono utilizando la herramienta Raster to Polygon.

Otra solución:

Si tu tienes Analista 3D entonces puede utilizar la licencia Herramienta de dominio rasterizado en 3D Analyst>Conversión>de raster (sin embargo, no tiene la función de suavizado)

Answer 2

¿Y si simplemente se digitaliza? Clic-clic-clic-hecho.

Answer 3

Yo uso este simple python script, crea un featureclass de polígono con la extensión de todos los raster presentes en una carpeta. Para usarlo, puede crear una herramienta en ArcToolbox o simplemente cambiar las variables InFolder y Dest (Destino).

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

Answer 4

Puede utilizar la herramienta Construir huella en la caja de herramientas de gestión de datos. Tienes que crear un conjunto de datos de mosaico a partir de tu raster.

Cree un nuevo mosaico en una geodatabase de archivos, añada su raster.

Abra su mosaico en ArcMap y extraiga la huella.

Eso es lo que yo haría.

Answer 5

Sé que ha pasado un tiempo, pero personalmente creo que reclasificar y convertir la trama en polígono puede ser waaaaay demasiado tiempo. En mi opinión, la forma más fácil de hacerlo es:

Obtener maxX, maxY, minX, minY (esto se puede hacer en Python con el comando describe) Ejecuta Create Fishnet con la extensión anterior y el número de filas/columnas como 1

Si tiene muchos rásteres, esto debe hacerse en Python y puede crear cientos de huellas en cuestión de segundos.