Plot.plot({
projection: { type: "mercator", domain: hydrocountiesEPSG4326 },
width,
marks: [
Plot.geo(basinsEPSG4326, Plot.pointer(Plot.centroid({fill:"lightblue"}))),
Plot.geo(riversEPSG4326, {stroke: "blue"}),
Plot.geo(basinsEPSG4326, {stroke: "lightgray"}),
Plot.geo(geoGRDams, {stroke: "darkblue", fill: "lightblue"}),
Plot.geo(hydrocountiesEPSG4326, {stroke: "black"}),
Plot.density(geoGRDamsGeometries, {x: (d) => d.coordinates[0], y: (d) => d.coordinates[1], bandwidth: 10, thresholds: 10, fill: "density", fillOpacity: 0.25, mixBlendMode: "multiply"}),
Plot.tip(riversEPSG4326.features, Plot.pointer(Plot.centroid({title: (d) => `${d.properties.name},${d.properties.shape_leng}`})))
],
})