PolygonXY.h

/*
 * (C) Copyright 2013 ECMWF.
 *
 * This software is licensed under the terms of the Apache Licence Version 2.0
 * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
 * In applying this licence, ECMWF does not waive the privileges and immunities
 * granted to it by virtue of its status as an intergovernmental organisation
 * nor does it submit to any jurisdiction.
 */

#pragma once

#include <type_traits>

#include "atlas/util/Polygon.h"

namespace atlas {
namespace util {

class PartitionPolygon;
class PolygonXY;
class ListPolygonXY;

//------------------------------------------------------------------------------------------------------

class PolygonXY : public PolygonCoordinates {
public:
    using Vector = ListPolygonXY;

    PolygonXY( const PartitionPolygon& );

    bool contains( const Point2& Pxy ) const override;

private:
    PointLonLat centroid_;
    double inner_radius_squared_{0};
    PointLonLat inner_coordinatesMin_;
    PointLonLat inner_coordinatesMax_;
};


//------------------------------------------------------------------------------------------------------

class ListPolygonXY : public PolygonCoordinates::Vector {
public:
    ListPolygonXY( const PartitionPolygons& partition_polygons ) {
        reserve( partition_polygons.size() );
        for ( auto& partition_polygon : partition_polygons ) {
            emplace_back( new PolygonXY( partition_polygon ) );
        }
    }
};

//------------------------------------------------------------------------------------------------------

}  // namespace util
}  // namespace atlas