[postgis-users] PostGIS Raster Decimal?

[Andreas Forø Tollefsen]

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
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/postgis-users/attachments/20110314/0c30cddb/attachment.html>