<div dir="ltr">Check out the last example (Multi-band versions) of ST_AddBand().<div><br></div><div><a href="http://postgis.refractions.net/documentation/manual-svn/RT_ST_AddBand.html">http://postgis.refractions.net/documentation/manual-svn/RT_ST_AddBand.html</a><br>
</div><div><br></div><div style>-bborie</div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Sun, Jul 7, 2013 at 8:43 PM, Mark Wynter <span dir="ltr"><<a href="mailto:mark@dimensionaledge.com" target="_blank">mark@dimensionaledge.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I'm hoping someone can suggest a way of efficiently combining 3 individual (but concordant) single band rasters into single raster having 3 bands?  We're working with PostgreSQL9.1, PostGIS2.1Beta.<br>

<br>
The starting point is a single tiled raster with 3 bands representing RGB values.<br>
<br>
Let's say we perform some color smoothing on each of the RGB bands using mapalgebra.<br>
<br>
CREATE TABLE resample_myrastertable_b1 AS<br>
SELECT rid, st_mapalgebrafctngb(rast, 1, '8BUI', 5, 5, 'st_mean4ma(float[][],text,text[])'::regprocedure, 'ignore', NULL) as rast FROM original_rastertable;<br>
<br>
CREATE TABLE resample_myrastertable_b2 AS<br>
SELECT rid, st_mapalgebrafctngb(rast, 2, '8BUI', 5, 5, 'st_mean4ma(float[][],text,text[])'::regprocedure, 'ignore', NULL) as rast FROM original_rastertable;<br>
<br>
CREATE TABLE resample_myrastertable_b3 AS<br>
SELECT rid, st_mapalgebrafctngb(rast, 3, '8BUI', 5, 5, 'st_mean4ma(float[][],text,text[])'::regprocedure, 'ignore', NULL) as rast FROM original_rastertable;<br>
<br>
This gives us 3 single-band results rasters.<br>
<br>
How do we now join the individual smoothed RGB "results" rasters into a single multiband "results" raster?   We'd also like to add a 4th band to the "results" raster based on say a calculation which references the smoothed values in bands 1-3.<br>

<br>
I'm getting reasonably comfortable writing MapAlgebra callback expressions.   Its more a case of how do we efficiently handle table operations (e.g. appending single band raster outputs) each time we perform some MapAlgebra?<br>

<br>
Any suggestions and or a simple worked example would be most appreciated.<br>
<br>
Many thanks<br>
<br>
Mark<br>
<br>
<br>
<br>
<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org">postgis-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users</a><br>
</blockquote></div><br></div>