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

Coordinate rounding that ripples everywhere

Browse Source
This commit is contained in:
Erica Fischer
2023-07-12 14:17:58 -07:00
parent 4b662171d0
commit a5441ac21b
72 changed files with 4544 additions and 4506 deletions
Download Patch File Download Diff File Expand all files Collapse all files
projection.cpp
+6 -6
View File
@@ -57,8 +57,8 @@ void lonlat2tile(double lon, double lat, int zoom, long long *x, long long *y) {
double lat_rad = lat * M_PI / 180;
unsigned long long n = 1LL << zoom;
long long llx = n * ((lon + 180) / 360);
long long lly = n * (1 - (log(tan(lat_rad) + 1 / cos(lat_rad)) / M_PI)) / 2;
long long llx = (long long) std::round(n * ((lon + 180) / 360));
long long lly = (long long) std::round(n * (1 - (log(tan(lat_rad) + 1 / cos(lat_rad)) / M_PI)) / 2);
*x = llx;
*y = lly;
@@ -84,12 +84,12 @@ void epsg3857totile(double ix, double iy, int zoom, long long *x, long long *y)
ix = 40000000.0;
}
*x = ix * (1LL << 31) / 6378137.0 / M_PI + (1LL << 31);
*y = ((1LL << 32) - 1) - (iy * (1LL << 31) / 6378137.0 / M_PI + (1LL << 31));
*x = (long long) std::round(ix * (1LL << 31) / 6378137.0 / M_PI + (1LL << 31));
*y = (long long) std::round(((1LL << 32) - 1) - (iy * (1LL << 31) / 6378137.0 / M_PI + (1LL << 31)));
if (zoom != 0) {
*x = std::round((double) *x / (1LL << (32 - zoom)));
*y = std::round((double) *y / (1LL << (32 - zoom)));
*x = (long long) std::round((double) *x / (1LL << (32 - zoom)));
*y = (long long) std::round((double) *y / (1LL << (32 - zoom)));
}
}