Changes between Version 34 and Version 35 of WKTRaster/SpecificationWorking03

Timestamp:

03/14/11 11:15:03 (14 years ago)

Author:

pracine

Comment:

—

WKTRaster/SpecificationWorking03

@@ -197 +197 @@
 
 ----
-== '''Objective FV.03 - Being able to use "group by" to accumulate tiles to form a new raster.''' ==
+== '''Objective FV.03 - Being able to use neighbour pixels in MapAlgebra expressions.''' ==
+
+For now ST_MapAlgebra expressions refer only to the pixel being computed. e.g. "rast * 100". The original plan was to allow refering to neighbour pixels using two coordinated relative to the pixel being computed. e.g. "rast[-1,0] * 100" where rast[-1,0] refer to the value of the pixel one pixel to the left of the pixel being computed. However this syntax might prove to be hard to use when many neighbours are to be used.
+
+An alternative syntax would involve another function name (e.g. ST_MapAlgebraNgb) and a way to define a neighbour rectangular region around the computed pixel (e.g.: "2,2" meaning a rectangle encompassing the two neighbour pixels in each direction) and a function to call with this matrix of pixel values. A complete example might look like:
+
+SELECT ST_MapAlgebraNgb(rast, "ST_Avg", 2, 2, "ignore") 
+
+So this would mean "for each pixel, compute the average of all the 1 + 8 + 16 = 25 pixels surrounding the current pixel and "ignore" pixels with nodata values."
+
+The "ST_Avg" summarizing function should accept three parameters: an array of float8 values, a X and a Y dimension, and optionnally a "what to do with nodata values". The possible value for this parameter could be:
+
+ *"NULL": If any value is a nodata value, return NULL.
+ *"ignore": Ignore any nodata value so that if 4 pixels on 25 are nodata values, do the computation with the 21 remaining.
+ *a value: Replace any nodata value with this value and compute.
+
+Users could write their own map algebra summarizing functions. 
+
+A more sophisticated version would pass a georeferenced raster instead of just a value matrix so that summarizing function could use this geoereference (e.g. to determine a ).
+
+Any remaining parameters to ST_MapAlgebraNgb could be passed to the summarizing functions for its own need.
+
+
+----
+== '''Objective FV.04 - Being able to use "group by" to accumulate tiles to form a new raster.''' ==