Changes between Version 41 and Version 42 of import layers postgis

Timestamp:

Aug 23, 2018 1:25:53 PM (7 years ago)

Author:

maria

Comment:

--

import layers postgis

@@ -595 +595 @@
 ||75270||Lis e Ribeiras Costeiras||t||0101000020B30E000000E78C289886F1C00070F085892BD740||8.23462515206324||
 
-The target 75270 is the largest flow of this basin
+The target 75270 is the largest flow of this basin[[BR]]
 Be careful! The target and source attributes are calculating randomly (so this example can be a different number for you) 
 
-[[Image(source:eda/data/Docs/trac/sudoang/lis_ribeiras_example.png, 300px)]][[BR]]
-
-We select everything in basin 75270
+[[Image(source:eda/data/Docs/trac/sudoang/lis_ribeiras_example.png, 200px)]][[BR]]
+
+We select everything in basin 75270[[BR]]
 The following function (recusrsive) selects everyting in the basin by a recursive query, using source and target created by the pg_routing algorythm (it is much quicker than a spatial recursive)
 
@@ -656 +656 @@
 }}}
 
-Different options to calculate topology
+Different options to calculate topology:
+ -- For our example "Lis e Ribeiras Costeiras", in the largest flow "75270"
 
 {{{#!sql
@@ -665 +666 @@
 UPDATE portugal.rivers set chemin=get_path(75270,u.source) from portugal.upstream_segments(75270) u
 where rivers.source = u.source; -- 5:13
+}}}
+
+ -- Since this works, we use the second option to calculate topology in each basin (that are ordered by north/south)
+
+{{{#!sql
+--select portugal.get_path(target, source) from portugal.rivers where nome = 'Minho';--Rivers in Spain
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Ancora e Ribeiras Costeiras';--37.4 secs
+--select portugal.get_path(target, source) from portugal.rivers where nome = 'Lima';--Target outside the basin
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Neiva e Ribeiras Costeiras';--01:20 min
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Cávado e Ribeiras Costeiras';--07:46 min
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Ave';--06:29 min
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Leça e Ribeiras Costeiras';--01:30 min
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Douro';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Vouga e Ribeiras Costeiras';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Mondego';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Lis e Ribeiras Costeiras';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Ribeiras do Oeste';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Tejo';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Sado';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Ribeiras do Alentejo';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Mira';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Ribeiras do Algarve';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Arade';--
+select portugal.get_path(target, source) from portugal.rivers where nome = 'Guadiana';--
+}}}
+
+ -- May be we can use the third option to calculate topology at the same time, but we have not tried yet
+
+{{{#!sql
+UPDATE portugal.rivers set chemin=get_path(u.target,u.source) from portugal.upstream_segments(target) u
+where rivers.target = u.target, rivers.source = u.source; -- could it work for all basins at the same time?
 }}}
 
@@ -719 +751 @@
 select portugal.write_chemin('Lis e Ribeiras Costeiras');
 select * from portugal.downstream_points where nome = 'Lis e Ribeiras Costeiras';
-}}}
+
+-- For the rest of basins (except Minho and Lima for the moment)
+select portugal.write_chemin('Neiva e Ribeiras Costeiras' ,'Cávado e Ribeiras Costeiras','Ave',
+        'Leça e Ribeiras Costeiras','Douro','Vouga e Ribeiras Costeiras','Mondego','Lis e Ribeiras Costeiras'
+        'Ribeiras do Oeste','Tejo','Sado','Ribeiras do Alentejo','Mira','Ribeiras do Algarve','Arade');
+}}}