srs.go

// Copyright 2021 Airbus Defence and Space
//
// Licensed under the Apache License, Version 2.0 (the "License");
// you may not use this file except in compliance with the License.
// You may obtain a copy of the License at
//
//     http://www.apache.org/licenses/LICENSE-2.0
//
// Unless required by applicable law or agreed to in writing, software
// distributed under the License is distributed on an "AS IS" BASIS,
// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
// See the License for the specific language governing permissions and
// limitations under the License.

package godal

import "fmt"

type srWKTOpts struct {
	errorHandler ErrorHandler
}

//WKTExportOption is an option that can be passed to SpatialRef.WKT()
//
// Available WKTExportOptions are:
//  - ErrLogger
type WKTExportOption interface {
	setWKTExportOpt(sro *srWKTOpts)
}

type trnOpts struct {
	errorHandler ErrorHandler
}

// TransformOption is an option that can be passed to NewTransform
//
// Available TransformOptions are:
//  - ErrLogger
type TransformOption interface {
	setTransformOpt(o *trnOpts)
}

func (sr *SpatialRef) setBoundsOpt(o *boundsOpts) {
	o.sr = sr
}

type boundsOpts struct {
	sr           *SpatialRef
	errorHandler ErrorHandler
}

// BoundsOption is an option that can be passed to Dataset.Bounds or Geometry.Bounds
//
// Available options are:
//  - *SpatialRef
//  - ErrLogger
type BoundsOption interface {
	setBoundsOpt(o *boundsOpts)
}

type createSpatialRefOpts struct {
	errorHandler ErrorHandler
}

// CreateSpatialRefOption is an option that can be passed when creating a new spatial
// reference object
//
// Available options are:
//  - ErrLogger
type CreateSpatialRefOption interface {
	setCreateSpatialRefOpt(so *createSpatialRefOpts)
}

func reprojectBounds(bnds [4]float64, src, dst *SpatialRef) ([4]float64, error) {
	var ret [4]float64
	trn, err := NewTransform(src, dst)
	if err != nil {
		return ret, fmt.Errorf("create coordinate transform: %w", err)
	}
	defer trn.Close()
	x := []float64{bnds[0], bnds[0], bnds[2], bnds[2]}
	y := []float64{bnds[1], bnds[3], bnds[3], bnds[1]}
	err = trn.TransformEx(x, y, nil, nil)
	if err != nil {
		return ret, fmt.Errorf("reproject bounds: %w", err)
	}
	ret[0] = x[0]
	ret[1] = y[0]
	ret[2] = x[0]
	ret[3] = y[0]
	for i := 1; i < 4; i++ {
		if x[i] < ret[0] {
			ret[0] = x[i]
		}
		if x[i] > ret[2] {
			ret[2] = x[i]
		}
		if y[i] < ret[1] {
			ret[1] = y[i]
		}
		if y[i] > ret[3] {
			ret[3] = y[i]
		}
	}
	return ret, nil
}