Le portail francophone de la géomatique

-- intersection entre la grille et les communes avec calcul de la surface commune
DROP TABLE IF EXISTS fra_h3_res6_com2022;
CREATE TABLE fra_h3_res6_com2022 AS
SELECT a.*, b.codgeo, ST_AREA(ST_INTERSECTION(a.geometry, b.geometry)) AS surface
FROM fra_h3_res6 a, a_com2022 b
WHERE ST_INTERSECTS(a.geometry, b.geometry)
ORDER BY h3, surface DESC;
-- sélection de la commune avec la plus grande surface
DROP TABLE IF EXISTS fra_h3_res6_com2022_codgeo;
CREATE TABLE fra_h3_res6_com2022_codgeo AS
SELECT DISTINCT ON (h3) *
FROM fra_h3_res6_com2022
ORDER BY h3, surface DESC
;