Mercurial > gemma
comparison cmd/octreediff/main.go @ 2583:5c6a7e69b02d
Morphological differences: Store calculated iso-lines in database.
author | Sascha L. Teichmann <sascha.teichmann@intevation.de> |
---|---|
date | Mon, 11 Mar 2019 16:42:38 +0100 |
parents | 647a58ee9ae9 |
children |
comparison
equal
deleted
inserted
replaced
2582:ab819b1b748b | 2583:5c6a7e69b02d |
---|---|
63 ) | 63 ) |
64 ), | 64 ), |
65 4326 | 65 4326 |
66 ) | 66 ) |
67 ` | 67 ` |
68 insertContourSQLClipped = ` | |
69 WITH joined AS ( | |
70 SELECT | |
71 sr.area AS area, | |
72 sr.date_info AS date_info | |
73 FROM waterway.sounding_results sr JOIN | |
74 waterway.bottlenecks bn ON sr.bottleneck_id = bn.id | |
75 WHERE bn.bottleneck_id = $1 | |
76 ), | |
77 inter AS ( | |
78 SELECT | |
79 ST_Buffer( | |
80 ST_intersection( | |
81 (SELECT ST_Transform(area::geometry, $2::int) FROM joined WHERE date_info = $3::date), | |
82 (SELECT ST_Transform(area::geometry, $2::int) FROM joined WHERE date_info = $4::date) | |
83 ), | |
84 0.001) AS area | |
85 ) | |
86 INSERT INTO diff_contour_lines ( | |
87 height, | |
88 lines | |
89 ) | |
90 SELECT | |
91 $5, | |
92 ST_Transform( | |
93 ST_Multi( | |
94 ST_intersection( | |
95 ST_CollectionExtract( | |
96 ST_SimplifyPreserveTopology( | |
97 ST_Multi(ST_Collectionextract( | |
98 ST_MakeValid(ST_GeomFromWKB($6, $2::integer)), 2)), | |
99 $7 | |
100 ), | |
101 2 | |
102 ), | |
103 area | |
104 ) | |
105 ), | |
106 4326 | |
107 ) | |
108 FROM inter` | |
109 ) | 68 ) |
110 | 69 |
111 func check(err error) { | 70 func check(err error) { |
112 if err != nil { | 71 if err != nil { |
113 log.Fatalf("error: %v\n", err) | 72 log.Fatalf("error: %v\n", err) |