diegoripley/tippecanoe
1
0
Fork 0
mirror of https://github.com/felt/tippecanoe.git synced 2026-10-02 16:35:40 +02:00
Code Issues Packages Projects Releases Wiki Activity

Trying to cover too-small triangles with bigger triangles

Browse Source
This commit is contained in:
Erica Fischer
2023-10-26 08:57:28 -07:00
parent a3a0e36b9d
commit 546d2fd63d
2 changed files with 22 additions and 27 deletions
Download Patch File Download Diff File Expand all files Collapse all files
earcut.cpp
+14 -26
View File
@@ -74,47 +74,35 @@ drawvec fix_by_triangulation(drawvec const &dv, int z, int detail) {
std::vector<N> indices = mapbox::earcut<N>(polygon);
drawvec out2;
bool again = true;
while (again) {
again = false;
for (size_t i = 0; i + 2 < indices.size(); i += 3) {
long long cx = 0, cy = 0;
for (size_t j = 0; j < 3; j++) {
cx += out[indices[i + j]].x;
cy += out[indices[i + j]].y;
}
cx /= 3;
cy /= 3;
for (size_t j = 0; j < 3; j++) {
size_t v1 = i + j;
size_t v2 = i + ((j + 1) % 3);
size_t v3 = i + ((j + 2) % 3);
if (std::llround(out[indices[v1]].x / scale) == std::llround(out[indices[v2]].x / scale) &&
std::llround(out[indices[v1]].y / scale) == std::llround(out[indices[v2]].y / scale)) {
double ang = atan2(out[indices[v1]].y - cy, out[indices[v1]].x - cx);
double dx = scale * cos(ang) * sqrt(2) * 20;
double dy = scale * sin(ang) * sqrt(2) * 20;
if (can_adjust(out, indices, v1, &dx, &dy)) {
out[indices[v1]].x += dx;
out[indices[v1]].y += dy;
again = true;
}
if (std::llabs(std::llround(out[indices[v1]].x / scale) - std::llround(out[indices[v2]].x / scale)) < 2 &&
std::llabs(std::llround(out[indices[v1]].y / scale) - std::llround(out[indices[v2]].y / scale)) < 2) {
double ang = atan2(out[indices[v2]].y - out[indices[v1]].y, out[indices[v2]].x - out[indices[v1]].x);
ang = atan2(out[indices[v2]].y - cy, out[indices[v2]].x - cx);
dx = scale * cos(ang) * sqrt(2) * 20;
dy = scale * sin(ang) * sqrt(2) * 20;
if (can_adjust(out, indices, v2, &dx, &dy)) {
out[indices[v2]].x += dx;
out[indices[v2]].y += dy;
again = true;
drawvec tri;
tri.emplace_back(VT_MOVETO, (long long) out[indices[v3]].x, (long long) out[indices[v3]].y);
tri.emplace_back(VT_LINETO, out[indices[v1]].x - scale * cos(ang) * sqrt(2) * 2, out[indices[v1]].y - scale * sin(ang) * sqrt(2) * 2);
tri.emplace_back(VT_LINETO, out[indices[v2]].x + scale * cos(ang) * sqrt(2) * 2, out[indices[v2]].y + scale * sin(ang) * sqrt(2) * 2);
tri.emplace_back(VT_LINETO, (long long) out[indices[v3]].x, (long long) out[indices[v3]].y);
printf("%f\n", get_area(tri, 0, tri.size()));
for (auto const &d : tri) {
out2.push_back(d);
}
}
}
}
}
drawvec out2;
for (size_t i = 0; i < out.size(); i++) {
if (out[i].op == VT_MOVETO) {
size_t j;
tile.cpp
+8 -1
View File
@@ -613,6 +613,12 @@ void *partial_feature_worker(void *v) {
drawvec geom = (*partials)[i].geoms[0];
if (t == VT_POLYGON) {
geom = clean_or_clip_poly(geom, 0, 0, false, false);
geom = fix_by_triangulation(geom, z, out_detail);
}
#if 0
if (t == VT_POLYGON) {
drawvec cleaned = clean_or_clip_poly(geom, 0, 0, false, false);
drawvec fixed = fix_by_triangulation(geom, z, out_detail);
@@ -620,6 +626,7 @@ void *partial_feature_worker(void *v) {
geom = fixed;
}
}
#endif
to_tile_scale(geom, z, out_detail);
@@ -634,7 +641,7 @@ void *partial_feature_worker(void *v) {
check_polygon(geom);
}
if (geom.size() < 3) {
if (geom.size() < 4) {
if (area > 0) {
// area is in world coordinates, calculated before scaling down
geom = revive_polygon(before, area, z, out_detail);