查找点表的相应栅格值

Question

在我的PostgreSQL数据库中，有一个具有相同投影的栅格表和一个点​​表。点表具有列，例如id，地址和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