Il 29/05/2012 20:06, Bborie Park ha scritto: > WITH ref AS ( > SELECT ST_MakeEmptyRaster(1, 1, 0, 0, 1, -1, 0, 0, 0) AS rast > ) > SELECT > ST_AsRaster(g.geom, ref.rast, '32BUI', g.id, 0) > FROM mygeomtable g > CROSS JOIN ref > > That would convert your geometries into rasters, all of which have the > same alignment with pixel values set to the particular geometry's id. > So, for 20 input geometries, there would be 20 output rasters.
but: WITH ref AS ( SELECT ST_MakeEmptyRaster(1, 1, 0, 0, 1, -1, 0, 0, 3003) AS rast ) SELECT ST_AsRaster(g.geom, ref.rast, '32BUI', g.id, 0) FROM province g CROSS JOIN ref; ERROR: rt_raster_gdal_rasterize: Unable to add band to GDALDataset CONTEXT: PL/pgSQL function "st_asraster" line 26 at RETURN All the best. -- Paolo Cavallini - Faunalia www.faunalia.eu Full contact details at www.faunalia.eu/pc _______________________________________________ postgis-users mailing list postgis-users@postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users