<div dir="ltr"><div>You are right! after a long list of messages, it finally get a map that looks alright.<br></div><div>The two rasters do intersect because they come from the same file:<br></div><div><br></div><div>In case you'd like to see it here is a small example (1.3 MB):<br><a href="https://dl.dropboxusercontent.com/u/1571018/prova.tif" target="_blank">https://dl.dropboxusercontent.com/u/1571018/prova.tif</a><br><br></div><div>-- uploaded to postgress as:<br></div><div>raster2pgsql -C -s 4326 -t 10x10 -I -M -Y prova.tif rastertmp.prova | psql -d bd_cide<br><br></div><div>-- create the 1/null map (for positive/negative values)<br></div><div>CREATE TABLE rastertmp.provapositiu (rid serial, rast raster);<br>INSERT INTO rastertmp.provapositiu<br>    SELECT rid, ST_RECLASS (ST_Band(rast, 1), 1, '[1-10000]:1', '2BUI', -9) <br>    FROM rastertmp.prova;<br></div><br>-- looking OK<br><br>CREATE TABLE rastertmp.provanet (rid serial, rast raster);<br>INSERT INTO rastertmp.provanet<br>    SELECT t1.rid, ST_Union(ST_Intersection(t1.rast, 1, t2.rast, 1, 'BAND1'))<br>    FROM rastertmp.prova AS t1, rastertmp.provapositiu AS t2<br>    WHERE ST_Intersects(t1.rast, t2.rast)<br>    GROUP BY t1.rid;<br><div><br>NOTICE:  The two rasters provided have no intersection.  Returning no band raster<br>CONTEXT:  PL/pgSQL function st_intersection(raster,integer,raster,integer,text,double precision[]) line 20 at assignment<br>NOTICE:  Could not find raster band of index 1 when setting pixel value. Nodata value not set. Returning original raster<br>CONTEXT:  PL/pgSQL function st_intersection(raster,integer,raster,integer,text,double precision[]) line 21 at assignment<br>NOTICE:  The two rasters provided have no intersection.  Returning no band raster<br>CONTEXT:  PL/pgSQL function st_intersection(raster,integer,raster,integer,text,double precision[]) line 20 at assignment<br>NOTICE:  Could not find raster band of index 1 when setting pixel value. Nodata value not set. Returning original raster<br>CONTEXT:  PL/pgSQL function st_intersection(raster,integer,raster,integer,text,double precision[]) line 21 at assignment<br></div><div><br>etc... a long list of these messages ...<br></div><div><br></div><div>But at the end it produces the map that seems OK.<br></div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr">Juli<div><span style="color:rgb(51,51,153)">--</span><br><font color="#333399"><b>CIDE, CSIC</b>  |  <a href="http://www.uv.es/jgpausas" target="_blank">www.uv.es/jgpausas</a>  | <br><br></font></div></div></div></div></div></div>
<br><div class="gmail_quote">On Mon, Jul 6, 2015 at 4:12 PM, Pierre Racine <span dir="ltr"><<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">These are not error. Only notices. Are you sure the query returned nothing?<br>
<br>
Otherwise I would restart with something more simple to double check some prerequired things.<br>
<br>
1) Make sure your two raster table actually contains pixels with values<br>
<br>
2) Make sure some tiles actually intersects<br>
<br>
Do you have a way to display the footprints of your rasters tiles?<br>
<span class="im HOEnZb"><br>
Pierre<br>
<br>
> -----Original Message-----<br>
> From: juli g. pausas [mailto:<a href="mailto:juli.g.pausas@uv.es">juli.g.pausas@uv.es</a>]<br>
</span><span class="im HOEnZb">> Sent: Monday, July 06, 2015 9:58 AM<br>
> To: Pierre Racine<br>
> Cc: PostGIS Users Discussion<br>
> Subject: Re: [postgis-users] raster, stats conditioned to a set of values<br>
><br>
><br>
</span><div class="HOEnZb"><div class="h5">> INSERT INTO rastertmp.provanet<br>
>     SELECT t1.rid, ST_Union(ST_Intersection(t1.rast, 1, t2.rast, 1, 'BAND1'))<br>
>     FROM rastertmp.prova AS t1, rastertmp.provapositiu AS t2<br>
>     WHERE ST_Intersects(t1.rast, t2.rast)<br>
>     GROUP BY t1.rid;<br>
><br>
><br>
> NOTICE:  The two rasters provided have no intersection.  Returning no band<br>
> raster<br>
> CONTEXT:  PL/pgSQL function<br>
> st_intersection(raster,integer,raster,integer,text,double precision[]) line 20<br>
> at assignment<br>
> NOTICE:  Could not find raster band of index 1 when setting pixel value.<br>
> Nodata value not set. Returning original raster<br>
><br>
><br>
><br>
><br>
> Juli<br>
> --<br>
> CIDE, CSIC  |  <a href="http://www.uv.es/jgpausas" rel="noreferrer" target="_blank">www.uv.es/jgpausas</a>  |<br>
><br>
><br>
><br>
> On Mon, Jul 6, 2015 at 3:52 PM, Pierre Racine<br>
> <<a href="mailto:Pierre.Racine@sbf.ulaval.ca">Pierre.Racine@sbf.ulaval.ca</a>> wrote:<br>
><br>
><br>
>       > Adding a GROUP BY t1.rid, doesn't solve the problem, still get an<br>
> error.<br>
><br>
>       Which error?<br>
><br>
><br>
<br>
</div></div></blockquote></div><br></div>