From 1fcab351ae05ec0e3751ed044c9f0117f0f1eaa7 Mon Sep 17 00:00:00 2001
From: "Alexis A." <alexis.athlani@gmail.com>
Date: Thu, 19 Sep 2024 11:54:35 +0200
Subject: [PATCH] fix(project): make has_zonage_urbanisme use point_on_surface

---
 project/models/project_base.py | 27 ++++++++++++++++++++++++++-
 1 file changed, 26 insertions(+), 1 deletion(-)

diff --git a/project/models/project_base.py b/project/models/project_base.py
index f2bce0821..d7f83a05b 100644
--- a/project/models/project_base.py
+++ b/project/models/project_base.py
@@ -569,7 +569,32 @@ def has_no_ocsge_coverage(self) -> bool:
 
     @property
     def has_zonage_urbanisme(self) -> bool:
-        return ArtifAreaZoneUrba.objects.filter(zone_urba__mpoly__intersects=self.combined_emprise).exists()
+        """
+        Cette fonction vérifie si l'emprise du diagnostic intersecte
+        des zonages d'urbanisme. Pour effectuer cette vérification,
+        on commence par sélectionner les zonages avec ST_Intersects,
+        puis on vérifie que le point central de chaque zone est
+        bien dans l'emprise du diagnostic. Cette vérification est
+        nécessaire car la méthode ST_Intersects peut retourner des
+        zones qui intersectent très légèrement l'emprise du projet
+        sur un bord.
+        """
+        zone_urbas_intersecting = ArtifAreaZoneUrba.objects.filter(
+            zone_urba__mpoly__intersects=self.combined_emprise,
+        ).values_list("zone_urba")
+
+        return (
+            ZoneUrba.objects.filter(
+                checksum__in=zone_urbas_intersecting,
+            )
+            .annotate(
+                point_on_surface=PointOnSurface("mpoly"),
+            )
+            .filter(
+                point_on_surface__intersects=self.combined_emprise,
+            )
+            .exists()
+        )
 
     def get_ocsge_millesimes(self):
         """Return all OCS GE millésimes available within project cities and between