haskal
/
kicad
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Upgrade Clipper2 to 1.3.0 

Fixes a number of minor inflation issues including slivers when
overlapping points are inflated

Fixes https://gitlab.com/kicad/code/kicad/-/issues/16182
This commit is contained in:
Seth Hillbrand 2023-11-28 16:01:43 -08:00
parent 9e48b388fe
commit daf178b64f
12 changed files with 863 additions and 949 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
libs/kimath/src/geometry/shape_poly_set.cpp
View File

@ -1217,11 +1217,11 @@ void SHAPE_POLY_SET::inflate2( int aAmount, int aCircleSegCount, CORNER_STRATEGY
Paths64 paths; Paths64 paths;
c.Execute( aAmount, paths ); c.Execute( aAmount, paths );
Clipper2Lib::SimplifyPaths( paths, std::abs( aAmount ) * coeff, false ); Clipper2Lib::SimplifyPaths( paths, std::abs( aAmount ) * coeff, true );
Clipper64 c2; Clipper64 c2;
c2.PreserveCollinear = false; c2.PreserveCollinear( false );
c2.ReverseSolution = false; c2.ReverseSolution( false );
c2.AddSubject( paths ); c2.AddSubject( paths );
c2.Execute(ClipType::Union, FillRule::Positive, tree); c2.Execute(ClipType::Union, FillRule::Positive, tree);
} }
@ -1317,8 +1317,8 @@ void SHAPE_POLY_SET::inflateLine2( const SHAPE_LINE_CHAIN& aLine, int aAmount, i
Clipper2Lib::SimplifyPaths( paths2, std::abs( aAmount ) * coeff, false ); Clipper2Lib::SimplifyPaths( paths2, std::abs( aAmount ) * coeff, false );
Clipper64 c2; Clipper64 c2;
c2.PreserveCollinear = false; c2.PreserveCollinear( false );
c2.ReverseSolution = false; c2.ReverseSolution( false );
c2.AddSubject( paths2 ); c2.AddSubject( paths2 );
c2.Execute( ClipType::Union, FillRule::Positive, tree ); c2.Execute( ClipType::Union, FillRule::Positive, tree );
} }

110
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.core.h vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 26 July 2023 * * Date : 24 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : Core Clipper Library structures and functions * * Purpose : Core Clipper Library structures and functions *
@ -19,6 +19,7 @@
#include <algorithm> #include <algorithm>
#include <climits> #include <climits>
#include <numeric> #include <numeric>
#include "clipper2/clipper.version.h"
namespace Clipper2Lib namespace Clipper2Lib
{ {
@ -42,18 +43,27 @@ namespace Clipper2Lib
"Invalid scale (either 0 or too large)"; "Invalid scale (either 0 or too large)";
static const char* non_pair_error = static const char* non_pair_error =
"There must be 2 values for each coordinate"; "There must be 2 values for each coordinate";
static const char* undefined_error =
"There is an undefined error in Clipper2";
#endif #endif
// error codes (2^n) // error codes (2^n)
const int precision_error_i = 1; // non-fatal const int precision_error_i = 1; // non-fatal
const int scale_error_i = 2; // non-fatal const int scale_error_i = 2; // non-fatal
const int non_pair_error_i = 4; // non-fatal const int non_pair_error_i = 4; // non-fatal
const int range_error_i = 64; const int undefined_error_i = 32; // fatal
const int range_error_i = 64;
#ifndef PI #ifndef PI
static const double PI = 3.141592653589793238; static const double PI = 3.141592653589793238;
#endif #endif
static const int MAX_DECIMAL_PRECISION = 8; // see https://github.com/AngusJohnson/Clipper2/discussions/564
#ifdef CLIPPER2_MAX_PRECISION
const int MAX_DECIMAL_PRECISION = CLIPPER2_MAX_PRECISION;
#else
const int MAX_DECIMAL_PRECISION = 8; // see Discussions #564
#endif
static const int64_t MAX_COORD = INT64_MAX >> 2; static const int64_t MAX_COORD = INT64_MAX >> 2;
static const int64_t MIN_COORD = -MAX_COORD; static const int64_t MIN_COORD = -MAX_COORD;
static const int64_t INVALID = INT64_MAX; static const int64_t INVALID = INT64_MAX;
@ -73,6 +83,8 @@ namespace Clipper2Lib
throw Clipper2Exception(scale_error); throw Clipper2Exception(scale_error);
case non_pair_error_i: case non_pair_error_i:
throw Clipper2Exception(non_pair_error); throw Clipper2Exception(non_pair_error);
case undefined_error_i:
throw Clipper2Exception(undefined_error);
case range_error_i: case range_error_i:
throw Clipper2Exception(range_error); throw Clipper2Exception(range_error);
} }
@ -81,6 +93,7 @@ namespace Clipper2Lib
#endif #endif
} }
//By far the most widely used filling rules for polygons are EvenOdd //By far the most widely used filling rules for polygons are EvenOdd
//and NonZero, sometimes called Alternate and Winding respectively. //and NonZero, sometimes called Alternate and Winding respectively.
//https://en.wikipedia.org/wiki/Nonzero-rule //https://en.wikipedia.org/wiki/Nonzero-rule
@ -137,7 +150,7 @@ namespace Clipper2Lib
friend std::ostream& operator<<(std::ostream& os, const Point& point) friend std::ostream& operator<<(std::ostream& os, const Point& point)
{ {
os << point.x << "," << point.y << "," << point.z << " "; os << point.x << "," << point.y << "," << point.z;
return os; return os;
} }
@ -174,7 +187,7 @@ namespace Clipper2Lib
friend std::ostream& operator<<(std::ostream& os, const Point& point) friend std::ostream& operator<<(std::ostream& os, const Point& point)
{ {
os << point.x << "," << point.y << " "; os << point.x << "," << point.y;
return os; return os;
} }
#endif #endif
@ -222,6 +235,14 @@ namespace Clipper2Lib
using Paths64 = std::vector< Path64>; using Paths64 = std::vector< Path64>;
using PathsD = std::vector< PathD>; using PathsD = std::vector< PathD>;
static const Point64 InvalidPoint64 = Point64(
(std::numeric_limits<int64_t>::max)(),
(std::numeric_limits<int64_t>::max)());
static const PointD InvalidPointD = PointD(
(std::numeric_limits<double>::max)(),
(std::numeric_limits<double>::max)());
// Rect ------------------------------------------------------------------------ // Rect ------------------------------------------------------------------------
template <typename T> template <typename T>
@ -237,19 +258,13 @@ namespace Clipper2Lib
T right; T right;
T bottom; T bottom;
Rect() :
left(0),
top(0),
right(0),
bottom(0) {}
Rect(T l, T t, T r, T b) : Rect(T l, T t, T r, T b) :
left(l), left(l),
top(t), top(t),
right(r), right(r),
bottom(b) {} bottom(b) {}
Rect(bool is_valid) Rect(bool is_valid = true)
{ {
if (is_valid) if (is_valid)
{ {
@ -258,10 +273,12 @@ namespace Clipper2Lib
else else
{ {
left = top = (std::numeric_limits<T>::max)(); left = top = (std::numeric_limits<T>::max)();
right = bottom = -(std::numeric_limits<int64_t>::max)(); right = bottom = (std::numeric_limits<T>::lowest)();
} }
} }
bool IsValid() const { return left != (std::numeric_limits<T>::max)(); }
T Width() const { return right - left; } T Width() const { return right - left; }
T Height() const { return bottom - top; } T Height() const { return bottom - top; }
void Width(T width) { right = left + width; } void Width(T width) { right = left + width; }
@ -309,10 +326,13 @@ namespace Clipper2Lib
((std::max)(top, rec.top) <= (std::min)(bottom, rec.bottom)); ((std::max)(top, rec.top) <= (std::min)(bottom, rec.bottom));
}; };
bool operator==(const Rect<T>& other) const {
return left == other.left && right == other.right &&
top == other.top && bottom == other.bottom;
}
friend std::ostream& operator<<(std::ostream& os, const Rect<T>& rect) { friend std::ostream& operator<<(std::ostream& os, const Rect<T>& rect) {
os << "(" os << "(" << rect.left << "," << rect.top << "," << rect.right << "," << rect.bottom << ") ";
<< rect.left << "," << rect.top << "," << rect.right << "," << rect.bottom
<< ")";
return os; return os;
} }
}; };
@ -340,10 +360,16 @@ namespace Clipper2Lib
return result; return result;
} }
static const Rect64 MaxInvalidRect64 = Rect64( static const Rect64 InvalidRect64 = Rect64(
INT64_MAX, INT64_MAX, INT64_MIN, INT64_MIN); (std::numeric_limits<int64_t>::max)(),
static const RectD MaxInvalidRectD = RectD( (std::numeric_limits<int64_t>::max)(),
MAX_DBL, MAX_DBL, -MAX_DBL, -MAX_DBL); (std::numeric_limits<int64_t>::lowest)(),
(std::numeric_limits<int64_t>::lowest)());
static const RectD InvalidRectD = RectD(
(std::numeric_limits<double>::max)(),
(std::numeric_limits<double>::max)(),
(std::numeric_limits<double>::lowest)(),
(std::numeric_limits<double>::lowest)());
template <typename T> template <typename T>
Rect<T> GetBounds(const Path<T>& path) Rect<T> GetBounds(const Path<T>& path)
@ -490,26 +516,6 @@ namespace Clipper2Lib
return result; return result;
} }
inline PathD Path64ToPathD(const Path64& path)
{
return TransformPath<double, int64_t>(path);
}
inline PathsD Paths64ToPathsD(const Paths64& paths)
{
return TransformPaths<double, int64_t>(paths);
}
inline Path64 PathDToPath64(const PathD& path)
{
return TransformPath<int64_t, double>(path);
}
inline Paths64 PathsDToPaths64(const PathsD& paths)
{
return TransformPaths<int64_t, double>(paths);
}
template<typename T> template<typename T>
inline double Sqr(T val) inline double Sqr(T val)
{ {
@ -565,7 +571,7 @@ namespace Clipper2Lib
inline void StripDuplicates( Path<T>& path, bool is_closed_path) inline void StripDuplicates( Path<T>& path, bool is_closed_path)
{ {
//https://stackoverflow.com/questions/1041620/whats-the-most-efficient-way-to-erase-duplicates-and-sort-a-vector#:~:text=Let%27s%20compare%20three%20approaches%3A //https://stackoverflow.com/questions/1041620/whats-the-most-efficient-way-to-erase-duplicates-and-sort-a-vector#:~:text=Let%27s%20compare%20three%20approaches%3A
path.erase(std::unique(path.begin(), path.end()),path.end()); path.erase(std::unique(path.begin(), path.end()), path.end());
if (is_closed_path) if (is_closed_path)
while (path.size() > 1 && path.back() == path.front()) path.pop_back(); while (path.size() > 1 && path.back() == path.front()) path.pop_back();
} }
@ -723,8 +729,9 @@ namespace Clipper2Lib
} }
} }
inline Point64 GetClosestPointOnSegment(const Point64& offPt, template<typename T>
const Point64& seg1, const Point64& seg2) inline Point<T> GetClosestPointOnSegment(const Point<T>& offPt,
const Point<T>& seg1, const Point<T>& seg2)
{ {
if (seg1.x == seg2.x && seg1.y == seg2.y) return seg1; if (seg1.x == seg2.x && seg1.y == seg2.y) return seg1;
double dx = static_cast<double>(seg2.x - seg1.x); double dx = static_cast<double>(seg2.x - seg1.x);
@ -734,9 +741,14 @@ namespace Clipper2Lib
static_cast<double>(offPt.y - seg1.y) * dy) / static_cast<double>(offPt.y - seg1.y) * dy) /
(Sqr(dx) + Sqr(dy)); (Sqr(dx) + Sqr(dy));
if (q < 0) q = 0; else if (q > 1) q = 1; if (q < 0) q = 0; else if (q > 1) q = 1;
return Point64( if constexpr (std::numeric_limits<T>::is_integer)
seg1.x + static_cast<int64_t>(nearbyint(q * dx)), return Point<T>(
seg1.y + static_cast<int64_t>(nearbyint(q * dy))); seg1.x + static_cast<T>(nearbyint(q * dx)),
seg1.y + static_cast<T>(nearbyint(q * dy)));
else
return Point<T>(
seg1.x + static_cast<T>(q * dx),
seg1.y + static_cast<T>(q * dy));
} }
enum class PointInPolygonResult { IsOn, IsInside, IsOutside }; enum class PointInPolygonResult { IsOn, IsInside, IsOutside };

