使用ST_RPCRectify实现基于RPC的栅格校正-云原生数据库 PolarDB-阿里云-云原生数据库 PolarDB(PolarDB)-阿里云帮助中心

根据栅格影像的RPC（Rational Polynomial Coefficients）参数对栅格进行校正操作，返回校正后的栅格对象。

语法

raster ST_RPCRectify(raster rast,
 raster dem default NULL,
                    cstring rectifyOption default '{}',
                    cstring storageOption default '{}',
   cstring parallelOption default '{}')

参数

参数名称	描述
rast	需要校正的Raster对象。
dem	校正所参考的DEM栅格对象。如不需要请传入NULL。
rectifyOption	JSON字符串校正选项。
storageOption	返回结果的存储选项，格式为JSON字符串。
parallelOption	JSON字符串的并行选项。

描述

说明栅格对象必须具备RPC相关信息。

rectifyOption为JSON格式的字符串数组，每个子JSON对象指定参数如下：

参数名称	描述	类型	默认值	说明
resample	重采样方式	string	'Near'	栅格重采样方式，仅支持如下几种： Near Average Cubic Bilinear
outSrid	结果栅格空间参考	integer	与原Raster一致	无
nodataValue	按波段指定新的nodata值。	float8 float8[]	NULL	nodataValue可指定为单个值或数组。如果指定为单个值表示输出栅格对象的所有波段使用同一个nodata值。如果指定为数组[0, 0, 0]，则数组元素个数必须与栅格的波段数一致。如不指定则优先使用原栅格数据nodata，如原栅格无nodata值，使用0作为nodata。
RPC_DEM_MISSING_VALUE	当DEM值为nodata或栅格数据在DEM范围外时使用的高程值。	float8	0	无

storageOption参数如下：

参数名称	描述	类型	默认值	说明
chunking	是否使用分块存储。	boolean	与原Raster一致	无
chunkdim	分块的维度信息。	string	与原Raster一致	只有chunking参数为true时，该参数有效。
chunktable	分块表名称。	string	`''`	如果传入`''`值，则会在当前会话中创建一个随机表名的临时表用于存放转换后的Raster对象。该临时表只在当前会话中有效，会话结束后临时表也会随之删除。如果需要保存转换后的Raster对象，则需要指定具体表名。如果您使用的并行方式不支持匿名表，则需要提前创建表。
compression	压缩算法类型。	string	与原Raster一致	目前仅支持如下类型： NONE JPEG ZLIB PNG LZO LZ4 JP2K
quality	压缩质量。	integer	与原Raster一致	仅针对JPEG压缩算法。
interleaving	交错方式。	string	与原Raster一致	仅支持如下几种方式： bip：Band interleaved by pixel bil：Band interleaved by line bsq：Band Sequential
endian	字节序。	string	与原Raster一致	仅支持如下几种： NDR：Little endian XDR：Big endian

示例

CREATE TABLE if not exists raster_table(id integer, rast raster);

-- RPC栅格数据
INSERT INTO raster_table values(1, ST_ImportFrom('t_chunk','<RAST_DATA_DIR>/rpc.tif', '{}'));

-- DEM数据
INSERT INTO raster_table values(2, ST_ImportFrom('t_chunk','<RAST_DATA_DIR>/rpc_dem.tif', '{}'));

CREATE TABLE raster_result(id integer, rast raster);

--不指定dem，使用最邻近采样
INSERT INTO raster_result(id, rast) 
select 10, ST_RPCRectify(rast,
                       NULL,    
                         '{"resample":"Near"}',
                         '{"chunking":true,"chunkdim":"(256,256,1)","interleaving":"bsq","chunktable":"t_chunk"}') 
from raster_table 
where id =1;

--指定dem，使用最邻近采样，输出srid为32635，并使用并行计算
INSERT INTO raster_result(id, rast) 
select 11, ST_RPCRectify(rast,
                       （SELECT rast FROM raster_table WHERE id = 2),    
                         '{"resample":"Near", "outSrid": 32635}',
                         '{"chunking":true,"chunkdim":"(256,256,1)","interleaving":"bsq","chunktable":"t_chunk"}',
                         '{"parallel": 8}') 
from raster_table 
where id =1;