[postgis-users] PostGIS Raster Decimal?

[Andreas Forø Tollefsen]

Hi Pierre,
I got memory issues as well.
I tried with this query:
SELECT ST_MapAlgebra(rast, 'rast*10') INTO algtest FROM gpw90;


WARNING:  Pixel value for 32BF band got truncated from 47519.1 to 47519.1
WARNING:  Pixel value for 32BF band got truncated from 185494 to 185494
WARNING:  Pixel value for 32BF band got truncated from 28485.1 to 28485.1
WARNING:  Pixel value for 32BF band got truncated from 10927.1 to 10927.1


ERROR:  out of memory
DETAIL:  Failed on request of size 19.

********** Error **********

ERROR: out of memory
SQL state: 53200
Detail: Failed on request of size 19.

2011/3/14 Paragon Corporation <lr at pcorp.us>

>  Pierre,
>
> Actually the single raster version is documented.  Though I did run into
> memory issues -- so still think we need some work on that.
>
> http://www.postgis.org/documentation/manual-svn/RT_ST_MapAlgebra.html
>
> Regina
> http://www.postgis.us
>
>
>
>  ------------------------------
> *From:* postgis-users-bounces at postgis.refractions.net [mailto:
> postgis-users-bounces at postgis.refractions.net] *On Behalf Of *Andreas Forø
> Tollefsen
> *Sent:* Monday, March 14, 2011 2:03 PM
> *To:* PostGIS Users Discussion
> *Cc:* Pierre Racine
>
> *Subject:* Re: [postgis-users] PostGIS Raster Decimal?
>
> Hi Pierre,
>
> I will check it out tomorrow and give some feedback.
>
> Cheers,
> Andreas
>
> 2011/3/14 Pierre Racine <Pierre.Racine at sbf.ulaval.ca>
>
>>  I’m very aware of the problem. Will try to fix it the best we can.
>>
>>
>>
>> Why not testing ST_MapAlgebra and give us your feedback? It is still to
>> document but the specifications are in this page:
>>
>>
>>
>> http://trac.osgeo.org/postgis/wiki/WKTRaster/SpecificationWorking02
>>
>>
>>
>> Pierre
>>
>>
>>
>> *From:* postgis-users-bounces at postgis.refractions.net [mailto:
>> postgis-users-bounces at postgis.refractions.net] *On Behalf Of *Andreas
>> Forø Tollefsen
>> *Sent:* 14 mars 2011 11:24
>> *To:* PostGIS Users Discussion
>> *Subject:* Re: [postgis-users] PostGIS Raster Decimal?
>>
>>
>>
>> Hi Pierre.
>>
>>
>>
>> Yes. I think that would be an appropriate temporary fix.
>>
>> 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 :)
>>
>>
>>
>> Good luck.
>>
>>
>>
>> Best,
>>
>> Andreas
>>
>>
>>
>> 2011/3/14 Pierre Racine <Pierre.Racine at sbf.ulaval.ca>
>>
>> Andreas,
>>
>> This is a well known problem. See ticket
>> http://trac.osgeo.org/postgis/ticket/650
>>
>> I will discuss this problem with Frank this week at the code sprint in
>> Montreal.
>>
>> 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).
>>
>> Pierre
>>
>> >-----Original Message-----
>> >From: Pierre Racine
>> >Sent: 14 mars 2011 10:55
>> >To: Pierre Racine
>> >Subject: FW: [postgis-users] PostGIS Raster Decimal?
>> >
>> >
>> >
>> >From: postgis-users-bounces at postgis.refractions.net [mailto:
>> postgis-users-
>> >bounces at postgis.refractions.net] On Behalf Of Andreas Forø Tollefsen
>> >Sent: 14 mars 2011 10:53
>> >To: PostGIS Users Discussion
>> >Subject: [postgis-users] PostGIS Raster Decimal?
>>
>> >
>> >Hi. Again.
>> >
>> >I have been working on some raster data in PostGIS lately.
>> >However, i have some issues with the raster values become integer after
>> importing the sql to postgis.
>> >For instance this procedure:
>> >
>> >C:\prio_grid\source\gpw>c:\python26\python
>> c:\prio_grid\script\raster2pgsql.py -
>> >r c:\prio_grid\source\gpw\lrc30p90\glp90ag30\w001001.adf -t gpw90 -s 4326
>> -o gpw
>> >90.sql -I -M
>> >
>> >Then to database:
>> >C:\prio_grid\source\gpw>psql -h 192.168.1.55 -d priogrid -f gpw90.sql
>> >
>> >Querying this data:
>> >SELECT gpw90.gid, ((gpw90.gpw90val).val) AS gpw90
>> >INTO popgrid
>> >FROM (SELECT priogrid_land.gid, ST_Intersection(gpw90.rast,
>> priogrid_land.centroid) AS gpw90val FROM
>> >gpw90, priogrid_land WHERE ST_Intersects(priogrid_land.centroid,
>> gpw90.rast)) AS gpw90
>> >WHERE gpw90.gid = 139303
>> >GROUP BY gid,((gpw90.gpw90val))
>> >;
>> >
>> >Gives:
>> >gid; gpw90
>> >139303;39849
>> >
>> >The value in the original raster is:
>> >39849.2
>> >
>> >Question is then. How can i ensure that raster remain decimal and not
>> integer after this import and
>> >query process?
>> >
>> >Thanks.
>> >Andreas
>> _______________________________________________
>> postgis-users mailing list
>> postgis-users at postgis.refractions.net
>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>>
>>
>>
>> _______________________________________________
>> postgis-users mailing list
>> postgis-users at postgis.refractions.net
>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>>
>>
>
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20110315/4a05bce0/attachment.html>