// 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 }