2020-11-16 14:47:43 +08:00
|
|
|
/*
|
2022-04-29 15:36:16 +08:00
|
|
|
* Poly2Tri Copyright (c) 2009-2022, Poly2Tri Contributors
|
2020-11-16 14:47:43 +08:00
|
|
|
* https://github.com/jhasse/poly2tri
|
|
|
|
*
|
|
|
|
* All rights reserved.
|
|
|
|
*
|
|
|
|
* Redistribution and use in source and binary forms, with or without modification,
|
|
|
|
* are permitted provided that the following conditions are met:
|
|
|
|
*
|
|
|
|
* * Redistributions of source code must retain the above copyright notice,
|
|
|
|
* this list of conditions and the following disclaimer.
|
|
|
|
* * Redistributions in binary form must reproduce the above copyright notice,
|
|
|
|
* this list of conditions and the following disclaimer in the documentation
|
|
|
|
* and/or other materials provided with the distribution.
|
|
|
|
* * Neither the name of Poly2Tri nor the names of its contributors may be
|
|
|
|
* used to endorse or promote products derived from this software without specific
|
|
|
|
* prior written permission.
|
|
|
|
*
|
|
|
|
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
|
|
|
|
* "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
|
|
|
|
* LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR
|
|
|
|
* A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR
|
|
|
|
* CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL,
|
|
|
|
* EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO,
|
|
|
|
* PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR
|
|
|
|
* PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF
|
|
|
|
* LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
|
|
|
|
* NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
|
|
|
|
* SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
|
|
|
|
*/
|
|
|
|
#include "sweep_context.h"
|
|
|
|
#include <algorithm>
|
|
|
|
#include "advancing_front.h"
|
|
|
|
|
|
|
|
namespace p2t {
|
|
|
|
|
2022-04-29 15:36:16 +08:00
|
|
|
SweepContext::SweepContext(std::vector<Point*> polyline) : points_(std::move(polyline)),
|
2020-11-16 14:47:43 +08:00
|
|
|
front_(nullptr),
|
|
|
|
head_(nullptr),
|
|
|
|
tail_(nullptr),
|
|
|
|
af_head_(nullptr),
|
|
|
|
af_middle_(nullptr),
|
|
|
|
af_tail_(nullptr)
|
|
|
|
{
|
|
|
|
InitEdges(points_);
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::AddHole(const std::vector<Point*>& polyline)
|
|
|
|
{
|
|
|
|
InitEdges(polyline);
|
2022-04-29 15:36:16 +08:00
|
|
|
for (auto i : polyline) {
|
|
|
|
points_.push_back(i);
|
2020-11-16 14:47:43 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::AddPoint(Point* point) {
|
|
|
|
points_.push_back(point);
|
|
|
|
}
|
|
|
|
|
|
|
|
std::vector<Triangle*> &SweepContext::GetTriangles()
|
|
|
|
{
|
|
|
|
return triangles_;
|
|
|
|
}
|
|
|
|
|
|
|
|
std::list<Triangle*> &SweepContext::GetMap()
|
|
|
|
{
|
|
|
|
return map_;
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::InitTriangulation()
|
|
|
|
{
|
|
|
|
double xmax(points_[0]->x), xmin(points_[0]->x);
|
|
|
|
double ymax(points_[0]->y), ymin(points_[0]->y);
|
|
|
|
|
|
|
|
// Calculate bounds.
|
2022-04-29 15:36:16 +08:00
|
|
|
for (auto& point : points_) {
|
|
|
|
Point& p = *point;
|
2020-11-16 14:47:43 +08:00
|
|
|
if (p.x > xmax)
|
|
|
|
xmax = p.x;
|
|
|
|
if (p.x < xmin)
|
|
|
|
xmin = p.x;
|
|
|
|
if (p.y > ymax)
|
|
|
|
ymax = p.y;
|
|
|
|
if (p.y < ymin)
|
|
|
|
ymin = p.y;
|
|
|
|
}
|
|
|
|
|
|
|
|
double dx = kAlpha * (xmax - xmin);
|
|
|
|
double dy = kAlpha * (ymax - ymin);
|
2021-06-01 11:47:19 +08:00
|
|
|
head_ = new Point(xmin - dx, ymin - dy);
|
|
|
|
tail_ = new Point(xmax + dx, ymin - dy);
|
2020-11-16 14:47:43 +08:00
|
|
|
|
|
|
|
// Sort points along y-axis
|
|
|
|
std::sort(points_.begin(), points_.end(), cmp);
|
|
|
|
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::InitEdges(const std::vector<Point*>& polyline)
|
|
|
|
{
|
|
|
|
size_t num_points = polyline.size();
|
|
|
|
for (size_t i = 0; i < num_points; i++) {
|
|
|
|
size_t j = i < num_points - 1 ? i + 1 : 0;
|
|
|
|
edge_list.push_back(new Edge(*polyline[i], *polyline[j]));
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
Point* SweepContext::GetPoint(size_t index)
|
|
|
|
{
|
|
|
|
return points_[index];
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::AddToMap(Triangle* triangle)
|
|
|
|
{
|
|
|
|
map_.push_back(triangle);
|
|
|
|
}
|
|
|
|
|
2021-06-01 11:47:19 +08:00
|
|
|
Node* SweepContext::LocateNode(const Point& point)
|
2020-11-16 14:47:43 +08:00
|
|
|
{
|
|
|
|
// TODO implement search tree
|
2021-06-01 11:47:19 +08:00
|
|
|
return front_->LocateNode(point.x);
|
2020-11-16 14:47:43 +08:00
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::CreateAdvancingFront()
|
|
|
|
{
|
|
|
|
|
|
|
|
// Initial triangle
|
2021-06-01 11:47:19 +08:00
|
|
|
Triangle* triangle = new Triangle(*points_[0], *head_, *tail_);
|
2020-11-16 14:47:43 +08:00
|
|
|
|
|
|
|
map_.push_back(triangle);
|
|
|
|
|
|
|
|
af_head_ = new Node(*triangle->GetPoint(1), *triangle);
|
|
|
|
af_middle_ = new Node(*triangle->GetPoint(0), *triangle);
|
|
|
|
af_tail_ = new Node(*triangle->GetPoint(2));
|
|
|
|
front_ = new AdvancingFront(*af_head_, *af_tail_);
|
|
|
|
|
|
|
|
// TODO: More intuitive if head is middles next and not previous?
|
|
|
|
// so swap head and tail
|
|
|
|
af_head_->next = af_middle_;
|
|
|
|
af_middle_->next = af_tail_;
|
|
|
|
af_middle_->prev = af_head_;
|
|
|
|
af_tail_->prev = af_middle_;
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::RemoveNode(Node* node)
|
|
|
|
{
|
|
|
|
delete node;
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::MapTriangleToNodes(Triangle& t)
|
|
|
|
{
|
|
|
|
for (int i = 0; i < 3; i++) {
|
|
|
|
if (!t.GetNeighbor(i)) {
|
|
|
|
Node* n = front_->LocatePoint(t.PointCW(*t.GetPoint(i)));
|
|
|
|
if (n)
|
|
|
|
n->triangle = &t;
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::RemoveFromMap(Triangle* triangle)
|
|
|
|
{
|
|
|
|
map_.remove(triangle);
|
|
|
|
}
|
|
|
|
|
|
|
|
void SweepContext::MeshClean(Triangle& triangle)
|
|
|
|
{
|
|
|
|
std::vector<Triangle *> triangles;
|
|
|
|
triangles.push_back(&triangle);
|
|
|
|
|
|
|
|
while(!triangles.empty()){
|
|
|
|
Triangle *t = triangles.back();
|
|
|
|
triangles.pop_back();
|
|
|
|
|
|
|
|
if (t != nullptr && !t->IsInterior()) {
|
|
|
|
t->IsInterior(true);
|
|
|
|
triangles_.push_back(t);
|
|
|
|
for (int i = 0; i < 3; i++) {
|
|
|
|
if (!t->constrained_edge[i])
|
|
|
|
triangles.push_back(t->GetNeighbor(i));
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
}
|
|
|
|
|
|
|
|
SweepContext::~SweepContext()
|
|
|
|
{
|
|
|
|
|
|
|
|
// Clean up memory
|
|
|
|
|
|
|
|
delete head_;
|
|
|
|
delete tail_;
|
|
|
|
delete front_;
|
|
|
|
delete af_head_;
|
|
|
|
delete af_middle_;
|
|
|
|
delete af_tail_;
|
|
|
|
|
2022-04-29 15:36:16 +08:00
|
|
|
for (auto ptr : map_) {
|
|
|
|
delete ptr;
|
2020-11-16 14:47:43 +08:00
|
|
|
}
|
|
|
|
|
2022-04-29 15:36:16 +08:00
|
|
|
for (auto& i : edge_list) {
|
|
|
|
delete i;
|
2020-11-16 14:47:43 +08:00
|
|
|
}
|
|
|
|
}
|
|
|
|
|
2021-06-01 11:47:19 +08:00
|
|
|
} // namespace p2t
|