void Polygon::initializeFromBounds(const BOX3D& box)
{
GEOSCoordSequence* coords = GEOSCoordSeq_create_r(m_ctx, 5, 3);
auto set_coordinate = [coords, this](int pt_num, const double&x,
const double& y, const double& z)
{
if (!GEOSCoordSeq_setX_r(m_ctx, coords, pt_num, x))
throw pdal_error("unable to set x for coordinate sequence");
if (!GEOSCoordSeq_setY_r(m_ctx, coords, pt_num, y))
throw pdal_error("unable to set y for coordinate sequence");
if (!GEOSCoordSeq_setZ_r(m_ctx, coords, pt_num, z))
throw pdal_error("unable to set z for coordinate sequence");
};
set_coordinate(0, box.minx, box.miny, box.minz);
set_coordinate(1, box.minx, box.maxy, box.minz);
set_coordinate(2, box.maxx, box.maxy, box.maxz);
set_coordinate(3, box.maxx, box.miny, box.maxz);
set_coordinate(4, box.minx, box.miny, box.minz);
GEOSGeometry* ring = GEOSGeom_createLinearRing_r(m_ctx, coords);
if (!ring)
throw pdal_error("unable to create linear ring from BOX2D");
m_geom = GEOSGeom_createPolygon_r(m_ctx, ring, 0, 0);
if (!m_geom)
throw pdal_error("unable to create polygon from linear ring in "
"BOX2D constructor");
prepare();
}
static VALUE impl_copy_from(VALUE klass, VALUE factory, VALUE original, char subtype)
{
VALUE result;
const GEOSGeometry* original_geom;
GEOSContextHandle_t context;
const GEOSCoordSequence* original_coord_seq;
GEOSCoordSequence* coord_seq;
GEOSGeometry* geom;
result = Qnil;
original_geom = RGEO_GEOMETRY_DATA_PTR(original)->geom;
if (original_geom) {
context = RGEO_FACTORY_DATA_PTR(factory)->geos_context;
if (subtype == 1 && GEOSGetNumCoordinates_r(context, original_geom) != 2) {
original_geom = NULL;
}
if (original_geom) {
original_coord_seq = GEOSGeom_getCoordSeq_r(context, original_geom);
if (original_coord_seq) {
coord_seq = GEOSCoordSeq_clone_r(context, original_coord_seq);
if (coord_seq) {
geom = subtype == 2 ? GEOSGeom_createLinearRing_r(context, coord_seq) : GEOSGeom_createLineString_r(context, coord_seq);
if (geom) {
result = rgeo_wrap_geos_geometry(factory, geom, klass);
}
}
}
}
}
return result;
}
static VALUE cmethod_create_linear_ring(VALUE module, VALUE factory, VALUE array)
{
VALUE result = Qnil;
GEOSCoordSequence* coord_seq = coord_seq_from_array(factory, array, 1);
if (coord_seq) {
RGeo_FactoryData* factory_data = RGEO_FACTORY_DATA_PTR(factory);
GEOSGeometry* geom = GEOSGeom_createLinearRing_r(factory_data->geos_context, coord_seq);
if (geom) {
result = rgeo_wrap_geos_geometry(factory, geom, factory_data->globals->geos_linear_ring);
}
}
return result;
}
