You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
206 lines
7.4 KiB
206 lines
7.4 KiB
// Copyright (c) 2017, Oracle and/or its affiliates. All rights reserved.
|
|
//
|
|
// 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 also distributed 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 included with MySQL.
|
|
//
|
|
// 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 difference functor and function.
|
|
|
|
#include <memory> // std::unique_ptr
|
|
|
|
#include <boost/geometry.hpp>
|
|
|
|
#include "sql/gis/difference_functor.h"
|
|
#include "sql/gis/disjoint_functor.h"
|
|
#include "sql/gis/geometries.h"
|
|
#include "sql/gis/geometries_traits.h"
|
|
|
|
namespace bg = boost::geometry;
|
|
|
|
namespace gis {
|
|
|
|
Difference::Difference(double semi_major, double semi_minor)
|
|
: m_semi_major(semi_major),
|
|
m_semi_minor(semi_minor),
|
|
m_geographic_pl_pa_strategy(
|
|
bg::srs::spheroid<double>(semi_major, semi_minor)),
|
|
m_geographic_ll_la_aa_strategy(
|
|
bg::srs::spheroid<double>(semi_major, semi_minor)) {}
|
|
|
|
Geometry *Difference::operator()(const Geometry *g1, const Geometry *g2) const {
|
|
return apply(*this, g1, g2);
|
|
}
|
|
|
|
Geometry *Difference::eval(const Geometry *g1, const Geometry *g2) const {
|
|
DBUG_ASSERT(false);
|
|
throw not_implemented_exception::for_non_projected(*g1, *g2);
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Cartesian_linestring, *)
|
|
|
|
Geometry *Difference::eval(const Cartesian_linestring *g1,
|
|
const Cartesian_multilinestring *g2) const {
|
|
std::unique_ptr<Cartesian_multilinestring> result(
|
|
new Cartesian_multilinestring());
|
|
bg::difference(*g1, *g2, *result);
|
|
return result.release();
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Cartesian_multipoint, *)
|
|
|
|
Geometry *Difference::eval(const Cartesian_multipoint *g1,
|
|
const Cartesian_multipoint *g2) const {
|
|
std::unique_ptr<Cartesian_multipoint> result(new Cartesian_multipoint());
|
|
bg::difference(*g1, *g2, *result);
|
|
return result.release();
|
|
}
|
|
|
|
Geometry *Difference::eval(const Cartesian_multipoint *g1,
|
|
const Cartesian_multilinestring *g2) const {
|
|
std::unique_ptr<Cartesian_multipoint> result(new Cartesian_multipoint());
|
|
bg::difference(*g1, *g2, *result);
|
|
return result.release();
|
|
}
|
|
|
|
Geometry *Difference::eval(const Cartesian_multipoint *g1,
|
|
const Cartesian_multipolygon *g2) const {
|
|
Disjoint disjoint(m_semi_major, m_semi_minor);
|
|
std::set<Cartesian_point> points;
|
|
for (auto &pt : *g1)
|
|
if (disjoint(&pt, g2)) points.insert(pt);
|
|
|
|
std::unique_ptr<Cartesian_multipoint> result(new Cartesian_multipoint());
|
|
for (auto &pt : points) result->push_back(pt);
|
|
|
|
return result.release();
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Cartesian_multilinestring, *)
|
|
|
|
Geometry *Difference::eval(const Cartesian_multilinestring *g1,
|
|
const Cartesian_multilinestring *g2) const {
|
|
std::unique_ptr<Cartesian_multilinestring> result(
|
|
new Cartesian_multilinestring());
|
|
bg::difference(*g1, *g2, *result);
|
|
return result.release();
|
|
}
|
|
|
|
Geometry *Difference::eval(const Cartesian_multilinestring *g1,
|
|
const Cartesian_multipolygon *g2) const {
|
|
std::unique_ptr<Cartesian_multilinestring> result(
|
|
new Cartesian_multilinestring());
|
|
bg::difference(*g1, *g2, *result);
|
|
return result.release();
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Cartesian_multipolygon, *)
|
|
|
|
Geometry *Difference::eval(const Cartesian_multipolygon *g1,
|
|
const Cartesian_multipolygon *g2) const {
|
|
std::unique_ptr<Cartesian_multipolygon> result(new Cartesian_multipolygon());
|
|
bg::difference(*g1, *g2, *result);
|
|
return result.release();
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Geographic_linestring, *)
|
|
|
|
Geometry *Difference::eval(const Geographic_linestring *g1,
|
|
const Geographic_multilinestring *g2) const {
|
|
std::unique_ptr<Geographic_multilinestring> result(
|
|
new Geographic_multilinestring());
|
|
bg::difference(*g1, *g2, *result, m_geographic_ll_la_aa_strategy);
|
|
return result.release();
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Geographic_multipoint, *)
|
|
|
|
Geometry *Difference::eval(const Geographic_multipoint *g1,
|
|
const Geographic_multipoint *g2) const {
|
|
std::unique_ptr<Geographic_multipoint> result(new Geographic_multipoint());
|
|
bg::difference(*g1, *g2, *result); // Default strategy is OK.
|
|
return result.release();
|
|
}
|
|
|
|
Geometry *Difference::eval(const Geographic_multipoint *g1,
|
|
const Geographic_multilinestring *g2) const {
|
|
std::unique_ptr<Geographic_multipoint> result(new Geographic_multipoint());
|
|
bg::difference(*g1, *g2, *result, m_geographic_pl_pa_strategy);
|
|
return result.release();
|
|
}
|
|
|
|
Geometry *Difference::eval(const Geographic_multipoint *g1,
|
|
const Geographic_multipolygon *g2) const {
|
|
Disjoint disjoint(m_semi_major, m_semi_minor);
|
|
std::set<Geographic_point> points;
|
|
for (auto &pt : *g1)
|
|
if (disjoint(&pt, g2)) points.insert(pt);
|
|
|
|
std::unique_ptr<Geographic_multipoint> result(new Geographic_multipoint());
|
|
for (auto &pt : points) result->push_back(pt);
|
|
|
|
return result.release();
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Geographic_multilinestring, *)
|
|
|
|
Geometry *Difference::eval(const Geographic_multilinestring *g1,
|
|
const Geographic_multilinestring *g2) const {
|
|
std::unique_ptr<Geographic_multilinestring> result(
|
|
new Geographic_multilinestring());
|
|
bg::difference(*g1, *g2, *result, m_geographic_ll_la_aa_strategy);
|
|
return result.release();
|
|
}
|
|
|
|
Geometry *Difference::eval(const Geographic_multilinestring *g1,
|
|
const Geographic_multipolygon *g2) const {
|
|
std::unique_ptr<Geographic_multilinestring> result(
|
|
new Geographic_multilinestring());
|
|
bg::difference(*g1, *g2, *result, m_geographic_ll_la_aa_strategy);
|
|
return result.release();
|
|
}
|
|
|
|
//////////////////////////////////////////////////////////////////////////////
|
|
|
|
// difference(Geographic_multipolygon, *)
|
|
|
|
Geometry *Difference::eval(const Geographic_multipolygon *g1,
|
|
const Geographic_multipolygon *g2) const {
|
|
std::unique_ptr<Geographic_multipolygon> result(
|
|
new Geographic_multipolygon());
|
|
bg::difference(*g1, *g2, *result, m_geographic_ll_la_aa_strategy);
|
|
return result.release();
|
|
}
|
|
|
|
} // namespace gis
|
|
|