如何加快对栅格数据库的查询？

我在postgresql / postgis中有一个带有以下列的栅格数据库：

（ID，rast，data_of_data）。

“ rast”是具有WKT格式的栅格文件的列。以下示例查询用于查找WGS84系统（30.424，-1.66）中某个点的DN值（对于2002-01-09）：

SELECT 
     st_value(rast,(st_GeomFromText('POINT(30.424 -1.66)', 4326))) as val
FROM 
     my_table
WHERE
     date_of_data='2002-01-09'

是否有一种方法（例如空间索引）来加快此类查询的速度？

这是一个令人兴奋的问题！您要查询的栅格有多大？WKTRaster 作为BLOB存储在数据库中。为了找到特定点的值，使用（dx，dy）步骤和旋转从已知的（x_0，y_0）角坐标行/列索引（i，j）计算。使用（i，j）已知，ST_Value（）函数可以以正确的字节偏移量访问实际数据。

这意味着，在回答一个查询点时，DB必须平均读取至少一半的数据块（取决于实现方式，它实际上可能一直在读取所有数据）。因此，我想这WKTRaster性能会受到影响，当数据的BLOB得到太大。切片数据集应加快查询速度。看一下本教程中如何处理SRTM数据（以6000x6000像素块为单位）。他们实际上将数据平铺成非常小的50x50像素，这清楚地表明我的猜测可能与事实相差不远。

在空间上为栅格数据建立索引可能只会对边界框建立索引，这对您的问题没有真正的帮助。

我发现加快了PostGIS栅格计算的两个方面，就是在栅格中使用整数值，并在可能的情况下使用多波段栅格。在这种情况下，是否可以将DN值存储为整数？

另一个想法（我不确定这里是否有意义）是使用多波段栅格。例如，如果您要查看每月的数据切片，则每个月可能是一个栅格图层。然后，您可以通过查询分层栅格来检索不同时间片上某个点的多个值。我发现这种方法比查询单独的栅格要快得多。

最后，当您加载数据时，会有TILE_SIZE的-t标志。您可以探索正在使用的图块大小是否适合您的查询。

根据数据的分布，仅通过索引date_of_data列即可获得一些非常好的加速效果。

您可以使用EXPLAIN ANALYZE语法来确定是否正在使用索引。