// Copyright (c) 2017, 2024, Oracle and/or its affiliates. // // This program is free software; you can redistribute it and/or modify // it under the terms of the GNU General Public License, version 2.0, // as published by the Free Software Foundation. // // This program is designed to work with certain software (including // but not limited to OpenSSL) that is licensed under separate terms, // as designated in a particular file or component or in included license // documentation. The authors of MySQL hereby grant you an additional // permission to link the program and your derivative works with the // separately licensed software that they have either included with // the program or referenced in the documentation. // // This program is distributed in the hope that it will be useful, // but WITHOUT ANY WARRANTY; without even the implied warranty of // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the // GNU General Public License, version 2.0, for more details. // // You should have received a copy of the GNU General Public License // along with this program; if not, write to the Free Software // Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA 02110-1301 USA. /// @file /// /// This file implements the length functor and function. #include // std::isfinite #include #include "my_inttypes.h" // MYF #include "my_sys.h" // my_error #include "mysqld_error.h" // Error codes #include "sql/dd/types/spatial_reference_system.h" // dd::Spatial_reference_system #include "sql/gis/geometries.h" #include "sql/gis/geometries_traits.h" #include "sql/gis/length.h" #include "sql/gis/length_functor.h" #include "sql/sql_exception_handler.h" // handle_gis_exception namespace bg = boost::geometry; namespace bgs = boost::geometry::srs; namespace bgsd = boost::geometry::strategy::distance; namespace gis { Length::Length(double major, double minor) : m_geographic_strategy(new bgsd::andoyer>( bgs::spheroid(major, minor))) {} double Length::operator()(const Geometry &g1) const { return apply(*this, g1); } double Length::eval(const Geometry &) const { assert(false); /* purecov: inspected */ throw std::exception(); } double Length::eval(const Geographic_linestring &g1) const { return bg::length(g1, *m_geographic_strategy); } double Length::eval(const Cartesian_linestring &g1) const { return bg::length(g1); } ///////////////////////////////////////////////////////////////////////////// double Length::eval(const Geographic_multilinestring &g1) const { return bg::length(g1, *m_geographic_strategy); } double Length::eval(const Cartesian_multilinestring &g1) const { return bg::length(g1); } ///////////////////////////////////////////////////////////////////////////// bool length(const dd::Spatial_reference_system *srs, const Geometry *g1, double *length, bool *null) noexcept { try { assert(srs == nullptr || ((srs->is_cartesian() && g1->coordinate_system() == Coordinate_system::kCartesian) || (srs->is_geographic() && g1->coordinate_system() == Coordinate_system::kGeographic))); if ((*null = (g1->is_empty()))) return false; if ((g1->type() != Geometry_type::kLinestring) && (g1->type() != Geometry_type::kMultilinestring)) { *null = true; return false; } Length len(srs ? srs->semi_major_axis() : 0.0, srs ? srs->semi_minor_axis() : 0.0); *length = len(*g1); } catch (...) { handle_gis_exception("st_length"); return true; } if (!std::isfinite(*length) || *length < 0.0) { my_error(ER_DATA_OUT_OF_RANGE, MYF(0), "Length", "st_length"); return true; } return false; } } // namespace gis