Чи нормально групуватися за геометрією в PostGIS?

Чи добре згрупуватися за геометрією?

Ми робимо багато підрахунків точок за геометрією багатокутника, що передбачає перетин даних спочатку, щоб підрахувати екземпляри учнів. за межами школи, блокової групи тощо:

SELECT
  n.nbhd_id
  , count(*) AS count_burglaries
FROM denver.crime AS c
  JOIN denver.neighborhoods AS n
    ON ST_Intersects(c.geom, n.geom)

WHERE c.offense_ty LIKE 'BURG%'
GROUP BY n.nbhd_id

Другий крок - загортання в підзапит, щоб приєднати геометрію від таблиці багатокутника назад до запиту підрахунку:

SELECT

count.*
, n.nbhd_name
, n.geom
FROM (
SELECT
n.nbhd_id
, count(*) AS count_burglaries
FROM denver.crime AS c
JOIN denver.neighborhoods AS n
ON ST_Intersects(c.geom, n.geom)

WHERE c.offense_ty LIKE 'BURG%'
GROUP BY n.nbhd_id

) AS count

JOIN denver.neighborhoods AS n
ON count.nbhd_id = n.nbhd_id

Але, здається, ви також можете використовувати геометрію в групі BY:

SELECT
n.nbhd_id
, n.nbhd_name
, n.geom

, count(*) AS count_burglaries
FROM denver.crime AS c
JOIN denver.neighborhoods AS n
ON ST_Intersects(c.geom, n.geom)

WHERE c.offense_ty LIKE 'BURG%'
GROUP BY n.nbhd_id, nbhd_name, n.geom

Чи правильним є підхід геометрії до GROUP BY?

З групуванням за геометрією в цьому випадку майже напевно немає нічого поганого, оскільки ви вже групуєтесь за унікальним ідентифікатором ( nbhd_id). І, як ви зазначаєте, це заощаджує вам приєднання та робить запит чистішим.

Важливо знати, що GROUP BY geomзастереження в PostGIS 2.3 та попередніх версіях насправді групує рядки на основі рівності обмежувальної рамки, а не геометричної рівності. Що стосується реальних даних, що не піддаються обробці, це часто має той самий ефект, що і групування на геометричній рівності (і набагато швидше), але може не дати очікуваних результатів, коли кілька різних багатокутників мають одне і те ж обмежувальне поле. Починаючи з PostGIS 2.4, GROUP BYпрацює на жорсткому визначенні геометричної рівності з урахуванням початкової точки та орієнтації геометрій.

Оскільки ви вже групуєтесь nbhd_id, ви все одно побачите окремі рядки, навіть якщо два різних полігони сусідства розділяють обмежувальне поле.