view pkg/imports/polygon.go @ 972:17a03a84b0e8

Split out polygon code out of sounding result importer source file.
author Sascha L. Teichmann <sascha.teichmann@intevation.de>
date Thu, 18 Oct 2018 12:22:08 +0200
parents
children a244b18cb916
line wrap: on
line source

package imports

import (
	"bytes"
	"encoding/binary"
	"fmt"
	"math"

	shp "github.com/jonas-p/go-shp"
)

type (
	Point struct {
		X float64
		Y float64
	}
	LineString []Point
	Polygon    []LineString
)

const (
	wkbNDR     byte   = 1
	wkbPolygon uint32 = 3
)

func shapeToPolygon(s shp.Shape) (Polygon, error) {
	switch p := s.(type) {
	case *shp.Polygon:
		return toPolygon(p.NumParts, p.Parts, p.Points), nil
	case *shp.PolygonZ:
		return toPolygon(p.NumParts, p.Parts, p.Points), nil
	case *shp.PolygonM:
		return toPolygon(p.NumParts, p.Parts, p.Points), nil
	}
	return nil, fmt.Errorf("Unsupported shape type %T", s)
}

func toPolygon(numParts int32, parts []int32, points []shp.Point) Polygon {
	out := make(Polygon, numParts)
	pos := 0
	for i := range out {
		ps := parts[i]
		line := make(LineString, ps)
		for j := int32(0); j < ps; j, pos = j+1, pos+1 {
			p := &points[pos]
			line[j] = Point{p.X, p.Y}
		}
		out[i] = line
	}
	return out
}

func (p Polygon) AsWBK() []byte {
	if p == nil {
		return nil
	}
	// pre-calculate size to avoid reallocations.
	size := 1 + 4 + 4
	for _, ring := range p {
		size += 4 + len(ring)*2*8
	}

	buf := bytes.NewBuffer(make([]byte, 0, size))

	binary.Write(buf, binary.LittleEndian, wkbNDR)
	binary.Write(buf, binary.LittleEndian, wkbPolygon)
	binary.Write(buf, binary.LittleEndian, uint32(len(p)))

	for _, ring := range p {
		binary.Write(buf, binary.LittleEndian, uint32(len(ring)))
		for _, v := range ring {
			binary.Write(buf, binary.LittleEndian, math.Float64bits(v.X))
			binary.Write(buf, binary.LittleEndian, math.Float64bits(v.Y))
		}
	}

	return buf.Bytes()
}