@trackunit/geo-json-utils 1.14.42 → 1.14.43

package/index.cjs.js

@@ -1380,7 +1380,7 @@ const geoJsonPolygonDifference = (subject, clips) => {
 // infinity; clamping keeps the round-trip finite for shapes near the poles.
 const WEB_MERCATOR_MAX_LAT = 85.05112878;
 const DEG_TO_RAD$1 = Math.PI / 180;
-const RAD_TO_DEG = 180 / Math.PI;
+const RAD_TO_DEG$1 = 180 / Math.PI;
 /**
  * Project a `[lng, lat]` position to Web Mercator metres (EPSG:3857).
  *
@@ -1399,8 +1399,8 @@ const projectLngLatToWebMercator = (position) => {
 /** Inverse of `projectLngLatToWebMercator`: Web Mercator metres back to `[lng, lat]`. */
 const unprojectWebMercatorToLngLat = (position) => {
     const [x, y] = position;
-    const lng = (x / EARTH_RADIUS) * RAD_TO_DEG;
-    const lat = (2 * Math.atan(Math.exp(y / EARTH_RADIUS)) - Math.PI / 2) * RAD_TO_DEG;
+    const lng = (x / EARTH_RADIUS) * RAD_TO_DEG$1;
+    const lat = (2 * Math.atan(Math.exp(y / EARTH_RADIUS)) - Math.PI / 2) * RAD_TO_DEG$1;
     return [lng, lat];
 };
 const mapPolygonalCoordinates = (geometry, project) => {
@@ -1800,6 +1800,57 @@ const tuGeoJsonRectangularBoxPolygonSchema = zod.z
 });
 
 const DEG_TO_RAD = Math.PI / 180;
