<div>Hi Pierre, </div>I got memory issues as well.<div>I tried with this query:</div><div><div>SELECT ST_MapAlgebra(rast, 'rast*10') INTO algtest FROM gpw90;</div><div><br></div><div><br></div><div><div>WARNING: Pixel value for 32BF band got truncated from 47519.1 to 47519.1</div>
<div>WARNING: Pixel value for 32BF band got truncated from 185494 to 185494</div><div>WARNING: Pixel value for 32BF band got truncated from 28485.1 to 28485.1</div><div>WARNING: Pixel value for 32BF band got truncated from 10927.1 to 10927.1</div>
<div><br></div><div><br></div><div>ERROR: out of memory</div><div>DETAIL: Failed on request of size 19.</div><div><br></div><div>********** Error **********</div><div><br></div><div>ERROR: out of memory</div><div>SQL state: 53200</div>
<div>Detail: Failed on request of size 19.</div><div><br></div><div class="gmail_quote">2011/3/14 Paragon Corporation <span dir="ltr"><<a href="mailto:lr@pcorp.us">lr@pcorp.us</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
<div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial">Pierre,</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial">Actually the single raster version is documented. Though
I did run into memory issues -- so still think we need some work on
that.</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial"><a href="http://www.postgis.org/documentation/manual-svn/RT_ST_MapAlgebra.html" target="_blank">http://www.postgis.org/documentation/manual-svn/RT_ST_MapAlgebra.html</a></font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial">Regina</font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial"><a href="http://www.postgis.us" target="_blank">http://www.postgis.us</a></font></span></div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial"></font></span> </div>
<div dir="ltr" align="left"><span><font color="#0000ff" size="2" face="Arial"></font></span> </div><font color="#0000ff" size="2" face="Arial"></font><font color="#0000ff" size="2" face="Arial"></font><font color="#0000ff" size="2" face="Arial"></font><br>
<div dir="ltr" lang="en-us" align="left">
<hr>
<font size="2" face="Tahoma"><div class="im"><b>From:</b>
<a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a>
[mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a>] <b>On Behalf Of
</b>Andreas Forø Tollefsen<br></div><b>Sent:</b> Monday, March 14, 2011 2:03
PM<br><b>To:</b> PostGIS Users Discussion<br><b>Cc:</b> Pierre
Racine<div><div></div><div class="h5"><br><b>Subject:</b> Re: [postgis-users] PostGIS Raster
Decimal?<br></div></div></font><br></div><div><div></div><div class="h5">
<div></div>Hi Pierre,
<div><br></div>
<div>I will check it out tomorrow and give some feedback.</div>
<div><br></div>
<div>Cheers,</div>
<div>Andreas<br><br>
<div class="gmail_quote">2011/3/14 Pierre Racine <span dir="ltr"><<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>></span><br>
<blockquote style="border-left:#ccc 1px solid;margin:0px 0px 0px 0.8ex;padding-left:1ex" class="gmail_quote">
<div lang="EN-US" vlink="purple" link="blue">
<div>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt">I’m very
aware of the problem. Will try to fix it the best we can. </span></p>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt"></span> </p>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt">Why not
testing ST_MapAlgebra and give us your feedback? It is still to document but
the specifications are in this page:</span></p>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt"></span> </p>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt"><a href="http://trac.osgeo.org/postgis/wiki/WKTRaster/SpecificationWorking02" target="_blank">http://trac.osgeo.org/postgis/wiki/WKTRaster/SpecificationWorking02</a></span></p>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt"></span> </p>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt">Pierre</span></p>
<p class="MsoNormal"><span style="color:#1f497d;font-size:11pt"></span> </p>
<div style="border-bottom:medium none;border-left:blue 1.5pt solid;padding-bottom:0cm;padding-left:4pt;padding-right:0cm;border-top:medium none;border-right:medium none;padding-top:0cm">
<div>
<div style="border-bottom:medium none;border-left:medium none;padding-bottom:0cm;padding-left:0cm;padding-right:0cm;border-top:#b5c4df 1pt solid;border-right:medium none;padding-top:3pt">
<p class="MsoNormal"><b><span style="font-size:10pt">From:</span></b><span style="font-size:10pt"> <a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a> [mailto:<a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a>] <b>On Behalf
Of </b>Andreas Forø Tollefsen<br><b>Sent:</b> 14 mars 2011 11:24<br><b>To:</b>
PostGIS Users Discussion<br><b>Subject:</b> Re: [postgis-users] PostGIS Raster
Decimal?</span></p></div></div>
<div>
<div></div>
<div>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Hi Pierre.</p>
<div>
<p class="MsoNormal"> </p></div>
<div>
<p class="MsoNormal">Yes. I think that would be
an appropriate temporary fix.</p></div>
<div>
<p class="MsoNormal">Still, a bit problematic when working with climate data
which very often include decimal digits and often have weekly and monthly
rasters for 100's of years :)</p></div>
<div>
<p class="MsoNormal"> </p></div>
<div>
<p class="MsoNormal">Good luck.</p></div>
<div>
<p class="MsoNormal"> </p></div>
<div>
<p class="MsoNormal">Best,</p></div>
<div>
<p class="MsoNormal">Andreas</p>
<div>
<p class="MsoNormal"> </p>
<div>
<p class="MsoNormal">2011/3/14 Pierre Racine <<a href="mailto:Pierre.Racine@sbf.ulaval.ca" target="_blank">Pierre.Racine@sbf.ulaval.ca</a>></p>
<p class="MsoNormal">Andreas,<br><br>This is a well known problem. See ticket <a href="http://trac.osgeo.org/postgis/ticket/650" target="_blank">http://trac.osgeo.org/postgis/ticket/650</a><br><br>I will
discuss this problem with Frank this week at the code sprint in
Montreal.<br><br>For now you can use some GIS with map algebra capability to
multiply your raster by 10 before importing or test our brand new
ST_MapAlgebra function (still to
debug).<br><br>Pierre<br><br>>-----Original Message-----<br>>From:
Pierre Racine<br>>Sent: 14 mars 2011 10:55<br>>To: Pierre
Racine<br>>Subject: FW: [postgis-users] PostGIS Raster
Decimal?<br>><br>><br>><br>>From: <a href="mailto:postgis-users-bounces@postgis.refractions.net" target="_blank">postgis-users-bounces@postgis.refractions.net</a> [mailto:<a href="mailto:postgis-users-" target="_blank">postgis-users-</a><br>
><a href="mailto:bounces@postgis.refractions.net" target="_blank">bounces@postgis.refractions.net</a>] On Behalf Of Andreas Forø
Tollefsen<br>>Sent: 14 mars 2011 10:53<br>>To: PostGIS Users
Discussion<br>>Subject: [postgis-users] PostGIS Raster Decimal?</p>
<div>
<div>
<p class="MsoNormal">><br>>Hi. Again.<br>><br>>I have been working
on some raster data in PostGIS lately.<br>>However, i have some issues with
the raster values become integer after importing the sql to
postgis.<br>>For instance this
procedure:<br>><br>>C:\prio_grid\source\gpw>c:\python26\python
c:\prio_grid\script\raster2pgsql.py -<br>>r
c:\prio_grid\source\gpw\lrc30p90\glp90ag30\w001001.adf -t gpw90 -s 4326 -o
gpw<br>>90.sql -I -M<br>><br>>Then to
database:<br>>C:\prio_grid\source\gpw>psql -h 192.168.1.55 -d priogrid
-f gpw90.sql<br>><br>>Querying this data:<br>>SELECT gpw90.gid,
((gpw90.gpw90val).val) AS gpw90<br>>INTO popgrid<br>>FROM (SELECT
priogrid_land.gid, ST_Intersection(gpw90.rast, priogrid_land.centroid) AS
gpw90val FROM<br>>gpw90, priogrid_land WHERE
ST_Intersects(priogrid_land.centroid, gpw90.rast)) AS gpw90<br>>WHERE
gpw90.gid = 139303<br>>GROUP BY
gid,((gpw90.gpw90val))<br>>;<br>><br>>Gives:<br>>gid;
gpw90<br>>139303;39849<br>><br>>The value in the original raster
is:<br>>39849.2<br>><br>>Question is then. How can i ensure that
raster remain decimal and not integer after this import and<br>>query
process?<br>><br>>Thanks.<br>>Andreas<br>_______________________________________________<br>postgis-users
mailing list<br><a href="mailto:postgis-users@postgis.refractions.net" target="_blank">postgis-users@postgis.refractions.net</a><br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a></p>
</div></div></div>
<p class="MsoNormal"> </p></div></div></div></div></div></div></div><br>_______________________________________________<br>postgis-users
mailing list<br><a href="mailto:postgis-users@postgis.refractions.net" target="_blank">postgis-users@postgis.refractions.net</a><br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br></blockquote></div><br></div></div></div></div>
<br>_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a><br>
<a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_blank">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>
<br></blockquote></div><br></div></div>