¿Obtener una imagen rasterizada como matriz en Python con ArcGIS Desktop?

Question

Cuando empecé a trabajar con Python y ArcGIS 9.3, supuse que habría una forma sencilla de introducir una imagen rasterizada en una matriz de Python para poder manipularla antes de volver a almacenarla como otra imagen rasterizada. Sin embargo, parece que no puedo encontrar cómo hacer esto.

Si es posible, ¿cómo?

Answer 1

http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/RasterToNumPyArray/000v0000012z000000/

ArcGIS 10 tiene la capacidad de escribir y leer matrices numPy.

Answer 2

No creo que esto sea posible con ArcGIS <= 9.3.1

Utilizo el código abierto API GDAL para tareas como ésta.

Answer 3

¿Acceso a ArcObjects desde Python? discute la integración de arcobjects con python.

Tal vez el código en esta muestra podría ser adaptado para que pueda ser llamado desde python.

No estoy seguro de si hay una forma de pasar un array byref a python. Si la hay, entonces IPixelBlock.PixelDatabyRef valdría la pena intentarlo.

Answer 4

No estoy seguro de que puedas manipular el raster píxel por píxel, pero puedes usar los objetos de geoprocesamiento en conjunción con la API de python.

Puedes usar cualquier caja de herramientas para ese tipo de manipulación. Un ejemplo de script sería:

#import arcgisscripting

gp = arcgisscripting.create(9.3)

gp.AddToolbox("SA") # addint spatial analyst toolbox

rasterA = @"C:\rasterA.tif"
rasterB = @"C:\rasterB.tif"

rasterC = @"C:\rasterC.tif" # this raster does not yet exist
rasterD = @"C:\rasterD.tif" # this raster does not yet exist

gp.Minus_SA(rasterA,rasterB,rasterC)

gp.Times_SA(rasterA,rasterB,rasterD)

# lets try to use more complex functions

# lets build and expression first

expression1 = "slope( " + rasterC + ")"
expression2 = "(" + rasterC " + " rasterD + ") - " + rasterA 

gp.SingleOutputMapAlgebra_SA(expression1,@"C:\result_exp1.tif")
gp.SingleOutputMapAlgebra_SA(expression2,@"C:\result_exp2.tif")

A continuación se presenta un seguimiento de su pregunta . Todavía no es posible. No estoy seguro de la versión 10.0.

Answer 5

La forma más fácil sería convertir el raster a netCDF y luego abrirlo y recorrer la cuadrícula. Hice prácticamente lo mismo para un proyecto que implicaba convertir los rásteres en datos de características basados en los datos asignados a las celdas del ráster. Estuve mirando esto durante años, y llegué a la conclusión de que recorrer los datos de la cuadrícula sería más fácil desde netCDF.