19 |
20 |
Namespace: GIS
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 |
30 |
31 | GIS
32 |
33 |
34 |
35 |
36 |
37 |
38 |
39 |
40 |
41 |
42 |
43 |
44 |
45 |
46 |
47 |
48 |
49 |
50 |
51 |
52 |
53 |
54 |
55 |
56 |
57 |
58 |
59 |
60 |
61 |
62 |
63 |
64 |
65 |
66 |
67 |
68 |
69 |
70 |
71 | - Source:
72 |
75 |
76 |
77 |
78 |
79 |
80 | - See:
81 | -
82 |
83 | - the
union function
84 |
85 |
86 |
87 |
88 |
89 |
90 |
91 |
92 |
93 |
94 |
118 |
119 |
120 |
121 |
122 |
19 |
20 |
Source: deprecated/table.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | /**
30 | * @function <s>properties/remove</s>
31 | * @deprecated
32 | * @summary From now on, use directly `geojson.features.map(d => d.properties)`
33 | */
34 | export function table(geojson) {
35 | return JSON.parse(JSON.stringify(geojson.features.map((d) => d.properties)));
36 | }
37 |
38 |
39 |
40 |
41 |
42 |
43 |
44 |
19 |
20 |
Source: gis/densify.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | import initGeosJs from "geos-wasm";
30 | import { geojsonToGeosGeom } from "../helpers/geojsonToGeosGeom";
31 | import { geosGeomToGeojson } from "../helpers/geosGeomToGeojson";
32 | import { featurecollection } from "../featurecollection.js";
33 |
34 | /**
35 | * Densify a geoJSON with GEOS-WASM
36 | *
37 | * Example: {@link https://observablehq.com/@neocartocnrs/densify?collection=@neocartocnrs/geotoolbox Observable notebook}
38 | *
39 | * @param {object|array} x - The targeted FeatureCollection / Features / Geometries
40 | * @param {object} options - Optional parameters
41 | * @param {number} options.dist - The minimal distance between nodes
42 | *
43 | */
44 | export async function densify(x, options = { dist: 1 }) {
45 | // TODO: This will create a new GEOS instance with every call
46 | // to geosunion. Ideally, we should create a single instance
47 | // when the library is loaded and then just pass it around
48 | const geos = await initGeosJs();
49 | x = featurecollection(x);
50 |
51 | const geosGeom = geojsonToGeosGeom(x, geos);
52 | const newGeom = geos.GEOSDensify(geosGeom, options.dist);
53 | const densiygeom = geosGeomToGeojson(newGeom, geos).geometries;
54 |
55 | x.features.forEach((d, i) => (d.geometry = densiygeom[i]));
56 | return x;
57 | }
58 |
59 |
60 |
61 |
62 |
63 |
64 |
65 |
19 |
20 |
Source: gis/dissolve.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | import { geoArea } from "d3-geo";
30 | import { featurecollection } from "../featurecollection.js";
31 |
32 | /**
33 | * Dissolve multipart geometries to single part geometries.
34 | *
35 | * Example: {@link https://observablehq.com/@neocartocnrs/dissolve?collection=@neocartocnrs/geotoolbox Observable notebook}
36 | *
37 | * @param geojson - The target GeoJSON FeatureCollection / array of Features / array of Geometries
38 | *
39 | */
40 | export function dissolve(geojson) {
41 | geojson = featurecollection(geojson);
42 | let result = [];
43 | geojson.features.forEach((d) => {
44 | result.push(sp(d));
45 | });
46 |
47 | const keys = Object.keys(geojson).filter((e) => e != "features");
48 | const obj = {};
49 | keys.forEach((d) => {
50 | obj[d] = geojson[d];
51 | });
52 | obj.features = result.flat();
53 |
54 | return obj;
55 | }
56 |
57 | function sp(feature) {
58 | let result = [];
59 |
60 | if (feature.geometry.type.includes("Multi")) {
61 | feature.geometry.coordinates.forEach((d) => {
62 | result.push({
63 | type: "Feature",
64 | properties: feature.properties,
65 | geometry: {
66 | type: feature.geometry.type.replace("Multi", ""),
67 | coordinates: d,
68 | },
69 | });
70 | });
71 | } else {
72 | result.push({ ...feature });
73 | }
74 |
75 | const totalArea = geoArea(feature);
76 | result.forEach((d) => (d.__share = geoArea(d) / totalArea));
77 |
78 | return JSON.parse(JSON.stringify(result));
79 | }
80 |
81 |
82 |
83 |
84 |
85 |
86 |
87 |
19 |
20 |
Source: gis/geolines.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | /**
30 | * Returns a GeoJSON FeatureCollection of natural geographic lines such as
31 | * the equator, tropics & polar circles.
32 | *
33 | * Example: {@link https://observablehq.com/@neocartocnrs/geolines?collection=@neocartocnrs/geotoolbox Observable notebook}
34 | *
35 | *
36 | */
37 | export function geolines() {
38 | let features = [];
39 | let arr = [
40 | ["Equator", 0],
41 | ["Tropic of Cancer", 23.43656],
42 | ["Tropic of Capricorn", -23.43636],
43 | ["Arctic Circle", 66.56345],
44 | ["Antarctic Circle", -66.56364],
45 | ];
46 |
47 | arr.forEach((d) => {
48 | features.push({
49 | type: "Feature",
50 | properties: { name: d[0], latitude: d[1] },
51 | geometry: line(d[1]),
52 | });
53 | });
54 |
55 | return { type: "FeatureCollection", features: features };
56 | }
57 |
58 | function line(lat) {
59 | let arr = [];
60 | let i = -180;
61 | while (i <= 180) {
62 | arr.push([i, lat]);
63 | i += 2.5;
64 | }
65 | return { type: "MultiLineString", coordinates: [arr] };
66 | }
67 |
68 |
69 |
70 |
71 |
72 |
73 |
74 |
19 |
20 |
Source: gis/simplify.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | import { topology } from "topojson-server";
30 | import { feature, quantize as quantiz } from "topojson-client";
31 | import { presimplify, quantile, simplify as simple } from "topojson-simplify";
32 |
33 | const topojson = Object.assign(
34 | {},
35 | { topology, presimplify, quantiz, quantile, simple, feature }
36 | );
37 |
38 | /**
39 | * @function simplify
40 | * @export
41 | * @description Simplify geometries. The `simplify()` function allows to simplify a geometry using <code>topojson-simplify</code> library. The parameter k difine the The quantile of the simplification. By default, the generalization level is calculated automatically to ensure smooth map display.
42 | * <br/><br/>
43 | * @param {object} data - a GeoJSON FeatureCollection
44 | * @param {object} options - Optional parameters
45 | * @param {number} [options.k = undefined] - quantile of the simplification (from 0 to 1). If not defened, the generalization level is calculated automatically to ensure smooth map display.
46 | * @param {number} [options.quantize = undefined] - A smaller quantizeAmount produces a smaller file. Typical values are between 1e4 and 1e6, although it depends on the resolution you want in the output file.
47 | * @param {number} [options.arcs = 15000] - Instead of the k parameter, you can determine the level of generalization by targeting a specific number of arcs. The result will be an approximation.
48 | * @param {boolean} [options.deepcopy = true] - Use true to ensure that the input object is not modified and to create a new object.
49 | * @example
50 | * geotoolbox.simplify(*a geojson*, {k: 0.1})
51 | */
52 |
53 | export function simplify(
54 | data,
55 | { k = undefined, quantize = undefined, arcs = 15000, deepcopy = true } = {}
56 | ) {
57 | // deep copy ?
58 | let geojson;
59 | if (deepcopy) {
60 | geojson = JSON.parse(JSON.stringify(data));
61 | } else {
62 | geojson = data;
63 | }
64 |
65 | let topo = topojson.topology({ foo: geojson });
66 | let simpl = topojson.presimplify(topo);
67 |
68 | if (k == undefined) {
69 | k = arcs / simpl.arcs.flat().length;
70 | k = k > 1 ? 1 : k;
71 | }
72 |
73 | simpl = topojson.simple(simpl, topojson.quantile(simpl, k));
74 | if (quantize) {
75 | simpl = topojson.quantiz(simpl, quantize);
76 | }
77 |
78 | geojson.features = topojson.feature(
79 | simpl,
80 | Object.keys(simpl.objects)[0]
81 | ).features;
82 | return geojson;
83 | }
84 |
85 |
86 |
87 |
88 |
89 |
90 |
91 |
19 |
20 |
Source: gis/tissot.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | import { geoCircle } from "d3-geo";
30 |
31 | /**
32 | * Generate Tissot's indicatrix.
33 | *
34 | * Example: {@link https://observablehq.com/@neocartocnrs/tissot?collection=@neocartocnrs/geotoolbox Observable notebook}
35 | *
36 | * @param {number} step - The distance between each circle
37 | *
38 | */
39 | export function tissot(step) {
40 | const circle = geoCircle()
41 | .center((d) => d)
42 | .radius(step / 4)
43 | .precision(10);
44 | const features = [];
45 | for (let y = -80; y <= 80; y += step) {
46 | for (let x = -180; x < 180; x += step) {
47 | features.push({
48 | type: "Feature",
49 | properties: {},
50 | geometry: {
51 | type: "MultiPolygon",
52 | coordinates: [circle([x, y]).coordinates],
53 | },
54 | });
55 | }
56 | }
57 |
58 | return { type: "FeatureCollection", features: features };
59 | }
60 |
61 |
62 |
63 |
64 |
65 |
66 |
67 |
19 |
20 |
Source: properties/add.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | // Add
30 | import { str2fun } from "../helpers/str2fun.js";
31 |
32 | /**
33 | * @function properties/add
34 | * @deprecated
35 | * @summary From now on, use {@link addproperty}
36 | */
37 | export function add({ x, field, expression }) {
38 | let data = [...x.features.map((d) => ({ ...d.properties }))];
39 |
40 | // Get keys
41 | let keys = [];
42 | x.features
43 | .map((d) => d.properties)
44 | .forEach((d) => {
45 | keys.push(Object.keys(d));
46 | });
47 | keys = Array.from(new Set(keys.flat()));
48 |
49 | keys.forEach((d) => {
50 | expression = expression.replace(d, `d.${d}`);
51 | });
52 |
53 | expression = "d=> " + expression;
54 |
55 | let newfield = data.map(str2fun(expression));
56 | // let newfield = data.map((d) => d.pop / d.gdp);
57 |
58 | data.forEach((d, i) => {
59 | d = Object.assign(d, { [field]: newfield[i] });
60 | });
61 |
62 | let output = JSON.parse(JSON.stringify(x));
63 | output.features.map((d, i) => (d.properties = { ...data[i] }));
64 | return output;
65 | }
66 |
67 |
68 |
69 |
70 |
71 |
72 |
73 |
19 |
20 |
Source: properties/head.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | /**
30 | * @function properties/head
31 | * @deprecated
32 | * @summary From now on, use {@link head}
33 | */
34 | export function head({ x, field, nb = 10 }) {
35 | let features = [...x.features];
36 | features = features
37 | .filter((d) => d.properties[field] != "")
38 | .filter((d) => d.properties[field] != null)
39 | .filter((d) => d.properties[field] != undefined)
40 | .filter((d) => d.properties[field] != +Infinity)
41 | .filter((d) => d.properties[field] != -Infinity)
42 | .filter((d) => d.properties[field] != NaN);
43 |
44 | let head = features.sort(
45 | (a, b) => +b.properties[field] - +a.properties[field]
46 | );
47 | features = features.slice(0, nb);
48 | let output = JSON.parse(JSON.stringify(x));
49 | output.features = features;
50 | return output;
51 | }
52 |
53 |
54 |
55 |
56 |
57 |
58 |
59 |
19 |
20 |
Source: properties/keep.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | import { remove } from "./remove.js";
30 |
31 | /**
32 | * @function properties/keep
33 | * @deprecated
34 | * @summary From now on, use {@link columns}
35 | */
36 | export function keep({ x, fields }) {
37 | // Get all keys
38 | let keys = [];
39 | x.features
40 | .map((d) => d.properties)
41 | .forEach((d) => {
42 | keys.push(Object.keys(d));
43 | });
44 | keys = Array.from(new Set(keys.flat()));
45 |
46 | // Fields to be removed
47 | let diff = keys.filter((k) => !fields.includes(k));
48 | return remove({ x, field: diff });
49 | }
50 |
51 |
52 |
53 |
54 |
55 |
56 |
57 |
19 |
20 |
Source: properties/remove.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | /**
30 | * @function properties/remove
31 | * @deprecated
32 | * @summary From now on, use {@link columns}
33 | */
34 | export function remove({ x, field }) {
35 | let data = [...x.features.map((d) => ({ ...d.properties }))];
36 | data.forEach((d) => {
37 | if (Array.isArray(field)) {
38 | field.forEach((e) => delete d[e]);
39 | } else {
40 | delete d[field];
41 | }
42 | });
43 | let output = JSON.parse(JSON.stringify(x));
44 | output.features.map((d, i) => (d.properties = { ...data[i] }));
45 | return output;
46 | }
47 |
48 |
49 |
50 |
51 |
52 |
53 |
54 |
19 |
20 |
Source: properties/select.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | import { str2fun } from "../helpers/str2fun.js";
30 |
31 | /**
32 | * @function properties/select
33 | * @deprecated
34 | * @summary From now on, use {@link filter}
35 | */
36 | export function select({ x, expression }) {
37 | let features = [...x.features];
38 |
39 | // Get keys
40 | let keys = [];
41 | x.features
42 | .map((d) => d.properties)
43 | .forEach((d) => {
44 | keys.push(Object.keys(d));
45 | });
46 | keys = Array.from(new Set(keys.flat()));
47 |
48 | keys.forEach((d) => {
49 | expression = expression.replace(d, `d.properties.${d}`);
50 | });
51 |
52 | expression = "d => " + expression;
53 |
54 | let output = JSON.parse(JSON.stringify(x));
55 | output.features = features.filter(str2fun(expression));
56 | return output;
57 | }
58 |
59 |
60 |
61 |
62 |
63 |
64 |
65 |
19 |
20 |
Source: properties/subset.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | /**
30 | * @function properties/subset
31 | * @deprecated
32 | * @summary From now on, use {@link filter}
33 | */
34 | export function subset({ x, field, selection, inverse = false }) {
35 | let features = [...x.features];
36 | selection = !Array.isArray(selection) ? [selection] : selection;
37 |
38 | if (inverse) {
39 | selection = Array.from(
40 | new Set(features.map((d) => d.properties[field]))
41 | ).filter((d) => !selection.includes(d));
42 | }
43 | let result = [];
44 |
45 | selection.forEach((e) => {
46 | result.push(features.filter((d) => d.properties[field] == e));
47 | });
48 |
49 | let output = JSON.parse(JSON.stringify(x));
50 | output.features = result.flat();
51 | return output;
52 | }
53 |
54 |
55 |
56 |
57 |
58 |
59 |
60 |
19 |
20 |
Source: properties/table.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | /**
30 | * @function properties/remove
31 | * @deprecated
32 | * @summary From now on, use directly `geojson.features.map(d => d.properties)`
33 | */
34 | export function table(geojson) {
35 | return JSON.parse(JSON.stringify(geojson.features.map((d) => d.properties)));
36 | }
37 |
38 |
39 |
40 |
41 |
42 |
43 |
44 |
19 |
20 |
Source: properties/tail.js
21 |
22 |
23 |
24 |
25 |
26 |
27 |
28 |
29 | /**
30 | * @function properties/tail
31 | * @deprecated
32 | * @summary From now on, use {@link tail}
33 | */
34 | export function tail({ x, field, nb = 10 }) {
35 | let features = [...x.features];
36 | features = features
37 | .filter((d) => d.properties[field] != "")
38 | .filter((d) => d.properties[field] != null)
39 | .filter((d) => d.properties[field] != undefined)
40 | .filter((d) => d.properties[field] != +Infinity)
41 | .filter((d) => d.properties[field] != -Infinity)
42 | .filter((d) => d.properties[field] != NaN);
43 |
44 | let head = features.sort(
45 | (a, b) => +a.properties[field] - +b.properties[field]
46 | );
47 | features = features.slice(0, nb);
48 | let output = JSON.parse(JSON.stringify(x));
49 | output.features = features;
50 | return output;
51 | }
52 |
53 |
54 |
55 |
56 |
57 |
58 |
59 |