名称
ST_MapAlgebraFctNgb — 1 波段版本：使用用户定义的 PostgreSQL 函数的地图代数最近邻。返回一个栅格，其值是 PLPGSQL 用户函数的结果，该函数涉及输入栅格波段中值的邻域。
语法
raster ST_MapAlgebraFctNgb(raster rast, integer band, text pixeltype, integer ngbwidth, integer ngbheight, regprocedure onerastngbuserfunc, text nodatamode, text[] VARIADIC args);
描述
![[Warning]](../images/warning.png) | |
ST_MapAlgebraFctNgb 已在 2.1.0 中弃用。请改用 ST_MapAlgebra（回调函数版本）。 |
（一个栅格版本）返回一个栅格，其值是 PLPGSQL 用户函数的结果，该函数涉及输入栅格波段中值的邻域。用户函数将像素值邻域作为数字数组，对于每个像素，返回用户函数的结果，用函数结果替换当前检查像素的像素值。
- rast
评估用户函数的栅格。
- band
要评估的栅格的波段号。默认为 1。
- pixeltype
输出栅格的结果像素类型。必须是 ST_BandPixelType 中列出的类型，或者留空或设置为 NULL。如果没有传入或设置为 NULL，将默认为
rast的像素类型。如果结果大于像素类型允许的值，则会截断结果。- ngbwidth
邻域的宽度（以像元为单位）。
- ngbheight
邻域的高度（以像元为单位）。
- onerastngbuserfunc
PLPGSQL/psql 用户函数，用于应用于栅格单波段的邻域像素。第一个元素是表示矩形像素邻域的二维数字数组
- nodatamode
定义为邻域像素（为 nodata 或 NULL）传递给函数的值
'ignore': 计算将忽略邻域中遇到的任何 NODATA 值 -- 此标志必须发送到用户回调函数，并且用户函数决定如何忽略它。
'NULL': 邻域中遇到的任何 NODATA 值将导致结果像素为 NULL -- 在这种情况下，将跳过用户回调函数。
'value': 邻域中遇到的任何 NODATA 值都将替换为参考像素（邻域中心的那个）。请注意，如果此值为 NODATA，则行为与 'NULL'（对于受影响的邻域）相同
- args
传递到用户函数的参数。
可用性：2.0.0
示例
示例利用作为单个切片加载的 katrina 栅格，如 http://trac.osgeo.org/gdal/wiki/frmts_wtkraster.html 中所述，然后在 ST_Rescale 示例中准备
--
-- A simple 'callback' user function that averages up all the values in a neighborhood.
--
CREATE OR REPLACE FUNCTION rast_avg(matrix float[][], nodatamode text, variadic args text[])
RETURNS float AS
$$
DECLARE
_matrix float[][];
x1 integer;
x2 integer;
y1 integer;
y2 integer;
sum float;
BEGIN
_matrix := matrix;
sum := 0;
FOR x in array_lower(matrix, 1)..array_upper(matrix, 1) LOOP
FOR y in array_lower(matrix, 2)..array_upper(matrix, 2) LOOP
sum := sum + _matrix[x][y];
END LOOP;
END LOOP;
RETURN (sum*1.0/(array_upper(matrix,1)*array_upper(matrix,2) ))::integer ;
END;
$$
LANGUAGE 'plpgsql' IMMUTABLE COST 1000;
-- now we apply to our raster averaging pixels within 2 pixels of each other in X and Y direction --
SELECT ST_MapAlgebraFctNgb(rast, 1, '8BUI', 4,4,
'rast_avg(float[][], text, text[])'::regprocedure, 'NULL', NULL) As nn_with_border
FROM katrinas_rescaled
limit 1;
|
 我们栅格的第一波段
|
 在 4x4 像素内对像素进行平均后生成的新栅格
|