105 lines
3.4 KiB
SQL
105 lines
3.4 KiB
SQL
/*
|
|
==========================================================================
|
|
|
|
phase 2 : préparation des données
|
|
|
|
==========================================================================
|
|
*/
|
|
|
|
-- dans la base redadeg on a chargé des données provenant de umap
|
|
-- ces données sont en 3857. On va les passer en 2154 (Lambert 93).
|
|
|
|
|
|
|
|
|
|
-- ici on fait une grosse requête qui va recaler les PK secteurs sur les nœuds routables
|
|
-- ça insère les points de secteurs dans la table en 2154
|
|
|
|
TRUNCATE TABLE phase_2_pk_secteur ;
|
|
|
|
-- on fait une table de base qui contient les relations entre le point à recaler et les poins de référence
|
|
-- on limite à un buffer de 25
|
|
WITH candidates AS
|
|
(
|
|
SELECT
|
|
pk_org.id AS pk_id,
|
|
pk_org.name AS name,
|
|
pk_org.secteur_id AS secteur_id,
|
|
node.id AS node_id,
|
|
ST_Distance(pk_org.the_geom, ST_ClosestPoint(node.the_geom, pk_org.the_geom)) AS distance,
|
|
ST_Snap(
|
|
pk_org.the_geom, -- le point d'origine à recaler
|
|
ST_ClosestPoint(node.the_geom, pk_org.the_geom), -- le point le plus près dans la couche de nœuds
|
|
ST_Distance(pk_org.the_geom, ST_ClosestPoint(node.the_geom, pk_org.the_geom))* 1.01 -- dans la distance de ce plus proche point
|
|
) AS the_geom
|
|
FROM
|
|
(SELECT id::integer, name, secteur_id, ST_Transform(the_geom,2154) AS the_geom FROM phase_2_pk_secteur_3857) AS pk_org,
|
|
(SELECT id, the_geom FROM osm_roads_pgr_vertices_pgr) AS node
|
|
WHERE
|
|
ST_INTERSECTS(node.the_geom, ST_BUFFER(ST_Transform(pk_org.the_geom,2154) ,25) )
|
|
ORDER BY pk_org.id, ST_Distance(pk_org.the_geom, ST_ClosestPoint(node.the_geom, pk_org.the_geom))
|
|
)
|
|
-- à partir de cette table on va faire une jointure entre les PK org et les nœuds ramenés par la sous-requête
|
|
INSERT INTO phase_2_pk_secteur
|
|
SELECT
|
|
a.pk_id AS id,
|
|
a.name,
|
|
b.node_id AS pgr_node_id,
|
|
a.secteur_id::integer,
|
|
b.the_geom
|
|
FROM candidates a JOIN
|
|
(
|
|
-- on fait une table qui ordonne les points d'accroche
|
|
SELECT
|
|
pk_id, node_id,
|
|
-- le rang va permettre de donner le rang de chaque rapprochement et d'en faire un critère
|
|
row_number() over (partition by pk_id order by min(distance)) AS rang,
|
|
min(distance) AS min_distance,
|
|
the_geom
|
|
FROM candidates
|
|
GROUP BY pk_id, secteur_id, node_id, the_geom
|
|
) b ON a.pk_id = b.pk_id
|
|
WHERE b.rang = 1
|
|
GROUP BY a.pk_id, a.name, a.secteur_id, b.node_id, b.the_geom
|
|
ORDER BY a.pk_id ;
|
|
|
|
|
|
|
|
-- on recale également les points de nettoyage sur le tracé
|
|
TRUNCATE TABLE phase_2_point_nettoyage ;
|
|
|
|
WITH candidates AS
|
|
(
|
|
SELECT
|
|
pt_org.ogc_fid AS pt_id,
|
|
edge.id AS edge_id,
|
|
ST_Distance(pt_org.the_geom, ST_ClosestPoint(edge.the_geom, pt_org.the_geom)) AS distance,
|
|
ST_Snap(
|
|
pt_org.the_geom, -- le point d'origine à recaler
|
|
ST_ClosestPoint(edge.the_geom, pt_org.the_geom), -- le point le plus près dans la couche de nœuds
|
|
ST_Distance(pt_org.the_geom, ST_ClosestPoint(edge.the_geom, pt_org.the_geom))* 1.01 -- dans la distance de ce plus proche point
|
|
) AS the_geom
|
|
FROM
|
|
(SELECT ogc_fid::integer, ST_Transform(the_geom,2154) AS the_geom FROM phase_2_point_nettoyage_3857) AS pt_org,
|
|
(SELECT id, the_geom FROM osm_roads_pgr) AS edge
|
|
WHERE
|
|
ST_INTERSECTS(edge.the_geom, ST_BUFFER(ST_Transform(pt_org.the_geom,2154) ,2) )
|
|
ORDER BY pt_org.ogc_fid, ST_Distance(pt_org.the_geom, ST_ClosestPoint(edge.the_geom, pt_org.the_geom))
|
|
)
|
|
INSERT INTO phase_2_point_nettoyage
|
|
SELECT
|
|
nextval('phase_2_point_nettoyage_id_seq'::regclass),
|
|
pt_id,
|
|
edge_id,
|
|
distance,
|
|
the_geom
|
|
FROM candidates ;
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|