Mercurial > gemma
view pkg/imports/sec.go @ 3733:ec86a7155377 concave-hull
Estimated too large triangles as triangles which have an edge which is at least 3.5 times as long as the standard dev of the longest egde per inner triangle.
author | Sascha L. Teichmann <sascha.teichmann@intevation.de> |
---|---|
date | Mon, 24 Jun 2019 11:39:09 +0200 |
parents | 4c585b5d4fe8 |
children | 6c5c15b2fb64 |
line wrap: on
line source
// This is Free Software under GNU Affero General Public License v >= 3.0 // without warranty, see README.md and license for details. // // SPDX-License-Identifier: AGPL-3.0-or-later // License-Filename: LICENSES/AGPL-3.0.txt // // Copyright (C) 2018,2019 by via donau // – Österreichische Wasserstraßen-Gesellschaft mbH // Software engineering by Intevation GmbH // // Author(s): // * Sascha L. Teichmann <sascha.teichmann@intevation.de> package imports import ( "context" "database/sql" "time" "gemma.intevation.de/gemma/pkg/models" ) type Section struct { Name string `json:"name"` From models.Isrs `json:"from"` To models.Isrs `json:"to"` Tolerance float32 `json:"tolerance"` ObjNam string `json:"objnam"` NObjNam *string `json:"nobjnam"` Source string `json:"source-organization"` Date models.Date `json:"date-info"` } const SECJobKind JobKind = "sec" type secJobCreator struct{} func init() { RegisterJobCreator(SECJobKind, secJobCreator{}) } func (secJobCreator) Description() string { return "section" } func (secJobCreator) AutoAccept() bool { return false } func (secJobCreator) Create() Job { return new(Section) } func (secJobCreator) Depends() [2][]string { return [2][]string{ {"sections"}, {"distance_marks_virtual", "waterway_axis", "waterway_area"}, } } const ( secDeleteSQL = ` DELETE FROM waterway.sections WHERE staging_done AND name = ( SELECT name FROM waterway.sections WHERE id = ( SELECT key from import.track_imports WHERE import_id = $1 AND relation = 'waterway.sections'::regclass) AND NOT staging_done )` secStageDoneSQL = ` UPDATE waterway.sections SET staging_done = true WHERE id IN ( SELECT key from import.track_imports WHERE import_id = $1 AND relation = 'waterway.sections'::regclass)` secInsertSQL = ` WITH bounds (b) AS (VALUES ( ($1::char(2), $2::char(3), $3::char(5), $4::char(5), $5::int)::isrs ), ( ($6::char(2), $7::char(3), $8::char(5), $9::char(5), $10::int)::isrs)), r AS (SELECT isrsrange( (SELECT b FROM bounds ORDER BY b USING <~ FETCH FIRST ROW ONLY), (SELECT b FROM bounds ORDER BY b USING >~ FETCH FIRST ROW ONLY)) AS r), axs AS ( SELECT ISRSrange_axis((SELECT r FROM r), $16::double precision) AS axs) INSERT INTO waterway.sections ( name, section, area, objnam, nobjnam, date_info, source_organization ) VALUES ( $11, (SELECT r FROM r), ISRSrange_area( (SELECT axs FROM axs), (SELECT ST_Collect(CAST(area AS geometry)) FROM waterway.waterway_area)), $12, $13, $14, $15) RETURNING id` ) // StageDone moves the imported section out of the staging area. func (secJobCreator) StageDone( ctx context.Context, tx *sql.Tx, id int64, ) error { if _, err := tx.ExecContext(ctx, secDeleteSQL, id); err != nil { return err } _, err := tx.ExecContext(ctx, secStageDoneSQL, id) return err } // CleanUp of a section import is a NOP. func (*Section) CleanUp() error { return nil } // Do executes the actual section import. func (sec *Section) Do( ctx context.Context, importID int64, conn *sql.Conn, feedback Feedback, ) (interface{}, error) { start := time.Now() if sec.Date.Time.IsZero() { sec.Date = models.Date{Time: start} } feedback.Info("Storing section '%s'", sec.Name) tx, err := conn.BeginTx(ctx, nil) if err != nil { return nil, err } defer tx.Rollback() var nobjnm sql.NullString if sec.NObjNam != nil { nobjnm = sql.NullString{String: *sec.NObjNam, Valid: true} } feedback.Info("Section from %s to %s.", sec.From.String(), sec.To.String()) feedback.Info("Tolerance used to snap waterway axis: %g", sec.Tolerance) var id int64 if err := tx.QueryRowContext( ctx, secInsertSQL, sec.From.CountryCode, sec.From.LoCode, sec.From.FairwaySection, sec.From.Orc, sec.From.Hectometre, sec.To.CountryCode, sec.To.LoCode, sec.To.FairwaySection, sec.To.Orc, sec.To.Hectometre, sec.Name, sec.ObjNam, nobjnm, sec.Date.Time, sec.Source, sec.Tolerance, ).Scan(&id); err != nil { return nil, handleError(err) } if err := track(ctx, tx, importID, "waterway.sections", id); err != nil { return nil, err } feedback.Info("Storing section '%s' took %s", sec.Name, time.Since(start)) if err := tx.Commit(); err != nil { return nil, err } feedback.Info("Import of section was successful") summary := sec // to provide full data for review return summary, nil }