Changes between Version 28 and Version 29 of WKTRasterGSoC
- Timestamp:
- 04/02/13 15:40:41 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
WKTRasterGSoC
v28 v29 13 13 14 14 15 ~~ '''2) An optimized version of the two rasters ST_MapAlgebra SQL function''' -- Bborie has already done this.15 ~~ '''2) An optimized version of the two rasters ST_MapAlgebra SQL function''' ~~ -- Bborie has already done this. 16 16 17 17 Mentor: Pierre Racine (pierre.racine at sbf.ulaval.ca) 18 18 19 19 20 [http://www.quantdec.com/SYSEN597/GTKAV/section9/map_algebra.htm Map algebra] allows the creation of a new raster as a mathematical (or logical) function of one or two other rasters. e.g. raster3 = (raster1 + raster2) / 2. This is a very common raster operation in the GIS world. 20 ~~ [http://www.quantdec.com/SYSEN597/GTKAV/section9/map_algebra.htm Map algebra] allows the creation of a new raster as a mathematical (or logical) function of one or two other rasters. e.g. raster3 = (raster1 + raster2) / 2. This is a very common raster operation in the GIS world. ~~ 21 21 22 A [http://postgis.net/docs/manual-dev/RT_ST_MapAlgebraExpr two rasters version of the ST_MapAlgebra SQL function already exist]. It allows a PostGIS user to create a new raster coverage as a function of two existing one with a simple SQL query. e.g.: SELECT ST_MapAlgebra(rast1, rast2, "(rast1 + rast2) / 2") FROM coverage1, coverage2 WHERE ST_Intersects(rast1, rast2). The extent of the resulting raster can be equivalent to the extent of the FIRST, the SECOND, the INTERSECTION or the UNION of both raster. 22 ~~ A [http://postgis.net/docs/manual-dev/RT_ST_MapAlgebraExpr two rasters version of the ST_MapAlgebra SQL function already exist]. It allows a PostGIS user to create a new raster coverage as a function of two existing one with a simple SQL query. e.g.: SELECT ST_MapAlgebra(rast1, rast2, "(rast1 + rast2) / 2") FROM coverage1, coverage2 WHERE ST_Intersects(rast1, rast2). The extent of the resulting raster can be equivalent to the extent of the FIRST, the SECOND, the INTERSECTION or the UNION of both raster. ~~ 23 23 24 The current implementation compute the new raster one pixel at a time. It is possible, but complicated, to optimize the algorythm in order to fill many areas as a chunck of data instead of pixels by pixels. The project includes reviewing the proposed algorythm and implementing it. You will familiarize yourself with PostGIS development (the most used open source geospatial database), geospatial imagery concepts, open source development tools and code in C and PL/pgSQL. ~~~24 ~~ The current implementation compute the new raster one pixel at a time. It is possible, but complicated, to optimize the algorythm in order to fill many areas as a chunck of data instead of pixels by pixels. The project includes reviewing the proposed algorythm and implementing it. You will familiarize yourself with PostGIS development (the most used open source geospatial database), geospatial imagery concepts, open source development tools and code in C and PL/pgSQL. ~~ 25 25 26 26