Я загрузил функцию hex_grid, она отлично работает см. Здесь , но теперь у меня есть шестиугольная сетка, я хотел бы сложить все значения из перекрывающихся точек центроида из слоя многоугольника. Я использовал с утверждениями, чтобы держать вещи немного аккуратными, оба оператора select работают как положено, моя проблема возникает при присоединении, ясно, что я сделал что-то не так, но я просто не вижу этого, любая помощь будет оценена. Вот код Гекс и точки находятся в одном и том же crs.
with hex AS(
SELECT hex_grid(
1,
(ST_XMIN(sa.geom)),
(ST_YMIN(sa.geom)),
(ST_XMAX(sa.geom)),
(ST_YMAX(sa.geom)),
32635, 32635, 32635) as geom_h
FROM geodata.study_area sa),
points AS(
SELECT count(*) as count, ST_Centroid(su.geom) AS geom_p
FROM clas.survey_unit su
group by su.geom
)
SELECT hex.geom_h, points.count
FROM hex
JOIN points
ON ST_Contains(points.geom_p, hex.geom_h)
group by hex.geom_h, points.count;