view cmd/octree2contour/store.go @ 919:271561dce2e6

Store contour lines in 2D We have the height as an attribute value and use the lines only for display on a map. There's no need to store height with every vertex additionally.
author Tom Gottfried <tom@intevation.de>
date Thu, 04 Oct 2018 16:32:47 +0200
parents fb9faf2c4f60
children 862cc5ba4d1d
line wrap: on
line source

package main

import (
	"database/sql"
	"time"

	"gemma.intevation.de/gemma/pkg/octree"
)

type result struct {
	h     float64
	lines octree.MultiLineStringZ
}

const (
	deleteSQL = `
DELETE FROM waterway.sounding_results_contour_lines
WHERE sounding_result_id IN (SELECT id
FROM waterway.sounding_results
WHERE bottleneck_id = $1 AND date_info = $2)
`
	insertSQL = `
INSERT INTO waterway.sounding_results_contour_lines
  (sounding_result_id, height, lines)
SELECT
  sr.id,
  $1,
  ST_Transform(
    ST_CollectionExtract(
      ST_Force2D(
        ST_Intersection(
          ST_Transform(sr.area::geometry, $3::integer),
          ST_GeomFromWKB($2, $3::integer)
        )
      ),
      2
    ),
    4326
  )
FROM waterway.sounding_results sr
WHERE bottleneck_id = $4 AND date_info = $5
`
)

func store(
	all []result, epsg uint32,
	bottleneck string, date time.Time,
) error {

	return run(func(db *sql.DB) error {

		tx, err := db.Begin()
		if err != nil {
			return err
		}
		defer tx.Rollback()

		if _, err := tx.Exec(deleteSQL, bottleneck, date); err != nil {
			return err
		}

		stmt, err := tx.Prepare(insertSQL)
		if err != nil {
			return err
		}

		for _, r := range all {
			if _, err := stmt.Exec(
				r.h, r.lines.AsWKB(), epsg,
				bottleneck, date,
			); err != nil {
				return err
			}
		}

		return tx.Commit()
	})
}