+const RAD_TO_DEG = 180 / Math.PI;
+/**
+ * Web Mercator is undefined at the poles; tile services clamp latitude to this
+ * value (the latitude where the projection becomes square). Projecting beyond it
+ * produces unbounded y, so we clamp before projecting.
+ */
+const MAX_WEB_MERCATOR_LAT = 85.05112878;
+/**
+ * Inner projection that accepts a pre-computed `worldSize` (= tileSize * 2^zoom).
+ *
+ * Exposed so hot loops that project dozens of points per frame (e.g. edge-label
+ * layout) can hoist `Math.pow(2, zoom)` out of the per-vertex work. Most callers
+ * should prefer the zoom-based {@link lngLatToWebMercatorPx} wrapper.
+ */
+const lngLatToMercatorPxWS = (lng, lat, worldSize) => {
+    const clampedLat = Math.max(-MAX_WEB_MERCATOR_LAT, Math.min(MAX_WEB_MERCATOR_LAT, lat));
+    const x = ((lng + 180) / 360) * worldSize;
+    const latRad = clampedLat * DEG_TO_RAD;
+    const y = (0.5 - Math.log(Math.tan(Math.PI / 4 + latRad / 2)) / (2 * Math.PI)) * worldSize;
+    return [x, y];
+};
+/**
+ * Inverse of {@link lngLatToMercatorPxWS}: map a pixel coordinate back to
+ * [lng, lat] given a pre-computed `worldSize`.
+ */
+const mercatorPxToLngLatWS = (px, py, worldSize) => {
+    const lng = (px / worldSize) * 360 - 180;
+    const lat = Math.atan(Math.sinh(Math.PI * (1 - (2 * py) / worldSize))) * RAD_TO_DEG;
+    return [lng, lat];
+};
+/**
+ * Project a [lng, lat] coordinate to absolute Web Mercator pixel space at the
+ * given zoom. This is the same projection formula map adapters use internally,
+ * so pixel-space measurements derived here agree with the rendered map.
+ */
+const lngLatToWebMercatorPx = (lng, lat, zoom, tileSize = 256) => lngLatToMercatorPxWS(lng, lat, tileSize * Math.pow(2, zoom));
+/**
+ * Inverse of {@link lngLatToWebMercatorPx}: map an absolute Web Mercator pixel
+ * coordinate back to [lng, lat] at the given zoom.
+ */
+const webMercatorPxToLngLat = (px, py, zoom, tileSize = 256) => mercatorPxToLngLatWS(px, py, tileSize * Math.pow(2, zoom));
+/**
+ * Convert a pixel distance to latitude degrees at a given zoom level,
+ * accounting for Web Mercator latitude distortion (latitude pixels scale by
+ * sec(lat)).
+ */
+const pixelsToLatDegrees = (pixels, zoom, latDeg, tileSize = 256) => {
+    const pxPerDeg = (tileSize * Math.pow(2, zoom)) / 360;
+    const secLat = 1 / Math.cos(latDeg * DEG_TO_RAD);
+    return pixels / (pxPerDeg * secLat);
+};
 /**
  * Compute the pixel length of a line segment at a given zoom level,
  * accounting for Web Mercator latitude distortion.
@@ -1871,8 +1922,12 @@ exports.isFullyContainedInGeoJsonPolygon = isFullyContainedInGeoJsonPolygon;
 exports.isGeoJsonPointInPolygon = isGeoJsonPointInPolygon;
 exports.isGeoJsonPositionInLinearRing = isGeoJsonPositionInLinearRing;
 exports.isPositionInsideRing = isPositionInsideRing;
+exports.lngLatToMercatorPxWS = lngLatToMercatorPxWS;
+exports.lngLatToWebMercatorPx = lngLatToWebMercatorPx;
+exports.mercatorPxToLngLatWS = mercatorPxToLngLatWS;
 exports.normalizeLongitudes = normalizeLongitudes;
 exports.padGeoJsonBbox = padGeoJsonBbox;
+exports.pixelsToLatDegrees = pixelsToLatDegrees;
 exports.projectLngLatToWebMercator = projectLngLatToWebMercator;
 exports.projectPolygonalToWebMercator = projectPolygonalToWebMercator;
 exports.scaleGeoJsonBbox = scaleGeoJsonBbox;
@@ -1889,3 +1944,4 @@ exports.validateBbox = validateBbox;
 exports.validateBboxWithFallback = validateBboxWithFallback;
 exports.validateFeatureCollection = validateFeatureCollection;
 exports.validatePosition = validatePosition;
+exports.webMercatorPxToLngLat = webMercatorPxToLngLat;

package/index.esm.js

@@ -1378,7 +1378,7 @@ const geoJsonPolygonDifference = (subject, clips) => {
 // infinity; clamping keeps the round-trip finite for shapes near the poles.
 const WEB_MERCATOR_MAX_LAT = 85.05112878;
 const DEG_TO_RAD$1 = Math.PI / 180;
-const RAD_TO_DEG = 180 / Math.PI;
+const RAD_TO_DEG$1 = 180 / Math.PI;
 /**
  * Project a `[lng, lat]` position to Web Mercator metres (EPSG:3857).
  *
@@ -1397,8 +1397,8 @@ const projectLngLatToWebMercator = (position) => {
 /** Inverse of `projectLngLatToWebMercator`: Web Mercator metres back to `[lng, lat]`. */
 const unprojectWebMercatorToLngLat = (position) => {
     const [x, y] = position;
-    const lng = (x / EARTH_RADIUS) * RAD_TO_DEG;
-    const lat = (2 * Math.atan(Math.exp(y / EARTH_RADIUS)) - Math.PI / 2) * RAD_TO_DEG;
+    const lng = (x / EARTH_RADIUS) * RAD_TO_DEG$1;
+    const lat = (2 * Math.atan(Math.exp(y / EARTH_RADIUS)) - Math.PI / 2) * RAD_TO_DEG$1;
     return [lng, lat];
 };
 const mapPolygonalCoordinates = (geometry, project) => {
@@ -1798,6 +1798,57 @@ const tuGeoJsonRectangularBoxPolygonSchema = z
 });
 
 const DEG_TO_RAD = Math.PI / 180;
