PREFIX xsd: #+ summary: Elite density in Amsterdam Neighbourhoods 1884. #+ tags: #+ - Amsterdam #+ - GIS #+ - Elite #+ endpoint: ... PREFIX rdf: PREFIX qb: PREFIX cedar: PREFIX sdmx-dimension: PREFIX sdmx-code: PREFIX geo: PREFIX rdfs: #select * # ?totPop ?elitePop ?districtWKTTooltip SELECT ?districtWKTLabel ?districtWKT ?districtWKTHeight ("0"^^ as ?districtWKTZ) (CONCAT('viridis,', STR(?eliteDensityIndex)) AS ?districtWKTColor) WHERE { # determine population per district { SELECT ?district (SUM(?pop) AS ?totPop) WHERE { ?obs a qb:Observation . ?obs cedar:population ?pop . ?obs sdmx-dimension:refArea ?district . ?slice a qb:Slice. ?slice qb:observation ?obs. ?slice sdmx-dimension:refPeriod "1889"^^xsd:int . } GROUP BY ?district } # determine elite per district { SELECT ?district (count(?obs) AS ?elitePop) WHERE { ?obs "1884" . ?obs ?district . } GROUP BY ?district } # get data about districts { SELECT ?district ?districtWKT ?districtLabel WHERE { ?district a . # ?district geo:hasGeometry / geo:asWKT ?districtWKT . # 2020-05-31 by @rlzijdeman: line above ought to be correct, # but data have no asWKT. Therefore replaced with line below: ?district geo:hasGeometry ?districtWKT . ?district rdfs:label ?districtLabel . } } BIND(1000*(5*(xsd:float(?elitePop))/(xsd:float(?totPop))) AS ?eliteDensity) BIND(CONCAT( "Neighbourhood: ", ?districtLabel, "\n", "Year: 1899", "\n", "Elite density: ", STR(?eliteDensity), "\n", "Source: ", "https://github.com/CLARIAH/ATM" ) AS ?districtWKTLabel) BIND((?eliteDensity / 50) AS ?eliteDensityIndex) BIND(?eliteDensity * 100 as ?districtWKTHeight) } ORDER BY ?districtLabel

Documentation

Triply