Message posté par : Christophe Vergon ---------------------------------------- Bonjour, Problématique : j'ai deux référentiels topographiques de type vecteur, ici le parcellaire actuel (ensemble A) et le parcellaire napoléonien (Ensemble N). J'ai géoréférencé sans déformation (transfo d'Helmert) l'ensemble N en me servant de A comme référence. Je constate que des arcs qui devraient être strictement égaux sont légèrement différents, jusque là normal. Comment déformer l'ensemble N pour le faire caler exactement avec A il existe de nombreuses méthodes, mais quelle que soit la méthode, il faut connaître le vecteur constitué par deux sommets homologues de A et N. Etape 1 : Création d'un buffer à partir des linestring de A. On cherche ensuite toutes les linestring de N qui sont strictement incluses dans un élément de A, on retient le couple idA et idN dans la table nommée ici "testedgeid" Etape 2 : Utilisation de la distance de Hausdorff pour comparer les éléments de A et de N pour chaque ligne de testedgeid ----------------- Code : WITH p as (SELECT ida, idn, ST_HausdorffDistance(a.geom,n.geom) d, MIN(ST_HausdorffDistance(a.geom,n.geom)) OVER (PARTITION BY ida) mind, a.geom ag ,n.geom ng INTO teshausdorff FROM testedgeid, toponapo.edge_data n, topo.edge_data a WHERE a.edge_id=ida AND n.edge_id=idn ORDER BY a.edge_id,n.edge_id ) SELECT * INTO teshausdorff8 FROM p WHERE d=mind AND mind<=2 ORDER BY ida,idn ----------------- Ceci permet de déterminer pour chaque élément de A l'élément de N compris dans le buffer qui est le plus semblable à l'élément de A. On choisit ensuite de façon empirique la valeur "acceptable pour la comparaison des formes (ici 2) 3 Une fois les couples à comparer déterminés par les étapes précédentes, on transforme la LINESTRING en sommets, ce pour A, ce Pour N. Puis on appaire les sommets de A et ceux de N en fonction de la distance mini. Chaque sommet de A est comparé à chaque sommet de N on retient le couple a,n tel que la distance a,n soit minimale. Puis on calcul le vecteur : ----------------- Code : WITH p as (SELECT edge_id as ea, idn, geom FROM topo.edge_data,teshausdorff8 WHERE ida=edge_id ), p1 as (SELECT ea, st_dumppoints(geom) as d, geom as ag, idn FROM p), p2 as (SELECT ea,(d).path as ap,(d).geom as ga,ag, idn FROM p1), q as (SELECT edge_id en,geom ng, ida FROM toponapo.edge_data,teshausdorff8 WHERE idn=edge_id ), q1 as (SELECT en, st_dumppoints(ng) as d,ng,ida from q), q2 as (SELECT en,(d).path as np,(d).geom as gn,ng,ida FROM q1), K as (SELECT ida,idn,ap,np,st_distance(ga,gn) as d, min(st_distance(ga,gn)) OVER (PARTITION BY ida,idn,ga) mindp,ga,gn FROM p2,q2 WHERE idn=en AND ida=ea ORDER BY ea,idn,ap,np ) SELECT k.*, st_makeline(ga,gn) vg INTO testvec9 FROM k WHERE d=mindp ----------------- Réalisé à l'échelle d'une commune ça donne des résultats sympas. En jaune N en bleu A en rouge le vecteur de déformation. ---------------------------------------- Ce message est accompagné de fichiers, pour les télécharger, suivre le lien ci-dessous. ---------------------------------------- Le message est situé https://georezo.net/forum/viewtopic.php?pid=329636#p329636 Pour y répondre : geobd(a)ml.georezo.net ou reply de votre messagerie Pour vous désabonner connectez-vous sur le forum puis Profil / Abonnement -- Association GeoRezo - le portail géomatique https://georezo.net