查找点表的相应栅格值

Question

在我的PostgreSQL数据库中，有一个栅格表和一个具有相同投影的点表。点表有列，如id，address和geom。栅格表包含id，r_proj4和rast等列。栅格表的每一行描述一个栅格磁贴。如何检索每个点的栅格值？

我希望粗略的指导如何解决这个问题和PostGIS代码示例。

Answer 1

就在我的脑海中，我会说st_intersection()可能会让你拥有包含几何兴趣点的栅格（如果这就是你要找的东西）。您可以将它用作标准，同时内部连接两个表，如下所示：

select p1.*, r1.rast
from points p1
join rasters r1
on (st_intersection(r1.rast, p1.geom)).geom = p1.geom