mirror of
https://github.com/ClickHouse/ClickHouse.git
synced 2024-10-01 14:10:50 +00:00
fd85fffa0d
* Fix compile with boost 1.65.1+ and clang 3.8 ( https://bugs.freebsd.org/bugzilla/show_bug.cgi?id=222439 ) * Fix compile CallPointInPolygon on clang4 (Thanks to @vavrusa)
286 lines
10 KiB
C++
286 lines
10 KiB
C++
#include <Functions/FunctionFactory.h>
|
|
#include <Functions/FunctionsGeo.h>
|
|
#include <Functions/GeoUtils.h>
|
|
#include <Functions/ObjectPool.h>
|
|
|
|
#include <boost/geometry.hpp>
|
|
#include <boost/geometry/geometries/point_xy.hpp>
|
|
#include <boost/geometry/geometries/polygon.hpp>
|
|
|
|
#include <Interpreters/ExpressionActions.h>
|
|
#include <DataTypes/DataTypeTuple.h>
|
|
#include <Columns/ColumnTuple.h>
|
|
#include <IO/WriteHelpers.h>
|
|
#include <DataTypes/DataTypeArray.h>
|
|
#include <Columns/ColumnArray.h>
|
|
#include <Common/ProfileEvents.h>
|
|
|
|
namespace ProfileEvents
|
|
{
|
|
extern const Event PolygonsAddedToPool;
|
|
extern const Event PolygonsInPoolAllocatedBytes;
|
|
}
|
|
|
|
namespace DB
|
|
{
|
|
|
|
namespace ErrorCodes
|
|
{
|
|
extern const int TOO_LESS_ARGUMENTS_FOR_FUNCTION;
|
|
extern const int BAD_ARGUMENTS;
|
|
extern const int ILLEGAL_TYPE_OF_ARGUMENT;
|
|
}
|
|
|
|
namespace FunctionPointInPolygonDetail
|
|
{
|
|
|
|
template <typename Polygon, typename PointInPolygonImpl>
|
|
ColumnPtr callPointInPolygonImplWithPool(const IColumn & x, const IColumn & y, Polygon & polygon)
|
|
{
|
|
using Pool = ObjectPoolMap<PointInPolygonImpl, std::string>;
|
|
/// C++11 has thread-safe function-local statics on most modern compilers.
|
|
static Pool known_polygons;
|
|
|
|
auto factory = [& polygon]()
|
|
{
|
|
GeoUtils::normalizePolygon(polygon);
|
|
auto ptr = std::make_unique<PointInPolygonImpl>(polygon);
|
|
|
|
/// To allocate memory.
|
|
ptr->init();
|
|
|
|
ProfileEvents::increment(ProfileEvents::PolygonsAddedToPool);
|
|
ProfileEvents::increment(ProfileEvents::PolygonsInPoolAllocatedBytes, ptr->getAllocatedBytes());
|
|
|
|
return ptr.release();
|
|
};
|
|
|
|
std::string serialized_polygon = GeoUtils::serialize(polygon);
|
|
auto impl = known_polygons.get(serialized_polygon, factory);
|
|
|
|
return GeoUtils::pointInPolygon(x, y, *impl);
|
|
}
|
|
|
|
template <typename Polygon, typename PointInPolygonImpl>
|
|
ColumnPtr callPointInPolygonImpl(const IColumn & x, const IColumn & y, Polygon & polygon)
|
|
{
|
|
PointInPolygonImpl impl(polygon);
|
|
return GeoUtils::pointInPolygon(x, y, impl);
|
|
}
|
|
|
|
}
|
|
|
|
template <template <typename> typename PointInPolygonImpl, bool use_object_pool = false>
|
|
class FunctionPointInPolygon : public IFunction
|
|
{
|
|
public:
|
|
|
|
template <typename Type>
|
|
using Point = boost::geometry::model::d2::point_xy<Type>;
|
|
template <typename Type>
|
|
using Polygon = boost::geometry::model::polygon<Point<Type>, false>;
|
|
template <typename Type>
|
|
using Box = boost::geometry::model::box<Point<Type>>;
|
|
|
|
static const char * name;
|
|
|
|
static FunctionPtr create(const Context & context)
|
|
{
|
|
return std::make_shared<FunctionPointInPolygon<PointInPolygonImpl, use_object_pool>>();
|
|
}
|
|
|
|
String getName() const override
|
|
{
|
|
return name;
|
|
}
|
|
|
|
bool isVariadic() const override
|
|
{
|
|
return true;
|
|
}
|
|
|
|
size_t getNumberOfArguments() const override
|
|
{
|
|
return 0;
|
|
}
|
|
|
|
DataTypePtr getReturnTypeImpl(const DataTypes & arguments) const override
|
|
{
|
|
if (arguments.size() < 2)
|
|
{
|
|
throw Exception("Too few arguments", ErrorCodes::TOO_LESS_ARGUMENTS_FOR_FUNCTION);
|
|
}
|
|
|
|
auto getMsgPrefix = [this](size_t i) { return "Argument " + toString(i + 1) + " for function " + getName(); };
|
|
|
|
for (size_t i = 1; i < arguments.size(); ++i)
|
|
{
|
|
auto * array = checkAndGetDataType<DataTypeArray>(arguments[i].get());
|
|
if (array == nullptr && i != 1)
|
|
throw Exception(getMsgPrefix(i) + " must be array of tuples.", ErrorCodes::ILLEGAL_TYPE_OF_ARGUMENT);
|
|
|
|
auto * tuple = checkAndGetDataType<DataTypeTuple>(array ? array->getNestedType().get() : arguments[i].get());
|
|
if (tuple == nullptr)
|
|
throw Exception(getMsgPrefix(i) + " must contains tuple.", ErrorCodes::ILLEGAL_TYPE_OF_ARGUMENT);
|
|
|
|
const DataTypes & elements = tuple->getElements();
|
|
|
|
if (elements.size() != 2)
|
|
throw Exception(getMsgPrefix(i) + " must have exactly two elements.", ErrorCodes::BAD_ARGUMENTS);
|
|
|
|
for (auto j : ext::range(0, elements.size()))
|
|
{
|
|
if (!elements[j]->isNumeric())
|
|
{
|
|
throw Exception(getMsgPrefix(i) + " must contains numeric tuple at position " + toString(j + 1),
|
|
ErrorCodes::ILLEGAL_TYPE_OF_ARGUMENT);
|
|
}
|
|
}
|
|
}
|
|
|
|
return std::make_shared<DataTypeUInt8>();
|
|
}
|
|
|
|
void executeImpl(Block & block, const ColumnNumbers & arguments, size_t result) override
|
|
{
|
|
|
|
const IColumn * point_col = block.getByPosition(arguments[0]).column.get();
|
|
auto const_tuple_col = checkAndGetColumn<ColumnConst>(point_col);
|
|
if (const_tuple_col)
|
|
point_col = &const_tuple_col->getDataColumn();
|
|
auto tuple_col = checkAndGetColumn<ColumnTuple>(point_col);
|
|
|
|
if (!tuple_col)
|
|
{
|
|
throw Exception("First argument for function " + getName() + " must be constant array of tuples.",
|
|
ErrorCodes::ILLEGAL_COLUMN);
|
|
}
|
|
|
|
const auto & tuple_block = tuple_col->getData();
|
|
const auto & x = tuple_block.safeGetByPosition(0);
|
|
const auto & y = tuple_block.safeGetByPosition(1);
|
|
|
|
bool use_float64 = checkDataType<DataTypeFloat64>(x.type.get()) || checkDataType<DataTypeFloat64>(y.type.get());
|
|
|
|
auto & result_column = block.safeGetByPosition(result).column;
|
|
|
|
if (use_float64)
|
|
result_column = executeForType<Float64>(*x.column, *y.column, block, arguments);
|
|
else
|
|
result_column = executeForType<Float32>(*x.column, *y.column, block, arguments);
|
|
|
|
if (const_tuple_col)
|
|
result_column = std::make_shared<ColumnConst>(result_column, const_tuple_col->size());
|
|
}
|
|
|
|
private:
|
|
|
|
Float64 getCoordinateFromField(const Field & field)
|
|
{
|
|
switch (field.getType())
|
|
{
|
|
case Field::Types::Float64:
|
|
return field.get<Float64>();
|
|
case Field::Types::Int64:
|
|
return field.get<Int64>();
|
|
case Field::Types::UInt64:
|
|
return field.get<UInt64>();
|
|
default:
|
|
{
|
|
std::string msg = "Expected numeric field, but got ";
|
|
throw Exception(msg + Field::Types::toString(field.getType()), ErrorCodes::LOGICAL_ERROR);
|
|
}
|
|
}
|
|
}
|
|
|
|
template <typename Type>
|
|
ColumnPtr executeForType(const IColumn & x, const IColumn & y, Block & block, const ColumnNumbers & arguments)
|
|
{
|
|
Polygon<Type> polygon;
|
|
|
|
auto getMsgPrefix = [this](size_t i) { return "Argument " + toString(i + 1) + " for function " + getName(); };
|
|
|
|
for (size_t i = 1; i < arguments.size(); ++i)
|
|
{
|
|
auto const_col = checkAndGetColumn<ColumnConst>(block.getByPosition(arguments[i]).column.get());
|
|
auto array_col = const_col ? checkAndGetColumn<ColumnArray>(&const_col->getDataColumn()) : nullptr;
|
|
auto tuple_col = array_col ? checkAndGetColumn<ColumnTuple>(&array_col->getData()) : nullptr;
|
|
|
|
if (!tuple_col)
|
|
throw Exception(getMsgPrefix(i) + " must be constant array of tuples.", ErrorCodes::ILLEGAL_COLUMN);
|
|
|
|
const auto & tuple_block = tuple_col->getData();
|
|
const auto & column_x = tuple_block.safeGetByPosition(0).column;
|
|
const auto & column_y = tuple_block.safeGetByPosition(1).column;
|
|
|
|
if (!polygon.outer().empty())
|
|
polygon.inners().emplace_back();
|
|
|
|
auto & container = polygon.outer().empty() ? polygon.outer() : polygon.inners().back();
|
|
|
|
auto size = column_x->size();
|
|
|
|
if (size == 0)
|
|
throw Exception(getMsgPrefix(i) + " shouldn't be empty.", ErrorCodes::ILLEGAL_COLUMN);
|
|
|
|
for (auto j : ext::range(0, size))
|
|
{
|
|
Type x_coord = getCoordinateFromField((*column_x)[j]);
|
|
Type y_coord = getCoordinateFromField((*column_y)[j]);
|
|
container.push_back(Point<Type>(x_coord, y_coord));
|
|
}
|
|
|
|
/// Polygon assumed to be closed. Allow user to escape repeating of first point.
|
|
if (!boost::geometry::equals(container.front(), container.back()))
|
|
container.push_back(container.front());
|
|
}
|
|
|
|
auto callImpl = use_object_pool
|
|
? FunctionPointInPolygonDetail::callPointInPolygonImplWithPool<Polygon<Type>, PointInPolygonImpl<Type>>
|
|
: FunctionPointInPolygonDetail::callPointInPolygonImpl<Polygon<Type>, PointInPolygonImpl<Type>>;
|
|
|
|
return callImpl(x, y, polygon);
|
|
}
|
|
|
|
};
|
|
|
|
template <typename Type>
|
|
using Point = boost::geometry::model::d2::point_xy<Type>;
|
|
|
|
template <typename Type>
|
|
using PointInPolygonCrossingStrategy = boost::geometry::strategy::within::crossings_multiply<Point<Type>>;
|
|
template <typename Type>
|
|
using PointInPolygonWindingStrategy = boost::geometry::strategy::within::winding<Point<Type>>;
|
|
template <typename Type>
|
|
using PointInPolygonFranklinStrategy = boost::geometry::strategy::within::franklin<Point<Type>>;
|
|
|
|
template <typename Type>
|
|
using PointInPolygonCrossing = GeoUtils::PointInPolygon<PointInPolygonCrossingStrategy<Type>, Type>;
|
|
template <typename Type>
|
|
using PointInPolygonWinding = GeoUtils::PointInPolygon<PointInPolygonWindingStrategy<Type>, Type>;
|
|
template <typename Type>
|
|
using PointInPolygonFranklin = GeoUtils::PointInPolygon<PointInPolygonFranklinStrategy<Type>, Type>;
|
|
template <typename Type>
|
|
using PointInPolygonWithGrid = GeoUtils::PointInPolygonWithGrid<Type>;
|
|
|
|
template <>
|
|
const char * FunctionPointInPolygon<PointInPolygonCrossing>::name = "pointInPolygonCrossing";
|
|
template <>
|
|
const char * FunctionPointInPolygon<PointInPolygonWinding>::name = "pointInPolygonWinding";
|
|
template <>
|
|
const char * FunctionPointInPolygon<PointInPolygonFranklin>::name = "pointInPolygonFranklin";
|
|
template <>
|
|
const char * FunctionPointInPolygon<PointInPolygonWithGrid, true>::name = "pointInPolygon";
|
|
|
|
void registerFunctionsGeo(FunctionFactory & factory)
|
|
{
|
|
factory.registerFunction<FunctionGreatCircleDistance>();
|
|
factory.registerFunction<FunctionPointInEllipses>();
|
|
|
|
factory.registerFunction<FunctionPointInPolygon<PointInPolygonFranklin>>();
|
|
factory.registerFunction<FunctionPointInPolygon<PointInPolygonWinding>>();
|
|
factory.registerFunction<FunctionPointInPolygon<PointInPolygonCrossing>>();
|
|
factory.registerFunction<FunctionPointInPolygon<PointInPolygonWithGrid, true>>();
|
|
}
|
|
}
|