ST_MapAlgebra (expression version)
Prev  11.12. 栅格处理：地图代数  Next

Name

ST_MapAlgebra (expression version) — 表达式版本 - 在给定一两个输入栅格、波段索引和一个或多个用户指定的 SQL 表达式的情况下，返回单波段栅格。

Synopsis

raster ST_MapAlgebra(raster rast, integer nband, text pixeltype, text expression, double precision nodataval=NULL);

raster ST_MapAlgebra(raster rast, text pixeltype, text expression, double precision nodataval=NULL);

raster ST_MapAlgebra(raster rast1, integer nband1, raster rast2, integer nband2, text expression, text pixeltype=NULL, text extenttype=INTERSECTION, text nodata1expr=NULL, text nodata2expr=NULL, double precision nodatanodataval=NULL);

raster ST_MapAlgebra(raster rast1, raster rast2, text expression, text pixeltype=NULL, text extenttype=INTERSECTION, text nodata1expr=NULL, text nodata2expr=NULL, double precision nodatanodataval=NULL);

描述

表达式版本 - 在给定一两个输入栅格、波段索引和一个或多个用户指定的 SQL 表达式的情况下，返回单波段栅格。

可用性：2.1.0

描述：变体 1 和 2（一个栅格）

通过应用由输入栅格 (rast) 上的expression定义的有效 PostgreSQL 代数运算来创建新的单波段栅格。 如果未提供 nband，则假定为波段 1。 新栅格将具有与原始栅格相同的地理参考、宽度和高度，但只有一个波段。

如果传入像素类型，则新栅格将具有该像素类型的波段。 如果将像素类型传递为 NULL，则新的栅格带将具有与输入栅格带相同的像素类型。

  • 允许expression的关键字

    1. [rast] - 感兴趣像素的像素值

    2. [rast.val]- 感兴趣像素的像素值

    3. [rast.x] - 感兴趣像素的基于 1 的像素列

    4. [rast.y] - 感兴趣像素的从 1 开始的像素行

描述：变体 3 和 4（两个栅格）

通过对两个输入栅格波段 rast1， (rast2) 上的expression定义的两个波段应用有效的 PostgreSQL 代数运算，创建一个新的单波段栅格。 如果没有 band1，则指定 band2，则假定为 band 1。 生成的栅格将在第一个栅格定义的网格上对齐（比例、倾斜和像素角）。 生成的栅格将具有由extenttype 参数定义的范围。

expression

涉及两个栅格的 PostgreSQL 代数表达式和 PostgreSQL 定义的函数/运算符，用于定义像素相交时的像素值。 例如 (([rast1] [rast2])/2.0)::integer

pixeltype

输出栅格的结果像素类型。 必须是ST_BandPixelType 中列出的一项，可省略或设置为 NULL。 如果未传入或设置为 NULL，将默认为第一个栅格的像素类型。

extenttype

控制生成的栅格的范围

  1. INTERSECTION - 新栅格的范围是两个栅格的交集。 这是默认设置。

  2. UNION - 新栅格的范围是两个栅格的并集。

  3. FIRST - 新栅格的范围与第一个栅格的范围相同。

  4. SECOND - 新栅格的范围与第二个栅格的范围相同。

nodata1expr

仅涉及 rast2 的代数表达式或定义当 rast1 的像素为无数据值并且空间对应的 rast2 像素具有值时返回的内容的常量。

nodata2expr

仅涉及 rast1 常量的代数表达式，该常量定义当 rast2 的像素为无数据值且空间对应的 rast1 像素具有值时返回的内容。

nodatanodataval

当空间对应的 rast1 和 rast2 像素均为无数据值时要返回的数值常量。

  • expression、nodata1expr 和 nodata2expr 中允许使用的关键字

    1. [rast1] - rast1 中感兴趣的像素的像素值

    2. [rast1.val] - rast1 中感兴趣像素的像素值

    3. [rast1.x] - 来自 rast1 的感兴趣像素的基于 1 的像素列

    4. [rast1.y] - 来自 rast1 的感兴趣像素的基于 1 的像素行

    5. [rast2] - rast2 中感兴趣像素的像素值

    6. [rast2.val] - rast2 中感兴趣像素的像素值

    7. [rast2.x] - 来自 rast2 的感兴趣像素的基于 1 的像素列

    8. [rast2.y] - 来自 rast2 的感兴趣像素的基于 1 的像素行

示例：变体1和2

WITH foo AS (
    SELECT ST_AddBand(ST_MakeEmptyRaster(10, 10, 0, 0, 1, 1, 0, 0, 0), '32BF'::text, 1, -1) AS rast
)
SELECT
    ST_MapAlgebra(rast, 1, NULL, 'ceil([rast]*[rast.x]/[rast.y]+[rast.val])')
FROM foo;

示例：变体3和4

WITH foo AS (
    SELECT 1 AS rid, ST_AddBand(ST_AddBand(ST_AddBand(ST_MakeEmptyRaster(2, 2, 0, 0, 1, -1, 0, 0, 0), 1, '16BUI', 1, 0), 2, '8BUI', 10, 0), 3, '32BUI'::text, 100, 0) AS rast UNION ALL
    SELECT 2 AS rid, ST_AddBand(ST_AddBand(ST_AddBand(ST_MakeEmptyRaster(2, 2, 0, 1, 1, -1, 0, 0, 0), 1, '16BUI', 2, 0), 2, '8BUI', 20, 0), 3, '32BUI'::text, 300, 0) AS rast
)
SELECT
    ST_MapAlgebra(
        t1.rast, 2,
        t2.rast, 1,
        '([rast2] + [rast1.val]) / 2'
    ) AS rast
FROM foo t1
CROSS JOIN foo t2
WHERE t1.rid = 1
    AND t2.rid = 2;