31:Tacf, 32:T475,center = { const center = { p2006: [0, 0], p2011: [0, 0], p2018: [0, 0], p2021: [0, 0] }; let tot2006 = 0; let tot2011 = 0; let tot2018 = 0; let tot2021 = 0; for (let c of cells) { tot2006 += c.TOT_P_2006; tot2011 += c.TOT_P_2011; tot2018 += c.TOT_P_2018; tot2021 += c.TOT_P_2021; const pos = getPosition(c.GRD_ID); center.p2006[0] += c.TOT_P_2006 * pos.x; center.p2006[1] += c.TOT_P_2006 * pos.y; center.p2011[0] += c.TOT_P_2011 * pos.x; center.p2011[1] += c.TOT_P_2011 * pos.y; center.p2018[0] += c.TOT_P_2018 * pos.x; center.p2018[1] += c.TOT_P_2018 * pos.y; center.p2021[0] += c.TOT_P_2021 * pos.x; center.p2021[1] += c.TOT_P_2021 * pos.y; } center.p2006[0] /= tot2006; center.p2006[1] /= tot2006; center.p2011[0] /= tot2011; center.p2011[1] /= tot2011; center.p2018[0] /= tot2018; center.p2018[1] /= tot2018; center.p2021[0] /= tot2021; center.p2021[1] /= tot2021; //handle case when no 2006 data //if (isNaN(center.p2006[0]) || isNaN(center.p2006[1])) { //center.p2006[0] = -1000000; //center.p2006[1] = -1000000; //} return center; }33:[] 2f:{"title":"Population center movement","cells":[{"id":0,"value":"# Population center movement\n\nThis map shows the movement of the population center of various European countries. Source: [Eurostat](https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/grids)","pinned":false,"hidden":false,"mode":"md"},{"id":134,"value":"viewof geo = Inputs.select(cc, {\n value: \"AT\",\n label: \"Country / EU\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":96,"value":"$31","pinned":false,"hidden":false,"mode":"html"},{"id":188,"value":"## Annex","pinned":false,"hidden":false,"mode":"md"},{"id":7,"value":"Get population grids.","pinned":false,"hidden":false,"mode":"md"},{"id":5,"value":"import { grid5000 as db, getPosition } from \"@jgaffuri/eurpopgrid\"","pinned":true,"hidden":false,"mode":"ojs"},{"id":9,"value":"Select cells needed to compute the population centers.","pinned":false,"hidden":false,"mode":"md"},{"id":13,"value":"cells = db.query(\n \"SELECT GRD_ID, TOT_P_2006, TOT_P_2011, TOT_P_2018, TOT_P_2021 FROM data WHERE \" +\n (geo == \"EU\" ? euWhere : \"CNTR_ID='\" + geo + \"'\")\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":343,"value":"Compute the population centers.","pinned":false,"hidden":false,"mode":"md"},{"id":38,"value":"$32","pinned":true,"hidden":false,"mode":"ojs"},{"id":541,"value":"Another way to compute the total population.","pinned":false,"hidden":false,"mode":"md"},{"id":399,"value":"/*totals = db.query(\n \"SELECT sum(pop2006),sum(pop2011),sum(pop2018) FROM data WHERE \" +\n (geo == \"EU\" ? euWhere : \"CNTR_ID='\" + geo + \"'\")\n)*/","pinned":true,"hidden":false,"mode":"ojs"},{"id":90,"value":"## Map elements","pinned":false,"hidden":false,"mode":"md"},{"id":543,"value":"### Main map, with country","pinned":false,"hidden":false,"mode":"md"},{"id":313,"value":"w = width * 0.3","pinned":true,"hidden":false,"mode":"ojs"},{"id":108,"value":"h = 300","pinned":true,"hidden":false,"mode":"ojs"},{"id":111,"value":"projection = d3\n .geoIdentity()\n .reflectY(true)\n .fitExtent(\n [\n [30, 30],\n [w - 60, h - 60]\n ],\n geomsMap\n )","pinned":true,"hidden":false,"mode":"ojs"},{"id":106,"value":"path = d3.geoPath(projection)","pinned":true,"hidden":false,"mode":"ojs"},{"id":250,"value":"circleRadius = 4","pinned":true,"hidden":false,"mode":"ojs"},{"id":265,"value":"colors = {\n return { p2006: \"salmon\", p2011: \"salmon\", p2018: \"salmon\", p2021: \"salmon\" };\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":446,"value":"### Inset map, with trajectory","pinned":false,"hidden":false,"mode":"md"},{"id":439,"value":"wI = w * 0.8","pinned":true,"hidden":false,"mode":"ojs"},{"id":437,"value":"hI = h * 0.8","pinned":true,"hidden":false,"mode":"ojs"},{"id":346,"value":"Build the polyline linking the successive population centers.","pinned":false,"hidden":false,"mode":"md"},{"id":285,"value":"line = {\n const cs = [];\n\n if (center.p2006[0]) cs.push(center.p2006);\n if (center.p2011[0]) cs.push(center.p2011);\n if (center.p2018[0]) cs.push(center.p2018);\n if (center.p2021[0]) cs.push(center.p2021);\n\n return {\n type: \"LineString\",\n coordinates: cs\n };\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":225,"value":"projectionI = d3\n .geoIdentity()\n .reflectY(true)\n .fitExtent(\n [\n [40, 40],\n [wI - 80, hI - 80]\n ],\n line\n )","pinned":true,"hidden":false,"mode":"ojs"},{"id":223,"value":"pathI = d3.geoPath(projectionI)","pinned":true,"hidden":false,"mode":"ojs"},{"id":329,"value":"circleRadiusI = 8","pinned":true,"hidden":false,"mode":"ojs"},{"id":466,"value":"### Scale bar","pinned":false,"hidden":false,"mode":"md"},{"id":505,"value":"Compute data for a nice scale bar.","pinned":false,"hidden":false,"mode":"md"},{"id":459,"value":"scaleBar = {\n //const bn = d3.geoBounds(line);\n const p0 = projectionI.invert([0, 0]);\n const p1 = projectionI.invert([wI / 2, 0]);\n const d = p1[0] - p0[0];\n const d10 = Math.pow(10, Math.floor(Math.log10(d)));\n const d_ =\n d10 * 8 <= d\n ? d10 * 8\n : d10 * 6 <= d\n ? d10 * 6\n : d10 * 5 <= d\n ? d10 * 5\n : d10 * 3 <= d\n ? d10 * 3\n : d10 * 2 <= d\n ? d10 * 2\n : d10;\n return {\n size: ((wI / 2) * d_) / d,\n label: d_ >= 1000 ? d_ / 1000 + \" km\" : d_ + \" m\"\n };\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":452,"value":"### Geometries","pinned":false,"hidden":false,"mode":"md"},{"id":550,"value":"Get European Union countries.","pinned":false,"hidden":false,"mode":"md"},{"id":81,"value":"geomsTJ = await d3.json(\n //\"https://raw.githubusercontent.com/eurostat/Nuts2json/master/pub/v2/2021/3035/03M/0.json\"\n \"https://gisco-services.ec.europa.eu/distribution/v2/nuts/topojson/NUTS_RG_03M_2021_3035.json\"\n)","pinned":true,"hidden":false,"mode":"ojs"},{"id":84,"value":"geomsGJ = {\n let g = topojson.feature(geomsTJ, geomsTJ.objects[\"nutsrg\"]);\n g.features = g.features.filter((r) => r.properties.LEVL_CODE == 0);\n g.features = g.features.filter(\n (r) =>\n r.id != \"AL\" &&\n r.id != \"IS\" &&\n r.id != \"UK\" &&\n r.id != \"TR\" &&\n r.id != \"RS\" &&\n r.id != \"CH\" &&\n r.id != \"LI\" &&\n r.id != \"NO\" &&\n r.id != \"ME\" &&\n r.id != \"MK\"\n );\n return g;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":557,"value":"Map geometries.","pinned":false,"hidden":false,"mode":"md"},{"id":149,"value":"geomsMap = {\n return {\n type: \"FeatureCollection\",\n features: geomsGJ.features.filter((r) => (geo == \"EU\" ? true : r.id == geo))\n };\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":208,"value":"Prepare list of country codes + EU for the menu.","pinned":false,"hidden":false,"mode":"md"},{"id":132,"value":"cc = {\n const o = geomsGJ.features.map((r) => r.id);\n o.sort();\n o.push(\"EU\");\n return o;\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":389,"value":"Prepare SQL 'WHERE' clause, for EU case.","pinned":false,"hidden":false,"mode":"md"},{"id":393,"value":"euWhere = {\n const o = geomsGJ.features.map((r) => \"CNTR_ID='\" + r.id + \"'\");\n return o.join(\" OR \");\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":667,"value":"Country names from country codes","pinned":false,"hidden":false,"mode":"md"},{"id":670,"value":"cn = {\n const d = new Intl.DisplayNames([\"en\"], { type: \"region\" });\n const out = {};\n for (let c of cc) {\n out[c] = d.of(c);\n }\n\n out.EL = \"Greece\";\n return out; //d.of(\"FR\");\n}","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W33","stdlib":"1","titling":"auto"}