25:["$","title",null,{"children":"Leaflet POC before looping | mheard03 | Observable"}] 49:T561,map = { const height = width * mapSize.height / mapSize.width; const container = html`

`; yield container; // const corner1 = L.latLng(-90, -180) // const corner2 = L.latLng(90, 180) // const maxBounds = L.latLngBounds(corner1, corner2) const renderer = new L.SVG({ padding: 0 }); const map = L.map(container, { zoomAnimation: false, zoomSnap: 0, topoMove: true, // maxBounds: [[-90, -Infinity], [90, Infinity]], worldCopyJump: true, renderer }); const topo = new LeafletTopology(countries, map); const topoLayer = topo.createLayer(countries.objects.countries); topoLayer.addTo(map); mutable topoRef = topo; // const topoLayer = new TopoJsonLayer(countries, countries.objects.countries, { padding: 1 }); // topoLayer.addTo(map); // const bounds = L.latLngBounds(L.latLng(-90, -180), L.latLng(90, 180)); // map.fitBounds(bounds); map.setMaxBounds([[-90, -360], [90, 360]]) function fitWorld() { const zoomY = map.getBoundsZoom(L.latLngBounds(L.latLng(-90, 0), L.latLng(90, 0))); const zoomX = map.getBoundsZoom(L.latLngBounds(L.latLng(0, -180), L.latLng(0, 180))); map.setMinZoom(Math.max(zoomX, zoomY)); } fitWorld(); map.fitBounds([[-90, -360], [90, 360]]); invalidation.then(() => map.remove()); container.map = map; yield container; }4a:T825,ClippedSVG = { const { stamp, splitWords } = L.Util; const { create } = L.SVG; const parent = L.SVG.prototype; const uniqueLayers = new Map(); return L.SVG.extend({ _setPath(layer, path) { layer._path.setAttribute('d', path); layer._path.setAttribute('d', path); }, options: { padding: 0.1, continuous: true, }, _initPath(layer) { const layerStamp = stamp(layer); layer.options.noClip = true; layer.options.smoothFactor = 0; const clipPath = layer._clipPath = create('clipPath'); clipPath.setAttribute("id", `clip_${layerStamp}`); const path = layer._path = create('path'); path.setAttribute("id", `path_${layerStamp}`); clipPath.append(path); const instances = layer._instances = []; const use = layer._use = create('use'); use.setAttribute("href", "#" + layer._path.id); if (layer.options.className) { path.classList.add(...splitWords(layer.options.className)); } if (layer.options.interactive) { path.classList.add('leaflet-interactive'); } this._updateStyle(layer); this._layers[layerStamp] = layer; }, _addPath(layer) { if (!this._rootGroup || !this._defs) { this._initContainer(); } this._defs.appendChild(layer._clipPath); this._rootGroup.appendChild(layer._use); layer.addInteractiveTarget(layer._use); }, _removePath(layer) { layer._clipPath.remove(); layer._path.remove(); layer.removeInteractiveTarget(layer._use); delete this._layers[stamp(layer)]; }, _updateStyle(layer) { parent._updateStyle.call(this, layer); const use = layer._use, options = layer.options; if (options.clip || options.clip === undefined) { use.setAttribute('clip-path', `url(#${layer._clipPath.id})`); } else { use.removeAttribute('clip-path'); } }, _bringToFront(layer) { L.DomUtil.toFront(layer._use); }, _bringToBack(layer) { L.DomUtil.toBack(layer._use); } }); }4b:T8ac,_ClippedSVG = { const { stamp, splitWords } = L.Util; const { create } = L.SVG; const parent = L.SVG.prototype; const uniqueLayers = new Map(); return L.SVG.extend({ options: { padding: 0.1, continuous: true, }, _initContainer() { parent._initContainer.call(this); this._defs = create("defs"); this._container.prepend(this._defs); }, _destroyContainer() { parent._destroyContainer.call(this); delete this._defs; }, _initPath(layer) { const layerStamp = stamp(layer); layer.options.noClip = true; layer.options.smoothFactor = 0; const clipPath = layer._clipPath = create('clipPath'); clipPath.setAttribute("id", `clip_${layerStamp}`); const path = layer._path = create('path'); path.setAttribute("id", `path_${layerStamp}`); clipPath.append(path); const instances = layer._instances = []; const use = layer._use = create('use'); use.setAttribute("href", "#" + layer._path.id); if (layer.options.className) { path.classList.add(...splitWords(layer.options.className)); } if (layer.options.interactive) { path.classList.add('leaflet-interactive'); } this._updateStyle(layer); this._layers[layerStamp] = layer; }, _addPath(layer) { if (!this._rootGroup || !this._defs) { this._initContainer(); } this._defs.appendChild(layer._clipPath); this._rootGroup.appendChild(layer._use); layer.addInteractiveTarget(layer._use); }, _removePath(layer) { layer._clipPath.remove(); layer._path.remove(); layer.removeInteractiveTarget(layer._use); delete this._layers[stamp(layer)]; }, _updateStyle(layer) { parent._updateStyle.call(this, layer); const use = layer._use, options = layer.options; if (options.clip || options.clip === undefined) { use.setAttribute('clip-path', `url(#${layer._clipPath.id})`); } else { use.removeAttribute('clip-path'); } }, _bringToFront(layer) { L.DomUtil.toFront(layer._use); }, _bringToBack(layer) { L.DomUtil.toBack(layer._use); } }); }4c:T162d,{ // Fork of https://github.com/perliedman/leaflet-realtime L.Realtime = L.Layer.extend({ options: { getFeatureId: function(f) { return f.properties.id; }, updateFeature: function(feature, oldLayer) { if (!oldLayer) return; var type = feature.geometry && feature.geometry.type var coordinates = feature.geometry && feature.geometry.coordinates switch (type) { case 'Point': oldLayer.setLatLng(L.GeoJSON.coordsToLatLng(coordinates)); break; case 'LineString': case 'MultiLineString': oldLayer.setLatLngs(L.GeoJSON.coordsToLatLngs(coordinates, type === 'LineString' ? 0 : 1)); break; case 'Polygon': case 'MultiPolygon': oldLayer.setLatLngs(L.GeoJSON.coordsToLatLngs(coordinates, type === 'Polygon' ? 1 : 2)); break; default: return null; } return oldLayer; }, logErrors: true, removeMissing: false }, initialize: function(geojson, options) { options ??= {}; L.setOptions(this, options); this._container = options.container || L.geoJson(null, options); this._features = {}; this._featureLayers = {}; this.update(geojson); }, update: function(geojson) { this._onNewData(false, geojson); }, remove: function(geojson) { var features = L.Util.isArray(geojson) ? geojson : geojson.features ? geojson.features : [geojson], exit = {}, i, len, fId; for (i = 0, len = features.length; i < len; i++) { fId = this.options.getFeatureId(features[i]); this._container.removeLayer(this._featureLayers[fId]); exit[fId] = this._features[fId]; delete this._features[fId]; delete this._featureLayers[fId]; } this.fire('update', { features: this._features, enter: {}, update: {}, exit: exit }); return this; }, getLayer: function(featureId) { return this._featureLayers[featureId]; }, getFeature: function(featureId) { return this._features[featureId]; }, getBounds: function() { var container = this._container; if (container.getBounds) { return container.getBounds(); } throw new Error('Container has no getBounds method'); }, onAdd: function(map) { map.addLayer(this._container); }, onRemove: function(map) { map.removeLayer(this._container); }, _onNewData: function(removeMissing, geojson) { var layersToRemove = [], enter = {}, update = {}, exit = {}, seenFeatures = {}, i, len, feature; var handleData = L.bind(function(geojson) { var features = L.Util.isArray(geojson) ? geojson : geojson.features; if (features) { for (i = 0, len = features.length; i < len; i++) { // only add this if geometry or geometries are set and not null feature = features[i]; if (feature.geometries || feature.geometry || feature.features || feature.coordinates) { handleData(feature); } } return; } var container = this._container; var options = this.options; if (options.filter && !options.filter(geojson)) { return; } var f = L.GeoJSON.asFeature(geojson); var fId = options.getFeatureId(f); var oldLayer = this._featureLayers[fId]; var layer = this.options.updateFeature(f, oldLayer); if (!layer) { layer = L.GeoJSON.geometryToLayer(geojson, options); if (!layer) { return; } layer.defaultOptions = layer.options; layer.feature = f; if (options.onEachFeature) { options.onEachFeature(geojson, layer); } if (options.style && layer.setStyle) { layer.setStyle(options.style(geojson)); } } layer.feature = f; if (container.resetStyle) { container.resetStyle(layer); } if (oldLayer) { update[fId] = geojson; if (oldLayer != layer) { layersToRemove.push(oldLayer); container.addLayer(layer); } } else { enter[fId] = geojson; container.addLayer(layer); } this._featureLayers[fId] = layer; this._features[fId] = seenFeatures[fId] = f; }, this); handleData(geojson); if (removeMissing) { exit = this._removeUnknown(seenFeatures); } for (i = 0; i < layersToRemove.length; i++) { this._container.removeLayer(layersToRemove[i]); } this.fire('update', { features: this._features, enter: enter, update: update, exit: exit }); }, _onError: function(err, msg) { if (this.options.logErrors) { console.warn(err, msg); } this.fire('error', { error: err, message: msg }); }, _removeUnknown: function(known) { var fId, removed = {}; for (fId in this._featureLayers) { if (!known[fId]) { this._container.removeLayer(this._featureLayers[fId]); removed[fId] = this._features[fId]; delete this._featureLayers[fId]; delete this._features[fId]; } } return removed; } }); L.realtime = function(geojson, options) { return new L.Realtime(geojson, options); }; }4d:T4c1,PolylineArc = { const levelsDeep = { LineString: 0, MultiLineString: 1, Polygon: 1, MultiPolygon: 2 }; return L.Polyline.extend({ ...polyPathOverrides, options: { smoothFactor: 0, noClip: true }, initialize(topology, geometry, options) { L.Util.setOptions(this, options); this._setTopology(topology); this._setGeometry({ type: geometry?.type, arcs: geometry?.arcs, properties: geometry?.properties ?? {} }); this.__latlngs = false; Object.defineProperty(this, "_latlngs", { get() { if (this.__latlngs === false) { if (!this._topology || !this._geometry) return false; const feature = topojson.feature(this._topology, this._geometry); const latlngs = L.GeoJSON.coordsToLatLngs(feature.geometry.coordinates, levelsDeep[geometry.type]); this._setLatLngs(latlngs); } return this.__latlngs; }, set(value) { this.__latlngs = value; this.redraw(); }, enumerable: true, configurable: true }) }, /* _updatePath: function() { this.__latlngs = false; } */ }); }4e:T1a13,LeafletTopology = { function coordsEqual(a, b) { if (a === b) return true; return (a.length == b.length && a[0] === b[0] && a[1] === b[1]); } function arcsEqual(a, b) { if (a === b) return true; if (!a || !b) return false; // This isn't a heuristic, it's just how topojson works return (a.length == b.length && coordsEqual(a[0], b[0]) && coordsEqual(a[1], b[1])); } function getBasePxArea(crs, latLngBounds) { const pointNw = crs.latLngToPoint(latLngBounds.getNorthWest(), 1); const pointSe = crs.latLngToPoint(latLngBounds.getSouthEast(), 1); const bounds = L.bounds(pointNw, pointSe); const size = bounds.getSize(); return size.x * size.x + size.y * size.y; } class LeafletTopology { constructor(topology, map) { this._map = map; this._basePxArea = undefined; this._ringPxAreas = new Map(); this._arcLayerIds = new Map(); this._loadTopology(topology); const fnUpdateArcs = this._updateArcs.bind(this); const startUpdating = (function() { if (this._layerUpdateInterval) return; this._layerUpdateInterval = setInterval(fnUpdateArcs, 100); }).bind(this); const stopUpdating = (function() { clearInterval(this._layerUpdateInterval); this._layerUpdateInterval = undefined; this._updateArcs(); }).bind(this); map.on("movestart", startUpdating); map.on("moveend", stopUpdating); map.once("unload", function() { this._map.off("movestart", startUpdating); this._map.off("moveend", stopUpdating); this._map = undefined; stopUpdating(); }, this); } createLayer(geometry, options) { if (!geometry) return null; if (geometry.type == "GeometryCollection") { const layers = geometry.geometries.map(g => this.createLayer(g, options)).filter(g => g); return new L.FeatureGroup(layers); } let layer; switch (geometry.type) { case 'LineString': case 'MultiLineString': layer = new PolylineArc(this, geometry, options); break; case 'Polygon': case 'MultiPolygon': layer = new PolygonArc(this, geometry, options); break; default: throw new Error('Invalid/unsupported GeoJSON object.'); } const stamp = L.Util.stamp(layer); const arcIds = geometry.arcs.flat().flat().map(i => (i < 0) ? ~i : i); arcIds.forEach(arcId => this._arcLayerIds.get(arcId).add(stamp)); return layer; } _updateArcs() { const map = this._map; if (!map) return; const bounds = map.getBounds(); const visiblePxArea = getBasePxArea(this._getCrs(), bounds); const pctVisible = visiblePxArea / this._basePxArea; const nw = bounds.getNorthWest(); const se = bounds.getSouthEast(); const bbox = Bbox.fromPoints([[nw.lng, nw.lat], [se.lng, se.lat]]); const visibleArcIds = new Set(this._index.search(...bbox)); // TODO: options.smoothFactor const smoothFactor = 1; const tolerance = Math.min(this._lowResFloorPx, smoothFactor * pctVisible); const dirtyArcIds = new Set(); this._lowResArcs.forEach((arc, i) => { if (visibleArcIds.has(i)) { arc = this._sourceArcs[i].filter(c => c[2] >= tolerance).map(c => [c[0], c[1]]); } if (!arcsEqual(arc, this.arcs[i])) dirtyArcIds.add(i); this.arcs[i] = arc; }); const stampsToUpdate = new Set(); for (let arcId of dirtyArcIds) { for (let stamp of this._arcLayerIds.get(arcId)) { stampsToUpdate.add(stamp); } } console.log(`Redrawing ${stampsToUpdate.size} objects @ tolerance ${smoothFactor / pctVisible}; ${this.arcs.flat().length}/${this._sourceArcs.flat().length} points in topology.`); this._map.eachLayer(l => { if (l.markDirty && stampsToUpdate.has(L.Util.stamp(l))) { l.markDirty(); } }); } _getCrs() { return this._map?.options.crs; } _loadTopology(topology) { this.type = topology.type; this.objects = topology.objects; // Unquantize and create _sourceArcs const transformCoord = topojson.transform(topology.transform); const sourceArcs = this._sourceArcs = topology.arcs.map(arc => { return arc.map(c => transformCoord(c)); }); // Build index const index = this._index = new Flatbush(sourceArcs.length); sourceArcs.forEach(arc => index.add(...Bbox.fromPoints(arc))); index.finish(); // Presimplify the topology this.presimplifyTopology(); // Build _lowResArcs and arcs this._lowResArcs = sourceArcs.map(a => a.filter(c => c[2] >= this._lowResFloorPx).map(c => [c[0], c[1]])); this.arcs = [...this._lowResArcs]; // Initialize this._arcLayerIds for (let arcId = 0; arcId < sourceArcs.length; arcId++) { if (this._arcLayerIds.has(arcId)) continue; this._arcLayerIds.set(arcId, new Set()); } console.log(`Topology loaded. ${sourceArcs.length} arcs, ${sourceArcs.flat().length} points.`); } presimplifyTopology() { const sourceArcs = this._sourceArcs; // Presimplify a temporary copy of this topology, with lat/lng converted to screen coords at zoom level 1 const crs = this._getCrs(); const presimplified = topojson.presimplify({ type: "Topology", arcs: sourceArcs.map(arc => { return arc.map(([lng,lat,value]) => { const { x, y } = crs.latLngToPoint(L.latLng(lat, lng), 1); return [x, y]; }); }) }); // Calculate this._basePxArea, the area (in pixels) of the topology at zoom level 1 const [n, e, s, w] = Bbox.fromPoints(presimplified.arcs.flat()); this._basePxArea = getBasePxArea(crs, L.latLngBounds(L.latLng(n, w), L.latLng(s, e))); // Update sourceArcs to include the new presimplification scores, expressed as a fraction of this._basePxArea sourceArcs.forEach((arc, i) => { arc.forEach((coord, j) => { if (coord[2] == Infinity) return; coord[2] = presimplified.arcs[i][j][2]; }); }); const allScores = sourceArcs.flat().map(c => c[2]).filter(c => c > 0 && isFinite(c)); allScores.sort((a,b) => a-b); const lowResFloorIndex = Math.floor(allScores.length - simple.arcs.length / 3); this._lowResFloorPx = allScores[lowResFloorIndex]; } } return LeafletTopology; }4f:[] 29:{"id":"9d8f72354379a70d","slug":null,"title":"Leaflet POC before looping","thumbnail":"49e3b8e7b5c25c76853ac31b4ed41e895e2facb91573182f4e68e995d0ccb899","access_level":"public","create_time":"$D2023-10-25T15:31:28.224Z","update_time":"$D2023-10-25T15:31:39.144Z","view_time":null,"first_public_time":"$D2023-10-25T15:31:28.224Z","has_unpublished_edits":null,"fork":{"id":"7e72891a23b5d28c","slug":null,"title":"Leaflet POC","owner":{"id":"02d4ccf8951f2f02","login":"mheard03","name":""}},"stars":0,"starred":false,"owner":{"type":"team","id":"02d4ccf8951f2f02","login":"mheard03","name":"","avatar_url":"https://avatars-next.observableusercontent.com/ae42ada5dc924d7cd620a4ae3bb6279a3e82fb649668ee6c0729bdada1ebf016","tier":"starter_2024","bio":"","following":false},"version":2891,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"6afb6cac5f188ec0","login":"mheard03","name":"","avatar_url":"https://avatars-next.observableusercontent.com/ae42ada5dc924d7cd620a4ae3bb6279a3e82fb649668ee6c0729bdada1ebf016","tier":null,"bio":"","following":false}],"members":[{"type":"individual","id":"6afb6cac5f188ec0","login":"mheard03","name":"","avatar_url":"https://avatars-next.observableusercontent.com/ae42ada5dc924d7cd620a4ae3bb6279a3e82fb649668ee6c0729bdada1ebf016","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"titling":"auto","title":"Leaflet POC before looping","cells":[{"id":0,"value":"# Leaflet POC before looping","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":464,"value":"// worldCopyJump","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":354,"value":"{\n const layers = [];\n map.map.eachLayer(l => layers.push(l));\n return layers;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1043,"value":"m = new L.Map(DOM.element(\"div\"));","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1125,"value":"map.map.getCenter()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2435,"value":"map.map.getZoom()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1110,"value":"map.map.flyTo({lat: -0.6364206182756988, lng: -49.91245580521718}, 7.228212509167584);","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2451,"value":"topoRef","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2446,"value":"mutable topoRef = false;","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2479,"value":"crs.scale(map.map.getZoom())","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2613,"value":"layers.filter(l => l.markDirty)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2609,"value":"layers = {\n const layers = [];\n map.map.eachLayer(l => layers.push(l));\n return layers;\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2618,"value":"topoRef._sourceArcs.flat().filter(c => c[2] >= 1)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2616,"value":"{\n const path = [...document.querySelectorAll(\"path.leaflet-interactive\")].filter(p => p.getAttribute(\"d\").length > 150)[0]\n return path.getAttribute(\"d\")\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":49,"value":"$49","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1057,"value":"countries = {\n const geojson = await fetch(\"https://cdn.jsdelivr.net/gh/martynafford/natural-earth-geojson/50m/cultural/ne_50m_admin_0_countries_lakes.json\").then(r => r.json());\n for (let c of geojson.features) {\n delete c.properties;\n }\n let topology = topojson.topology({ countries: geojson });\n return topology; // topology.presimplify\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1139,"value":"[\n L.Map.prototype.options.crs.wrapLng,\n L.Map.prototype.options.crs.wrapLat\n]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1145,"value":"L.Map.prototype.options.crs.latLngToPoint(L.latLng(0,180), -8)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1149,"value":"L.Map.prototype.options.crs.zoom(1)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":528,"value":"// countries = await fetch(\"https://cdn.jsdelivr.net/npm/world-atlas@2/countries-10m.json\").then(r => r.json());","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":548,"value":"L.CRS.Earth","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":434,"value":"turf.featureCollection()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":512,"value":"sp = d3.scalePow()\n .exponent(2)\n .domain([0, 1, 3])\n .range([0.5, 0.5, 6]);","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":514,"value":"sp(3)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":731,"value":"map.getrend","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":397,"value":"rtLayer = undefined;\n/*\n{\n const scale = d3.scalePow()\n .exponent(2)\n .domain([0, 1, 5])\n .range([1, 1, 10]);\n \n return L.realtime(turf.featureCollection([]), {\n noClip: true,\n smoothFactor: 1,\n weight: 1,\n fnWeight(lineWidth, renderer) {\n return scale(renderer._zoom) * lineWidth;\n }\n });\n}\n*/","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":470,"value":"CanvasCustomStroke = {\n const parent = L.Canvas.prototype;\n \n return L.Canvas.extend({\n _fillStroke(ctx, layer) {\n \t\tconst options = layer.options;\n \n \t\tif (options.fill) {\n \t\t\tctx.globalAlpha = options.fillOpacity;\n \t\t\tctx.fillStyle = options.fillColor || options.color;\n \t\t\tctx.fill(options.fillRule || 'evenodd');\n \t\t}\n\n let lineWidth = options.weight;\n if (options.fnWeight) {\n lineWidth = options?.fnWeight(lineWidth, this) ?? lineWidth;\n }\n\n \t\tif (options.stroke && lineWidth !== 0) {\n \t\t\tif (ctx.setLineDash) {\n \t\t\t\tctx.setLineDash(layer.options && layer.options._dashArray || []);\n \t\t\t}\n \t\t\tctx.globalAlpha = options.opacity;\n \t\t\tctx.lineWidth = lineWidth;\n \t\t\tctx.strokeStyle = options.color;\n \t\t\tctx.lineCap = options.lineCap;\n \t\t\tctx.lineJoin = options.lineJoin;\n this._stroke(ctx, layer, lineWidth);\n \t\t}\n \t},\n _stroke(ctx, lineWidth) {\n ctx.stroke();\n }\n })\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":446,"value":"CanvasInnerStroke = {\n const parent = CanvasCustomStroke.prototype;\n \n return CanvasCustomStroke.extend({\n _stroke(ctx, layer) {\n // From https://github.com/Leaflet/Leaflet/issues/3665#issue-96911836\n ctx.save();\n ctx.clip();\n ctx.lineWidth *= 2;\n ctx.stroke();\n ctx.restore();\n }\n });\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":178,"value":"clippedSVG = new ClippedSVG()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":442,"value":"$4a","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":163,"value":"$4b","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":119,"value":"testFeature = {\n const featureCoords = geoJson.features.map(f => f.geometry)\n .map(g => g.type == \"Polygon\" ? g.coordinates : g.coordinates.flat())\n .map((fc,i) => ({ index: i, xCoords: fc.flat().map(c => c[0] == 0 ? Infinity : c[0]) }))\n const leftmost = _.minBy(featureCoords, fc => _.min(fc.xCoords));\n return geoJson.features[leftmost.index];\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":65,"value":"mapSize = ({ width: 2048, height: 904 });","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":26,"value":"wrappingCrs = {\n const simple = L.CRS.Simple;\n const wrappingCrs = L.extend({}, simple, {\n wrapLng: [0, mapSize.width],\n infinite: true,\n projection: {\n project: function(latLng) {\n const latLng2 = L.latLng(mapSize.height - latLng.lat, latLng.lng);\n return simple.projection.project.call(this, latLng2);\n },\n unproject: function(point) {\n const latLng = simple.projection.unproject.call(this, point);\n return L.latLng(mapSize.height - latLng.lat, latLng.lng);\n },\n bounds: {\n min: { x: 0, y: 0 },\n max: { x: mapSize.width, y: mapSize.height }\n }\n }\n })\n return wrappingCrs;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":317,"value":"topoJson","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":322,"value":"geoContiguous = {\n const layer = topoJson.objects.traced;\n const geoJson = topojson.feature(topoJson, layer);\n \n \n \n}\n//turf.coordAll(geoJson)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":331,"value":"geoJson","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":553,"value":"topoJson","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":563,"value":"contigTopology = \n return makeContiguous(st2, st2.objects.traced);\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":699,"value":"Bbox.toString()","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":19,"value":"---\n## Imports","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":28,"value":"import { smoothTopology, makeContiguous, sampleTopology as topoJson, topojson, Bbox, sampleGeoJson as geoJson } from \"8c9dbd3e4bfdb2e9\"","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":702,"value":"Flatbush = require('flatbush@4.2.0/flatbush.js');","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":874,"value":"topojson.feature(topoJson, topoJson.objects.traced)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":872,"value":"topoJson.objects","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":864,"value":"leafletBase = await require('leaflet@1.2.0');","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":878,"value":"topoJson.objects.traced.geometries[0]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":881,"value":"topojson.feature(topoJson, topoJson.objects.traced.geometries[0])","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":913,"value":"_.intersection(Object.keys(L.Polygon.prototype), Object.keys(L.Polyline.prototype))","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":898,"value":"// pa = new PolygonArc(topoJson, topoJson.objects.traced.geometries[0]);","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":959,"value":"topoJson.objects.traced","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":990,"value":"topoJson.objects.traced.geometries[0]","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":870,"value":"$4c","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1159,"value":"d3 = require(\"d3\", \"d3-geo\", \"d3-geo-projection\", \"d3-regression\");","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1162,"value":"simple = topojson.presimplify(countries);","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2349,"value":"L = require('leaflet@1.9.4/dist/leaflet-src.js')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":8,"value":"html``","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2797,"value":"{\n const allGeometries = Object.values(simple.objects).flatMap(o => o.geometries);\n return _.groupBy(allGeometries.map(g => g.arcs.flat().flat()).flat().map(id => (id < 0 ? ~id : id)));\n // .flatMap(o => o.arcs.flat().flat()))).filter(i => i.length > 1)\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2793,"value":"## Topology classes for Leaflet","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":918,"value":"polyPathOverrides = ({\n getTopology() {\n return this._topology;\n },\n setTopology(topology) {\n this._setTopology(topology);\n return this.redraw();\n },\n getGeometry() {\n return this._geometry;\n },\n setGeometry(geometry) {\n this._setGeometry({ type: geometry?.type, arcs: geometry?.arcs });\n return this.redraw(); \n },\n setLatLngs(latlngs) {\n throw \"Can't call setLatLngs on a topology path.\";\n },\n markDirty() {\n this._latlngs = false;\n this.redraw();\n },\n isEmpty() {\n return !this._geometry?.arcs?.length;\n },\n _setTopology(topology) {\n this._topology = topology;\n this.markDirty(); \n },\n _setGeometry(geometry) {\n this._geometry = geometry;\n this.markDirty(); \n },\n _simplifyPoints() {\n return;\n },\n _clipPoints() {\n this._parts = this._rings;\n },\n});","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":876,"value":"$4d","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":925,"value":"PolygonArc = L.Polygon.extend({\n ...polyPathOverrides,\n options: {\n smoothFactor: 0,\n noClip: true,\n fill: true\n },\n initialize() {\n PolylineArc.prototype.initialize.call(this, ...arguments);\n },\n /*\n _updatePath: function () {\n this.__latlngs = false;\n L.Polygon.prototype._updatePath.call(this);\n },\n */\n});","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2851,"value":"Math.sqrt(simple.arcs.length)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2839,"value":"{\n const allScores = simple.arcs.flat().map(c => c[2]).filter(c => c > 0 && isFinite(c));\n allScores.sort((a,b) => a-b);\n const index = Math.floor(allScores.length - simple.arcs.length / 3);\n return allScores[index];\n // return { asl: allScores.length, al: simple.arcs.length, index };\n return allScores.toReversed();\n}","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2874,"value":"allScores = simple.arcs.flat().map(c => c[2]).filter(c => c > 0 && isFinite(c));","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2876,"value":"stdev = d3.deviation(allScores)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2879,"value":"d3.mean(allScores) + 2 * stdev","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2870,"value":"lSimple._basePxArea","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":2868,"value":"lSimple = new LeafletTopology(simple, map.map)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":1816,"value":"$4e","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":84,"value":"turf = await require('@turf/turf@6.5.0/turf.min.js')","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":241,"value":"svg``","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","tool":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W4f","stdlib":"1"}} 24:null