Le portail francophone de la géomatique

Merci de votre réponse,
En faite, je ne suis pas trop à l'aise avec les scripts python et du coup je me demandais s'il n'était pas possible dans la calculatrice de champ de rajouter un attribut Dist Haussdorf qui calculerait automatiquement celle-ci sur chaque paire de surfaces (issue d'une jointure entre mes 2 couches).
Car en regardant la doc sur l'opérateur hausdorff_distance, il semble qu'on ne puisse l'appliquer que sur une seule paire de surface et non sur un ensemble de paires de surface
(exemple de code pour une paire :
hausdorff_distance( geom_from_wkt('LINESTRING (130 0, 0 0, 0 150)'),geom_from_wkt('LINESTRING (10 10, 10 150, 130 10)'),0.5))

Je pense que ce serait plus simple et plus efficace de laisser ce type de traitement à PostGIS, qui dispose d'une fonction de distance de Hausdorff.

Cela permettrait d'affiner le traitement exact à réaliser, en fonction de si vous voulez comparer les polygones spatialement ou juste leurs formes géométriques. À priori c'est spatial étant donné que vous avez deux millésimes de données.

Vous pourrez ainsi utiliser les indexes spatiaux sur les données afin de ne comparer que les géométries qui le méritent, dont les bounding box s'intersectent par exemple.

Exemple de requête :

select
  a.id
  , b.id
  , st_HausdorffDistance(a.geom, b.geom) as dist
from
  old_table as a
join
  new_table as b
on
  st_intersects(a.geom, b.geom);

Dans le cas où on a un id commun on peut directement jointurer sur cet ID :

select
  a.id, st_HausdorffDistance(a.geom, b.geom) as dist
from
  old_table as a
join
  new_table as b
on
  a.id = b.id

Pour la doc c'est par là : https://postgis.net/docs/ST_HausdorffDistance.html

Tout cela est exécutable directement dans le DB Manager de QGIS.

Merci beaucoup de vos conseils mais pour être franc, je n'ai jamais utilisé PostGis.
Je vais essayer de m'en sortir, pour l'instant je sèche sur l'import de mes couches (message d'erreur : Erreur 9 :le fournisseur de données virtual n'a aucune méthode createEmptyLayer) et ne sais pas où positionner le code que vous proposez  dans l'interface de Postgis.
La route est longue ...