@seasketch/geoprocessing/dist/src/toolbox/antimeridian.test.js

Geoprocessing and reporting framework for SeaSketch 2.0

import { describe, test, expect } from "vitest";
import { splitFeatureAntimeridian, splitSketchAntimeridian, splitBBoxAntimeridian, } from "./antimeridian.js";
import { bbox, feature, polygon } from "@turf/turf";
import { toFeaturePolygonArray } from "../helpers/index.js";
import { genSampleSketch, genSampleSketchCollectionFromSketches, isSketch, isSketchCollection, } from "../helpers/sketch.js";
import deepEqual from "fast-deep-equal";
describe("splitFeature", () => {
    test("splitFeature Polygon antimeridian", async () => {
        // This polygon crosses, it is to the left of -180, and to the right up to 175
        const poly = feature({
            type: "Polygon",
            coordinates: [
                [
                    [-175, -16],
                    [-188, -16],
                    [-188, -15],
                    [-175, -15],
                    [-175, -16],
                ],
            ],
        });
        const result = toFeaturePolygonArray(splitFeatureAntimeridian(poly));
        // Split into a multipolygon with two polygons
        expect(result.length).toBe(1);
        expect(result[0].geometry.type).toBe("MultiPolygon");
        expect(result[0].bbox).toEqual([-175, -16, 172, -15]); // new bounding box should have been cleaned to be within -180 to 180
        expect(result[0].geometry.coordinates).toEqual([
            [
                [
                    [-180, -15],
                    [-175, -15],
                    [-175, -16],
                    [-180, -16],
                    [-180, -15],
                ],
            ],
            [
                [
                    [180, -16],
                    [172, -16],
                    [172, -15],
                    [180, -15],
                    [180, -16],
                ],
            ],
        ]);
    });
});
describe("splitSketch", () => {
    test("splitSketch Sketch not crossing should be untouched", async () => {
        const poly = genSampleSketch({
            type: "Polygon",
            coordinates: [
                [
                    [2, 2],
                    [2, 8],
                    [8, 8],
                    [8, 2],
                    [2, 2],
                ],
            ],
        });
        const splitPoly = splitSketchAntimeridian(poly);
        expect(deepEqual(splitPoly, poly)).toBe(true); // no change
    });
    test("splitSketch SketchCollection not crossing should be untouched", async () => {
        const sketch = genSampleSketch({
            type: "Polygon",
            coordinates: [
                [
                    [2, 2],
                    [2, 8],
                    [8, 8],
                    [8, 2],
                    [2, 2],
                ],
            ],
        });
        const sc = genSampleSketchCollectionFromSketches([sketch]);
        const splitSc = splitSketchAntimeridian(sc);
        expect(deepEqual(splitSc, sc)).toBe(true); // no change
    });
    test("splitSketch Sketch antimeridian", async () => {
        const sketch = genSampleSketch({
            type: "Polygon",
            coordinates: [
                [
                    [175, -16],
                    [-178, -16],
                    [-178, -15],
                    [175, -15],
                    [175, -16],
                ],
            ],
        });
        const result = splitSketchAntimeridian(sketch);
        // Split into a multipolygon with two polygons
        if (isSketch(result)) {
            expect(result.geometry.type).toBe("MultiPolygon");
            expect(result.geometry.coordinates).toEqual([
                [
                    [
                        [180, -15],
                        [175, -15],
                        [175, -16],
                        [180, -16],
                        [180, -15],
                    ],
                ],
                [
                    [
                        [-180, -16],
                        [-178, -16],
                        [-178, -15],
                        [-180, -15],
                        [-180, -16],
                    ],
                ],
            ]);
        }
        else {
            throw new Error("Should not get here");
        }
    });
    test("splitSketch SketchCollection antimeridian", async () => {
        const sketch = genSampleSketch({
            type: "Polygon",
            coordinates: [
                [
                    [175, -16],
                    [-178, -16],
                    [-178, -15],
                    [175, -15],
                    [175, -16],
                ],
            ],
        });
        const sc = genSampleSketchCollectionFromSketches([sketch]);
        const splitSc = splitSketchAntimeridian(sc);
        // Split into a multipolygon with two polygons
        if (isSketchCollection(splitSc)) {
            expect(splitSc.features[0].geometry.type).toBe("MultiPolygon");
            expect(deepEqual(sc.properties, splitSc.properties)).toBe(true);
            expect(deepEqual(sc.bbox, splitSc.bbox)).toBe(true);
            expect(splitSc.features[0].geometry.coordinates).toEqual([
                [
                    [
                        [180, -15],
                        [175, -15],
                        [175, -16],
                        [180, -16],
                        [180, -15],
                    ],
                ],
                [
                    [
                        [-180, -16],
                        [-178, -16],
                        [-178, -15],
                        [-180, -15],
                        [-180, -16],
                    ],
                ],
            ]);
        }
        else {
            throw new Error("Should not reach here");
        }
    });
});
describe("splitBBoxAntimeridian", () => {
    test("splitBBoxAntimeridian bbox crosses antimeridian at -180", async () => {
        // This poly extends 10 degrees on each side of the antimeridian at -180
        const poly = polygon([
            [
                [-170, 10],
                [-190, 10],
                [-190, 0],
                [-170, 0],
                [-170, 10],
            ],
        ]);
        const polyBox = bbox(poly);
        const boxes = splitBBoxAntimeridian(polyBox);
        expect(boxes.length).toBe(2);
        expect(boxes[0]).toEqual([170, 0, 180, 10]); // portion extending below -180 shifted to be less than 180
        expect(boxes[1]).toEqual([-180, 0, -170, 10]); // portion within -180 left as-is
    });
    test("splitBBoxAntimeridian bbox crosses antimeridian at +180", async () => {
        // This poly extends 10 degrees on each side of the antimeridian at -180
        const poly = polygon([
            [
                [170, 10],
                [190, 10],
                [190, 0],
                [170, 0],
                [170, 10],
            ],
        ]);
        const polyBox = bbox(poly);
        const boxes = splitBBoxAntimeridian(polyBox);
        expect(boxes.length).toBe(2);
        expect(boxes[0]).toEqual([170, 0, 180, 10]); // portion within 180 left as-is
        expect(boxes[1]).toEqual([-180, 0, -170, 10]); // portion extending beyond 180 shifted to be greater than -180
    });
});
//# sourceMappingURL=antimeridian.test.js.map