Le total des surfaces des zonages d'urbanisme ne correspond pas à la surface du territoire

[alexisig]

Exemple du problème pour Auch (dont la surface est égale à 7 341ha):

Le problème vient de cette requête :

ArtifAreaZoneUrba.objects.filter(zone_urba__in=ZoneUrba.objects.intersect(self.combined_emprise))

Intersect selectionne les objets qui croisent l'emprise du diagnostic (combined_emprise) mais également ceux qui touchent l'emprise du diagnostic, sans forcement le croiser. Cela a pour conséquence de retourner tous les zonages du territoire PLUS ceux voisins, qui touchent uniquement le territoire sur ses frontières.

Solution courante : combiner ST_Intersectavec un ST_touches== false
Cependant cette solution ne fonctionne pas toujours dans notre cas car les géométries des zonages sont parfois "mal dessinés", c'est à dire que les surfaces de zonages voisins se croisent parfois.

Solution choisie : pour chaque zonage, récupérer un PointOnSurface et vérifier si celui-ci intersecte avec l'emprise du diagnostic.

Résultat pour auch :