3 votos

Python lazo Polígono a trama no funciona

Quiero convertir un gran número de diferentes shp de polígonos a ráster en ArcGIS 10.1, el uso de un campo existente columna "clase" para la nueva trama clases. No estoy muy fluido en python, pero pensé que esto iba a funcionar, y que no:


 # Import arcpy module
 import arcpy
... # Local variables:
... project_input = "D:\\RESEARCH_PROJECTS\\NORD-STAR\\AMAZONIA\\1-RAW_DATA\\INPE\\PRODES DIGITAL\\2004\\"
... raster_output = "D:\\RESEARCH_PROJECTS\\NORD-STAR\\AMAZONIA\\1-RAW_DATA\\INPE\\PRODES DIGITAL\\2004\\convert\\"
... # List features
... arcpy.env.workspace = project_input
... fcList = arcpy.ListFeatureClasses() 
... # Loop
... for featureClass in fcList:
... # Output
...     raster_outputpath = raster_output + featureClass
... 
... # Process: Polygon to Raster
...     arcpy.PolygonToRaster_conversion(featureClass, "CLASS", raster_output, "CELL_CENTER", "NONE", "30")
...     print "finished polygon to raster" 
...     
Runtime error  Traceback (most recent call last):   File "", line 20, in    File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\conversion.py", line 2436, in PolygonToRaster     raise e ExecuteError: ERROR 999999: Error executing function. Failed to execute (PolygonToRaster). 

El script se llama realmente el primer polígono, pero que es esto:

Ejecución: PolygonToRaster "D:\RESEARCH_PROJECTS\NORD-STAR\AMAZONIA\1-RAW_DATA\INPE\PRODES DIGITAL\2004\PDigital2004_00157_pol.shp" CLASE "D:\RESEARCH_PROJECTS\NORD-STAR\AMAZONIA\1-RAW_DATA\INPE\PRODES DIGITAL\2004\convertir\" CELL_CENTER NINGUNO 30 Hora de Inicio: Wed Nov 21 15:17:11 2012 de ERROR 999999: Error al ejecutar la función. No pudo ejecutar (PolygonToRaster). Error en el Mié Nov 21 15:17:11 2012 (Transcurrido Tiempo: 0,00 segundos)

Gracias por su ayuda

Javier

6voto

MattB Puntos 76

El problema es que el código define raster_outputpath pero intenta guardar el raster en raster_output.

Cambiar esta línea:

arcpy.PolygonToRaster_conversion(featureClass, "CLASS", raster_output, "CELL_CENTER", "NONE", "30")

A:

arcpy.PolygonToRaster_conversion(featureClass, "CLASS", raster_outputpath, "CELL_CENTER", "NONE", "30")

3voto

Aaron Puntos 25882

Pruebe las siguientes sugerencias:

Las líneas 4 y 5 archivos han de personajes tales como "-" y "de", que puede conducir a errores. He reorganizado el script, espero que estos cambios ayudarán.

# Import arcpy module
import arcpy
from arcpy import env

#Set working environment
env.workspace = "D:\\RESEARCH_PROJECTS\\NORD-STAR\\AMAZONIA\\1-RAW_DATA\\INPE\\PRODES DIGITAL\\2004"
Dir = env.workspace

#List FCs
fcList = arcpy.ListFeatureClasses() 

# Loop
for fc in fcList:

     output = Dir + "\\" + fc + "_ras"

     # Process: Polygon to Raster
     arcpy.PolygonToRaster_conversion(fc, "CLASS", output, "CELL_CENTER", "", "30")

print "finished polygon to raster" 

1voto

Prachur Puntos 111

EDITAR

Me corrió una prueba utilizando el siguiente en un polígono que representa la CA de límite:

 arcpy.PolygonToRaster_conversion(featureClass, "CLASS", raster_output, "CELL_CENTER", "NONE", 30)

Funcionó bien, por lo que el problema puede estar en el polígono o tu estructura de directorios. Me gustaría tratar de ejecutar la Reparación de la Geometría en el polígono de las características y asegúrese de que los directorios no tienen espacios en ellos (parece una tontería, pero tener espacios en ALGUNOS directorios, puede causar problemas). Además, si usted está funcionando esta en ArcMap como opuesto a un IDE, deshabilitar el procesamiento en segundo plano.

Respuesta inicial:

No pude encontrar nada sobre la forma correcta de los valores de entrada para el priority_filed parámetro (el quinto parámetro). Trate de reemplazar "NINGUNO" con vacío entre comillas dobles ( "" ), ya que es un parámetro opcional y puede así por defecto si usted no tiene un uso específico para ello:

arcpy.PolygonToRaster_conversion(featureClass, "CLASS", raster_output, "CELL_CENTER", "", "30")

1voto

Eric G Puntos 131

Cuando miro ESRI ayuda en línea parece arcpy.PolygonToRaster_conversion() espera un número para el cell_size y parece que está pasando una cadena. Probar esto y ver si funciona. arcpy.PolygonToRaster_conversion(featureClass, "CLASE",raster_output, "CELL_CENTER", "NINGUNO", 30)

0voto

Amir Rachum Puntos 118

Yo tuve un problema similar con la misma herramienta.

A menudo utilizo una .formato tif, como la salida y que resuelve los problemas que surgen con el formato de CUADRÍCULA.

JJH

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X