Mercurial > gemma

comparison cmd/octreediff/main.go @ 2583:5c6a7e69b02d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
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)