radialColour = (context, alpha=1) => {
let interpolateColScheme = d3.interpolateRgbBasis(d3.schemeBlues[9].slice(2))
let colour = d3.scaleDiverging(interpolateColScheme)
.domain([1, 15, 30])
let grd = context.createRadialGradient(0,0, r(colour.domain()[0]), 0, 0, r(colour.domain()[2]))
let stops = 50
for (let i = 0; i <= stops; i ++){
let stop = i/stops;
grd.addColorStop(stop, d3.color(interpolateColScheme(stop)).copy({opacity:alpha}))
}
return grd
}