cómo superponer shapefile y de trama?

Question

Tengo un shapefile con polígonos. Y tengo un global de archivo raster. Quiero superposición de la shapefile de polígonos en la trama de cuadrícula y calcular la media de trama valor para cada polígono.

¿Cómo puedo hacer esto con GDAL, la escritura de los resultados para el archivo de forma?

Answer 1

En R se puede hacer

library(raster)
library(rgdal)
r <- raster('raster_filename')
p <- readOGR('shp_path', 'shp_file')
e <- extract(r, p, fun=mean)

e es un vector con la media de los valores de celda ráster para cada polígono.

Answer 2

Siguiendo los consejos que me dieron en el gdal-dev lista de correo, he utilizado StarSpan:

starspan --vector V --raster R1 R2 ... --stats mystats.csv avg mode

Los resultados se guardan en formato CSV. En ese momento, que ya era suficiente para mí, pero debería ser posible de alguna manera a forjar un Shapefile de esa información.

Answer 3

Carga tu archivo shapefile y su trama en PostGIS 2.0 y hacer:

SELECT (ST_SummaryStats(ST_Clip(rast, geom))).*
FROM rastertable, geomtable

Answer 4

No creo que GDAL es la mejor herramienta para ello, pero se puede usar gdal_rasterize "limpiar" todos los valores fuera del polígono.

Algo así como:

gdal_translate -a_nodata 0 original.tif work.tif
gdal_rasterize -burn 0 -b 1 -i work.tif yourpolygon.shp -l yourpolygon
gdalinfo -stats work.tif
rm work.tif

El gdal_rasterize programa modifica el archivo, así que hacer una copia para trabajar. También nos marca a un determinado valor (cero en este caso) nodata. La "quema 0 -b 1" significa quemar un valor de cero en la banda 1 del archivo de destino (en el trabajo.tif). La "-i" significa invertir la rasterización de modo de quema de valores fuera de la poligonal en lugar de dentro de ella. El gdalinfo comando con -estadísticas de los informes sobre la banda de estadísticas. Creo que va a excluir el valor nodata (que hemos marcado anteriormente con -a_nodata).

Answer 5

Transformar la forma de archivo de mapa de bits por gdal_rasterize y el uso del código en http://www.spatial-ecology.net/dokuwiki/doku.php?id=wiki:geo_tools para calcular zonal estadística para cada polígono. Puede ejecutar http://km.fao.org/OFwiki/index.php/Oft-reclass si desea obtener una tif con su rásteres de estadística. Disfrutar el código Ciao Giuseppe