back to first page[..][[BR]]
back to ["RHT"] ["Noeuds mer RHT"][[BR]]


Le parcours est fait dans R par la méthode distance_sea.R voir main_RHT_France.R
{{{
#!sql
--pour des raisons d'efficacité on sauve un texte depuis R et on recharge en sql
--Sauvegarde sous R
--write.table(rht@data[,c("id_drain","ltree_from_sea","dmer")],file=str_c(datawd,"/dataEDArht/rht_ltree_distance.csv"),sep=";",row.names=FALSE)
--Restauration en sql
create table rht.temp_ltree_distance(
	id_drain integer primary key,
	ltreetext text,
	distance bigint);
copy rht.temp_ltree_distance from 'E:/workspace/EDAdata/dataEDArht/rht_ltree_distance.csv' with csv header delimiter as ';' NULL AS 'NA';--114600

-- lancer le script ltree.sql qui se trouve dans postgres/share
alter table rht.rht_topology add column chemin ltree;
update rht.rht_topology t set chemin=sub.chemin from (
select text2ltree(ltreetext) as chemin,id_drain from rht.temp_ltree_distance) as sub

where t.id_drain=sub.id_drain; --114600

-- mise à jour des distances mer

update rht.rht_topology t set dmer=sub.distance from (
select distance,id_drain from rht.temp_ltree_distance) as sub
where t.id_drain=sub.id_drain
and dmer is null;  --113517
}}}


[[Image(source:data/Docs/trac/rht/distances_mer.JPG)]] 
== Récupération des id_drains mer des tronçons ==
{{{
#!sql
-- quels sont les id_drain mer des tronçons
select cast(ltree2text(subltree(chemin,0,1)) as integer) from rht.rht_topology ;
}}}