2b:null 2c:["$","title",null,{"children":"World Atlas Maker v5 (now make it Airocean) | Vivek Rao | Observable"}] 31:T6e7,viewof topo2 = { let topo2 = topojson.topology( Object.fromEntries( Object.entries(topo.objects) .filter(([key]) => topo_layers.includes(key)) .map(([key, o]) => [ rename_alternatives && key === "alternatives" ? "countries" : key, preproject ? d3.geoProject(topojson.feature(topo, o), projection) : topojson.feature(topo, o) ]) ) ); topo2 = topojson.presimplify(topo2, topojson.sphericalTriangleArea); // x 15e-5 for spherical topo2 = topojson.simplify(topo2, 15e-5 * 10 ** -line_simplification); topo2 = topojson.quantize(topo2, 10 ** topo_quantization); topo2 = JSON.parse(JSON.stringify(topo2)); Object.entries(topo2.objects).forEach(([key, o]) => o.geometries.forEach((g) => { if ("iso_n3" in g.properties) { g.id = g.properties.iso_n3; delete g.properties.iso_n3; } if ("adm0_a3" in g.properties) { g.properties.a3 = g.properties.adm0_a3; delete g.properties.adm0_a3; } if ("fill" in g.properties) { delete g.properties.fill; } delete g.properties.pieces; if (g.type === "MultiPolygon") { g.arcs = d3 .sort(g.arcs, (p) => d3.geoArea(topojson.feature(topo2, { type: "Polygon", arcs: p })) ) .reverse(); } }) ); const dl = JSON.stringify(topo2); const size = (dl.length / 1024).toFixed(1); const compressed = (evaluate_compressed_size(dl) / 1024).toFixed(1); mutable topto2 = topo2; return Object.assign( DOM.download( () => new Blob([dl], { type: "application/json" }), "world.topo.json", `Download topojson (~${size}kB, zip: ~${compressed}kB)` ), { value: topo2 } ); }32:T707,viewof clean_properties = { const fields = "name fill featurecla scalerank labelrank continent region_un region_wb region_who name_ar name_bn name_de name_en name_es name_fa name_fr name_el name_he name_hi name_hu name_id name_it name_ja name_ko name_nl name_pl name_pt name_ru name_sv name_tr name_uk name_ur name_vi name_zh name_zht".split( /\s/ ); function largestPolyCentroid(geom) { if (geom.type === "MultiPolygon") geom = { type: "Polygon", coordinates: d3.greatest(geom.coordinates, (d) => d3.geoArea({ type: "Polygon", coordinates: d }) ) || [[]] }; if (geom.type === "Polygon") return d3.geoCentroid(geom); } const centroids = new Map( [ ...pieces.features, ...topojson.feature(topo, topo.objects.alternatives).features ].map((d) => [d.id, largestPolyCentroid(d.geometry)]) ); const A = d3 .sort( Array.from(properties, ([id, props]) => { const c = centroids.get(id); if (!c) mutable debug = id; return { id, ...(c && { lon: c[0], lat: c[1] }), ...Object.fromEntries(fields.map((k) => [k, props[k]])) }; }), ({ id }) => id ) .map((d) => Object.assign(d, { region_who: who_regions.get(d.id), fill: countries.features.find((e) => e.id === d.id)?.properties?.fill }) ); const ids = new Set(d3.map(A, (d) => d.id)); for (const [i, d] of islands.entries()) { if (ids.has(d.id)) continue; A.push({ id: d.id, name: d.name, fill: colors[i % colors.length], featurecla: "small nation", scalerank: 7, labelrank: 7, lon: d.lon, lat: d.lat }); } return Inputs.table(A); }33:T65c,{ // see skipped point/India if (true) return md`*topological error fixed*`; return Plot.plot({ projection: { type: "mercator", domain: d3.geoCircle().center([78.9167, 33.3867]).radius(2)() }, marks: [ Plot.geo(topojson.feature(topo, topo.objects.alternatives), { fill: (d) => d.id, fillOpacity: 0.1, title: (d) => d.id }), void Plot.geo(topojson.feature(topo, topo.objects.pieces), { fill: "#333", fillOpacity: 0.1 // title: (d) => d.id }), Plot.geo(topojson.merge(topo, topo.objects.alternatives.geometries), { strokeWidth: 4 }), Plot.line( [ [75.772705078125, 32.879587173066305], [76.79443359375, 33.128351191631566], [78.46435546875, 32.62087018318113], [79.178466796875, 32.491230287947594], [79.2333984375, 32.96258644191747], [78.804931640625, 33.52307880890422], [78.90380859375, 34.32529192442733], [77.838134765625, 35.523285179107816], [76.88232421875, 34.65128519895413], [75.706787109375, 34.49750272138159], [74.24560546875, 34.75966612466248], [73.751220703125, 34.334364487026306], [74.10003662109374, 33.43831750748322], [74.44335937499999, 32.75494243654723], [75.26184082031249, 32.26855544621476], [75.772705078125, 32.879587173066305] ], { stroke: "green", r: 7 } ), Plot.dot([[78.9167, 33.3867]], { fill: "red" }) ] }); // Oops, there is a topological error in Kashmir… }34:Tf0e,# ANNEX ## Full list of disputed areas *based on Natural Earth—disputed areas* ### ASIA #### Kashmir India, Pakistan and China have several disputes around Kashmir: - China controls *Aksai Chin*, claimed by India - China controls *Shaksam Valley*, claimed by India - India controls *Demchok*, claimed by China (small) - India controls *Jammu & Kashmir*, claimed by Pakistan - Pakistan controls *Gilgit-Baltistan*, claimed by India - Pakistan controls *Azad-Kashmir*, claimed by India - The *Siachen Glacier* is claimed both by India and Pakistan (small) - A good representation includes the *Line of Control* as a dashed line. On top of that: - In the East, India controls *3 valleys* and *Arunachal-Pradesh* (South Tibet), claimed by India - In the South, India controls *Junagadh and Manadavar*, claimed by Pakistan (In the same region, there is a dispute between *Bhutan* and *China* for two small sections.) 👉 ACTION: show Aksai Chin, Jammu & Kashmir, Gilgit-Baltistan and Azad-Kashmir (+ the glacier) as a greyed out/dashed territory + show a dashed LOC. #### Kurils The *Kurils Islands* are controlled by Russia, but the closest to Japan are claimed by Japan. IGNORE. #### Taiwan The Island of Taiwan is the independent Republic of China, claimed by Beijing. Taiwan controls a stream of very small islands very near the continent. 👉 ACTION: Taiwan greyed out / or as an independent country (switch). ### EUROPE #### Nagorno-Karabagh Claimed by Azerbaidjan (and generally recognized as part of Azerbaidjan), it is a self-governing entity. IGNORE. #### Georgia South Ossetia an Abkhasia are self-governing breakaway territories of Georgia, under the influence of Russia. IGNORE. #### Ukraine Two territories have broken away from Ukraine, under Russia’s influence: Donbass, and Crimea, which has been annexed unilaterally and illegally by Russia. In Natural Earth, *Donbass* is still part of Ukraine. But *Crimea* has been “given” to Russia. The United Nations recognized Crimea as Ukrainian in 2014, the same year Russia annexed it. 👉 ACTION: reattach Crimea to UKR. *Transnistria* is a small territory disputed between Ukraine and Moldova. IGNORE. #### Kosovo *Kosovo* broke away from Serbia; but still not fully included in the UN system. It has an ISO code which is partly used… current position in the file is as an independent country. We want its ISO code to be XKX rather than KOS. #### Northern Cyprus Controlled by Turkey, self-governing… shown as separate on most maps. IGNORE. #### Israel/Palestine Lots of details we won't cover here. 👉 ACTION: We added Gaza to Palestine, which was absent from NE. ### AFRICA #### Egypt/Sudan The Halayib triangle is disputed; a small triangle connected to Halayib is, on the opposite, claimed by nobody! IGNORE. #### Sudan/South Sudan *Abyei* is disputed. 👉 ACTION: show as a greyed out/dashed territory. #### South Sudan / Kenya The *Ilemi triangle* is disputed. IGNORE. #### Western Sahara Occupied and claimed by Morocco, its international borders according to the UN are still on the 27.66] parallel. Natural Earth gives a great chunk of SAH to MOR. 👉 ACTION: Map according to the United Nations guidelines. #### Somaliland not fully recognized, Somaliland is sometimes difficult to attach to a region (ex. for the WHO Somaliland is a part of Somalia). 👉 ACTION: what solution? #### Canary Islands Possibly detach from Spain. ### AMERICAS #### Guyanes - *Guyana Esequiba* - Almost half the territory of Guyana is claimed by Venezuela. IGNORE. - *Courantyne* - South of Guyana, claimed by Suriname. IGNORE. - *Lawa* - South of French Guyane, claimed by Suriname. IGNORE. #### Belize The southern half of Belize is claimed by Guatemala. IGNORE. And there are more (like Bolivia claiming access to the Ocean, etc).35:T1b95,projections = [ { name: "Equal Earth", value: () => d3.geoEqualEarth().rotate([-10, 0]) }, { name: "Airocean", value: d3.geoAirocean }, // { name: "Airy’s minimum error", value: d3.geoAiry }, // { name: "Aitoff", value: d3.geoAitoff }, // { name: "American polyconic", value: d3.geoPolyconic }, { name: "Armadillo", value: () => d3.geoArmadillo().rotate([-10, 0]), options: { clip: { type: "Sphere" } } }, { name: "August", value: () => d3.geoAugust().rotate([-10, 0]) }, { name: "azimuthal equal-area", value: d3.geoAzimuthalEqualArea }, { name: "azimuthal equidistant", value: d3.geoAzimuthalEquidistant }, // { name: "Baker dinomic", value: d3.geoBaker }, // { name: "Berghaus’ star", value: d3.geoBerghaus, options: { clip: { type: "Sphere" } } }, { name: "Bertin’s 1953", value: d3.geoBertin1953 }, // { name: "Boggs’ eumorphic", value: d3.geoBoggs }, // { name: "Boggs’ eumorphic (interrupted)", value: d3.geoInterruptedBoggs, options: { clip: { type: "Sphere" } } }, { name: "Bonne", value: d3.geoBonne }, { name: "Bottomley", value: d3.geoBottomley }, { name: "Bromley", value: d3.geoBromley }, // { name: "Butterfly (gnomonic)", value: d3.geoPolyhedralButterfly }, // { name: "Butterfly (Collignon)", value: d3.geoPolyhedralCollignon }, // { name: "Butterfly (Waterman)", value: d3.geoPolyhedralWaterman }, { name: "Cahill-Keyes", value: d3.geoCahillKeyes }, // { name: "Collignon", value: d3.geoCollignon }, // {name: "conic conformal", value: d3.geoConicConformal}, // Not suitable for world maps. { name: "conic equal-area", value: d3.geoConicEqualArea }, { name: "conic equidistant", value: d3.geoConicEquidistant }, // { name: "Craig retroazimuthal", value: d3.geoCraig }, // { name: "Craster parabolic", value: d3.geoCraster }, { name: "Cox", value: d3.geoCox }, // { name: "cubic", value: d3.geoCubic }, // { name: "cylindrical equal-area", value: d3.geoCylindricalEqualArea }, // { name: "cylindrical stereographic", value: d3.geoCylindricalStereographic }, { name: "dodecahedral", value: d3.geoDodecahedral }, { name: "Eckert I", value: d3.geoEckert1 }, { name: "Eckert II", value: d3.geoEckert2 }, { name: "Eckert III", value: d3.geoEckert3 }, { name: "Eckert IV", value: d3.geoEckert4 }, { name: "Eckert V", value: d3.geoEckert5 }, { name: "Eckert VI", value: d3.geoEckert6 }, // { name: "Eisenlohr conformal", value: d3.geoEisenlohr }, { name: "Equirectangular (plate carrée)", value: d3.geoEquirectangular }, // { name: "Fahey pseudocylindrical", value: d3.geoFahey }, // { name: "flat-polar parabolic", value: d3.geoMtFlatPolarParabolic }, // { name: "flat-polar quartic", value: d3.geoMtFlatPolarQuartic }, // { name: "flat-polar sinusoidal", value: d3.geoMtFlatPolarSinusoidal }, // { name: "Foucaut’s stereographic equivalent", value: d3.geoFoucaut }, // { name: "Foucaut’s sinusoidal", value: d3.geoFoucautSinusoidal }, // { name: "general perspective", value: d3.geoSatellite }, // {name: "Gilbert’s two-world", value: d3.geoGilbert}, // https://github.com/d3/d3-geo-projection/issues/165 // { name: "Gingery", value: d3.geoGingery, options: { clip: { type: "Sphere" } } }, // { name: "Ginzburg V", value: d3.geoGinzburg5 }, // { name: "Ginzburg VI", value: d3.geoGinzburg6 }, // { name: "Ginzburg VIII", value: d3.geoGinzburg8 }, // { name: "Ginzburg IX", value: d3.geoGinzburg9 }, { name: "Goode’s homolosine", value: d3.geoHomolosine }, { name: "Goode’s homolosine (interrupted)", value: d3.geoInterruptedHomolosine, options: { clip: { type: "Sphere" } } }, // { name: "gnomonic", value: d3.geoGnomonic }, // { name: "Gringorten square", value: d3.geoGringorten }, // { name: "Gringorten quincuncial", value: d3.geoGringortenQuincuncial }, // { name: "Guyou square", value: d3.geoGuyou }, { name: "Hammer", value: d3.geoHammer }, // { name: "Hammer retroazimuthal", value: d3.geoHammerRetroazimuthal, options: { clip: { type: "Sphere" } } }, // { name: "HEALPix", value: d3.geoHealpix, options: { clip: { type: "Sphere" } } }, { name: "Hill eucyclic", value: d3.geoHill }, // { name: "Hufnagel pseudocylindrical", value: d3.geoHufnagel }, // { name: "icosahedral", value: d3.geoIcosahedral }, { name: "Imago", value: d3.geoImago }, // { name: "Kavrayskiy VII", value: d3.geoKavrayskiy7 }, // { name: "Lagrange conformal", value: d3.geoLagrange }, { name: "Larrivée", value: d3.geoLarrivee }, // { name: "Laskowski tri-optimal", value: d3.geoLaskowski }, // {name: "Littrow retroazimuthal", value: d3.geoLittrow}, // Not suitable for world maps. // { name: "Loximuthal", value: d3.geoLoximuthal }, { name: "Mercator", value: d3.geoMercator }, { name: "Miller cylindrical", value: d3.geoMiller }, { name: "Mollweide", value: d3.geoMollweide }, // { // name: "Mollweide (Goode’s interrupted)", // value: d3.geoInterruptedMollweide, // options: { clip: { type: "Sphere" } } // }, // { // name: "Mollweide (interrupted hemispheres)", // value: d3.geoInterruptedMollweideHemispheres, // options: { clip: { type: "Sphere" } } // }, // { name: "Natural Earth", value: d3.geoNaturalEarth1 }, // { name: "Natural Earth II", value: d3.geoNaturalEarth2 }, // { name: "Nell–Hammer", value: d3.geoNellHammer }, // { name: "Nicolosi globular", value: d3.geoNicolosi }, { name: "orthographic", value: d3.geoOrthographic }, // { name: "Patterson cylindrical", value: d3.geoPatterson }, { name: "Peirce quincuncial", value: d3.geoPeirceQuincuncial }, // { name: "rectangular polyconic", value: d3.geoRectangularPolyconic }, { name: "Robinson", value: d3.geoRobinson }, // { name: "sinusoidal", value: d3.geoSinusoidal }, // { name: "sinusoidal (interrupted)", value: d3.geoInterruptedSinusoidal, options: { clip: { type: "Sphere" } } }, { name: "sinu-Mollweide", value: d3.geoSinuMollweide }, // { name: "sinu-Mollweide (interrupted)", value: d3.geoInterruptedSinuMollweide, options: { clip: { type: "Sphere" } } }, { name: "stereographic", value: d3.geoStereographic }, { name: "Lee’s tetrahedal", value: d3.geoTetrahedralLee }, { name: "Times", value: d3.geoTimes }, // { name: "Tobler hyperelliptical", value: d3.geoHyperelliptical }, // { name: "transverse Mercator", value: d3.geoTransverseMercator }, // { name: "Van der Grinten", value: d3.geoVanDerGrinten }, // { name: "Van der Grinten II", value: d3.geoVanDerGrinten2 }, // { name: "Van der Grinten III", value: d3.geoVanDerGrinten3 }, // { name: "Van der Grinten IV", value: d3.geoVanDerGrinten4 }, // { name: "Wagner IV", value: d3.geoWagner4 }, // { name: "Wagner VI", value: d3.geoWagner6 }, // { name: "Wagner VII", value: d3.geoWagner7 }, // { name: "Werner", value: d3.geoBonne ? () => d3.geoBonne().parallel(90) : null }, // { name: "Wiechel", value: d3.geoWiechel }, { name: "Winkel tripel", value: d3.geoWinkel3 } ].filter((p) => p.value)36:T9b0,controls = ({ geostitch: () => Inputs.toggle({ label: "geostitch", value: true }), guf: () => Inputs.select( new Map([ ["France (FRA)", "FRA"], ["Stand-alone (GUF)", "GUF"] // although useful, this creates a real issue with ids ]), { label: "Guyane" } ), abyi: () => Inputs.select( new Map([ ["Sudan (SDN)", "SDN"], ["South Sudan (SSD)", "SSD"], ["Abyei (stand-alone, ABYI)", "ABYI"] ]), { label: "Abyei" } ), crma: () => Inputs.select( new Map([ ["Ukraine (UKR)", "UKR"], ["Russia (RUS)", "RUS"], ["Crimea (stand-alone, CRMA)", "CRMA"] ]), { label: "Crimea" } ), esh: () => Inputs.form({ esh: Inputs.select( new Map([ ["United Nations", "ESH"], ["Morocco", "MAR"], ["Disputed (stand-alone piece)", "UN"] ]), { label: "Western Sahara" } ), cany: Inputs.select( new Map([ ["Spain", "ESP"], ["Stand-alone (CANY)", "CANY"] ]), { label: "Canary Islands - 50m" } ) }), kash: () => Inputs.form({ jamm: Inputs.select( new Map([ ["India (IND)", "IND"], ["Kashmir (KASH)", "KASH"], ["Jammu & Kashmir (stand-alone, JAMM)", "JAMM"] ]), { label: "Jammu & Kashmir" } ), gilg: Inputs.select( new Map([ ["Pakistan (PAK)", "PAK"], ["Kashmir (KASH)", "KASH"], ["Gilgit and Azad (stand-alone, GILG)", "GILG"] ]), { label: "Gilgit and Azad" } ), aksa: Inputs.select( new Map([ ["China (CHN)", "CHN"], ["Kashmir (KASH)", "KASH"], ["Aksai Chin (stand-alone, AKSA)", "AKSA"] ]), { label: "Aksai Chin" } ), ...(apply_sources === "50m" && { kas: Inputs.select( new Map([ ["Siachen Glacier (stand-alone, KAS)", "KAS"], ["Kashmir (KASH)", "KASH"] ]), { label: "Siachen Glacier" } ) }) }), twn: () => Inputs.select( new Map([ ["Taïwan (TWN)", "TWN"], ["China (CHN)", "CHN"] ]), { label: "Taïwan" } ), sol: () => Inputs.select( new Map([ ["Somaliland (SOL)", "SOL"] // ["Somalia (SOM)", "SOM"] ]), { label: "Somaliland" } ) })37:T746,PIPELINE_PROPERTIES = [ { name: "unzero_props", description: "removes the unicode-zero characters from the zip", process: (a) => { for (const d of a.features) { Object.keys(d.properties).forEach( (prop) => (d.properties[prop] = JSON.parse( JSON.stringify(d.properties[prop]).replace(/\\u0000/g, "") )) ); } return a; } }, { name: "lowercase_props", description: "lowercase all properties", process: (a) => { for (const d of a.features) { Object.keys(d.properties).forEach((PROP) => { const prop = PROP.toLowerCase(); if (prop !== PROP) { d.properties[prop] = d.properties[PROP]; delete d.properties[PROP]; } }); } return a; } }, { name: "fix_missing_iso_n3", description: "Fixes missing ISO_N3 codes for Norway, North Cyprus, Kosovo and Somaliland", comment: "still present in 2021—only Norway is 'fixed' in ISO_N3_EH", process: (a) => { const c = { Norway: "578", "N. Cyprus": "-1", Kosovo: "-2", Somaliland: "-3", France: "250" }; // for (const d of a.features) { if (d.properties.iso_n3 === "-99") d.properties.iso_n3 = c[d.properties.name]; } return a; } }, { name: "fix_missing_iso_a3", description: "Fixes missing ISO_A3 codes for South Sudan, Western Sahara, Somaliland", process: (a) => { const c = { "South Sudan": "SSD", // not "SDS" "Western Sahara": "ESH", // not "SAH" Kosovo: "XKX" // not "KOS" }; for (const d of a.features) { if (d.properties.admin in c) d.properties.adm0_a3 = c[d.properties.admin]; } return a; } } ]38:T4b37,PIPELINE = [ { name: "snap", description: "Snap all coordinates", process: (a) => { const projection = d3.geoEquirectangular().scale(1e5); const snap = new Map(); for (const arc of a.arcs) { for (const p of arc) { if (Math.abs(p[0]) > 179.5) continue; const q = projection(p); let snapped = false; for (const x of [ Math.round(q[0]), Math.round(q[0] - 0.5), Math.round(q[0] + 0.5) ]) { for (const y of [ Math.round(q[1]), Math.round(q[1] - 0.5), Math.round(q[1] + 0.5) ]) { const q0 = `${x},${y}`; if (!snapped && snap.has(q0)) { p[0] = snap.get(q0)[0]; p[1] = snap.get(q0)[1]; snapped = true; } } } if (!snapped) snap.set(`${Math.round(q[0])},${Math.round(q[1])}`, p); } } // remove duplicate points on an arc a.arcs = a.arcs.map((arc) => { let prev = arc[0]; const a0 = [prev]; for (let i = 1; i < arc.length; i++) { if (String(prev) === String(arc[i])) continue; let dd; if ((dd = d3.geoDistance(arc[i], [78.9167, 33.3867])) < 0.00001) { console.warn(arc[i], dd); continue; // skip a point in India, to fix a topological error } a0.push((prev = arc[i])); } return a0; }); return a; } }, { name: "Add disputed", description: "Add disputed as pieces", disabled: false, process: (a) => { // we're going to add the disputed features to the topojson, // and snap their coordinates to the countries' nearest points, // so that the collection of arcs is complete const countries = topojson.feature(a, a.objects.countries); // in 50m we read the disputed shapefile, and snap its coordinates // to the countries if (a.objects.disputed) { const index = d3.index( a.objects.disputed.geometries, (d) => (d.id = d.properties.BRK_A3) ); index.get("B05").id = "JAMM"; index.get("B07").id = "AKSA"; index.get("B08").id = "GILG"; index.get("B45").id = "KAS"; // Siachen Glacier index.get("B13").id = "ABYI"; const disputed = topojson.feature(a, a.objects.disputed); for (const f of disputed.features.filter((d) => [ "GILG", "JAMM", "AKSA", "ABYI", "B03", "B05", "B06", "B07", "B08", "B13", "B45", "B17", "B16" // shebaa farms ].includes(d.id) )) { f.geometry.coordinates = f.geometry.coordinates.map( coordsMatcher(countries, { precision: 0.001, _snap: undefined }) ); } const abyi = disputed.features.find((d) => d.id === "ABYI"); const ssd = countries.features.find((d) => d.id === "SSD"); // remove a point from ssd that is missing in ABYI (and is useless) a.arcs = a.arcs.map((d) => d.filter((p) => d3.geoDistance(p, [28.9323, 9.5495]) > 0.00001) ); a = topojson.topology( Object.assign( Object.fromEntries( Object.entries(a.objects).map(([key, o]) => [ key, topojson.feature(a, o) ]) ), { disputed } ) ); } else { const disputed = { type: "FeatureCollection", features: [abyei, ...kashmir.features].map((d, i) => { if (i === 0) /* Abyei */ d.geometry.coordinates = d.geometry.coordinates.map( coordsMatcher(countries, { precision: 0.0003, snap: 4 }) ); if (i === 1) d.geometry.coordinates = d.geometry.coordinates.map( coordsMatcher(countries, { precision: 0.001, snap: 13 }) ); if (i === 2) d.geometry.coordinates = d.geometry.coordinates.map( coordsMatcher(countries, { precision: 0.001, snap: 13 }) ); if (i === 3) d.geometry.coordinates = d.geometry.coordinates.map( coordsMatcher(countries, { precision: 0.001, snap: 5 }) ); return d; }) }; a = topojson.topology( Object.assign( Object.fromEntries( Object.entries(a.objects).map(([key, o]) => [ key, topojson.feature(a, o) ]) ), { disputed } ) ); } // we can now start the collection of pieces, by detaching existing polygons // (Crimea from RUS; French Guyane from FRA) and adding Abyei. const p = JSON.parse(JSON.stringify(a)); a.objects.pieces = { type: "GeometryCollection", geometries: detachCrma(p, detachGuyane(p)) }; // A specific technique here: to compute the difference of PAK with GILG, we // topo-merge PAK with the *inverse* of GILG; same for the two other shapes of KASH, // as well as Abyei from Sudan mergeInverse(a, "PAK", "GILG"); mergeInverse(a, "IND", "JAMM"); mergeInverse(a, "CHN", "AKSA"); mergeInverse(a, "SDN", "ABYI"); return a; function mergeInverse(a, main, piece) { const M = a.objects.pieces.geometries.find((d) => d.id === main); const P = a.objects.disputed.geometries.find((d) => d.id === piece); if (!M || !P) return; const merged = topojson.mergeArcs(a, [ M, { type: "Polygon", arcs: P.arcs.map((ring) => ring.map((i) => -1 - i)) } ]); if (merged.arcs.length < 2) { M.type = "Polygon"; M.arcs = merged.arcs[0]; } else { M.type = merged.type; M.arcs = merged.arcs; } a.objects.pieces.geometries.push(P); } } }, { name: "egy/isr/jor", description: "fix tripoint Egypt/Israel/Jordan", comment: "Still present in 2021 (110m)", process: (a) => JSON.parse( JSON.stringify(a).replace( /\[34\.9226\d+,29\.5013[23]\d+\]/gm, "[34.9226,29.50133]" ) ) }, { name: "gaza/psx", description: "Add Gaza to PSX; fix tripoint with Egypt", comment: "Gaza still missing in 2021", process: (a) => { a.arcs = JSON.parse( JSON.stringify(a.arcs) .replace(/(\[33\.77\d+,30\.96\d+\])/g, "$1,[34.20,31.28]") .replace(/\[34\.55\d+,31\.5\d+\]/g, "[34.65,31.40]") .replace(/\[34\.26\d+,31\.2[12]\d+\]/g, "[34.34,31.03]") ); const gaza = [ [ [34.2, 31.28], [34.48810713068136, 31.60553884533732], [34.65, 31.4], [34.34, 31.03], [34.2, 31.28] ] ]; a = featureTransform(a, (geo) => { const d = geo.features.find((d) => d.id === "PSX"); if (d && d.geometry.type === "Polygon") { d.geometry.type = "MultiPolygon"; d.geometry.coordinates = [d.geometry.coordinates, gaza]; } return geo; }); return a; } }, { name: "ESH/MAR", description: "UN border between Western Sahara and Morocco: 27°40'N", // disabled: true, process: (a) => { a = featureTransform(a, (geo, key) => { if (!["countries", "countries_lakes", "pieces"].includes(key)) return geo; // identify the tripoint on Algeria—we'll reuse it. const alg = geo.features.find((d) => d.id === "DZA"); if (!alg) return geo; const tripoint = d3.least(alg.geometry.coordinates.flat(), (p) => d3.geoDistance(p, [-8.6656, 27.66667]) ); function fix_tripoint(a) { return a.map((p) => d3.geoDistance(p, tripoint) > 1e-3 ? p : tripoint ); } const mar = geo.features.find((d) => d.id === "MAR"); const esh = geo.features.find((d) => d.id === "SAH" || d.id === "ESH"); const topo = topojson.topology({ mar, esh }); const joined = topojson.merge(topo, [ topo.objects.mar, topo.objects.esh ]); const cut = d3 .geoCircle() .center([0, 90]) .precision(0.5) .radius(90 - 27.66667)().coordinates; const top = pclip.intersect(joined.coordinates[0], cut, pclip.geo); const bottom = pclip.difference(joined.coordinates[0], cut, pclip.geo); const disputedWest = pclip.difference( mar.geometry.coordinates[0], cut, pclip.geo ); const disputedEast = pclip.difference( esh.geometry.coordinates[0], cut, pclip.geo ); mar.geometry.coordinates = top[0] .map(fix_tripoint) .map(remove_duplicate_points); esh.geometry.coordinates = bottom[0] .map(fix_tripoint) .map((d) => d.slice().reverse()) .map(remove_duplicate_points); if (key === "pieces") { esh.geometry.coordinates = disputedEast[0] .map(fix_tripoint) .map((d) => d.slice().reverse()) .map(remove_duplicate_points); const west = JSON.parse(JSON.stringify(esh)); west.id = "WEST"; west.properties.name = "WEST"; west.geometry.coordinates = disputedWest[0] .map(fix_tripoint) .map(remove_duplicate_points); geo.features.push(west); mutable berm = { berm: topojson.feature(topo, { type: "LineString", id: "berm", arcs: [0], properties: { name: "berm" } }) }; } return geo; }); return a; } }, { name: "sudan-caf", description: "fix SDN+SDS+CAF tripoint", comment: "Issue still valid in 2021. Not for 50m.", disabled: apply_sources.endsWith("50m"), process: (a) => featureTransform(a, (geo) => { const B = [23.886979580860668, 8.619729712933065], C = [24.567, 8.229]; // SDN: remove C ; CAF: add B before C { const d = geo.features.find((d) => d.id === "SDN"); if (d) { const remove = C; const geom = d.geometry.coordinates[0].slice(0, -1).filter((p) => { return !( Math.abs(p[0] - remove[0]) < 0.001 && Math.abs(p[1] - remove[1]) < 0.001 ); }); geom.push(geom[0]); d.geometry.coordinates[0] = geom; } } { const d = geo.features.find((d) => d.id === "CAF"); if (d) { let geom = d.geometry.coordinates[0].slice(0, -1); let f = -1; const find = C; geom.forEach((p, i) => { if ( Math.abs(p[0] - find[0]) < 0.001 && Math.abs(p[1] - find[1]) < 0.001 ) f = i; }); geom = Array.from( new Set( [...geom.slice(0, f), B, ...geom.slice(f, Infinity)].map((d) => d.join("/") ) ), (d) => d.split("/").map(Number) ); geom.push(geom[0]); d.geometry.coordinates[0] = geom; } } return geo; }), test: (a) => { const B = [23.886979580860668, 8.619729712933065]; return ( a.features.filter( (d) => d.geometry.coordinates.flat(Infinity).includes(B[0]) && d.geometry.coordinates.flat(Infinity).includes(B[1]) ).length === 3 ); } }, /* { name: "sudan-ethiopia", description: "fix SDN+SDS+ETH tripoint", comment: "Not useful anymore in 2021.", disabled: apply_sources.endsWith("50m") || apply_sources.startsWith("tiny"), process: (a) => featureTransform(a, (geo) => { // SDN: remove C ; ETH: add B after C const B = [33.963392794971185, 9.464285229420625], C = [33.975, 8.68456]; { const d = geo.features.find((d) => d.id === "SDN"); if (d) { const remove = C; const geom = d.geometry.coordinates[0].slice(0, -1).filter((p) => { return !( Math.abs(p[0] - remove[0]) < 0.001 && Math.abs(p[1] - remove[1]) < 0.001 ); }); geom.push(geom[0]); d.geometry.coordinates[0] = geom; } } { const d = geo.features.find((d) => d.id === "ETH"); if (d) { var geom = d.geometry.coordinates[0].slice(0, -1); let f = -1; const find = C; geom.forEach((p, i) => { if ( Math.abs(p[0] - find[0]) < 0.001 && Math.abs(p[1] - find[1]) < 0.001 ) f = i; }); geom = Array.from( new Set( [...geom.slice(0, f), B, ...geom.slice(f, Infinity)].map((d) => d.join("/") ) ), (d) => d.split("/").map(Number) ); geom.push(geom[0]); d.geometry.coordinates[0] = geom; } } return geo; }), test: (a) => { const B = [33.963392794971185, 9.464285229420625]; return ( a.features.filter( (d) => d.geometry.coordinates.flat(Infinity).includes(B[0]) && d.geometry.coordinates.flat(Infinity).includes(B[1]) ).length === 3 ); } },*/ { name: "add_kurils", description: "Adds three South Kuril Islands in Russia, for aesthetic reasons.", comment: "only in 110m", process: (a) => featureTransform(a, (geo) => { const d = geo.features.find((d) => d.id === "RUS"); // only for 110m! if (d && d.geometry.coordinates.length < 50) { d.geometry.coordinates.push([ [ [146.2, 44.52], [148.44, 44.9], [149.15, 45.64], [147.52, 45.39], [146.2, 44.52] ].reverse() ]), d.geometry.coordinates.push([ [ [149.63, 45.85], [149.85, 45.64], [151.12, 46.31], [150.42, 46.4], [149.63, 45.85] ].reverse() ]), d.geometry.coordinates.push([ [ [155.74, 50.76], [154.33, 49.09], [156.66, 50.54], [155.74, 50.76] ].reverse() ]), d.geometry.coordinates.push([ [ [151.56, 46.98], [151.87, 46.7], [152.66, 47.1], [152.4, 47.33], [151.56, 46.98] ].reverse() ]); } return geo; }) }, { name: "detach_canary", description: "Detach Canary Islands from Spain mainland.", comment: "only in 50m", process: (a) => featureTransform(a, (geo, key) => { if (ctr_esh.cany !== "CANY") return geo; if (!["pieces"].includes(key)) return geo; const d = geo.features.find((d) => d.id === "ESP"); const cany = d3.group(d.geometry.coordinates, (poly) => { const [lon, lat] = d3.geoCentroid({ type: "Polygon", coordinates: poly }); // Canary islands are the 7 polygons in this bounding box return 26 < lat && lat < 30 && -22 < lon && lon < -13; }); d.geometry.coordinates = cany.get(false); geo.features.push({ type: "Feature", id: "CANY", properties: { ...d.properties, name: "Canary Islands" }, geometry: { type: "MultiPolygon", coordinates: cany.get(true) } }); return geo; }) }, { name: "crimea", description: "Reattach Crimea to countries_lakes.", process: (a) => featureTransform(a, (geo, key) => { if (ctr_crma === "RUS") return geo; if (!["countries_lakes"].includes(key)) return geo; const r = geo.features.find((d) => d.id === "RUS"); const ukr = d3.least(r.geometry.coordinates, (poly) => d3.geoDistance( d3.geoCentroid({ type: "Polygon", coordinates: poly }), [34.5, 45.4] ) ); r.geometry.coordinates = r.geometry.coordinates.filter( (d) => d !== ukr ); const u = geo.features.find((d) => d.id === "UKR"); const t = topojson.topology({ ukr: { type: "GeometryCollection", geometries: [u.geometry, { type: "Polygon", coordinates: ukr }] } }); u.geometry = topojson.merge(t, t.objects.ukr.geometries); return geo; }) }, { name: "color", description: "Assign a color to each country, respecting neighbors", process: (topo) => { for (const [key, o] of Object.entries(topo.objects)) { // key = countries, countries_lakes, disputed, pieces… //if (key === "pieces") { // console.warn(key, o); //} const gm = d3.sort( o.geometries, // by starting with these countries we can reach 5 colors! (d) => !["CHN", "SAU", "SSD", "SRB"].includes(d.id) ); const neighbors = topojson.neighbors(gm); { // because of the intermediate pieces, these neighbors must be added by hand… const pak = gm.map((d) => d.id).indexOf("PAK"); const chn = gm.map((d) => d.id).indexOf("CHN"); if (pak > -1 && chn > -1) { neighbors[pak].push(chn); neighbors[chn].push(pak); } } const index = new Uint8Array(o.geometries.length).fill(255); gm.forEach((d, i) => { const u = [ ...d3.difference( d3.range(colors.length).map((j) => (j + i) % colors.length), Array.from(neighbors[i], (j) => index[j]) ) ][0]; index[i] = u ?? 255; // d.properties.colorIndex = u; d.properties.fill = colors[index[i]] || "lime"; // d.properties.neighbors = neighbors[i].map((j) => gm[j].id).join(","); }); } return topo; } }, { name: "simplify", description: "apply topojson.simplify", disabled: true, // no clear win process: (a) => { // compute triangles’ weights a = topojson.presimplify(a); a = topojson.simplify(a, 0.00001); return a; } } ]39:Td0f,worlds = { let a = {}; const steps = []; for (const { name, process } of SOURCES) { if (name === apply_sources) { a = await process(a); steps.push(a); } } // geo to topo if (!a.objects) { a = topojson.topology({ countries: a }); } // compute and transform the properties // then index them and add the relevant IDs to the countries objects let props = topojson.feature(a, a.objects.countries); for (const { name, process } of PIPELINE_PROPERTIES) { if (apply_properties.includes(name)) { props = await process(props); steps.push(props); } } const properties = d3.index( props.features.map((d) => d.properties), (d) => d.adm0_a3 ); a.objects.countries.geometries.forEach((d) => { d.id = d.properties.adm0_a3; d.properties = { adm0_a3: d.id, iso_n3: d.properties.iso_n3, // https://en.wikipedia.org/wiki/ISO_3166-1_numeric // note: we can't always use this as an ID since Guyane has the same as Metro France name: d.properties.name }; }); if (a.objects.countries_lakes) { a.objects.countries_lakes.geometries.forEach((d) => { if (d.properties.ADM0_A3 === "SDS") d.properties.ADM0_A3 = "SSD"; if (d.properties.ADM0_A3 === "SAH") d.properties.ADM0_A3 = "ESH"; d.id = d.properties.adm0_a3 || d.properties.ADM0_A3; d.properties = { adm0_a3: d.id }; }); } // geostitch if (geostitch) { a = featureTransform(a, d3.geoStitch); } for (const { name, process } of PIPELINE) { if (apply_steps.includes(name)) { a = await process(a, properties); steps.push(a); } } steps.push(a); const piece = { featurecla: "PIECE", scalerank: 3, labelrank: 6, level: 3, type: "Piece" }; properties.set("WEST", { ...piece, name: "Western part of Western Sahara", continent: "Africa", region_un: "Africa", region_wb: "Middle East & North Africa" }); properties.set("CRMA", { ...piece, name: "Crimea", continent: "Europe", region_un: "Europe", region_wb: "Europe & Central Asia" }); properties.set("CANY", { ...piece, name: "Canary Islands", continent: "Africa", region_un: "Europe", // ? region_wb: "Europe & Central Asia" // ? }); properties.set("GUF", { ...piece, name: "Guyane", continent: "South America", region_un: "Americas", region_wb: "Europe & Central Asia" }); properties.set("ABYI", { ...piece, name: "Abyei", continent: "Africa", region_un: "Africa", region_wb: "Sub-Saharan Africa" }); properties.set("AKSA", { ...piece, name: "Aksai Chin", continent: "Asia", region_un: "Asia", region_wb: "East Asia & Pacific" // like CHN }); properties.set("GILG", { ...piece, name: "Gilgit and Azad", continent: "Asia", region_un: "Asia", region_wb: "South Asia" // like PAK }); properties.set("JAMM", { ...piece, name: "Jammu & Kashmir", continent: "Asia", region_un: "Asia", region_wb: "South Asia" // like IND }); properties.set("KASH", { ...piece, name: "Kashmir Region", continent: "Asia", region_un: "Asia", region_wb: "Asia" }); properties.get("PSX").region_who = "EMRO"; steps.properties = properties; return steps; }3a:T768,// Base topojson has all the countries and lakes and disputed areas // Each zone is marked as a (multi)polygon // We'll build countries up from multiple zones. baseTopo = { const topo = topojson.topology({ ...(countries && { countries }), ...(countries_lakes && { countries_lakes }), ...(pieces && { pieces }), ...(disputed && { disputed }), ...(berm && { berm: berm.berm }) }); // merge optional pieces to alternatives if (disputed) { topo.objects.alternatives = { type: "GeometryCollection", geometries: d3 .rollups( topo.objects.pieces.geometries, (v) => [ Object.assign( v.find((d) => d.id === d.properties.adm0_a3)?.properties || { ...v[0].properties, name: v.map((d) => d.id).join(" + ") }, { pieces: v.map((d) => d.id).join(" + "), ...(v.find((d) => d.id === "FRA") && { name: "France" }) } ), topojson.mergeArcs(topo, v) ], ({ id }) => { if (id === "WEST") return ctr_esh.esh; // "MAR", "SAH" if (id === "ABYI") return ctr_abyi; // "SDN", "ABYI", "SSD" if (id === "AKSA") return ctr_kash.aksa; // "CHN", "KASH" if (id === "GILG") return ctr_kash.gilg; // "PAK", "KASH" if (id === "JAMM") return ctr_kash.jamm; // "IND", "KASH" if (id === "KAS") return ctr_kash.kas; // Siachen glacier if (id === "GUF") return ctr_guf; // "GUF", "FRA" if (id === "CRMA") return ctr_crma; // "UKR", "RUS", "CRMA" if (id === "TWN") return ctr_twn; // "TWN", "CHN" return id; } ) .map(([id, [properties, shape]]) => ({ id, ...shape, properties })) }; } return topo; }3b:T4c0,topo = { const topo = { ...baseTopo, objects: { ocean, land, land_lakes, ...(continents && { continents }), ...(region_un && { region_un }), ...(region_who && { region_who }), //income_grp, ...(region_wb && { region_wb }), ...(baseTopo.objects.pieces && { pieces: baseTopo.objects.pieces }), countries: baseTopo.objects.countries, ...(baseTopo.objects.countries_lakes && { countries_lakes: baseTopo.objects.countries_lakes }), ...(baseTopo.objects.alternatives && { alternatives: baseTopo.objects.alternatives }), ...(baseTopo.objects.disputed && { disputed: baseTopo.objects.disputed, disputed_abyei, disputed_kashmir }), taggedMesh, ...(disputedMesh && { disputedMesh }), arcsMesh: { type: "GeometryCollection", geometries: d3.range(baseTopo.arcs.length).map((i) => ({ id: `arc${i}`, properties: { name: `arc${i}` }, type: "LineString", arcs: [i] })) } } }; for (let key in topo.objects) { removeSelfNeighbors(topo, topo.objects[key].geometries); } return topo; }3c:T411,function removeSelfNeighbors(topo, geometries) { for (const geom of geometries) { switch (geom.type) { case "Polygon": geom.arcs = geom.arcs.map(arcClean).filter((d) => d.length > 0); break; case "MultiPolygon": geom.arcs = geom.arcs .map((arcs) => arcs.map(arcClean).filter((d) => d.length > 0)) .filter((d) => d.length > 0); break; } } return geometries; function arcClean(a) { let changed = false; // [a, a] for (let i = 1; i < a.length; ++i) { if (a[i] === a[i - 1]) { a = a.slice(); delete a[i - 1]; changed = true; } } // [a, -1-a] for (let i = 1; i < a.length; ++i) { if (a[i] === -1 - a[i - 1]) { a = a.slice(); delete a[i - 1]; delete a[i]; changed = true; } } if (changed) { a = a.filter((d) => d != null); if (a.length === 1 && topo.arcs[a[0] < 0 ? -1 - a[0] : a[0]].length === 2) a = []; } return a; } }3d:T9af,async function map({ center, scale, projection_ = d3.geoMercator(), displayObjects, t = topo, sphere, visibility } = {}) { const topo = t; // projection_.rotate([-90, 0]); visibility && (await visibility()); const height = (width * 0.58) | 0; const svg = d3 .create("svg") .attr("viewBox", [0, 0, width, height]) .attr("width", width); const hh = DOM.uid("hatched"); svg.append("defs") .html(` `); if (center) { projection_.rotate(d3.geoCentroid(center).map((d) => -d)); } else { projection_.fitExtent( [ [3, 3], [width - 3, height - 3] ], { type: "Sphere" } ); } if (scale) projection_.scale(scale); const path = d3.geoPath(projection_); const g = svg.append("g").attr("fill-opacity", 0.1); if (sphere) { g.append("g") .attr("id", "sphere") .selectAll() .data([{ type: "Sphere" }]) .join("path") .attr("d", path) .attr("fill", "lightblue") .attr("stroke", "black"); } displayObjects.forEach((layer) => { g.append("g") .attr("id", layer) .attr("fill-opacity", layer.startsWith("disputed") ? 1 : 0.1) .selectAll() .data(topojson.feature(topo, topo.objects[layer]).features) .join("path") .attr("d", path) .attr( "stroke", layer === "land" || layer.endsWith("Mesh") ? d3.scaleOrdinal(["black", ...d3.schemeTableau10]) : "none" ) .attr("fill", (d) => layer === "land" || layer.endsWith("Mesh") ? "none" : layer.startsWith("disputed") ? `${hh}` //"url(#disputedHatch)" : d.properties.fill ) .append("title") .text((d) => `${d.properties.name || d.id}`); }); displayObjects.forEach((layer) => { g.append("g") .attr("id", `${layer}-mesh`) .attr("fill", "none") .attr("stroke-width", 0.5) .selectAll() .data( topojson.feature(topo, mesh(topo, topo.objects[layer], `${layer}-mesh`)) .features ) .join("path") .attr("d", path) .attr("stroke", d3.scaleOrdinal(["black", ...d3.schemeTableau10])); }); return svg.node(); }3e:Td4c,kashmir = ({ type: "FeatureCollection", id: "GILG", features: [ { type: "Feature", id: "GILG", properties: { name: "Gilgit and Azad" }, geometry: { type: "Polygon", coordinates: [ [ [73.839111328125, 32.95336814579932], [74.454345703125, 32.75494243654723], [74.102783203125, 33.44060944370356], [73.76220703125, 34.334364487026306], [74.234619140625, 34.75966612466248], [75.7177734375, 34.50655662164561], [76.893310546875, 34.66935854524543], [77.80517578125, 35.53222622770337], [76.1956787109375, 35.902399875143615], [75.8880615234375, 36.672824886786564], [75.146484375, 37.142803443716836], [74.57244873046875, 37.02229110771148], [74.05609130859375, 36.83566824724438], [72.92724609375, 36.721273880045004], [72.50976562499999, 36.12900165569652], [72.7734375, 35.746512259918504], [73.17993164062499, 35.782170703266075], [74.02587890625, 35.02999636902566], [73.377685546875, 34.488447837809304], [73.839111328125, 32.95336814579932] ].reverse() ] } }, { type: "Feature", id: "JAMM", properties: { stroke: "#555555", "stroke-width": 2, "stroke-opacity": 1, fill: "#555555", "fill-opacity": 0.5, name: "Jammu Kashmir" }, geometry: { type: "Polygon", coordinates: [ [ [75.772705078125, 32.879587173066305], [76.79443359375, 33.128351191631566], [78.46435546875, 32.62087018318113], [79.178466796875, 32.491230287947594], [79.2333984375, 32.96258644191747], [78.804931640625, 33.52307880890422], [78.90380859375, 34.32529192442733], [77.838134765625, 35.523285179107816], [76.88232421875, 34.65128519895413], [75.706787109375, 34.49750272138159], [74.24560546875, 34.75966612466248], [73.751220703125, 34.334364487026306], [74.10003662109374, 33.43831750748322], [74.44335937499999, 32.75494243654723], [75.26184082031249, 32.26855544621476], [75.772705078125, 32.879587173066305] ].reverse() ] } }, { type: "Feature", id: "AKSA", properties: { stroke: "#555555", "stroke-width": 2, "stroke-opacity": 1, fill: "#555555", "fill-opacity": 0.5, name: "Aksai-Chin" }, geometry: { type: "Polygon", coordinates: [ [ [80.37597656249999, 35.505400093441324], [79.332275390625, 36.00467348670187], [77.838134765625, 35.523285179107816], [78.914794921875, 34.34343606848294], [78.826904296875, 33.52307880890422], [79.22241210937499, 32.98102014898148], [79.178466796875, 32.491230287947594], [79.617919921875, 32.7872745269555], [79.024658203125, 33.779147331286474], [80.068359375, 34.66935854524543], [80.37597656249999, 35.505400093441324] ].reverse() ] } } ] })3f:T475,// https://observablehq.com/@fil/geoformat function geoformat(n, exclude = [], align = true) { if (typeof n === "undefined") n = 14; function q(o) { if ( align && typeof o === "object" && typeof o[0] === "object" && typeof o[0][0] === "number" ) return "°°°[" + o.map(q).join(",") + "]°°°"; if (align && typeof o === "object" && typeof o[0] === "number") return "°°°[" + o.map(q).join(",") + "]°°°"; if (typeof o === "object") return o.map(q); if (typeof o === "number") return +o.toFixed(n); return o; } return function (x) { let st = JSON.stringify( x, function (key, val) { if (key === "coordinates") { return q(val); } if (exclude.indexOf(key) > -1) return undefined; return val; }, 2 ); return align ? st .replace(/"°°°/g, "") .replace(/°°°",\n */g, ", ") .replace(/°°°"/g, "") .replace(/°°°/g, "") // no repeats .replace(/(\[[\.\-0123456789]+,[\.\-0123456789]+\])(,\1)+,/g, "$1,") : st; }; }40:T456,/* * A tagged border is an arc that belongs to two polygons with a non-null tag */ function taggedBorders(topo, geometries, tag = (a, b) => a !== b) { const x = []; for (const d of geometries) { for (const i of d.type === "MultiPolygon" ? d.arcs.flat(2) : d.type === "Polygon" ? d.arcs.flat(1) : []) x[i] = d; } const borders = new Map(); for (let i = 0; i < x.length; i++) { const t = x[-i - 1] ? tag(x[i], x[-i - 1]) : null; if (t != null && t !== false) { borders.set(t, borders.get(t) || []); borders.get(t).push(i); } } return { type: "GeometryCollection", geometries: Array.from(borders, ([name, arcs]) => { const lines = []; for (const i of arcs) { lines.push([i]); } return lines.length === 1 ? { id: name, type: "LineString", arcs: lines[0], properties: { name } } : { id: name, type: "MultiLineString", arcs: lines, properties: { name } }; }) }; }41:[{"id":"d85ca74667f1e37292804e1090dafce1a6f14782abfe620f416d63c81411506b3591b753fe1973dfe4b7e851111b0620b12e20210bc9da35e84a0af20c06b2c4","href":"https://static.observableusercontent.com/files/d85ca74667f1e37292804e1090dafce1a6f14782abfe620f416d63c81411506b3591b753fe1973dfe4b7e851111b0620b12e20210bc9da35e84a0af20c06b2c4","name":"ne_50m_admin_0_countries.zip","type":"application/zip","lastModified":1679576954411.366,"size":792663},{"id":"f7a14b40ac728f7cb8019d9a9b710dca69725ff615fe9b09fb9223a8d907d375dda53a3229b4ec818b9780c2da1101d6e041e379e79cccac66caa06d904ed592","href":"https://static.observableusercontent.com/files/f7a14b40ac728f7cb8019d9a9b710dca69725ff615fe9b09fb9223a8d907d375dda53a3229b4ec818b9780c2da1101d6e041e379e79cccac66caa06d904ed592","name":"201707-ne_110m_admin_0_countries.zip","type":"application/zip","lastModified":1679576954411.366,"size":188561},{"id":"8333dafb5d06857977ae466ee9f11e85e8de4c2917d73cc82df9a5c43d93b027e8fd3103460f0b66935cdc55c656e7be450eb8444494ff46280d1f9c7dda0f9a","href":"https://static.observableusercontent.com/files/8333dafb5d06857977ae466ee9f11e85e8de4c2917d73cc82df9a5c43d93b027e8fd3103460f0b66935cdc55c656e7be450eb8444494ff46280d1f9c7dda0f9a","name":"202112-ne_110m_admin_0_countries.zip","type":"application/zip","lastModified":1679576954411.366,"size":209851},{"id":"31ae86f2a7086b1e9fbae34dbcc6169b2e14596289869803fc4d39a7bb86266620f588199d2a6aeaed6ca0884f4b7c48ae4944e64fc9e3d2153f4ec9421b6044","href":"https://static.observableusercontent.com/files/31ae86f2a7086b1e9fbae34dbcc6169b2e14596289869803fc4d39a7bb86266620f588199d2a6aeaed6ca0884f4b7c48ae4944e64fc9e3d2153f4ec9421b6044","name":"who_regions.csv","type":"text/csv","lastModified":1679576954411.366,"size":2287},{"id":"5b3c98512c167aa6aa7ada5502dc1bc56ec3ce50616be4987609f775c8395d16daa18fe68165ae9fb54a33539dc9d4f28ef6912cb8d3b60875a793b62e80f072","href":"https://static.observableusercontent.com/files/5b3c98512c167aa6aa7ada5502dc1bc56ec3ce50616be4987609f775c8395d16daa18fe68165ae9fb54a33539dc9d4f28ef6912cb8d3b60875a793b62e80f072","name":"2019xx-ne_110m_admin_0_countries.zip","type":"application/zip","lastModified":1679576954411.366,"size":204507},{"id":"7e6a2bdac4832a7d38f1ca42e3ccadc2fcf38db42058344df719a8d15a7d3764ae96a0bd41f590fc859cba0e136867f96495fd1cdd016606d5ebaa32f5f355cd","href":"https://static.observableusercontent.com/files/7e6a2bdac4832a7d38f1ca42e3ccadc2fcf38db42058344df719a8d15a7d3764ae96a0bd41f590fc859cba0e136867f96495fd1cdd016606d5ebaa32f5f355cd","name":"201808-ne_110m_admin_0_countries.zip","type":"application/zip","lastModified":1679576954411.366,"size":196764},{"id":"b7c6a57efbdab6c2267868095a3e7511a0f50ee6e2bccba4b45d5534951a81a29d45e94c1effb63545990191618a528b17ec7aea6c17a5b275131a0493e772b8","href":"https://static.observableusercontent.com/files/b7c6a57efbdab6c2267868095a3e7511a0f50ee6e2bccba4b45d5534951a81a29d45e94c1effb63545990191618a528b17ec7aea6c17a5b275131a0493e772b8","name":"tiny_countries@2.json","type":"application/json","lastModified":1679576954411.366,"size":3013609},{"id":"8369f6cfffe606b58c48f5682d88f63d7baf261d514dd6ce2122d59e29e312f4014c25dcc3a9f457685da1d335bd55aa714f9a3caf41e567b20061ddd366b7c5","href":"https://static.observableusercontent.com/files/8369f6cfffe606b58c48f5682d88f63d7baf261d514dd6ce2122d59e29e312f4014c25dcc3a9f457685da1d335bd55aa714f9a3caf41e567b20061ddd366b7c5","name":"tiny_countries_lakes@2.json","type":"application/json","lastModified":1679576954411.366,"size":3018291},{"id":"0d4770057a3c1d8a66b32768bcee1c02fd5dbc8793d6b2a8c54322311f529a243cb9e45983a88dd2656898b65637a3e4968356bb16327b4d381870f93114e8f4","href":"https://static.observableusercontent.com/files/0d4770057a3c1d8a66b32768bcee1c02fd5dbc8793d6b2a8c54322311f529a243cb9e45983a88dd2656898b65637a3e4968356bb16327b4d381870f93114e8f4","name":"ne_50m_admin_0_breakaway_disputed_areas.zip","type":"application/zip","lastModified":1679576954411.366,"size":46787},{"id":"37d42e865735ea49327358c1324d6d065dae1e95ced289026aee6c450ca0cc3a37e0fcc3039d6beebad166ba0c90a923f6b6fbea8dba90614cc7b3153fed276c","href":"https://static.observableusercontent.com/files/37d42e865735ea49327358c1324d6d065dae1e95ced289026aee6c450ca0cc3a37e0fcc3039d6beebad166ba0c90a923f6b6fbea8dba90614cc7b3153fed276c","name":"ne_50m_admin_0_countries_lakes.zip","type":"application/zip","lastModified":1679576954411.366,"size":810508},{"id":"146605beb718548e21a3bb5461e190c73a4b0590fa1a836b8709917821713b49d694d8a2d6c5de508eabbb3bf084246f09865e913daf4815dc82adc7f98d7a1f","href":"https://static.observableusercontent.com/files/146605beb718548e21a3bb5461e190c73a4b0590fa1a836b8709917821713b49d694d8a2d6c5de508eabbb3bf084246f09865e913daf4815dc82adc7f98d7a1f","name":"island-nations@1.csv","type":"text/csv","lastModified":1679576954411.366,"size":1106}] 2f:{"id":"bef97b29cfc092c2","slug":null,"title":"World Atlas Maker v5 (now make it Airocean)","thumbnail":"8fc6c20bc5dad03e4f05187df0a86f991e96dfb3fbd831a98caa99ff7b763f94","access_level":"public","create_time":"$D2023-03-23T13:09:14.411Z","update_time":"$D2023-03-23T14:19:15.401Z","view_time":null,"first_public_time":"$D2023-03-23T13:09:14.411Z","has_unpublished_edits":null,"fork":true,"stars":0,"starred":false,"owner":{"type":"team","id":"0b85d005387c1032","login":"vivrao9","name":"Vivek Rao","avatar_url":"https://avatars-next.observableusercontent.com/bef3b5fc8c76d305b837832196006da09e14e5ff3812149e979a7f004178f2bf","tier":"starter_2024"},"version":4877,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"16d1a97215fda001","login":"vivrao9","name":"Vivek Rao","avatar_url":"https://avatars-next.observableusercontent.com/bef3b5fc8c76d305b837832196006da09e14e5ff3812149e979a7f004178f2bf","tier":null}],"members":[{"type":"individual","id":"16d1a97215fda001","login":"vivrao9","name":"Vivek Rao","avatar_url":"https://avatars-next.observableusercontent.com/bef3b5fc8c76d305b837832196006da09e14e5ff3812149e979a7f004178f2bf","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"title":"World Atlas Maker v5 (now make it Airocean)","cells":[{"id":1643,"value":"# World Atlas Maker v5 (now make it Airocean)","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3909,"value":"

👉 User Guide | Choropleth","pinned":false,"hidden":false,"mode":"html","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":1620,"value":"viewof apply_sources = Inputs.select(\n SOURCES.map((d) => d.name),\n {\n label: \"sources\",\n value: this ? this.value : SOURCES.find((d) => !d.disabled)?.name\n }\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2239,"value":"viewof geostitch = controls.geostitch()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2930,"value":"Alternative attributions:","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":2944,"value":"viewof ctr_esh = controls.esh()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2798,"value":"viewof ctr_abyi = controls.abyi()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2827,"value":"viewof ctr_crma = controls.crma()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2940,"value":"viewof ctr_guf = controls.guf()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2831,"value":"viewof ctr_twn = controls.twn()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2804,"value":"viewof ctr_kash = controls.kash()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4468,"value":"viewof ctr_sol = controls.sol() // TODO(?)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3037,"value":"---\nDownloads:","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3024,"value":"viewof geoJsonLayer = Inputs.select(\n [\n \"alternatives\",\n \"countries\",\n \"countries_lakes\",\n \"land\",\n \"land_lakes\",\n \"continents\",\n \"region_un\",\n \"region_un\",\n \"region_who\",\n \"region_wb\"\n // \"disputed\",\n // \"pieces\"\n ].filter((d) => topo.objects[d]),\n { label: \"geojson layer\" }\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3022,"value":"{\n const dl = geoformat(4)(topojson.feature(topo, topo.objects[geoJsonLayer]));\n const size = (dl.length / 1024).toFixed(0);\n const compressed = (evaluate_compressed_size(dl) / 1024).toFixed(0);\n return DOM.download(\n () => new Blob([dl], { type: \"application/json\" }),\n `${geoJsonLayer}.geo.json`,\n `Download ${geoJsonLayer}.geo.json (~${size}kB, zip: ~${compressed}kB)`\n );\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3658,"value":"---\ntopojson","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3097,"value":"viewof topo_quantization = Inputs.range([3, 6], {\n step: 0.1,\n value: 6,\n label: \"← quantize\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4196,"value":"viewof line_simplification = Inputs.range([1, 3], {\n step: 0.05,\n value: 3,\n label: \"← simplify\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3817,"value":"viewof topo_layers = Inputs.checkbox(Object.keys(topo.objects), {\n label: \"layers\",\n value: [\"land\", \"alternatives\"]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4189,"value":"viewof preproject = Inputs.toggle({ label: \"preproject\" })","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4610,"value":"viewof rename_alternatives = Inputs.toggle({\n label: \"rename alternatives to countries\",\n value: true\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3003,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4618,"value":"mutable topto2 = null","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4621,"value":"topto2.objects.countries.geometries[0]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4379,"value":"topojson.mesh(\n topo,\n {\n type: \"GeometryCollection\",\n geometries: topo.objects.countries.geometries\n },\n (a, b) => a === b\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4574,"value":"Inputs.table(\n features.map(({ id, properties }) => ({\n id,\n a3: properties.a3,\n name: properties.name,\n ori: ori.get(id),\n DIFF: ori.get(id) !== properties.name\n // pieces: properties.pieces\n }))\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4577,"value":"original = d3.json(\n \"https://cdn.jsdelivr.net/npm/world-atlas@2/countries-110m.json\"\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4588,"value":"original.objects.countries.geometries.find(\n ({ id, properties }) => properties.name === \"Somaliland\"\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4579,"value":"ori = new Map(\n original.objects.countries.geometries.map(({ id, properties }) => [\n id,\n properties.name\n ])\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4387,"value":"features = topojson.feature(topo2, {\n type: \"GeometryCollection\",\n geometries: topo2.objects.countries.geometries\n}).features","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4526,"value":"features.filter((d) => d.properties.a3 === \"FRA\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3662,"value":"---","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3201,"value":"{\n const dl = d3.csvFormat(clean_properties);\n const size = (dl.length / 1024).toFixed(0);\n const compressed = (evaluate_compressed_size(dl) / 1024).toFixed(0);\n return DOM.download(\n () => new Blob([dl], { type: \"text/csv\" }),\n \"properties.csv\",\n `Download properties.csv (~${size}kB, zip: ~${compressed}kB)`\n );\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3130,"value":"$32","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2949,"value":"---\nMaps:","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":57,"value":"map({ projection_, visibility, center: false, displayObjects, sphere: false })","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4835,"value":"import { projectionInput2 } from \"@fil/d3-projections\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1293,"value":"// viewof projection // Input2\nviewof projection_ = projectionInput2({\n name: \"projection_\",\n value:\n new URLSearchParams(html``.search).get(\"projection_\") || \"Mercator\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":476,"value":"viewof displayObjects = Inputs.checkbox(Object.keys(topo.objects), {\n label: \"objects to display\",\n value: this?.value || [\"land\", \"alternatives\"]\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":850,"value":"### Abyei","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3705,"value":"Inputs.bind(controls.abyi(), viewof ctr_abyi)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":849,"value":"map({\n center: abyei,\n scale: 3000,\n displayObjects: [...displayObjects, \"disputed_abyei\", \"disputedMesh\"],\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":957,"value":"### Kashmir","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3717,"value":"Inputs.bind(controls.kash(), viewof ctr_kash)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":993,"value":"map({\n center: kashmir,\n scale: 3000,\n displayObjects: [...displayObjects, \"disputed_kashmir\", \"disputedMesh\"],\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4661,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1310,"value":"### Crimea","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3722,"value":"Inputs.bind(controls.crma(), viewof ctr_crma)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1308,"value":"map({\n center: { type: \"Point\", coordinates: [34.5, 46.2] },\n scale: 3000,\n displayObjects: [...displayObjects, \"disputed\", \"disputedMesh\"],\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1664,"value":"### Morocco / Western Sahara","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3712,"value":"Inputs.bind(controls.esh(), viewof ctr_esh)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1667,"value":"map({\n center: { type: \"Point\", coordinates: [-8, 27] },\n scale: 2000,\n displayObjects: [...displayObjects, \"disputed\", \"disputedMesh\"],\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1629,"value":"### Taiwan","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3727,"value":"Inputs.bind(controls.twn(), viewof ctr_twn)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1627,"value":"map({\n center: { type: \"Point\", coordinates: [120.9605, 23.6978] },\n scale: 1500,\n displayObjects,\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3366,"value":"### Guyane","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3697,"value":"Inputs.bind(controls.guf(), viewof ctr_guf)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3364,"value":"map({\n center: { type: \"Point\", coordinates: [-53.1258, 3.9339] },\n scale: 2500,\n displayObjects,\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3375,"value":"### Lakes (US/Canada)","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3380,"value":"map({\n center: { type: \"Point\", coordinates: [-85, 55] },\n scale: 700,\n displayObjects: [\"land\", \"countries\", \"countries_lakes\"],\n t: topo,\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3390,"value":"### Lakes (Africa)","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3394,"value":"map({\n center: { type: \"Point\", coordinates: [40, -10] },\n scale: 1200,\n displayObjects: [\"land\", \"countries\", \"countries_lakes\"],\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3445,"value":"### Kurils","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3443,"value":"map({\n center: { type: \"Point\", coordinates: [151.5, 46] },\n scale: 2000,\n displayObjects,\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3638,"value":"### Antemeridian cut","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3733,"value":"Inputs.bind(controls.geostitch(), viewof geostitch)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3636,"value":"map({\n center: { type: \"Point\", coordinates: [-180, 68] },\n scale: 3000,\n displayObjects: [\"land\", \"countries\", \"countries_lakes\"],\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":855,"value":"### Palestine & OPT","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":313,"value":"map({\n center: { type: \"Point\", coordinates: [34.3088, 31.3547] },\n scale: 5000,\n displayObjects,\n visibility\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4454,"value":"### Balkans (w/ Kosovo)","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":4452,"value":"map({\n center: { type: \"Point\", coordinates: [20, 44] },\n scale: 2500,\n displayObjects,\n visibility\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3882,"value":"Other options\n\nThe following options are here for development purposes. Experiment as you like, but they are probably all fine in their default values:","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":3884,"value":"viewof apply_properties = Inputs.checkbox(\n PIPELINE_PROPERTIES.map((d) => d.name),\n {\n label: \"properties\",\n value: PIPELINE_PROPERTIES.filter((d) => !d.disabled).map((d) => d.name)\n }\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3886,"value":"viewof apply_steps = Inputs.checkbox(\n PIPELINE.map((d) => d.name),\n {\n label: \"topological steps\",\n value: PIPELINE.filter((d) => !d.disabled).map((d) => d.name)\n }\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":409,"value":"viewof encoding = Inputs.select([\"utf-8\", \"windows-1252\"], { label: \"charset\" })","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":413,"value":"charset_control = properties.get(\"LBN\").name_ar","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":614,"value":"$34","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":111,"value":"----","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":518,"value":"o2017_countries110m = ({\n name: \"201707_countries110m\",\n zip: () => FileAttachment(\"201707-ne_110m_admin_0_countries.zip\").zip(),\n shpsource: \"ne_110m_admin_0_countries\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":522,"value":"o2018_countries110m = ({\n name: \"201808_countries110m\",\n zip: () => FileAttachment(\"201808-ne_110m_admin_0_countries.zip\").zip(),\n shpsource: \"ne_110m_admin_0_countries\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":549,"value":"o2019_countries110m = ({\n name: \"2019xx_countries110m\",\n zip: () => FileAttachment(\"2019xx-ne_110m_admin_0_countries.zip\").zip(),\n shpsource: \"ne_110m_admin_0_countries\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":552,"value":"o2021_countries110m = ({\n name: \"202112_countries110m\",\n zip: () => FileAttachment(\"202112-ne_110m_admin_0_countries.zip\").zip(),\n shpsource: \"ne_110m_admin_0_countries\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":557,"value":"countries110m = ({\n name: \"countries110m\",\n zipsource:\n \"https://naturalearth.s3.amazonaws.com/110m_cultural/ne_110m_admin_0_countries.zip\",\n shpsource: \"ne_110m_admin_0_countries\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":482,"value":"countries50m = ({\n name: \"countries50m\",\n zipsource:\n \"https://naturalearth.s3.amazonaws.com/50m_cultural/ne_50m_admin_0_countries.zip\",\n zip: () => FileAttachment(\"ne_50m_admin_0_countries.zip\").zip(),\n shpsource: \"ne_50m_admin_0_countries\"\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1290,"value":"import { viewof projection } with { projections } from \"@fil/d3-projections\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1519,"value":"$35","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3800,"value":"$36","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":45,"value":"localGeojson = () =>\n o2021_countries110m\n .zip()\n .then((z) =>\n Promise.all([\n z.file(`${o2021_countries110m.shpsource}.shp`).arrayBuffer(),\n z.file(`${o2021_countries110m.shpsource}.dbf`).arrayBuffer()\n ])\n )\n .then(([shp, dbf]) => shapefile.read(shp, dbf, { encoding }))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":515,"value":"o2017Geojson = () =>\n o2017_countries110m\n .zip()\n .then((z) =>\n Promise.all([\n z.file(`${o2017_countries110m.shpsource}.shp`).arrayBuffer(),\n z.file(`${o2017_countries110m.shpsource}.dbf`).arrayBuffer()\n ])\n )\n .then(([shp, dbf]) => shapefile.read(shp, dbf, { encoding }))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":524,"value":"o2018Geojson = () =>\n o2018_countries110m\n .zip()\n .then((z) =>\n Promise.all([\n z.file(`${o2018_countries110m.shpsource}.shp`).arrayBuffer(),\n z.file(`${o2018_countries110m.shpsource}.dbf`).arrayBuffer()\n ])\n )\n .then(([shp, dbf]) => shapefile.read(shp, dbf, { encoding }))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":562,"value":"o2019Geojson = () =>\n o2019_countries110m\n .zip()\n .then((z) =>\n Promise.all([\n z.file(`${o2019_countries110m.shpsource}.shp`).arrayBuffer(),\n z.file(`${o2019_countries110m.shpsource}.dbf`).arrayBuffer()\n ])\n )\n .then(([shp, dbf]) => shapefile.read(shp, dbf, { encoding }))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2429,"value":"tiny_both = async () => {\n const geojsons = {\n countries: await FileAttachment(\"tiny_countries@2.json\").json(),\n countries_lakes: await FileAttachment(\"tiny_countries_lakes@2.json\").json()\n };\n\n return topojson.topology(geojsons);\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3340,"value":"localGeojson50m = () =>\n Promise.all(\n [\n FileAttachment(\"ne_50m_admin_0_countries.zip\"),\n FileAttachment(\"ne_50m_admin_0_countries_lakes.zip\"),\n FileAttachment(\"ne_50m_admin_0_breakaway_disputed_areas.zip\")\n ].map((d) =>\n d\n .zip()\n .then((z) =>\n Promise.all([\n z.file(z.filenames.find((d) => d.endsWith(\".shp\"))).arrayBuffer(),\n z.file(z.filenames.find((d) => d.endsWith(\".dbf\"))).arrayBuffer()\n ])\n )\n .then(([shp, dbf]) => shapefile.read(shp, dbf, { encoding }))\n )\n ).then(([countries, countries_lakes, disputed]) =>\n topojson.topology({\n countries,\n countries_lakes,\n disputed\n })\n )","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1524,"value":"SOURCES = [\n {\n name: \"110m\",\n description: \"retrieves both 'tiny' countries with and without lakes\",\n process: tiny_both,\n disabled: false\n },\n {\n name: \"50m\",\n description: \"official zip from Natural Earth\",\n process: localGeojson50m,\n disabled: true\n }\n]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2199,"value":"$37","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":216,"value":"$38","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4099,"value":"mutable berm = null // will be used to add to the disputedMesh","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3966,"value":"mutable debug = null","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2260,"value":"function featureTransform(a, transform) {\n return topojson.topology(\n Object.fromEntries(\n Object.entries(a.objects).map(([key, o]) => [\n key,\n transform(topojson.feature(a, o), key)\n ])\n )\n );\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":712,"value":"who_regions = FileAttachment(\"who_regions.csv\")\n .csv()\n .then(\n (data) =>\n new Map([\n ...data.map((d) => [d.ISO_3_CODE, d.WHO_REGION]),\n [\"WEST\", \"AFRO\"], // like Western Sahara\n [\"XKX\", \"EURO\"], // Kosovo…\n [\"SOL\", \"EMRO\"], // Somaliland = Same WHO region as SOMALIA\n [\"PSX\", \"EMRO\"], // Palestine\n [\"GUF\", \"AMRO\"], // French Guyane is a member of PAHO\n [\"ABYI\", \"EMRO\"], // like Sudan\n [\"CRMA\", \"EURO\"], // Crimea\n [\"AKSA\", \"WPRO\"], // like China\n [\"GILG\", \"EMRO\"], // like Pakistan\n [\"JAMM\", \"SEARO\"], // like India\n [\"TWN\", \"WPRO\"] // Taiwan is either WPRO or UNKNOWN…\n ])\n )","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4008,"value":"colors = [].concat(d3.schemeTableau10).slice(0, 9)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":218,"value":"$39","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2352,"value":"outTopo = worlds[worlds.length - 1]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2354,"value":"properties = worlds.properties","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2444,"value":"countries = topojson.feature(outTopo, outTopo.objects.countries)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2461,"value":"countries_lakes = outTopo.objects.countries_lakes &&\n topojson.feature(outTopo, outTopo.objects.countries_lakes)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2786,"value":"pieces = outTopo.objects.pieces && topojson.feature(outTopo, outTopo.objects.pieces)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":792,"value":"disputed = outTopo.objects.disputed &&\n topojson.feature(outTopo, outTopo.objects.disputed)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":425,"value":"$3a","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1985,"value":"land = ({\n type: \"GeometryCollection\",\n geometries: [\n Object.assign(\n topojson.mergeArcs(baseTopo, baseTopo.objects.pieces.geometries),\n { properties: { name: \"land\" }, id: \"land\" }\n )\n ]\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4550,"value":"land_lakes = ({\n type: \"GeometryCollection\",\n geometries: [\n Object.assign(\n topojson.mergeArcs(baseTopo, baseTopo.objects.countries_lakes.geometries),\n { properties: { name: \"land_lakes\" }, id: \"land_lakes\" }\n )\n ]\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4202,"value":"ocean = {\n const ocean = JSON.parse(JSON.stringify(land.geometries[0]));\n ocean.id = \"ocean\";\n ocean.properties = { name: \"ocean\", fill: \"#add8e6\" };\n\n ocean.arcs = [\n ocean.arcs\n .sort(([a], [b]) => b.length - a.length) // Eurasia, Americas, …\n .map(([d]) =>\n d\n .slice()\n .reverse()\n .map((d) => -d - 1)\n ),\n ocean.arcs.flatMap((\n [, ...holes] // Caspian Sea\n ) =>\n holes.map((d) =>\n d\n .slice()\n .reverse()\n .map((d) => -d - 1)\n )\n )\n ];\n return { type: \"GeometryCollection\", geometries: [ocean] };\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1989,"value":"continents = ({\n type: \"GeometryCollection\",\n geometries: [\n ...d3.rollup(\n baseTopo.objects.pieces.geometries,\n (v) => topojson.mergeArcs(baseTopo, v),\n (d) => properties.get(d.id)?.continent || `UNKNOWN: ${d.id}`\n )\n ].map(([continent, o], i) =>\n Object.assign(o, {\n id: continent,\n properties: { name: continent, fill: d3.schemeAccent[i % 10] }\n })\n )\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1992,"value":"region_un = ({\n type: \"GeometryCollection\",\n geometries: [\n ...d3.rollup(\n baseTopo.objects.pieces.geometries,\n (v) => topojson.mergeArcs(baseTopo, v),\n (d) => properties.get(d.id)?.region_un\n )\n ].map(([region_un, o], i) =>\n Object.assign(o, {\n id: region_un,\n properties: { name: region_un, fill: d3.schemeAccent[i % 10] }\n })\n )\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1995,"value":"region_who = ({\n type: \"GeometryCollection\",\n geometries: [\n ...d3.rollup(\n baseTopo.objects.pieces.geometries,\n (v) => topojson.mergeArcs(baseTopo, v),\n (d) => who_regions.get(d.id) || `UNKNOWN: ${d.id}`\n )\n ].map(([region_who, o], i) =>\n Object.assign(o, {\n id: region_who,\n properties: { name: region_who, fill: d3.schemeAccent[i % 10] }\n })\n )\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2001,"value":"region_wb = ({\n type: \"GeometryCollection\",\n geometries: [\n ...d3.rollup(\n baseTopo.objects.pieces.geometries,\n (v) => topojson.mergeArcs(baseTopo, v),\n (d) => properties.get(d.id)?.region_wb || `UNKNOWN: ${d.id}`\n )\n ].map(([region_wb, o], i) =>\n Object.assign(o, {\n id: region_wb,\n properties: { name: region_wb, fill: d3.schemeAccent[i % 10] }\n })\n )\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2008,"value":"disputed_kashmir = baseTopo.objects.disputed && {\n type: \"GeometryCollection\",\n geometries: [\n Object.assign(\n topojson.mergeArcs(\n baseTopo,\n baseTopo.objects.disputed.geometries.filter((d) =>\n [\"AKSA\", \"JAMM\", \"GILG\", \"KASH\", \"KAS\"].includes(d.id)\n )\n ),\n { id: \"KASH\", properties: { name: \"KASH\" } }\n )\n ]\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2614,"value":"disputed_abyei = baseTopo.objects.disputed && {\n type: \"GeometryCollection\",\n geometries: [\n Object.assign(\n topojson.mergeArcs(\n baseTopo,\n baseTopo.objects.disputed.geometries.filter((d) => d.id === \"ABYI\")\n ),\n { id: \"ABYI\", properties: { name: \"Abyei\" } }\n )\n ]\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2384,"value":"mesh = (topo, objects, name = \"mesh\") => ({\n type: \"GeometryCollection\",\n geometries: [\n Object.assign(\n topojson.meshArcs(topo, objects, (a, b) => a !== b),\n { id: name, properties: { name } }\n )\n ]\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2016,"value":"taggedMesh = taggedBorders(baseTopo, baseTopo.objects.countries.geometries, (a, b) =>\n a && b && a !== b ? [a.id, b.id].sort().join(\".\") : null\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2020,"value":"// TODO: for ESH.MAR, this could include the opposite shapes (?)\ndisputedMesh = {\n if (!baseTopo.objects.pieces) return;\n\n const borders = taggedBorders(\n baseTopo,\n baseTopo.objects.pieces.geometries,\n (a, b) => (a && b && a !== b ? [a.id, b.id].sort().join(\".\") : null)\n );\n const disput = taggedBorders(\n baseTopo,\n baseTopo.objects.disputed.geometries,\n (a, b) => (a && b && a !== b ? [a.id, b.id].sort().join(\".\") : null)\n );\n return {\n type: \"GeometryCollection\",\n geometries: [\n ...[baseTopo.objects.berm] /\"ESH.MAR\"/,\n ...borders.geometries\n .filter((d) =>\n [\"CRMA.UKR\", /\"ESH.MAR\",/ \"ABYI.SDN\"].includes(d.properties.name)\n )\n .map((d) => ((d.arcs = [d.arcs[0]]), d)),\n ...disput.geometries.filter((d) => d.properties.name === \"GILG.JAMM\")\n ]\n };\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":451,"value":"$3b","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":4314,"value":"$3c","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":581,"value":"tests = {\n let a = countries;\n const tests = [];\n for (const { name, test, disabled } of PIPELINE) {\n if (test && !disabled) {\n const result = await test(a);\n tests.push({ name, result });\n }\n }\n return tests;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":689,"value":"// detach Guyane from France when it makes sense,\nfunction detachGuyane(topo) {\n return detachPolygon(\n topo,\n topo.objects.countries.geometries,\n d => d.id === \"FRA\",\n [-53.1258, 3.9339],\n { region_un: \"Americas\", continent: \"South America\" },\n \"GUF\"\n );\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1382,"value":"// option to detach Crimea from Russia\nfunction detachCrma(topo, geometries, reason) {\n return detachPolygon(\n topo,\n geometries,\n (d) => d.id === \"RUS\",\n [34.4997, 45.3453],\n reason,\n \"CRMA\"\n );\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1324,"value":"function detachPolygons(geometries, test, indices, reason = {}, id) {\n return geometries.flatMap((d) => {\n if (test(d)) {\n return [\n {\n id,\n type: \"MultiPolygon\",\n arcs: d.arcs.filter((d, i) => indices.includes(i)),\n properties: { ...d.properties, ...reason }\n },\n {\n id: d.id,\n type: \"MultiPolygon\",\n arcs: d.arcs.filter((d, i) => !indices.includes(i)),\n properties: d.properties\n }\n ].map((d) =>\n d.arcs.length === 1 ? { ...d, type: \"Polygon\", arcs: d.arcs[0] } : d\n );\n } else return [d];\n });\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2292,"value":"function detachPolygon(topo, geometries, test, coords, reason = {}, id) {\n return geometries.flatMap((d) => {\n if (test(d)) {\n const centroids = topojson\n .feature(topo, {\n type: \"GeometryCollection\",\n geometries: d.arcs.map((arc) => ({\n type: \"Polygon\",\n arcs: arc\n }))\n })\n .features.map(d3.geoCentroid);\n const closest = d3.leastIndex(centroids, (c) =>\n d3.geoDistance(c, coords)\n );\n\n return [\n {\n id,\n type: \"MultiPolygon\",\n arcs: d.arcs.filter((d, i) => i === closest),\n properties: { ...d.properties, ...reason }\n },\n {\n id: d.id,\n type: \"MultiPolygon\",\n arcs: d.arcs.filter((d, i) => i !== closest),\n properties: d.properties\n }\n ].map((d) =>\n d.arcs.length === 1 ? { ...d, type: \"Polygon\", arcs: d.arcs[0] } : d\n );\n } else return [d];\n });\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":955,"value":"$3d","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1449,"value":"// https://bl.ocks.org/mbostock/8ba407a7a53d80be62a8d54774663c2f\n// Usage:\n// topojson.mesh(topo, topo.objects.countries, border(\"Ukraine\", \"Crimea\")))\nfunction border(id0, id1) {\n return function (a, b) {\n return (\n (a.properties.name === id0 && b.properties.name === id1) ||\n (a.properties.name === id1 && b.properties.name === id0)\n );\n };\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":821,"value":"// find the closest coordinates of point x in the target geojson\ncoordsMatcher = (geo, { precision = 0.0003, snap } = {}) => {\n const c = geoVertices(geo);\n const voronoi = d3.geoVoronoi(c);\n return function (points) {\n let snapped = 0;\n for (let i = 0; i < points.length; i++) {\n const point = points[i];\n const candidate = c[voronoi.find(...point)];\n const distance = d3.geoDistance(candidate, point);\n const ok = distance < precision;\n snapped += ok;\n if (ok) points[i] = candidate;\n }\n if (false && snapped !== snap) {\n if (snap != null) throw new Error(`points snapped: ${snapped}/${snap}`);\n else console.info(`points snapped: ${snapped}`);\n }\n return points;\n };\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":941,"value":"$3e","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":777,"value":"abyei = ({\n type: \"Feature\",\n id: \"ABYI\",\n properties: { name: \"Abyei\" },\n geometry: {\n type: \"Polygon\",\n coordinates: [\n [\n [27.850341796875, 10.17437402751379],\n [29.003906249999996, 10.163560279490476],\n [29.0093994140625, 9.606166114941981],\n [28.965454101562504, 9.400290685848121],\n [27.9656982421875, 9.400290685848121],\n [27.833862304687496, 9.606166114941981],\n [27.850341796875, 10.17437402751379]\n ]\n ]\n }\n})","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":206,"value":"---\nLibs","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":61,"value":"d3 = require(\"d3@7\", \"d3-geo-projection@3\", \"d3-geo-voronoi@1\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":210,"value":"topojson = require(\"topojson-server@3\", \"topojson-simplify@3\", \"topojson-client@3\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":19,"value":"shapefile = require(\"shapefile@0.6\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":171,"value":"import { remote } from \"@observablehq/remote\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":335,"value":"import { diff } from \"@jobleonard/diff\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":396,"value":"$3f","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":827,"value":"geoVertices = function (feature, tag) {\n const vertices = [];\n const stream = {\n point:\n tag === undefined\n ? (lambda, phi) => vertices.push([lambda, phi])\n : (lambda, phi) => vertices.push([lambda, phi, tag]),\n lineStart: function () {},\n lineEnd: function () {},\n polygonStart: function () {},\n polygonEnd: function () {}\n };\n d3.geoStream(feature, stream);\n return vertices;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3927,"value":"$40","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1679,"value":"// See https://observablehq.com/@fil/hello-pclip\n// https://github.com/substack/pclip\npclip = Promise.all([\n import(\"https://cdn.skypack.dev/pclip\"),\n import(\"https://cdn.skypack.dev/pclip/geo\").then((d) => d.default),\n import(\"https://cdn.skypack.dev/pclip/xy\").then((d) => d.default)\n]).then(([pclip, geo, xy]) =>\n Object.assign(pclip.default, { ...pclip, geo, xy })\n)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1868,"value":"function remove_duplicate_points(ring) {\n const r = [ring[0]];\n for (let i = 1; i < ring.length; i++) {\n if (d3.geoDistance(ring[i], ring[i - 1]) > 0.000001) r.push(ring[i]);\n }\n r.push(ring[0].slice());\n return r;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3826,"value":"fflate = import(\"https://cdn.skypack.dev/fflate?min\")","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3846,"value":"function evaluate_compressed_size(text) {\n if (typeof text !== \"string\") text = JSON.stringify(text);\n return fflate.zipSync({ test: fflate.strToU8(text) }).length;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3992,"value":"islands = FileAttachment(\"island-nations@1.csv\").csv({ typed: true })","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W41","stdlib":"1","titling":"auto"}} 30:"$2f:body"