+const RAD_TO_DEG = 180 / Math.PI;
+/**
+ * Web Mercator is undefined at the poles; tile services clamp latitude to this
+ * value (the latitude where the projection becomes square). Projecting beyond it
+ * produces unbounded y, so we clamp before projecting.
+ */
+const MAX_WEB_MERCATOR_LAT = 85.05112878;
+/**
+ * Inner projection that accepts a pre-computed `worldSize` (= tileSize * 2^zoom).
+ *
+ * Exposed so hot loops that project dozens of points per frame (e.g. edge-label
+ * layout) can hoist `Math.pow(2, zoom)` out of the per-vertex work. Most callers
+ * should prefer the zoom-based {@link lngLatToWebMercatorPx} wrapper.
+ */
+const lngLatToMercatorPxWS = (lng, lat, worldSize) => {
+    const clampedLat = Math.max(-MAX_WEB_MERCATOR_LAT, Math.min(MAX_WEB_MERCATOR_LAT, lat));
+    const x = ((lng + 180) / 360) * worldSize;
+    const latRad = clampedLat * DEG_TO_RAD;
+    const y = (0.5 - Math.log(Math.tan(Math.PI / 4 + latRad / 2)) / (2 * Math.PI)) * worldSize;
+    return [x, y];
+};
+/**
+ * Inverse of {@link lngLatToMercatorPxWS}: map a pixel coordinate back to
+ * [lng, lat] given a pre-computed `worldSize`.
+ */
+const mercatorPxToLngLatWS = (px, py, worldSize) => {
+    const lng = (px / worldSize) * 360 - 180;
+    const lat = Math.atan(Math.sinh(Math.PI * (1 - (2 * py) / worldSize))) * RAD_TO_DEG;
+    return [lng, lat];
+};
+/**
+ * Project a [lng, lat] coordinate to absolute Web Mercator pixel space at the
+ * given zoom. This is the same projection formula map adapters use internally,
+ * so pixel-space measurements derived here agree with the rendered map.
+ */
+const lngLatToWebMercatorPx = (lng, lat, zoom, tileSize = 256) => lngLatToMercatorPxWS(lng, lat, tileSize * Math.pow(2, zoom));
+/**
+ * Inverse of {@link lngLatToWebMercatorPx}: map an absolute Web Mercator pixel
+ * coordinate back to [lng, lat] at the given zoom.
+ */
+const webMercatorPxToLngLat = (px, py, zoom, tileSize = 256) => mercatorPxToLngLatWS(px, py, tileSize * Math.pow(2, zoom));
+/**
+ * Convert a pixel distance to latitude degrees at a given zoom level,
+ * accounting for Web Mercator latitude distortion (latitude pixels scale by
+ * sec(lat)).
+ */
+const pixelsToLatDegrees = (pixels, zoom, latDeg, tileSize = 256) => {
+    const pxPerDeg = (tileSize * Math.pow(2, zoom)) / 360;
+    const secLat = 1 / Math.cos(latDeg * DEG_TO_RAD);
+    return pixels / (pxPerDeg * secLat);
+};
 /**
  * Compute the pixel length of a line segment at a given zoom level,
  * accounting for Web Mercator latitude distortion.
@@ -1824,4 +1875,4 @@ const edgePixelLength = (x0, y0, x1, y1, zoom, midLatDeg, tileSize = 256) => {
     return Math.sqrt(dxPx * dxPx + dyPx * dyPx);
 };
 
-export { EARTH_RADIUS, EMPTY_FEATURE_COLLECTION, boundingBoxCrossesMeridian, checkCrossesMeridian, computeGeometryCentroid, coordinatesToStandardFormat, denormalizeLongitude, distanceToGeoJsonPolygonBoundary, edgePixelLength, extractEdges, extractFirstPointCoordinate, extractPositionsFromGeometry, geoJsonBboxSchema, geoJsonFeatureCollectionSchema, geoJsonFeatureSchema, geoJsonGeometryCollectionSchema, geoJsonGeometrySchema, geoJsonLineStringSchema, geoJsonLinearRingSchema, geoJsonMultiLineStringSchema, geoJsonMultiPointSchema, geoJsonMultiPolygonSchema, geoJsonPointSchema, geoJsonPolygonDifference, geoJsonPolygonSchema, geoJsonPosition2dSchema, geoJsonPositionSchema, getBboxFromGeoJsonPolygon, getBoundingBoxFromGeoJsonBbox, getBoundingBoxFromGeoJsonPolygon, getExtremeGeoJsonPointFromPolygon, getGeoJsonPolygonFromBoundingBox, getGeoJsonPolygonIntersection, getMinMaxLongitudes, getMultipleCoordinatesFromGeoJsonObject, getPointCoordinateFromGeoJsonObject, getPointCoordinateFromGeoJsonPoint, getPolygonFromBbox, getPolygonFromPointAndRadius, isBboxInsideFeatureCollection, isFullyContainedInGeoJsonGeometry, isFullyContainedInGeoJsonPolygon, isGeoJsonPointInPolygon, isGeoJsonPositionInLinearRing, isPositionInsideRing, normalizeLongitudes, padGeoJsonBbox, projectLngLatToWebMercator, projectPolygonalToWebMercator, scaleGeoJsonBbox, splitPolygonAtAntimeridian, splitPolygonWithHolesAtAntimeridian, toFeatureCollection, toPosition2d, tuGeoJsonPointRadiusSchema, tuGeoJsonPolygonNoHolesSchema, tuGeoJsonRectangularBoxPolygonSchema, unprojectPolygonalFromWebMercator, unprojectWebMercatorToLngLat, validateBbox, validateBboxWithFallback, validateFeatureCollection, validatePosition };
+export { EARTH_RADIUS, EMPTY_FEATURE_COLLECTION, boundingBoxCrossesMeridian, checkCrossesMeridian, computeGeometryCentroid, coordinatesToStandardFormat, denormalizeLongitude, distanceToGeoJsonPolygonBoundary, edgePixelLength, extractEdges, extractFirstPointCoordinate, extractPositionsFromGeometry, geoJsonBboxSchema, geoJsonFeatureCollectionSchema, geoJsonFeatureSchema, geoJsonGeometryCollectionSchema, geoJsonGeometrySchema, geoJsonLineStringSchema, geoJsonLinearRingSchema, geoJsonMultiLineStringSchema, geoJsonMultiPointSchema, geoJsonMultiPolygonSchema, geoJsonPointSchema, geoJsonPolygonDifference, geoJsonPolygonSchema, geoJsonPosition2dSchema, geoJsonPositionSchema, getBboxFromGeoJsonPolygon, getBoundingBoxFromGeoJsonBbox, getBoundingBoxFromGeoJsonPolygon, getExtremeGeoJsonPointFromPolygon, getGeoJsonPolygonFromBoundingBox, getGeoJsonPolygonIntersection, getMinMaxLongitudes, getMultipleCoordinatesFromGeoJsonObject, getPointCoordinateFromGeoJsonObject, getPointCoordinateFromGeoJsonPoint, getPolygonFromBbox, getPolygonFromPointAndRadius, isBboxInsideFeatureCollection, isFullyContainedInGeoJsonGeometry, isFullyContainedInGeoJsonPolygon, isGeoJsonPointInPolygon, isGeoJsonPositionInLinearRing, isPositionInsideRing, lngLatToMercatorPxWS, lngLatToWebMercatorPx, mercatorPxToLngLatWS, normalizeLongitudes, padGeoJsonBbox, pixelsToLatDegrees, projectLngLatToWebMercator, projectPolygonalToWebMercator, scaleGeoJsonBbox, splitPolygonAtAntimeridian, splitPolygonWithHolesAtAntimeridian, toFeatureCollection, toPosition2d, tuGeoJsonPointRadiusSchema, tuGeoJsonPolygonNoHolesSchema, tuGeoJsonRectangularBoxPolygonSchema, unprojectPolygonalFromWebMercator, unprojectWebMercatorToLngLat, validateBbox, validateBboxWithFallback, validateFeatureCollection, validatePosition, webMercatorPxToLngLat };

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@trackunit/geo-json-utils",
-  "version": "1.14.42",
+  "version": "1.14.43",
   "repository": "https://github.com/Trackunit/manager",
   "license": "SEE LICENSE IN LICENSE.txt",
   "engines": {

package/src/WebMercatorUtils.d.ts

@@ -1,3 +1,33 @@
+/**
+ * Inner projection that accepts a pre-computed `worldSize` (= tileSize * 2^zoom).
+ *
+ * Exposed so hot loops that project dozens of points per frame (e.g. edge-label
+ * layout) can hoist `Math.pow(2, zoom)` out of the per-vertex work. Most callers
+ * should prefer the zoom-based {@link lngLatToWebMercatorPx} wrapper.
+ */
+export declare const lngLatToMercatorPxWS: (lng: number, lat: number, worldSize: number) => readonly [number, number];
+/**
+ * Inverse of {@link lngLatToMercatorPxWS}: map a pixel coordinate back to
+ * [lng, lat] given a pre-computed `worldSize`.
+ */
+export declare const mercatorPxToLngLatWS: (px: number, py: number, worldSize: number) => readonly [number, number];
+/**
+ * Project a [lng, lat] coordinate to absolute Web Mercator pixel space at the
+ * given zoom. This is the same projection formula map adapters use internally,
+ * so pixel-space measurements derived here agree with the rendered map.
+ */
+export declare const lngLatToWebMercatorPx: (lng: number, lat: number, zoom: number, tileSize?: number) => readonly [number, number];
+/**
+ * Inverse of {@link lngLatToWebMercatorPx}: map an absolute Web Mercator pixel
+ * coordinate back to [lng, lat] at the given zoom.
+ */
+export declare const webMercatorPxToLngLat: (px: number, py: number, zoom: number, tileSize?: number) => readonly [number, number];
+/**
+ * Convert a pixel distance to latitude degrees at a given zoom level,
+ * accounting for Web Mercator latitude distortion (latitude pixels scale by
+ * sec(lat)).
+ */
+export declare const pixelsToLatDegrees: (pixels: number, zoom: number, latDeg: number, tileSize?: number) => number;
 /**
  * Compute the pixel length of a line segment at a given zoom level,
  * accounting for Web Mercator latitude distortion.