36
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.engine.h vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 26 July 2023 * * Date : 22 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : This is the main polygon clipping module * * Purpose : This is the main polygon clipping module *
@ -10,8 +10,6 @@
#ifndef CLIPPER_ENGINE_H #ifndef CLIPPER_ENGINE_H
#define CLIPPER_ENGINE_H #define CLIPPER_ENGINE_H
constexpr auto CLIPPER2_VERSION = "1.2.2";
#include <cstdlib> #include <cstdlib>
#include <stdint.h> //#541 #include <stdint.h> //#541
#include <iostream> #include <iostream>
@ -21,12 +19,11 @@ constexpr auto CLIPPER2_VERSION = "1.2.2";
#include <numeric> #include <numeric>
#include <memory> #include <memory>
#include "clipper.core.h" #include "clipper2/clipper.core.h"
#ifdef None #ifdef None
#undef None #undef None
#endif #endif
namespace Clipper2Lib { namespace Clipper2Lib {
struct Scanline; struct Scanline;
@ -268,6 +265,8 @@ namespace Clipper2Lib {
inline void CheckJoinRight(Active& e, inline void CheckJoinRight(Active& e,
const Point64& pt, bool check_curr_x = false); const Point64& pt, bool check_curr_x = false);
protected: protected:
bool preserve_collinear_ = true;
bool reverse_solution_ = false;
int error_code_ = 0; int error_code_ = 0;
bool has_open_paths_ = false; bool has_open_paths_ = false;
bool succeeded_ = true; bool succeeded_ = true;
@ -286,9 +285,11 @@ namespace Clipper2Lib {
void AddPaths(const Paths64& paths, PathType polytype, bool is_open); void AddPaths(const Paths64& paths, PathType polytype, bool is_open);
public: public:
virtual ~ClipperBase(); virtual ~ClipperBase();
int ErrorCode() { return error_code_; }; int ErrorCode() const { return error_code_; };
bool PreserveCollinear = true; void PreserveCollinear(bool val) { preserve_collinear_ = val; };
bool ReverseSolution = false; bool PreserveCollinear() const { return preserve_collinear_;};
void ReverseSolution(bool val) { reverse_solution_ = val; };
bool ReverseSolution() const { return reverse_solution_; };
void Clear(); void Clear();
void AddReuseableData(const ReuseableDataContainer64& reuseable_data); void AddReuseableData(const ReuseableDataContainer64& reuseable_data);
#ifdef USINGZ #ifdef USINGZ
@ -350,12 +351,12 @@ namespace Clipper2Lib {
childs_.resize(0); childs_.resize(0);
} }
const PolyPath64* operator [] (size_t index) const PolyPath64* operator [] (size_t index) const
{ {
return childs_[index].get(); //std::unique_ptr return childs_[index].get(); //std::unique_ptr
} }
const PolyPath64* Child(size_t index) const PolyPath64* Child(size_t index) const
{ {
return childs_[index].get(); return childs_[index].get();
} }
@ -407,12 +408,12 @@ namespace Clipper2Lib {
childs_.resize(0); childs_.resize(0);
} }
const PolyPathD* operator [] (size_t index) const PolyPathD* operator [] (size_t index) const
{ {
return childs_[index].get(); return childs_[index].get();
} }
const PolyPathD* Child(size_t index) const PolyPathD* Child(size_t index) const
{ {
return childs_[index].get(); return childs_[index].get();
} }
@ -421,7 +422,8 @@ namespace Clipper2Lib {
PolyPathDList::const_iterator end() const { return childs_.cend(); } PolyPathDList::const_iterator end() const { return childs_.cend(); }
void SetScale(double value) { scale_ = value; } void SetScale(double value) { scale_ = value; }
double Scale() { return scale_; } double Scale() const { return scale_; }
PolyPathD* AddChild(const Path64& path) override PolyPathD* AddChild(const Path64& path) override
{ {
int error_code = 0; int error_code = 0;
@ -431,6 +433,14 @@ namespace Clipper2Lib {
return result; return result;
} }
PolyPathD* AddChild(const PathD& path)
{
auto p = std::make_unique<PolyPathD>(this);
PolyPathD* result = childs_.emplace_back(std::move(p)).get();
result->polygon_ = path;
return result;
}
void Clear() override void Clear() override
{ {
childs_.resize(0); childs_.resize(0);

835
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.export.h vendored
View File

File diff suppressed because it is too large Load Diff

80
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.h vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 16 July 2023 * * Date : 18 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : This module provides a simple interface to the Clipper Library * * Purpose : This module provides a simple interface to the Clipper Library *
@ -14,11 +14,11 @@
#include <type_traits> #include <type_traits>
#include <vector> #include <vector>
#include "clipper.core.h" #include "clipper2/clipper.core.h"
#include "clipper.engine.h" #include "clipper2/clipper.engine.h"
#include "clipper.offset.h" #include "clipper2/clipper.offset.h"
#include "clipper.minkowski.h" #include "clipper2/clipper.minkowski.h"
#include "clipper.rectclip.h" #include "clipper2/clipper.rectclip.h"
namespace Clipper2Lib { namespace Clipper2Lib {
@ -341,6 +341,19 @@ namespace Clipper2Lib {
details::OutlinePolyPathD(os, *pp.Child(i), i, preamble + " "); details::OutlinePolyPathD(os, *pp.Child(i), i, preamble + " ");
} }
template<typename T, typename U>
inline constexpr void MakePathGeneric(const T an_array,
size_t array_size, std::vector<U>& result)
{
result.reserve(array_size / 2);
for (size_t i = 0; i < array_size; i +=2)
#ifdef USINGZ
result.push_back( U{ an_array[i], an_array[i +1], 0} );
#else
result.push_back( U{ an_array[i], an_array[i + 1]} );
#endif
}
} // end details namespace } // end details namespace
inline std::ostream& operator<< (std::ostream& os, const PolyTree64& pp) inline std::ostream& operator<< (std::ostream& os, const PolyTree64& pp)
@ -391,22 +404,6 @@ namespace Clipper2Lib {
return true; return true;
} }
namespace details {
template<typename T, typename U>
inline constexpr void MakePathGeneric(const T list, size_t size,
std::vector<U>& result)
{
for (size_t i = 0; i < size; ++i)
#ifdef USINGZ
result[i / 2] = U{list[i], list[++i], 0};
#else
result[i / 2] = U{list[i], list[++i]};
#endif
}
} // end details namespace
template<typename T, template<typename T,
typename std::enable_if< typename std::enable_if<
std::is_integral<T>::value && std::is_integral<T>::value &&
@ -417,7 +414,7 @@ namespace Clipper2Lib {
const auto size = list.size() - list.size() % 2; const auto size = list.size() - list.size() % 2;
if (list.size() != size) if (list.size() != size)
DoError(non_pair_error_i); // non-fatal without exception handling DoError(non_pair_error_i); // non-fatal without exception handling
Path64 result(size / 2); // else ignores unpaired value Path64 result;
details::MakePathGeneric(list, size, result); details::MakePathGeneric(list, size, result);
return result; return result;
} }
@ -431,7 +428,7 @@ namespace Clipper2Lib {
{ {
// Make the compiler error on unpaired value (i.e. no runtime effects). // Make the compiler error on unpaired value (i.e. no runtime effects).
static_assert(N % 2 == 0, "MakePath requires an even number of arguments"); static_assert(N % 2 == 0, "MakePath requires an even number of arguments");
Path64 result(N / 2); Path64 result;
details::MakePathGeneric(list, N, result); details::MakePathGeneric(list, N, result);
return result; return result;
} }
@ -446,7 +443,7 @@ namespace Clipper2Lib {
const auto size = list.size() - list.size() % 2; const auto size = list.size() - list.size() % 2;
if (list.size() != size) if (list.size() != size)
DoError(non_pair_error_i); // non-fatal without exception handling DoError(non_pair_error_i); // non-fatal without exception handling
PathD result(size / 2); // else ignores unpaired value PathD result;
details::MakePathGeneric(list, size, result); details::MakePathGeneric(list, size, result);
return result; return result;
} }
@ -460,7 +457,7 @@ namespace Clipper2Lib {
{ {
// Make the compiler error on unpaired value (i.e. no runtime effects). // Make the compiler error on unpaired value (i.e. no runtime effects).
static_assert(N % 2 == 0, "MakePath requires an even number of arguments"); static_assert(N % 2 == 0, "MakePath requires an even number of arguments");
PathD result(N / 2); PathD result;
details::MakePathGeneric(list, N, result); details::MakePathGeneric(list, N, result);
return result; return result;
} }
@ -653,7 +650,7 @@ namespace Clipper2Lib {
} }
template <typename T> template <typename T>
inline Path<T> SimplifyPath(const Path<T> path, inline Path<T> SimplifyPath(const Path<T> &path,
double epsilon, bool isClosedPath = true) double epsilon, bool isClosedPath = true)
{ {
const size_t len = path.size(), high = len -1; const size_t len = path.size(), high = len -1;
@ -662,7 +659,7 @@ namespace Clipper2Lib {
std::vector<bool> flags(len); std::vector<bool> flags(len);
std::vector<double> distSqr(len); std::vector<double> distSqr(len);
size_t prior = high, curr = 0, start, next, prior2, next2; size_t prior = high, curr = 0, start, next, prior2;
if (isClosedPath) if (isClosedPath)
{ {
distSqr[0] = PerpendicDistFromLineSqrd(path[0], path[high], path[1]); distSqr[0] = PerpendicDistFromLineSqrd(path[0], path[high], path[1]);
@ -692,26 +689,25 @@ namespace Clipper2Lib {
next = GetNext(curr, high, flags); next = GetNext(curr, high, flags);
if (next == prior) break; if (next == prior) break;
// flag for removal the smaller of adjacent 'distances'
if (distSqr[next] < distSqr[curr]) if (distSqr[next] < distSqr[curr])
{ {
flags[next] = true; prior2 = prior;
next = GetNext(next, high, flags); prior = curr;
next2 = GetNext(next, high, flags);
distSqr[curr] = PerpendicDistFromLineSqrd(path[curr], path[prior], path[next]);
if (next != high || isClosedPath)
distSqr[next] = PerpendicDistFromLineSqrd(path[next], path[curr], path[next2]);
curr = next; curr = next;
next = GetNext(next, high, flags);
} }
else else
{
flags[curr] = true;
curr = next;
next = GetNext(next, high, flags);
prior2 = GetPrior(prior, high, flags); prior2 = GetPrior(prior, high, flags);
flags[curr] = true;
curr = next;
next = GetNext(next, high, flags);
if (isClosedPath || ((curr != high) && (curr != 0)))
distSqr[curr] = PerpendicDistFromLineSqrd(path[curr], path[prior], path[next]); distSqr[curr] = PerpendicDistFromLineSqrd(path[curr], path[prior], path[next]);
if (prior != 0 || isClosedPath) if (isClosedPath || ((prior != 0) && (prior != high)))
distSqr[prior] = PerpendicDistFromLineSqrd(path[prior], path[prior2], path[curr]); distSqr[prior] = PerpendicDistFromLineSqrd(path[prior], path[prior2], path[curr]);
}
} }
Path<T> result; Path<T> result;
result.reserve(len); result.reserve(len);
@ -721,7 +717,7 @@ namespace Clipper2Lib {
} }
template <typename T> template <typename T>
inline Paths<T> SimplifyPaths(const Paths<T> paths, inline Paths<T> SimplifyPaths(const Paths<T> &paths,
double epsilon, bool isClosedPath = true) double epsilon, bool isClosedPath = true)
{ {
Paths<T> result; Paths<T> result;

4
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.minkowski.h vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 28 January 2023 * * Date : 1 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : Minkowski Sum and Difference * * Purpose : Minkowski Sum and Difference *
@ -13,7 +13,7 @@
#include <cstdlib> #include <cstdlib>
#include <vector> #include <vector>
#include <string> #include <string>
#include "clipper.core.h" #include "clipper2/clipper.core.h"
namespace Clipper2Lib namespace Clipper2Lib
{ {

34
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.offset.h vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 15 May 2023 * * Date : 19 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : Path Offset (Inflate/Shrink) * * Purpose : Path Offset (Inflate/Shrink) *
@ -15,7 +15,9 @@
namespace Clipper2Lib { namespace Clipper2Lib {
enum class JoinType { Square, Round, Miter }; enum class JoinType { Square, Bevel, Round, Miter };
//Square : Joins are 'squared' at exactly the offset distance (more complex code)
//Bevel : Similar to Square, but the offset distance varies with angle (simple code & faster)
enum class EndType {Polygon, Joined, Butt, Square, Round}; enum class EndType {Polygon, Joined, Butt, Square, Round};
//Butt : offsets both sides of a path, with square blunt ends //Butt : offsets both sides of a path, with square blunt ends
@ -32,13 +34,13 @@ private:
class Group { class Group {
public: public:
Paths64 paths_in; Paths64 paths_in;
Paths64 paths_out; std::vector<bool> is_hole_list;
Path64 path; std::vector<Rect64> bounds_list;
int lowest_path_idx = -1;
bool is_reversed = false; bool is_reversed = false;
JoinType join_type; JoinType join_type;
EndType end_type; EndType end_type;
Group(const Paths64& _paths, JoinType _join_type, EndType _end_type) : Group(const Paths64& _paths, JoinType _join_type, EndType _end_type);
paths_in(_paths), join_type(_join_type), end_type(_end_type) {}
}; };
int error_code_ = 0; int error_code_ = 0;
@ -49,9 +51,10 @@ private:
double step_sin_ = 0.0; double step_sin_ = 0.0;
double step_cos_ = 0.0; double step_cos_ = 0.0;
PathD norms; PathD norms;
Path64 path_out;
Paths64 solution; Paths64 solution;
std::vector<Group> groups_; std::vector<Group> groups_;
JoinType join_type_ = JoinType::Square; JoinType join_type_ = JoinType::Bevel;
EndType end_type_ = EndType::Polygon; EndType end_type_ = EndType::Polygon;
double miter_limit_ = 0.0; double miter_limit_ = 0.0;
@ -64,14 +67,17 @@ private:
#endif #endif
DeltaCallback64 deltaCallback64_ = nullptr; DeltaCallback64 deltaCallback64_ = nullptr;
void DoSquare(Group& group, const Path64& path, size_t j, size_t k); size_t CalcSolutionCapacity();
void DoMiter(Group& group, const Path64& path, size_t j, size_t k, double cos_a); bool CheckReverseOrientation();
void DoRound(Group& group, const Path64& path, size_t j, size_t k, double angle); void DoBevel(const Path64& path, size_t j, size_t k);
void DoSquare(const Path64& path, size_t j, size_t k);
void DoMiter(const Path64& path, size_t j, size_t k, double cos_a);
void DoRound(const Path64& path, size_t j, size_t k, double angle);
void BuildNormals(const Path64& path); void BuildNormals(const Path64& path);
void OffsetPolygon(Group& group, Path64& path); void OffsetPolygon(Group& group, const Path64& path);
void OffsetOpenJoined(Group& group, Path64& path); void OffsetOpenJoined(Group& group, const Path64& path);
void OffsetOpenPath(Group& group, Path64& path); void OffsetOpenPath(Group& group, const Path64& path);
void OffsetPoint(Group& group, Path64& path, size_t j, size_t k); void OffsetPoint(Group& group, const Path64& path, size_t j, size_t k);
void DoGroupOffset(Group &group); void DoGroupOffset(Group &group);
void ExecuteInternal(double delta); void ExecuteInternal(double delta);
public: public:

5
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.rectclip.h vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 30 May 2023 * * Date : 1 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : FAST rectangular clipping * * Purpose : FAST rectangular clipping *
@ -13,8 +13,7 @@
#include <cstdlib> #include <cstdlib>
#include <vector> #include <vector>
#include <queue> #include <queue>
#include "clipper.h" #include "clipper2/clipper.core.h"
#include "clipper.core.h"
namespace Clipper2Lib namespace Clipper2Lib
{ {

6
thirdparty/clipper2/Clipper2Lib/include/clipper2/clipper.version.h vendored Normal file
View File

@ -0,0 +1,6 @@
#ifndef CLIPPER_VERSION_H
#define CLIPPER_VERSION_H
constexpr auto CLIPPER2_VERSION = "1.3.0";
#endif // CLIPPER_VERSION_H

99
thirdparty/clipper2/Clipper2Lib/src/clipper.engine.cpp vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 27 August 2023 * * Date : 22 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : This is the main polygon clipping module * * Purpose : This is the main polygon clipping module *
@ -15,6 +15,7 @@
#include <algorithm> #include <algorithm>
#include "clipper2/clipper.engine.h" #include "clipper2/clipper.engine.h"
#include "clipper2/clipper.h"
// https://github.com/AngusJohnson/Clipper2/discussions/334 // https://github.com/AngusJohnson/Clipper2/discussions/334
// #discussioncomment-4248602 // #discussioncomment-4248602
@ -1469,13 +1470,14 @@ namespace Clipper2Lib {
e2.outrec->front_edge = nullptr; e2.outrec->front_edge = nullptr;
e2.outrec->back_edge = nullptr; e2.outrec->back_edge = nullptr;
e2.outrec->pts = nullptr; e2.outrec->pts = nullptr;
SetOwner(e2.outrec, e1.outrec);
if (IsOpenEnd(e1)) if (IsOpenEnd(e1))
{ {
e2.outrec->pts = e1.outrec->pts; e2.outrec->pts = e1.outrec->pts;
e1.outrec->pts = nullptr; e1.outrec->pts = nullptr;
} }
else
SetOwner(e2.outrec, e1.outrec);
//and e1 and e2 are maxima and are about to be dropped from the Actives list. //and e1 and e2 are maxima and are about to be dropped from the Actives list.
e1.outrec = nullptr; e1.outrec = nullptr;
@ -1541,7 +1543,7 @@ namespace Clipper2Lib {
//NB if preserveCollinear == true, then only remove 180 deg. spikes //NB if preserveCollinear == true, then only remove 180 deg. spikes
if ((CrossProduct(op2->prev->pt, op2->pt, op2->next->pt) == 0) && if ((CrossProduct(op2->prev->pt, op2->pt, op2->next->pt) == 0) &&
(op2->pt == op2->prev->pt || (op2->pt == op2->prev->pt ||
op2->pt == op2->next->pt || !PreserveCollinear || op2->pt == op2->next->pt || !preserve_collinear_ ||
DotProduct(op2->prev->pt, op2->pt, op2->next->pt) < 0)) DotProduct(op2->prev->pt, op2->pt, op2->next->pt) < 0))
{ {
@ -1706,6 +1708,28 @@ namespace Clipper2Lib {
return op; return op;
} }
inline void TrimHorz(Active& horzEdge, bool preserveCollinear)
{
bool wasTrimmed = false;
Point64 pt = NextVertex(horzEdge)->pt;
while (pt.y == horzEdge.top.y)
{
//always trim 180 deg. spikes (in closed paths)
//but otherwise break if preserveCollinear = true
if (preserveCollinear &&
((pt.x < horzEdge.top.x) != (horzEdge.bot.x < horzEdge.top.x)))
break;
horzEdge.vertex_top = NextVertex(horzEdge);
horzEdge.top = pt;
wasTrimmed = true;
if (IsMaxima(horzEdge)) break;
pt = NextVertex(horzEdge)->pt;
}
if (wasTrimmed) SetDx(horzEdge); // +/-infinity
}
inline void ClipperBase::UpdateEdgeIntoAEL(Active* e) inline void ClipperBase::UpdateEdgeIntoAEL(Active* e)
{ {
@ -1717,9 +1741,13 @@ namespace Clipper2Lib {
if (IsJoined(*e)) Split(*e, e->bot); if (IsJoined(*e)) Split(*e, e->bot);
if (IsHorizontal(*e)) return; if (IsHorizontal(*e))
InsertScanline(e->top.y); {
if (!IsOpen(*e)) TrimHorz(*e, preserve_collinear_);
return;
}
InsertScanline(e->top.y);
CheckJoinLeft(*e, e->bot); CheckJoinLeft(*e, e->bot);
CheckJoinRight(*e, e->bot, true); // (#500) CheckJoinRight(*e, e->bot, true); // (#500)
} }
@ -2139,7 +2167,7 @@ namespace Clipper2Lib {
horz_seg_list_.end(), horz_seg_list_.end(),
[](HorzSegment& hs) { return UpdateHorzSegment(hs); }); [](HorzSegment& hs) { return UpdateHorzSegment(hs); });
if (j < 2) return; if (j < 2) return;
std::sort(horz_seg_list_.begin(), horz_seg_list_.end(), HorzSegSorter()); std::stable_sort(horz_seg_list_.begin(), horz_seg_list_.end(), HorzSegSorter());
HorzSegmentList::iterator hs1 = horz_seg_list_.begin(), hs2; HorzSegmentList::iterator hs1 = horz_seg_list_.begin(), hs2;
HorzSegmentList::iterator hs_end = hs1 +j; HorzSegmentList::iterator hs_end = hs1 +j;
@ -2451,35 +2479,6 @@ namespace Clipper2Lib {
} }
} }
inline bool HorzIsSpike(const Active& horzEdge)
{
Point64 nextPt = NextVertex(horzEdge)->pt;
return (nextPt.y == horzEdge.bot.y) &&
(horzEdge.bot.x < horzEdge.top.x) != (horzEdge.top.x < nextPt.x);
}
inline void TrimHorz(Active& horzEdge, bool preserveCollinear)
{
bool wasTrimmed = false;
Point64 pt = NextVertex(horzEdge)->pt;
while (pt.y == horzEdge.top.y)
{
//always trim 180 deg. spikes (in closed paths)
//but otherwise break if preserveCollinear = true
if (preserveCollinear &&
((pt.x < horzEdge.top.x) != (horzEdge.bot.x < horzEdge.top.x)))
break;
horzEdge.vertex_top = NextVertex(horzEdge);
horzEdge.top = pt;
wasTrimmed = true;
if (IsMaxima(horzEdge)) break;
pt = NextVertex(horzEdge)->pt;
}
if (wasTrimmed) SetDx(horzEdge); // +/-infinity
}
void ClipperBase::DoHorizontal(Active& horz) void ClipperBase::DoHorizontal(Active& horz)
/******************************************************************************* /*******************************************************************************
* Notes: Horizontal edges (HEs) at scanline intersections (ie at the top or * * Notes: Horizontal edges (HEs) at scanline intersections (ie at the top or *
@ -2505,10 +2504,10 @@ namespace Clipper2Lib {
else else
vertex_max = GetCurrYMaximaVertex(horz); vertex_max = GetCurrYMaximaVertex(horz);
// remove 180 deg.spikes and also simplify //// remove 180 deg.spikes and also simplify
// consecutive horizontals when PreserveCollinear = true //// consecutive horizontals when PreserveCollinear = true
if (vertex_max && !horzIsOpen && vertex_max != horz.vertex_top) //if (!horzIsOpen && vertex_max != horz.vertex_top)
TrimHorz(horz, PreserveCollinear); // TrimHorz(horz, PreserveCollinear);
int64_t horz_left, horz_right; int64_t horz_left, horz_right;
bool is_left_to_right = bool is_left_to_right =
@ -2537,6 +2536,9 @@ namespace Clipper2Lib {
if (IsHotEdge(horz) && IsJoined(*e)) if (IsHotEdge(horz) && IsJoined(*e))
Split(*e, e->top); Split(*e, e->top);
//if (IsHotEdge(horz) != IsHotEdge(*e))
// DoError(undefined_error_i);
if (IsHotEdge(horz)) if (IsHotEdge(horz))
{ {
while (horz.vertex_top != vertex_max) while (horz.vertex_top != vertex_max)
@ -2591,6 +2593,7 @@ namespace Clipper2Lib {
{ {
IntersectEdges(horz, *e, pt); IntersectEdges(horz, *e, pt);
SwapPositionsInAEL(horz, *e); SwapPositionsInAEL(horz, *e);
CheckJoinLeft(*e, pt);
horz.curr_x = e->curr_x; horz.curr_x = e->curr_x;
e = horz.next_in_ael; e = horz.next_in_ael;
} }
@ -2598,6 +2601,7 @@ namespace Clipper2Lib {
{ {
IntersectEdges(*e, horz, pt); IntersectEdges(*e, horz, pt);
SwapPositionsInAEL(*e, horz); SwapPositionsInAEL(*e, horz);
CheckJoinRight(*e, pt);
horz.curr_x = e->curr_x; horz.curr_x = e->curr_x;
e = horz.prev_in_ael; e = horz.prev_in_ael;
} }
@ -2633,9 +2637,6 @@ namespace Clipper2Lib {
AddOutPt(horz, horz.top); AddOutPt(horz, horz.top);
UpdateEdgeIntoAEL(&horz); UpdateEdgeIntoAEL(&horz);
if (PreserveCollinear && !horzIsOpen && HorzIsSpike(horz))
TrimHorz(horz, true);
is_left_to_right = is_left_to_right =
ResetHorzDirection(horz, vertex_max, horz_left, horz_right); ResetHorzDirection(horz, vertex_max, horz_left, horz_right);
} }
@ -2872,7 +2873,7 @@ namespace Clipper2Lib {
if (!outrec->bounds.IsEmpty()) return true; if (!outrec->bounds.IsEmpty()) return true;
CleanCollinear(outrec); CleanCollinear(outrec);
if (!outrec->pts || if (!outrec->pts ||
!BuildPath64(outrec->pts, ReverseSolution, false, outrec->path)){ !BuildPath64(outrec->pts, reverse_solution_, false, outrec->path)){
return false;} return false;}
outrec->bounds = GetBounds(outrec->path); outrec->bounds = GetBounds(outrec->path);
return true; return true;
@ -2947,7 +2948,7 @@ namespace Clipper2Lib {
Path64 path; Path64 path;
if (solutionOpen && outrec->is_open) if (solutionOpen && outrec->is_open)
{ {
if (BuildPath64(outrec->pts, ReverseSolution, true, path)) if (BuildPath64(outrec->pts, reverse_solution_, true, path))
solutionOpen->emplace_back(std::move(path)); solutionOpen->emplace_back(std::move(path));
} }
else else
@ -2955,7 +2956,7 @@ namespace Clipper2Lib {
// nb: CleanCollinear can add to outrec_list_ // nb: CleanCollinear can add to outrec_list_
CleanCollinear(outrec); CleanCollinear(outrec);
//closed paths should always return a Positive orientation //closed paths should always return a Positive orientation
if (BuildPath64(outrec->pts, ReverseSolution, false, path)) if (BuildPath64(outrec->pts, reverse_solution_, false, path))
solutionClosed.emplace_back(std::move(path)); solutionClosed.emplace_back(std::move(path));
} }
} }
@ -2978,7 +2979,7 @@ namespace Clipper2Lib {
if (outrec->is_open) if (outrec->is_open)
{ {
Path64 path; Path64 path;
if (BuildPath64(outrec->pts, ReverseSolution, true, path)) if (BuildPath64(outrec->pts, reverse_solution_, true, path))
open_paths.push_back(path); open_paths.push_back(path);
continue; continue;
} }
@ -3055,14 +3056,14 @@ namespace Clipper2Lib {
PathD path; PathD path;
if (solutionOpen && outrec->is_open) if (solutionOpen && outrec->is_open)
{ {
if (BuildPathD(outrec->pts, ReverseSolution, true, path, invScale_)) if (BuildPathD(outrec->pts, reverse_solution_, true, path, invScale_))
solutionOpen->emplace_back(std::move(path)); solutionOpen->emplace_back(std::move(path));
} }
else else
{ {
CleanCollinear(outrec); CleanCollinear(outrec);
//closed paths should always return a Positive orientation //closed paths should always return a Positive orientation
if (BuildPathD(outrec->pts, ReverseSolution, false, path, invScale_)) if (BuildPathD(outrec->pts, reverse_solution_, false, path, invScale_))
solutionClosed.emplace_back(std::move(path)); solutionClosed.emplace_back(std::move(path));
} }
} }
@ -3084,7 +3085,7 @@ namespace Clipper2Lib {
if (outrec->is_open) if (outrec->is_open)
{ {
PathD path; PathD path;
if (BuildPathD(outrec->pts, ReverseSolution, true, path, invScale_)) if (BuildPathD(outrec->pts, reverse_solution_, true, path, invScale_))
open_paths.push_back(path); open_paths.push_back(path);
continue; continue;
} }

463
thirdparty/clipper2/Clipper2Lib/src/clipper.offset.cpp vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 7 August 2023 * * Date : 28 November 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : Path Offset (Inflate/Shrink) * * Purpose : Path Offset (Inflate/Shrink) *
@ -20,38 +20,63 @@ const double floating_point_tolerance = 1e-12;
// Miscellaneous methods // Miscellaneous methods
//------------------------------------------------------------------------------ //------------------------------------------------------------------------------
void GetBoundsAndLowestPolyIdx(const Paths64& paths, Rect64& r, int & idx) inline bool ToggleBoolIf(bool val, bool condition)
{ {
idx = -1; return condition ? !val : val;
r = MaxInvalidRect64;
int64_t lpx = 0;
for (int i = 0; i < static_cast<int>(paths.size()); ++i)
for (const Point64& p : paths[i])
{
if (p.y >= r.bottom)
{
if (p.y > r.bottom || p.x < lpx)
{
idx = i;
lpx = p.x;
r.bottom = p.y;
}
}
else if (p.y < r.top) r.top = p.y;
if (p.x > r.right) r.right = p.x;
else if (p.x < r.left) r.left = p.x;
}
//if (idx < 0) r = Rect64(0, 0, 0, 0);
//if (r.top == INT64_MIN) r.bottom = r.top;
//if (r.left == INT64_MIN) r.left = r.right;
} }
bool IsSafeOffset(const Rect64& r, double abs_delta) void GetMultiBounds(const Paths64& paths, std::vector<Rect64>& recList)
{ {
return r.left > min_coord + abs_delta && recList.reserve(paths.size());
r.right < max_coord - abs_delta && for (const Path64& path : paths)
r.top > min_coord + abs_delta && {
r.bottom < max_coord - abs_delta; if (path.size() < 1)
{
recList.push_back(InvalidRect64);
continue;
}
int64_t x = path[0].x, y = path[0].y;
Rect64 r = Rect64(x, y, x, y);
for (const Point64& pt : path)
{
if (pt.y > r.bottom) r.bottom = pt.y;
else if (pt.y < r.top) r.top = pt.y;
if (pt.x > r.right) r.right = pt.x;
else if (pt.x < r.left) r.left = pt.x;
}
recList.push_back(r);
}
}
bool ValidateBounds(std::vector<Rect64>& recList, double delta)
{
int64_t int_delta = static_cast<int64_t>(delta);
int64_t big = MAX_COORD - int_delta;
int64_t small = MIN_COORD + int_delta;
for (const Rect64& r : recList)
{
if (!r.IsValid()) continue; // ignore invalid paths
else if (r.left < small || r.right > big ||
r.top < small || r.bottom > big) return false;
}
return true;
}
int GetLowestClosedPathIdx(std::vector<Rect64>& boundsList)
{
int i = -1, result = -1;
Point64 botPt = Point64(INT64_MAX, INT64_MIN);
for (const Rect64& r : boundsList)
{
++i;
if (!r.IsValid()) continue; // ignore invalid paths
else if (r.bottom > botPt.y || (r.bottom == botPt.y && r.left < botPt.x))
{
botPt = Point64(r.left, r.bottom);
result = static_cast<int>(i);
}
}
return result;
} }
PointD GetUnitNormal(const Point64& pt1, const Point64& pt2) PointD GetUnitNormal(const Point64& pt1, const Point64& pt2)
@ -125,6 +150,44 @@ inline void NegatePath(PathD& path)
} }
} }
//------------------------------------------------------------------------------
// ClipperOffset::Group methods
//------------------------------------------------------------------------------
ClipperOffset::Group::Group(const Paths64& _paths, JoinType _join_type, EndType _end_type):
paths_in(_paths), join_type(_join_type), end_type(_end_type)
{
bool is_joined =
(end_type == EndType::Polygon) ||
(end_type == EndType::Joined);
for (Path64& p: paths_in)
StripDuplicates(p, is_joined);
// get bounds of each path --> bounds_list
GetMultiBounds(paths_in, bounds_list);
if (end_type == EndType::Polygon)
{
is_hole_list.reserve(paths_in.size());
for (const Path64& path : paths_in)
is_hole_list.push_back(Area(path) < 0);
lowest_path_idx = GetLowestClosedPathIdx(bounds_list);
// the lowermost path must be an outer path, so if its orientation is negative,
// then flag the whole group is 'reversed' (will negate delta etc.)
// as this is much more efficient than reversing every path.
is_reversed = (lowest_path_idx >= 0) && is_hole_list[lowest_path_idx];
if (is_reversed) is_hole_list.flip();
}
else
{
lowest_path_idx = -1;
is_reversed = false;
is_hole_list.resize(paths_in.size());
}
}
//------------------------------------------------------------------------------ //------------------------------------------------------------------------------
// ClipperOffset methods // ClipperOffset methods
//------------------------------------------------------------------------------ //------------------------------------------------------------------------------
@ -147,10 +210,10 @@ void ClipperOffset::BuildNormals(const Path64& path)
norms.clear(); norms.clear();
norms.reserve(path.size()); norms.reserve(path.size());
if (path.size() == 0) return; if (path.size() == 0) return;
Path64::const_iterator path_iter, path_last_iter = --path.cend(); Path64::const_iterator path_iter, path_stop_iter = --path.cend();
for (path_iter = path.cbegin(); path_iter != path_last_iter; ++path_iter) for (path_iter = path.cbegin(); path_iter != path_stop_iter; ++path_iter)
norms.push_back(GetUnitNormal(*path_iter,*(path_iter +1))); norms.push_back(GetUnitNormal(*path_iter,*(path_iter +1)));
norms.push_back(GetUnitNormal(*path_last_iter, *(path.cbegin()))); norms.push_back(GetUnitNormal(*path_stop_iter, *(path.cbegin())));
} }
inline PointD TranslatePoint(const PointD& pt, double dx, double dy) inline PointD TranslatePoint(const PointD& pt, double dx, double dy)
@ -200,7 +263,25 @@ PointD IntersectPoint(const PointD& pt1a, const PointD& pt1b,
} }
} }
void ClipperOffset::DoSquare(Group& group, const Path64& path, size_t j, size_t k) void ClipperOffset::DoBevel(const Path64& path, size_t j, size_t k)
{
PointD pt1, pt2;
if (j == k)
{
double abs_delta = std::abs(group_delta_);
pt1 = PointD(path[j].x - abs_delta * norms[j].x, path[j].y - abs_delta * norms[j].y);
pt2 = PointD(path[j].x + abs_delta * norms[j].x, path[j].y + abs_delta * norms[j].y);
}
else
{
pt1 = PointD(path[j].x + group_delta_ * norms[k].x, path[j].y + group_delta_ * norms[k].y);
pt2 = PointD(path[j].x + group_delta_ * norms[j].x, path[j].y + group_delta_ * norms[j].y);
}
path_out.push_back(Point64(pt1));
path_out.push_back(Point64(pt2));
}
void ClipperOffset::DoSquare(const Path64& path, size_t j, size_t k)
{ {
PointD vec; PointD vec;
if (j == k) if (j == k)
@ -228,8 +309,8 @@ void ClipperOffset::DoSquare(Group& group, const Path64& path, size_t j, size_t
pt.z = ptQ.z; pt.z = ptQ.z;
#endif #endif
//get the second intersect point through reflecion //get the second intersect point through reflecion
group.path.push_back(Point64(ReflectPoint(pt, ptQ))); path_out.push_back(Point64(ReflectPoint(pt, ptQ)));
group.path.push_back(Point64(pt)); path_out.push_back(Point64(pt));
} }
else else
{ {
@ -238,28 +319,28 @@ void ClipperOffset::DoSquare(Group& group, const Path64& path, size_t j, size_t
#ifdef USINGZ #ifdef USINGZ
pt.z = ptQ.z; pt.z = ptQ.z;
#endif #endif
group.path.push_back(Point64(pt)); path_out.push_back(Point64(pt));
//get the second intersect point through reflecion //get the second intersect point through reflecion
group.path.push_back(Point64(ReflectPoint(pt, ptQ))); path_out.push_back(Point64(ReflectPoint(pt, ptQ)));
} }
} }
void ClipperOffset::DoMiter(Group& group, const Path64& path, size_t j, size_t k, double cos_a) void ClipperOffset::DoMiter(const Path64& path, size_t j, size_t k, double cos_a)
{ {
double q = group_delta_ / (cos_a + 1); double q = group_delta_ / (cos_a + 1);
#ifdef USINGZ #ifdef USINGZ
group.path.push_back(Point64( path_out.push_back(Point64(
path[j].x + (norms[k].x + norms[j].x) * q, path[j].x + (norms[k].x + norms[j].x) * q,
path[j].y + (norms[k].y + norms[j].y) * q, path[j].y + (norms[k].y + norms[j].y) * q,
path[j].z)); path[j].z));
#else #else
group.path.push_back(Point64( path_out.push_back(Point64(
path[j].x + (norms[k].x + norms[j].x) * q, path[j].x + (norms[k].x + norms[j].x) * q,
path[j].y + (norms[k].y + norms[j].y) * q)); path[j].y + (norms[k].y + norms[j].y) * q));
#endif #endif
} }
void ClipperOffset::DoRound(Group& group, const Path64& path, size_t j, size_t k, double angle) void ClipperOffset::DoRound(const Path64& path, size_t j, size_t k, double angle)
{ {
if (deltaCallback64_) { if (deltaCallback64_) {
// when deltaCallback64_ is assigned, group_delta_ won't be constant, // when deltaCallback64_ is assigned, group_delta_ won't be constant,
@ -280,29 +361,25 @@ void ClipperOffset::DoRound(Group& group, const Path64& path, size_t j, size_t k
if (j == k) offsetVec.Negate(); if (j == k) offsetVec.Negate();
#ifdef USINGZ #ifdef USINGZ
group.path.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y, pt.z)); path_out.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y, pt.z));
#else #else
group.path.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y)); path_out.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y));
#endif #endif
if (angle > -PI + 0.01) // avoid 180deg concave int steps = static_cast<int>(std::ceil(steps_per_rad_ * std::abs(angle))); // #448, #456
for (int i = 1; i < steps; ++i) // ie 1 less than steps
{ {
int steps = static_cast<int>(std::ceil(steps_per_rad_ * std::abs(angle))); // #448, #456 offsetVec = PointD(offsetVec.x * step_cos_ - step_sin_ * offsetVec.y,
for (int i = 1; i < steps; ++i) // ie 1 less than steps offsetVec.x * step_sin_ + offsetVec.y * step_cos_);
{
offsetVec = PointD(offsetVec.x * step_cos_ - step_sin_ * offsetVec.y,
offsetVec.x * step_sin_ + offsetVec.y * step_cos_);
#ifdef USINGZ #ifdef USINGZ
group.path.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y, pt.z)); path_out.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y, pt.z));
#else #else
group.path.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y)); path_out.push_back(Point64(pt.x + offsetVec.x, pt.y + offsetVec.y));
#endif #endif
}
} }
group.path.push_back(GetPerpendic(path[j], norms[j], group_delta_)); path_out.push_back(GetPerpendic(path[j], norms[j], group_delta_));
} }
void ClipperOffset::OffsetPoint(Group& group, Path64& path, size_t j, size_t k) void ClipperOffset::OffsetPoint(Group& group, const Path64& path, size_t j, size_t k)
{ {
// Let A = change in angle where edges join // Let A = change in angle where edges join
// A == 0: ie no change in angle (flat join) // A == 0: ie no change in angle (flat join)
@ -323,57 +400,51 @@ void ClipperOffset::OffsetPoint(Group& group, Path64& path, size_t j, size_t k)
} }
if (std::fabs(group_delta_) <= floating_point_tolerance) if (std::fabs(group_delta_) <= floating_point_tolerance)
{ {
group.path.push_back(path[j]); path_out.push_back(path[j]);
return; return;
} }
if (cos_a > 0.999) // almost straight - less than 2.5 degree (#424, #526) if (cos_a > -0.99 && (sin_a * group_delta_ < 0)) // test for concavity first (#593)
{
DoMiter(group, path, j, k, cos_a);
}
else if (cos_a > -0.99 && (sin_a * group_delta_ < 0))
{ {
// is concave // is concave
group.path.push_back(GetPerpendic(path[j], norms[k], group_delta_)); path_out.push_back(GetPerpendic(path[j], norms[k], group_delta_));
// this extra point is the only (simple) way to ensure that // this extra point is the only (simple) way to ensure that
// path reversals are fully cleaned with the trailing clipper // path reversals are fully cleaned with the trailing clipper
group.path.push_back(path[j]); // (#405) path_out.push_back(path[j]); // (#405)
group.path.push_back(GetPerpendic(path[j], norms[j], group_delta_)); path_out.push_back(GetPerpendic(path[j], norms[j], group_delta_));
}
else if (cos_a > 0.999 && join_type_ != JoinType::Round)
{
// almost straight - less than 2.5 degree (#424, #482, #526 & #724)
DoMiter(path, j, k, cos_a);
} }
else if (join_type_ == JoinType::Miter) else if (join_type_ == JoinType::Miter)
{ {
// miter unless the angle is so acute the miter would exceeds ML // miter unless the angle is sufficiently acute to exceed ML
if (cos_a > temp_lim_ - 1) DoMiter(group, path, j, k, cos_a); if (cos_a > temp_lim_ - 1) DoMiter(path, j, k, cos_a);
else DoSquare(group, path, j, k); else DoSquare(path, j, k);
} }
else if (cos_a > 0.99 || join_type_ == JoinType::Square) // 0.99 ~= 8.1 deg. else if (join_type_ == JoinType::Round)
DoSquare(group, path, j, k); DoRound(path, j, k, std::atan2(sin_a, cos_a));
else if ( join_type_ == JoinType::Bevel)
DoBevel(path, j, k);
else else
DoRound(group, path, j, k, std::atan2(sin_a, cos_a)); DoSquare(path, j, k);
} }
void ClipperOffset::OffsetPolygon(Group& group, Path64& path) void ClipperOffset::OffsetPolygon(Group& group, const Path64& path)
{ {
// when the path is contracting, make sure path_out.clear();
// there is sufficient space to do so. //#593
// nb: this will have a small impact on performance
double a = Area(path);
// contracting when orientation is opposite offset direction
if ((a < 0) != (group_delta_ < 0))
{
Rect64 rec = GetBounds(path);
if (std::fabs(group_delta_) * 2 > rec.Width()) return;
}
for (Path64::size_type j = 0, k = path.size() -1; j < path.size(); k = j, ++j) for (Path64::size_type j = 0, k = path.size() -1; j < path.size(); k = j, ++j)
OffsetPoint(group, path, j, k); OffsetPoint(group, path, j, k);
group.paths_out.push_back(group.path); solution.push_back(path_out);
} }
void ClipperOffset::OffsetOpenJoined(Group& group, Path64& path) void ClipperOffset::OffsetOpenJoined(Group& group, const Path64& path)
{ {
OffsetPolygon(group, path); OffsetPolygon(group, path);
std::reverse(path.begin(), path.end()); Path64 reverse_path(path);
std::reverse(reverse_path.begin(), reverse_path.end());
//rebuild normals // BuildNormals(path); //rebuild normals // BuildNormals(path);
std::reverse(norms.begin(), norms.end()); std::reverse(norms.begin(), norms.end());
@ -381,39 +452,28 @@ void ClipperOffset::OffsetOpenJoined(Group& group, Path64& path)
norms.erase(norms.begin()); norms.erase(norms.begin());
NegatePath(norms); NegatePath(norms);
group.path.clear(); OffsetPolygon(group, reverse_path);
OffsetPolygon(group, path);
} }
void ClipperOffset::OffsetOpenPath(Group& group, Path64& path) void ClipperOffset::OffsetOpenPath(Group& group, const Path64& path)
{ {
// do the line start cap // do the line start cap
if (deltaCallback64_) group_delta_ = deltaCallback64_(path, norms, 0, 0); if (deltaCallback64_) group_delta_ = deltaCallback64_(path, norms, 0, 0);
if (std::fabs(group_delta_) <= floating_point_tolerance) if (std::fabs(group_delta_) <= floating_point_tolerance)
group.path.push_back(path[0]); path_out.push_back(path[0]);
else else
{ {
switch (end_type_) switch (end_type_)
{ {
case EndType::Butt: case EndType::Butt:
#ifdef USINGZ DoBevel(path, 0, 0);
group.path.push_back(Point64(
path[0].x - norms[0].x * group_delta_,
path[0].y - norms[0].y * group_delta_,
path[0].z));
#else
group.path.push_back(Point64(
path[0].x - norms[0].x * group_delta_,
path[0].y - norms[0].y * group_delta_));
#endif
group.path.push_back(GetPerpendic(path[0], norms[0], group_delta_));
break; break;
case EndType::Round: case EndType::Round:
DoRound(group, path, 0, 0, PI); DoRound(path, 0, 0, PI);
break; break;
default: default:
DoSquare(group, path, 0, 0); DoSquare(path, 0, 0);
break; break;
} }
} }
@ -433,64 +493,42 @@ void ClipperOffset::OffsetOpenPath(Group& group, Path64& path)
group_delta_ = deltaCallback64_(path, norms, highI, highI); group_delta_ = deltaCallback64_(path, norms, highI, highI);
if (std::fabs(group_delta_) <= floating_point_tolerance) if (std::fabs(group_delta_) <= floating_point_tolerance)
group.path.push_back(path[highI]); path_out.push_back(path[highI]);
else else
{ {
switch (end_type_) switch (end_type_)
{ {
case EndType::Butt: case EndType::Butt:
#ifdef USINGZ DoBevel(path, highI, highI);
group.path.push_back(Point64(
path[highI].x - norms[highI].x * group_delta_,
path[highI].y - norms[highI].y * group_delta_,
path[highI].z));
#else
group.path.push_back(Point64(
path[highI].x - norms[highI].x * group_delta_,
path[highI].y - norms[highI].y * group_delta_));
#endif
group.path.push_back(GetPerpendic(path[highI], norms[highI], group_delta_));
break; break;
case EndType::Round: case EndType::Round:
DoRound(group, path, highI, highI, PI); DoRound(path, highI, highI, PI);
break; break;
default: default:
DoSquare(group, path, highI, highI); DoSquare(path, highI, highI);
break; break;
} }
} }
for (size_t j = highI, k = 0; j > 0; k = j, --j) for (size_t j = highI, k = 0; j > 0; k = j, --j)
OffsetPoint(group, path, j, k); OffsetPoint(group, path, j, k);
group.paths_out.push_back(group.path); solution.push_back(path_out);
} }
void ClipperOffset::DoGroupOffset(Group& group) void ClipperOffset::DoGroupOffset(Group& group)
{ {
Rect64 r;
int idx = -1;
//the lowermost polygon must be an outer polygon. So we can use that as the
//designated orientation for outer polygons (needed for tidy-up clipping)
GetBoundsAndLowestPolyIdx(group.paths_in, r, idx);
if (idx < 0) return;
if (group.end_type == EndType::Polygon) if (group.end_type == EndType::Polygon)
{ {
double area = Area(group.paths_in[idx]); // a straight path (2 points) can now also be 'polygon' offset
//if (area == 0) return; // probably unhelpful (#430) // where the ends will be treated as (180 deg.) joins
group.is_reversed = (area < 0); if (group.lowest_path_idx < 0) delta_ = std::abs(delta_);
if (group.is_reversed) group_delta_ = -delta_; group_delta_ = (group.is_reversed) ? -delta_ : delta_;
else group_delta_ = delta_;
}
else
{
group.is_reversed = false;
group_delta_ = std::abs(delta_) * 0.5;
} }
else
group_delta_ = std::abs(delta_);// *0.5;
double abs_delta = std::fabs(group_delta_); double abs_delta = std::fabs(group_delta_);
// do range checking if (!ValidateBounds(group.bounds_list, abs_delta))
if (!IsSafeOffset(r, abs_delta))
{ {
DoError(range_error_i); DoError(range_error_i);
error_code_ |= range_error_i; error_code_ |= range_error_i;
@ -500,10 +538,9 @@ void ClipperOffset::DoGroupOffset(Group& group)
join_type_ = group.join_type; join_type_ = group.join_type;
end_type_ = group.end_type; end_type_ = group.end_type;
if (!deltaCallback64_ && if (group.join_type == JoinType::Round || group.end_type == EndType::Round)
(group.join_type == JoinType::Round || group.end_type == EndType::Round))
{ {
//calculate a sensible number of steps (for 360 deg for the given offset) // calculate a sensible number of steps (for 360 deg for the given offset)
// arcTol - when arc_tolerance_ is undefined (0), the amount of // arcTol - when arc_tolerance_ is undefined (0), the amount of
// curve imprecision that's allowed is based on the size of the // curve imprecision that's allowed is based on the size of the
// offset (delta). Obviously very large offsets will almost always // offset (delta). Obviously very large offsets will almost always
@ -519,61 +556,80 @@ void ClipperOffset::DoGroupOffset(Group& group)
steps_per_rad_ = steps_per_360 / (2 * PI); steps_per_rad_ = steps_per_360 / (2 * PI);
} }
bool is_joined = std::vector<Rect64>::const_iterator path_rect_it = group.bounds_list.cbegin();
(end_type_ == EndType::Polygon) || std::vector<bool>::const_iterator is_hole_it = group.is_hole_list.cbegin();
(end_type_ == EndType::Joined); Paths64::const_iterator path_in_it = group.paths_in.cbegin();
Paths64::iterator path_iter; for ( ; path_in_it != group.paths_in.cend(); ++path_in_it, ++path_rect_it, ++is_hole_it)
for(path_iter = group.paths_in.begin(); path_iter != group.paths_in.end(); ++path_iter)
{ {
Path64 &path = *path_iter; if (!path_rect_it->IsValid()) continue;
StripDuplicates(path, is_joined); Path64::size_type pathLen = path_in_it->size();
Path64::size_type cnt = path.size(); path_out.clear();
if (cnt == 0 || ((cnt < 3) && group.end_type == EndType::Polygon))
continue;
group.path.clear(); if (pathLen == 1) // single point
if (cnt == 1) // single point - only valid with open paths
{ {
if (group_delta_ < 1) continue; if (group_delta_ < 1) continue;
const Point64& pt = (*path_in_it)[0];
//single vertex so build a circle or square ... //single vertex so build a circle or square ...
if (group.join_type == JoinType::Round) if (group.join_type == JoinType::Round)
{ {
double radius = abs_delta; double radius = abs_delta;
group.path = Ellipse(path[0], radius, radius); int steps = static_cast<int>(std::ceil(steps_per_rad_ * 2 * PI)); //#617
path_out = Ellipse(pt, radius, radius, steps);
#ifdef USINGZ #ifdef USINGZ
for (auto& p : group.path) p.z = path[0].z; for (auto& p : path_out) p.z = pt.z;
#endif #endif
} }
else else
{ {
int d = (int)std::ceil(abs_delta); int d = (int)std::ceil(abs_delta);
r = Rect64(path[0].x - d, path[0].y - d, path[0].x + d, path[0].y + d); Rect64 r = Rect64(pt.x - d, pt.y - d, pt.x + d, pt.y + d);
group.path = r.AsPath(); path_out = r.AsPath();
#ifdef USINGZ #ifdef USINGZ
for (auto& p : group.path) p.z = path[0].z; for (auto& p : path_out) p.z = pt.z;
#endif #endif
} }
group.paths_out.push_back(group.path); solution.push_back(path_out);
} continue;
else } // end of offsetting a single point
{
if ((cnt == 2) && (group.end_type == EndType::Joined))
{
if (group.join_type == JoinType::Round)
end_type_ = EndType::Round;
else
end_type_ = EndType::Square;
}
BuildNormals(path); // when shrinking outer paths, make sure they can shrink this far (#593)
if (end_type_ == EndType::Polygon) OffsetPolygon(group, path); // also when shrinking holes, make sure they too can shrink this far (#715)
else if (end_type_ == EndType::Joined) OffsetOpenJoined(group, path); if ((group_delta_ > 0) == ToggleBoolIf(*is_hole_it, group.is_reversed) &&
else OffsetOpenPath(group, path); (std::min(path_rect_it->Width(), path_rect_it->Height()) <= -group_delta_ * 2) )
} continue;
if ((pathLen == 2) && (group.end_type == EndType::Joined))
end_type_ = (group.join_type == JoinType::Round) ?
EndType::Round :
EndType::Square;
BuildNormals(*path_in_it);
if (end_type_ == EndType::Polygon) OffsetPolygon(group, *path_in_it);
else if (end_type_ == EndType::Joined) OffsetOpenJoined(group, *path_in_it);
else OffsetOpenPath(group, *path_in_it);
} }
solution.reserve(solution.size() + group.paths_out.size()); }
copy(group.paths_out.begin(), group.paths_out.end(), back_inserter(solution));
group.paths_out.clear();
size_t ClipperOffset::CalcSolutionCapacity()
{
size_t result = 0;
for (const Group& g : groups_)
result += (g.end_type == EndType::Joined) ? g.paths_in.size() * 2 : g.paths_in.size();
return result;
}
bool ClipperOffset::CheckReverseOrientation()
{
// nb: this assumes there's consistency in orientation between groups
bool is_reversed_orientation = false;
for (const Group& g : groups_)
if (g.end_type == EndType::Polygon)
{
is_reversed_orientation = g.is_reversed;
break;
}
return is_reversed_orientation;
} }
void ClipperOffset::ExecuteInternal(double delta) void ClipperOffset::ExecuteInternal(double delta)
@ -581,29 +637,29 @@ void ClipperOffset::ExecuteInternal(double delta)
error_code_ = 0; error_code_ = 0;
solution.clear(); solution.clear();
if (groups_.size() == 0) return; if (groups_.size() == 0) return;
solution.reserve(CalcSolutionCapacity());
if (std::abs(delta) < 0.5) if (std::abs(delta) < 0.5) // ie: offset is insignificant
{ {
Paths64::size_type sol_size = 0;
for (const Group& group : groups_) sol_size += group.paths_in.size();
solution.reserve(sol_size);
for (const Group& group : groups_) for (const Group& group : groups_)
{
solution.reserve(solution.size() + group.paths_in.size());
copy(group.paths_in.begin(), group.paths_in.end(), back_inserter(solution)); copy(group.paths_in.begin(), group.paths_in.end(), back_inserter(solution));
} return;
} }
else
{
temp_lim_ = (miter_limit_ <= 1) ?
2.0 :
2.0 / (miter_limit_ * miter_limit_);
delta_ = delta; temp_lim_ = (miter_limit_ <= 1) ?
std::vector<Group>::iterator git; 2.0 :
for (git = groups_.begin(); git != groups_.end(); ++git) 2.0 / (miter_limit_ * miter_limit_);
{
DoGroupOffset(*git); delta_ = delta;
if (!error_code_) continue; // all OK std::vector<Group>::iterator git;
solution.clear(); for (git = groups_.begin(); git != groups_.end(); ++git)
} {
DoGroupOffset(*git);
if (!error_code_) continue; // all OK
solution.clear();
} }
} }
@ -614,19 +670,17 @@ void ClipperOffset::Execute(double delta, Paths64& paths)
ExecuteInternal(delta); ExecuteInternal(delta);
if (!solution.size()) return; if (!solution.size()) return;
paths = solution; bool paths_reversed = CheckReverseOrientation();
//clean up self-intersections ... //clean up self-intersections ...
Clipper64 c; Clipper64 c;
c.PreserveCollinear = false; c.PreserveCollinear(false);
//the solution should retain the orientation of the input //the solution should retain the orientation of the input
c.ReverseSolution = reverse_solution_ != groups_[0].is_reversed; c.ReverseSolution(reverse_solution_ != paths_reversed);
#ifdef USINGZ #ifdef USINGZ
if (zCallback64_) { if (zCallback64_) { c.SetZCallback(zCallback64_); }
c.SetZCallback(zCallback64_);
}
#endif #endif
c.AddSubject(solution); c.AddSubject(solution);
if (groups_[0].is_reversed) if (paths_reversed)
c.Execute(ClipType::Union, FillRule::Negative, paths); c.Execute(ClipType::Union, FillRule::Negative, paths);
else else
c.Execute(ClipType::Union, FillRule::Positive, paths); c.Execute(ClipType::Union, FillRule::Positive, paths);
@ -640,18 +694,21 @@ void ClipperOffset::Execute(double delta, PolyTree64& polytree)
ExecuteInternal(delta); ExecuteInternal(delta);
if (!solution.size()) return; if (!solution.size()) return;
bool paths_reversed = CheckReverseOrientation();
//clean up self-intersections ... //clean up self-intersections ...
Clipper64 c; Clipper64 c;
c.PreserveCollinear = false; c.PreserveCollinear(false);
//the solution should retain the orientation of the input //the solution should retain the orientation of the input
c.ReverseSolution = reverse_solution_ != groups_[0].is_reversed; c.ReverseSolution (reverse_solution_ != paths_reversed);
#ifdef USINGZ #ifdef USINGZ
if (zCallback64_) { if (zCallback64_) {
c.SetZCallback(zCallback64_); c.SetZCallback(zCallback64_);
} }
#endif #endif
c.AddSubject(solution); c.AddSubject(solution);
if (groups_[0].is_reversed)
if (paths_reversed)
c.Execute(ClipType::Union, FillRule::Negative, polytree); c.Execute(ClipType::Union, FillRule::Negative, polytree);
else else
c.Execute(ClipType::Union, FillRule::Positive, polytree); c.Execute(ClipType::Union, FillRule::Positive, polytree);

130
thirdparty/clipper2/Clipper2Lib/src/clipper.rectclip.cpp vendored
View File

@ -1,6 +1,6 @@
/******************************************************************************* /*******************************************************************************
* Author : Angus Johnson * * Author : Angus Johnson *
* Date : 6 August 2023 * * Date : 8 September 2023 *
* Website : http://www.angusj.com * * Website : http://www.angusj.com *
* Copyright : Angus Johnson 2010-2023 * * Copyright : Angus Johnson 2010-2023 *
* Purpose : FAST rectangular clipping * * Purpose : FAST rectangular clipping *
@ -24,11 +24,11 @@ namespace Clipper2Lib {
for (const Point64& pt : path2) for (const Point64& pt : path2)
{ {
PointInPolygonResult pip = PointInPolygon(pt, path1); PointInPolygonResult pip = PointInPolygon(pt, path1);
switch (pip) switch (pip)
{ {
case PointInPolygonResult::IsOutside: ++io_count; break; case PointInPolygonResult::IsOutside: ++io_count; break;
case PointInPolygonResult::IsInside: --io_count; break; case PointInPolygonResult::IsInside: --io_count; break;
default: continue; default: continue;
} }
if (std::abs(io_count) > 1) break; if (std::abs(io_count) > 1) break;
} }
@ -66,6 +66,56 @@ namespace Clipper2Lib {
return true; return true;
} }
inline bool IsHorizontal(const Point64& pt1, const Point64& pt2)
{
return pt1.y == pt2.y;
}
inline bool GetSegmentIntersection(const Point64& p1,
const Point64& p2, const Point64& p3, const Point64& p4, Point64& ip)
{
double res1 = CrossProduct(p1, p3, p4);
double res2 = CrossProduct(p2, p3, p4);
if (res1 == 0)
{
ip = p1;
if (res2 == 0) return false; // segments are collinear
else if (p1 == p3 || p1 == p4) return true;
//else if (p2 == p3 || p2 == p4) { ip = p2; return true; }
else if (IsHorizontal(p3, p4)) return ((p1.x > p3.x) == (p1.x < p4.x));
else return ((p1.y > p3.y) == (p1.y < p4.y));
}
else if (res2 == 0)
{
ip = p2;
if (p2 == p3 || p2 == p4) return true;
else if (IsHorizontal(p3, p4)) return ((p2.x > p3.x) == (p2.x < p4.x));
else return ((p2.y > p3.y) == (p2.y < p4.y));
}
if ((res1 > 0) == (res2 > 0)) return false;
double res3 = CrossProduct(p3, p1, p2);
double res4 = CrossProduct(p4, p1, p2);
if (res3 == 0)
{
ip = p3;
if (p3 == p1 || p3 == p2) return true;
else if (IsHorizontal(p1, p2)) return ((p3.x > p1.x) == (p3.x < p2.x));
else return ((p3.y > p1.y) == (p3.y < p2.y));
}
else if (res4 == 0)
{
ip = p4;
if (p4 == p1 || p4 == p2) return true;
else if (IsHorizontal(p1, p2)) return ((p4.x > p1.x) == (p4.x < p2.x));
else return ((p4.y > p1.y) == (p4.y < p2.y));
}
if ((res3 > 0) == (res4 > 0)) return false;
// segments must intersect to get here
return GetIntersectPoint(p1, p2, p3, p4, ip);
}
inline bool GetIntersection(const Path64& rectPath, inline bool GetIntersection(const Path64& rectPath,
const Point64& p, const Point64& p2, Location& loc, Point64& ip) const Point64& p, const Point64& p2, Location& loc, Point64& ip)
{ {
@ -74,100 +124,84 @@ namespace Clipper2Lib {
switch (loc) switch (loc)
{ {
case Location::Left: case Location::Left:
if (SegmentsIntersect(p, p2, rectPath[0], rectPath[3], true)) if (GetSegmentIntersection(p, p2, rectPath[0], rectPath[3], ip)) return true;
GetIntersectPoint(p, p2, rectPath[0], rectPath[3], ip); else if ((p.y < rectPath[0].y) && GetSegmentIntersection(p, p2, rectPath[0], rectPath[1], ip))
else if (p.y < rectPath[0].y &&
SegmentsIntersect(p, p2, rectPath[0], rectPath[1], true))
{ {
GetIntersectPoint(p, p2, rectPath[0], rectPath[1], ip);
loc = Location::Top; loc = Location::Top;
return true;
} }
else if (SegmentsIntersect(p, p2, rectPath[2], rectPath[3], true)) else if (GetSegmentIntersection(p, p2, rectPath[2], rectPath[3], ip))
{ {
GetIntersectPoint(p, p2, rectPath[2], rectPath[3], ip);
loc = Location::Bottom; loc = Location::Bottom;
return true;
} }
else return false; else return false;
break;
case Location::Top: case Location::Top:
if (SegmentsIntersect(p, p2, rectPath[0], rectPath[1], true)) if (GetSegmentIntersection(p, p2, rectPath[0], rectPath[1], ip)) return true;
GetIntersectPoint(p, p2, rectPath[0], rectPath[1], ip); else if ((p.x < rectPath[0].x) && GetSegmentIntersection(p, p2, rectPath[0], rectPath[3], ip))
else if (p.x < rectPath[0].x &&
SegmentsIntersect(p, p2, rectPath[0], rectPath[3], true))
{ {
GetIntersectPoint(p, p2, rectPath[0], rectPath[3], ip);
loc = Location::Left; loc = Location::Left;
return true;
} }
else if (p.x > rectPath[1].x && else if (GetSegmentIntersection(p, p2, rectPath[1], rectPath[2], ip))
SegmentsIntersect(p, p2, rectPath[1], rectPath[2], true))
{ {
GetIntersectPoint(p, p2, rectPath[1], rectPath[2], ip);
loc = Location::Right; loc = Location::Right;
return true;
} }
else return false; else return false;
break;
case Location::Right: case Location::Right:
if (SegmentsIntersect(p, p2, rectPath[1], rectPath[2], true)) if (GetSegmentIntersection(p, p2, rectPath[1], rectPath[2], ip)) return true;
GetIntersectPoint(p, p2, rectPath[1], rectPath[2], ip); else if ((p.y < rectPath[1].y) && GetSegmentIntersection(p, p2, rectPath[0], rectPath[1], ip))
else if (p.y < rectPath[0].y &&
SegmentsIntersect(p, p2, rectPath[0], rectPath[1], true))
{ {
GetIntersectPoint(p, p2, rectPath[0], rectPath[1], ip);
loc = Location::Top; loc = Location::Top;
return true;
} }
else if (SegmentsIntersect(p, p2, rectPath[2], rectPath[3], true)) else if (GetSegmentIntersection(p, p2, rectPath[2], rectPath[3], ip))
{ {
GetIntersectPoint(p, p2, rectPath[2], rectPath[3], ip);
loc = Location::Bottom; loc = Location::Bottom;
return true;
} }
else return false; else return false;
break;
case Location::Bottom: case Location::Bottom:
if (SegmentsIntersect(p, p2, rectPath[2], rectPath[3], true)) if (GetSegmentIntersection(p, p2, rectPath[2], rectPath[3], ip)) return true;
GetIntersectPoint(p, p2, rectPath[2], rectPath[3], ip); else if ((p.x < rectPath[3].x) && GetSegmentIntersection(p, p2, rectPath[0], rectPath[3], ip))
else if (p.x < rectPath[3].x &&
SegmentsIntersect(p, p2, rectPath[0], rectPath[3], true))
{ {
GetIntersectPoint(p, p2, rectPath[0], rectPath[3], ip);
loc = Location::Left; loc = Location::Left;
return true;
} }
else if (p.x > rectPath[2].x && else if (GetSegmentIntersection(p, p2, rectPath[1], rectPath[2], ip))
SegmentsIntersect(p, p2, rectPath[1], rectPath[2], true))
{ {
GetIntersectPoint(p, p2, rectPath[1], rectPath[2], ip);
loc = Location::Right; loc = Location::Right;
return true;
} }
else return false; else return false;
break;
default: // loc == rInside default: // loc == rInside
if (SegmentsIntersect(p, p2, rectPath[0], rectPath[3], true)) if (GetSegmentIntersection(p, p2, rectPath[0], rectPath[3], ip))
{ {
GetIntersectPoint(p, p2, rectPath[0], rectPath[3], ip);
loc = Location::Left; loc = Location::Left;
return true;
} }
else if (SegmentsIntersect(p, p2, rectPath[0], rectPath[1], true)) else if (GetSegmentIntersection(p, p2, rectPath[0], rectPath[1], ip))
{ {
GetIntersectPoint(p, p2, rectPath[0], rectPath[1], ip);
loc = Location::Top; loc = Location::Top;
return true;
} }
else if (SegmentsIntersect(p, p2, rectPath[1], rectPath[2], true)) else if (GetSegmentIntersection(p, p2, rectPath[1], rectPath[2], ip))
{ {
GetIntersectPoint(p, p2, rectPath[1], rectPath[2], ip);
loc = Location::Right; loc = Location::Right;
return true;
} }
else if (SegmentsIntersect(p, p2, rectPath[2], rectPath[3], true)) else if (GetSegmentIntersection(p, p2, rectPath[2], rectPath[3], ip))
{ {
GetIntersectPoint(p, p2, rectPath[2], rectPath[3], ip);
loc = Location::Bottom; loc = Location::Bottom;
return true;
} }
else return false; else return false;
break;
} }
return true;
} }
inline Location GetAdjacentLocation(Location loc, bool isClockwise) inline Location GetAdjacentLocation(Location loc, bool isClockwise)