Re: [gdal-dev] Raster/vector intersections in PostGIS WKT Raster

Ivan Lucena ivan.lucena at pmldnet.com
Tue Jun 29 14:30:39 EDT 2010


Christian,

How did you added support for Oracle GeoRaster on GeoTools/GeoServer? Did you use the Oracle GeoRaster Java API or GDAL? Is that available?

I am looking forward to use it.

Regards,

Ivan


>  -------Original Message-------
>  From: christian.mueller at nvoe.at
>  To: Pierre Racine <Pierre.Racine at sbf.ulaval.ca>
>  Cc: gdal-dev at lists.osgeo.org <gdal-dev at lists.osgeo.org>
>  Subject: Re: [gdal-dev] Raster/vector intersections in PostGIS WKT Raster
>  Sent: Jun 29 '10 11:46
>  
>  Hi, superb job. Congratulations.
>  
>  I am a geotools/geoserver developer and one of my jobs is fetching  
>  image tiles  or the corresponding pyramid tiles from jdbc data  
>  sources. Some weeks ago, I added support for Oracle Georaster and I  
>  want to support WKT Georaster too.
>  
>  wktraster makes it quite simple to achieve this. Is there any database  
>  layout you recommend. (A table for the image and one for each pyramid,  
>  or put all the tiles into one table having an attribute "level",.....).
>  
>  I have seen "gdal2wktraster.py", is this the recommended way to import  
>  tiles into postgis ?
>  What is the table layout created/needed by this utility ?
>  
>  Cheers
>  Christian
>  
>  
>  
>  Quoting Pierre Racine <Pierre.Racine at sbf.ulaval.ca>:
>  
>  > Hi everybody,
>  >
>  > I'm pleased to announce that you can now do intersections between  
>  > rasters and geometries very much like you used to do  
>  > geometry/geometry intersections in PostGIS. For this, PostGIS WKT  
>  > Raster introduces two new functions to PostGIS:  
>  > ST_Intersects(geometry, raster, band) and ST_Intersection(geometry,  
>  > raster, band).
>  >
>  > As its geometry/geometry sister, ST_Intersects(geometry, raster,  
>  > band) returns TRUE if the withvalue area of a raster or a raster  
>  > tile (nodata value are ignored) intersects a geometry and  
>  > ST_Intersection(geometry, raster, band) returns the geometry/value  
>  > set of geometries representing the intersection between the geometry  
>  >  and each polygonized group of pixel sharing a same value from the  
>  > raster and its associated value. Example:
>  >
>  > SELECT polyid,
>  >        (ST_Intersection(the_geom, rast)).geom,
>  >        (ST_Intersection(the_geom, rast)).val
>  > FROM my_polygons, srtm_tiled
>  > WHERE ST_Intersects(rast, the_geom)
>  >
>  > or, a bit more complex, but much faster and returning the same result:
>  >
>  > SELECT polyid, (gv).geom, (gv).val
>  > FROM (SELECT polyid, ST_Intersection(the_geom, rast, 1) AS gv
>  >      FROM my_polygons, srtm_tiled
>  >      WHERE ST_Intersects(rast, the_geom)) foo
>  >
>  > These functions works with any kind of geometry so you can intersect  
>  >  any road, river, polygons or point layer with your favorite  
>  > elevation or land cover raster of any resolution and any size. You  
>  > can also load any number of raster in the database with one command  
>  > line to constitute a unique table raster coverage of any shape (no  
>  > matter if your group of raster do not form a rectangular coverage).
>  >
>  > You will find a complete tutorial on how to use those two new  
>  > functions here:  
>  > http://trac.osgeo.org/postgis/wiki/WKTRasterTutorial01. A reference  
>  > is available here:  
>  > http://postgis.refractions.net/documentation/manual-svn/RT_reference.html  
>  > and chapter 13 of "PostGIS in Action" is also a very good  
>  > introduction to WKT Raster.
>  >
>  > Compilation and installation info are available in the project home  
>  > page: http://trac.osgeo.org/postgis/wiki/WKTRaster
>  >
>  > These features are the result of two years of collaboration between  
>  > many people and companies interested in bringing raster into  
>  > PostGIS. I would like to thanks particularly Steve Cumming who  
>  > initially made all this possible, Sandro Santilli who wrote the base  
>  >  code, Mateusz Loskot who wrote the Python loader, Jorge Arevalo who  
>  >  wrote the GDAL driver and ST_DumpAsPolygons(), Regina and Leo Obe  
>  > who are doing so many things and believed in the project from the  
>  > beginning by adding a chapter about WKT Raster in "PostGIS in  
>  > Action" and also David Zwarg who wrote most of the setter functions.
>  >
>  > This is the first step toward a first complete raster/vector  
>  > analysis and manipulation SQL API. We hope to make of PostGIS with  
>  > WKT Raster the most powerful and complete GIS analysis and  
>  > manipulation high level language ever (!!!). WKT Raster should  
>  > normally be totally integrated into PostGIS in PostGIS release 2.0.  
>  > Upcoming functions should include ST_Reclass, ST_Clip, ST_AsRaster,  
>  > ST_Resample.
>  >
>  > Let us know your need and your experience with WKT Raster.
>  >
>  > Thanks,
>  >
>  > Pierre Racine
>  >
>  > _______________________________________________
>  > gdal-dev mailing list
>  > gdal-dev at lists.osgeo.org
>  > http://lists.osgeo.org/mailman/listinfo/gdal-dev
>  >
>  
>  
>  
>  ----------------------------------------------------------------
>  This message was sent using IMP, the Internet Messaging Program.
>  
>  
>  _______________________________________________
>  gdal-dev mailing list
>  gdal-dev at lists.osgeo.org
>  http://lists.osgeo.org/mailman/listinfo/gdal-dev
>  


More information about the gdal-dev mailing list