Comparison between distances from the ccm and our calculated distances

[cedric]

Nos distances sont écrites dans la base CCM table riversegments
Ci dessous une requete pour les comparer

• Requête 1 : Permet de comparer la distance calculée "cum_len_sea" et non calculée len_tom (les 100ères lignes sont affichées)

  select cum_len_sea,len_tom from rivernodes
   inner join riversegments on riversegments.tonode=rivernodes.id 
  limit 100
  

• Requête 2 : Compte le nombre de différences entre les deux distances (soit 30384)

  select count(*) from(
  select cum_len_sea,len_tom from rivernodes
   inner join riversegments on riversegments.tonode=rivernodes.id 
  where cum_len_sea != len_tom) as sub
  -- 30384
  

• Requête 3 : Compte le nombre de distances identiques (soit 86380)

  select count(*) from(
  select cum_len_sea,len_tom from rivernodes
   inner join riversegments on riversegments.tonode=rivernodes.id 
  where cum_len_sea != len_tom) as sub
  -- 86380
  

C'est vraiment curieux quand on regarde du plus près

• Requête 4 : Calcul de la différence entre les deux distances cum_len_sea et len_tom

  select cum_len_sea -len_tom as diff from rivernodes
   inner join riversegments on riversegments.tonode=rivernodes.id 
  

• Requête 5 : Compte le nombre de lignes avec des différences identiques

  select count(*),diff from(
  select cum_len_sea -len_tom as diff from rivernodes
   inner join riversegments on riversegments.tonode=rivernodes.id ) as sub
  group by diff
  

Résultats :

count	diff (m)
4	-87532
6	-57127
30338	-26407
36	-11560
86380	0

Table pour carte de comparaison

drop table if exists tabledistance;
CREATE TABLE tabledistance AS(
select cum_len_sea-len_tom as diff, rivernodes.the_geom from rivernodes
 inner join riversegments on riversegments.tonode=rivernodes.id 
where cum_len_sea != len_tom);
alter table tabledistance add column id serial;
alter table tabledistance ADD CONSTRAINT pktabledistance_id PRIMARY KEY(id);
--  
drop table tabledistance;

J'ai trouvé d'où vient le problème !! :

• Pour le calcul de la distance à la mer (script R), on utilise la variable "length", or pour les distances à la mer qui posent problèmes, on observe une différence entre "length" et "shape.length", le plus souvent c'est la première distance "length" sur laquelle se base le calcul qui pose problème. Pour les distances à la mer qui sont bonnes il n'y a pas de différence entre les variables "length" et "shape.length".

• Pour le calcul des distances à la mer il faudrait donc utiliser "shape.length" au lieu de "length", qui dep plus est plus précis (valeur non arrondie).

• En relançant main_CCM.R (EDA/EDACCM/main_CCM.R) :

Plus de différences observées pour le calcul des distances à la mer entre les deux méthodes.

[celine]

Je ne comprend toujours pas ni ai trouvé d'où venait cette différence, mais ça m'aura permis de me pencher plus sur les données.
Cédric, as tu trouvé d'où cela pouvait venir ?

[cedric]

Milestone Data integration deleted