J'ai commencé avec un fichier raster NetCDF .gri et .grd du Royaume-Uni fourni par un collègue. Je l'ai découpé dans R pour être uniquement Londres, exporté et converti en fichier ASC, puis importé dans PostGIS à l'aide des commandes suivantes dans R:
library(raster)
uk_raster <- raster("AnnMean2011.grd")
london_area <- extent(-720000.0,-630000.0,-50000.0,25000)
london_raster <- crop(uk_raster, london_area)
writeRaster(london_raster, filename="AnnMean2011.asc", format="ascii")
Et puis sur la ligne de commande Ubuntu:
raster2pgsql -I -C -s 10001 -t 20x20 AnnMean2011.asc annualmean | psql -d james_traffic
J'ai maintenant une table raster dans PostGIS. Le SRID de 10001 est d'ailleurs le suivant, fourni par un collègue:
INSERT INTO spatial_ref_sys(srid, auth_name, auth_srid, proj4text)
VALUES (10001,'CMAQ_Urban',10001,'+proj=lcc +a=6370000 +b=6370000 +lat_1=35 +lat_2=65 +lat_0=52 +lon_0=10 +x_0=000000 +y_0=00000');
Dans la même base de données, j'ai un fichier polygonal, SRID 27700, qui couvre Londres. Je voudrais calculer la valeur moyenne à l'intérieur de chaque polygone, à partir du raster.
J'essaie quelque chose comme ça, mais ce n'est pas bien:
select polygons.postcode, avg(st_value(joined_data.rast))
from (
select (ST_Intersection(raster.rast, 1, polygons.geom)).*
from raster, polygons
where ST_Intersects(raster.rast, 1, polygons.geom)
) joined_data
group by polygons.postcode
Comment pourrais-je m'y prendre s'il vous plaît?
Il semble que le polygone et le raster soient correctement alignés, je dois les convertir tous les deux en WGS84 je pense.