Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-25 Thread JamesH
"However, if I wanted to create a view or a new dataset which has the raster clipped by the polygon's geometry in addition to the stats, can I use just ST_Clip?" I think I understand my problem better now. What I ideally want is the composite raster (i.e. all three bands as one RGB raster) to be

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-25 Thread JamesH
>http://trac.osgeo.org/postgis/browser/trunk/raster/scripts/plpgsql/st_summarystatsagg.sql> >I copied the code and ran it, returned as successful but when I run the full query it says that function >ST_SummaryStatsagg(raster) does not exist. >Have I copied the function incorrect> This was my mista

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread JamesH
"What is the type of the geometries you are intersecting with the raster?" I am intersecting a shapefile of polygons. -If they are polygons, what proportion of your raster does all your geometry layer represent? They do cover a large proportion of the raster. ST_SummaryStatsAgg is not installed

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread Pierre Racine
: postgis-users@postgis.refractions.net > Subject: Re: [postgis-users] Geom/Rast Intersection St9bad_alloc > > "I am attempting to intersect a vector dataset of building polygons with a > 3-Band Thermal tiff to produce a dataset that shows the thermal readings for > each polygon.

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread JamesH
"I am attempting to intersect a vector dataset of building polygons with a 3-Band Thermal tiff to produce a dataset that shows the thermal readings for each polygon." If I took the 3 band jpeg of this raster and used ST_AsTiff to convert it to a single band tiff would this calculation be valid? Ho

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread JamesH
I am attempting to intersect a vector dataset of building polygons with a 3-Band Thermal tiff to produce a dataset that shows the thermal readings for each polygon. "What were you trying to achieve ? Did you read the reference manual ?" I hoped from the example in the reference "-- Return bands 2

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread Sandro Santilli
On Tue, Apr 24, 2012 at 03:51:27AM -0700, JamesH wrote: > Have tried ST_Intersection(geom, rast, ST_Band(rast, ARRAY[1,2,3])) in query What were you trying to achieve ? Did you read the reference manual ? --strk; ,--o-. | __/ |Delivering high quality PostGIS 2.0 ! | / 2.0 |

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread JamesH
"ST_Band ?" Thank you. Have tried ST_Intersection(geom, rast, ST_Band(rast, ARRAY[1,2,3])) in query above Returns: function st_intersection(geometry, raster, raster) does not exist. Have also attempted ST_Intersection(geom, rast, ST_Band(rast, '1,2,3')) but this returns the same error. Basic sy

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread Sandro Santilli
On Tue, Apr 24, 2012 at 03:06:15AM -0700, JamesH wrote: > I've checked and seen that my raster is actually three bands when I thought > it was just single band. > Is there a way to alter the query to account for this? ST_Band ? http://www.postgis.org/documentation/manual-svn/RT_ST_Band.html --s

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread JamesH
"St9Bad_alloc means the process run out of memory. There may be a memory leak or uncareful use in ST_Intersection(raster,geometry)." Thought it could be something like that. "Try restricting the scope of your queries, using subsets from clip_build and nclheat and storing the results in a temporar

Re: [postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread Sandro Santilli
On Tue, Apr 24, 2012 at 02:50:37AM -0700, JamesH wrote: > Hi all, > > I am attempting to intersect a thermal raster TIFF with a MasterMap > shapefile. > > Both sets of data are stored in my database successfully and the query I am > running is: > SELECT gid, ST_Intersection(geom, rast, 1) > FROM

[postgis-users] Geom/Rast Intersection St9bad_alloc

2012-04-24 Thread JamesH
Hi all, I am attempting to intersect a thermal raster TIFF with a MasterMap shapefile. Both sets of data are stored in my database successfully and the query I am running is: SELECT gid, ST_Intersection(geom, rast, 1) FROM clip_build, nclheat WHERE ST_Intersects(rast, geom) The raster file is qu