30:Tcf6,canvasSkyMap = function(shapes, options={}) { options = Object.assign(defaultOptions, options); const {height, width, projection, visible} = initMap(options); const context = DOM.context2d(width, height); const path = d3.geoPath(projection, context); path.pointRadius(options.pointRadius); function render() { context.save(); context.fillStyle = "#fff", context.fillRect(0, 0, width, height); if(options["Grid"] && visible["Grid"]) { context.beginPath(), path(grid), context.strokeStyle = options["Grid"].stroke, context.stroke(); } context.save(); if(options["Ecliptic"] && visible["Ecliptic"]) { context.beginPath(), path(Ecliptic), context.strokeStyle = options["Ecliptic"].stroke, context.setLineDash([4, 4]), context.stroke(); } context.restore(); if(options["MilkyWay"] && visible["MilkyWay"]) { context.beginPath(), path(MilkyWay), context.strokeStyle = options["MilkyWay"].stroke, context.stroke(); } if(options["DESI_N"] && visible["DESI_N"]) { context.beginPath(), path(DESInorth), context.fillStyle = options["DESI_N"].fill, context.globalAlpha = options["DESI_N"].fillOpacity, context.fill(), context.strokeStyle = options["DESI_N"].stroke, context.globalAlpha = options["DESI_N"].strokeOpacity, context.stroke(); } if(options["DESI_S"] && visible["DESI_S"]) { context.beginPath(), path(DESIsouth), context.fillStyle = options["DESI_S"].fill, context.globalAlpha = options["DESI_S"].fillOpacity, context.fill(), context.strokeStyle = options["DESI_S"].stroke, context.globalAlpha = options["DESI_S"].strokeOpacity, context.stroke(); } if(options["EBOSS_N"] && visible["EBOSS_N"]) { context.beginPath(), path(EBOSSnorth), context.fillStyle = options["EBOSS_N"].fill, context.globalAlpha = options["EBOSS_N"].fillOpacity, context.fill(), context.strokeStyle = options["EBOSS_N"].stroke, context.globalAlpha = options["EBOSS_N"].strokeOpacity, context.stroke(); } if(options["EBOSS_S"] && visible["EBOSS_S"]) { context.beginPath(), path(EBOSSsouth), context.fillStyle = options["EBOSS_S"].fill, context.globalAlpha = options["EBOSS_S"].fillOpacity, context.fill(), context.strokeStyle = options["EBOSS_S"].stroke, context.globalAlpha = options["EBOSS_S"].strokeOpacity, context.stroke(); } context.globalAlpha = 1; if(options["Shapes"] && visible["Shapes"]) { for(const {shape,stroke,fill,mode} of shapes) { context.beginPath(), path(shape); //if(mode) context.mixBlendMode= mode; if(stroke) context.strokeStyle = stroke, context.stroke(); if(fill) context.fillStyle= fill, context.fill(); } } context.restore(); context.beginPath(), path(outline), context.strokeStyle = "#000", context.stroke(); } const zoomdrag = (options.zoomdrag == "auto") ? (options.projectionName == "Orthographic") : options.zoomdrag; if(zoomdrag) { return d3.select(context.canvas) .call(zoom(projection) .on("zoom.render", () => render()) .on("end.render", () => render())) .call(() => render()) .node(); } else { render(); return context.canvas; } }31:Tbe6,svgSkyMap = function(shapes, options={}) { options = Object.assign(defaultOptions, options); const {height, width, projection, visible} = initMap(options); const path = d3.geoPath(projection); path.pointRadius(options.pointRadius); // ${options["Shapes"] && visible["Shapes"] ? const shapesSVG = shapes.map( ({shape,stroke,fill,mode}) => svg``); return html``; }32:T4b4,function initMap({width, height, visibleFeatures, projectionName, centerRA, centerDEC} = {}) { if(!height && !width) throw new Error("Must specify either height or width."); if(height && width) throw new Error("Cannot specify both height and width."); const visible = Object.fromEntries( mapFeatures.map(feature => [feature, visibleFeatures.includes(feature)])); // Lookup the projection by name, first in d3-geo-projection, then d3-geo (lifted into d3). const symbol = "geo" + projectionName; const projection = (d3gp?.[symbol] ?? d3[symbol])(); projection.rotate([centerRA, centerDEC]); if(height) { const [[x0, y0], [x1, y1]] = d3.geoPath(projection.fitHeight(height, outline)).bounds(outline); width = Math.ceil(x1 - x0); const l = Math.min(Math.ceil(y1 - y0), width); projection.scale(projection.scale() * (l - 1) / l).precision(0.2); } else { const [[x0, y0], [x1, y1]] = d3.geoPath(projection.fitWidth(width, outline)).bounds(outline); height = Math.ceil(y1 - y0); const l = Math.min(Math.ceil(x1 - x0), height); projection.scale(projection.scale() * (l - 1) / l).precision(0.2); } return {height, width, projection, visible}; }33:T4dc,circlePath = function(RA0, DEC0, size, {nphi=36, dtype=Float32Array}={}) { const deg2rad = Math.PI/180; const phi = dtype.from({length:nphi+1}, (_,i) => -i*2*Math.PI/nphi); const cosPhi = phi.map(Math.cos), sinPhi = phi.map(Math.sin); const angles = dtype.from([RA0, DEC0, size], a => a * deg2rad); const sinAngles = angles.map(Math.sin), cosAngles = angles.map(Math.cos); // Calculate the 3D coordinates (x,y,z) of unit vectors along the circle. const cosAngles12 = cosAngles[1]*cosAngles[2]; const sinAngles12 = sinAngles[1]*sinAngles[2]; const A = cosAngles[0]*cosAngles12; const B = cosAngles[0]*sinAngles12; const C = sinAngles[0]*sinAngles[2]; const D = sinAngles[0]*cosAngles12; const E = cosAngles[0]*sinAngles[2]; const F = sinAngles[0]*sinAngles12; const G = sinAngles[1]*cosAngles[2]; const H = cosAngles[1]*sinAngles[2]; const coordinates = Array.from({length:nphi}, (_,i) => dtype.from([ // RA = atan2(y, x) Math.atan2( D + (E*cosPhi[i] - F*sinPhi[i]), // y[i] A - (B*sinPhi[i] + C*cosPhi[i]) // x[i] )/deg2rad, // DEC = asin(z) Math.asin( G + H*sinPhi[i] // z[i] )/deg2rad ])) return {type: "Polygon", coordinates: [coordinates]}; }34:[{"id":"2f59c749853cd6689f4d7b0889f9cbc591da80e4a867e11b9f7b9adf208f532e811c9fcda57ff1e6675be16ffd7c7df47574acc25c477358deb49000a4a78876","href":"https://static.observableusercontent.com/files/2f59c749853cd6689f4d7b0889f9cbc591da80e4a867e11b9f7b9adf208f532e811c9fcda57ff1e6675be16ffd7c7df47574acc25c477358deb49000a4a78876","name":"eboss-south.json","type":"application/json","lastModified":1708715609865.944,"size":1627},{"id":"be33db4056ff168fe696e4d952b0fbe603f486aabd71adf0911d487d3599ee255231bf8f6c7b9bcfcf5cc4d693d38aa110659e2c0f81b9f5fb62b54ae1033c67","href":"https://static.observableusercontent.com/files/be33db4056ff168fe696e4d952b0fbe603f486aabd71adf0911d487d3599ee255231bf8f6c7b9bcfcf5cc4d693d38aa110659e2c0f81b9f5fb62b54ae1033c67","name":"eboss-north@1.json","type":"application/json","lastModified":1708716231033.02,"size":2228},{"id":"8c6f8af215a4048754bb96aa4f71c81b7d5c94bcd77fdf6787ca879e22217a84bb6da3f8ec5cf7873e1dfacccd136a4455bd1531e5a6f497bb491710e11807de","href":"https://static.observableusercontent.com/files/8c6f8af215a4048754bb96aa4f71c81b7d5c94bcd77fdf6787ca879e22217a84bb6da3f8ec5cf7873e1dfacccd136a4455bd1531e5a6f497bb491710e11807de","name":"Ecliptic@1.json","type":"application/json","lastModified":1621555285964.44,"size":3927},{"id":"aa99e5b71bd9c1387b279027586bdcc57443fcdb5723072e24115d1e66b3d47c3d0cc1019726d050bb7148703ca43e5c149e00cd29ae7d3e318bae1ad51fe020","href":"https://static.observableusercontent.com/files/aa99e5b71bd9c1387b279027586bdcc57443fcdb5723072e24115d1e66b3d47c3d0cc1019726d050bb7148703ca43e5c149e00cd29ae7d3e318bae1ad51fe020","name":"desi-north.json","type":"application/json","lastModified":1621960294380.624,"size":2166},{"id":"07220c735761a79ae175aabde11ddac07244fa642df291ffa726e5cf999b91f771de500a765c0da310e8553987909e5c8347110f5b6c1d50d5d414d2d1d36960","href":"https://static.observableusercontent.com/files/07220c735761a79ae175aabde11ddac07244fa642df291ffa726e5cf999b91f771de500a765c0da310e8553987909e5c8347110f5b6c1d50d5d414d2d1d36960","name":"desi-south.json","type":"application/json","lastModified":1621960147079.365,"size":1886},{"id":"e2fdc0e5fdb3a5dca3537497711241dd726eeb8d020981c7a0b05a5657c286ae36f8fc906f46a48b8aa6372a48732a4c1402974fb680081a8b6b92260ce73c7b","href":"https://static.observableusercontent.com/files/e2fdc0e5fdb3a5dca3537497711241dd726eeb8d020981c7a0b05a5657c286ae36f8fc906f46a48b8aa6372a48732a4c1402974fb680081a8b6b92260ce73c7b","name":"MilkyWay@2.json","type":"application/json","lastModified":1621553470728.542,"size":3906}] 2e:{"title":"Sky Maps","cells":[{"id":0,"value":"# Sky Maps","pinned":false,"hidden":false,"mode":"md"},{"id":299,"value":"viewof ra0 = Inputs.range([-180,180], {value:0, step:1, label:\"ra0\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":301,"value":"viewof dec0 = Inputs.range([-90,90], {value:0, step:1, label:\"dec0\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":305,"value":"viewof radius = Inputs.range([1,120], {value:30, step:1, label:\"radius\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":419,"value":"canvasSky = canvasSkyMap(skyShapes, {centerRA, centerDEC, projectionName, visibleFeatures})","pinned":false,"hidden":false,"mode":"ojs"},{"id":852,"value":"viewof projectionName = Inputs.select(projectionNames, {label: \"Projection\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":157,"value":"viewof centerRA = Inputs.range([-180,180], {value:-125, step:1, label:\"Center RA [deg]\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":990,"value":"viewof centerDEC = Inputs.range([-80,80], {value:0, step:1, label:\"Center DEC [deg]\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":234,"value":"viewof mapSize = Inputs.range([10,1990], {value:500, step:1, label:\"Map Size [pix]\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":251,"value":"viewof visibleFeatures = {\n return Inputs.checkbox(mapFeatures, {value: mapFeatures, label: \"Map Features\"});\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":413,"value":"svgSky = svgSkyMap(skyShapes, {centerRA, centerDEC, projectionName, visibleFeatures})","pinned":false,"hidden":false,"mode":"ojs"},{"id":25,"value":"md`----`","pinned":false,"hidden":false,"mode":"ojs"},{"id":456,"value":"defaultOptions = ({\n visibleFeatures: mapFeatures,\n Grid: {stroke:\"#ddd\"},\n GridLabels: {fill:\"#999\"},\n Ecliptic: {stroke:\"#999\"},\n MilkyWay: {stroke:\"#555\"},\n DESI_N: {stroke:\"#dc143c\", strokeOpacity:0.4, fill:\"#dc143c\", fillOpacity:0.1},\n DESI_S: {stroke:\"#dc143c\", strokeOpacity:0.4, fill:\"#dc143c\", fillOpacity:0.1},\n EBOSS_N: {stroke:\"#4682B4\", strokeOpacity:0.4, fill:\"#4682B4\", fillOpacity:0.2},\n EBOSS_S: {stroke:\"#4682B4\", strokeOpacity:0.4, fill:\"#4682B4\", fillOpacity:0.2},\n Shapes: { },\n pointRadius: 2.5,\n projectionName: \"Eckert IV\",\n centerRA: -125,\n width: 800,\n height: null,\n render: \"canvas\",\n zoomdrag: \"auto\", // true / false / auto = (projectionName==\"Orthographic\")\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":384,"value":"skyShapes = [\n {shape:circlePath(ra0, dec0, radius), stroke:\"#000\", fill:\"#0f06\", mode:\"multiply\"},\n {shape:{type: \"Point\", coordinates: [ra0,dec0]}, fill:\"#00f4\"}\n]","pinned":false,"hidden":false,"mode":"ojs"},{"id":900,"value":"function getColorScale(colorMap) {\n const capitalized = colorMap[0].toUpperCase() + colorMap.slice(1).toLowerCase();\n const interpolate = \"interpolate\" + capitalized;\n const scheme = \"scheme\" + capitalized;\n if(interpolate in d3) {\n return d3.scaleSequential(d3[interpolate]);\n }\n else if(scheme in d3) {\n return d3.scaleOrdinal(d3[scheme]);\n }\n else {\n throw new Error(`Invalid color map name: \"${colorMap}\"`);\n }\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":903,"value":"function skyPoints(data, {\n RA=\"RA\", DEC=\"DEC\", value=d => 1,\n colorMap=\"plasma\",\n clipFraction = 0.025,\n}) {\n const z = data.map(value);\n const zlo = d3.quantile(z, clipFraction);\n const zhi = d3.quantile(z, 1 - clipFraction);\n const scale = getColorScale(colorMap).domain([zlo, zhi]).clamp(true);\n return data.map(d => ({\n shape:{ type: \"Point\", coordinates: [d[RA], d[DEC]] },\n fill: d3.color(scale(value(d))).hex(), // SVG needs hex\n }));\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":684,"value":"dataRangeOriginal = function(data, {vmin=null, vmax=null, pmin=0, pmax=100, niceCount=0, symmetric=false}={}) {\n let vlo = (vmin!=null) ? vmin : d3.quantile(data, 100*pmin);\n let vhi = (vmax!=null) ? vmax : d3.quantile(data, 100*pmax);\n if(symmetric) {\n const vlim = Math.max(Math.abs(vlo), Math.abs(vhi));\n [vlo, vhi] = [-vlim,+vlim];\n }\n if(niceCount) [vlo, vhi] = d3.nice(vlo, vhi, niceCount);\n return [vlo, vhi];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":888,"value":"TODO: implement GridLabels option for canvas map","pinned":false,"hidden":false,"mode":"md"},{"id":891,"value":"function skyMap(shapes, options={}) {\n options = Object.assign(defaultOptions, options);\n return options.render === \"svg\" ? svgSkyMap(shapes, options) : canvasSkyMap(shapes, options);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":28,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":31,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":551,"value":"gridLabels = {\n const raLabelDec = -(45+60)/2;\n const raLabels = d3.range(0, 360, 60).map(ra=>({label:ra+\"°\", radec:[ra,raLabelDec]}));\n const decLabelRa = 80;\n const decLabels = d3.range(-30, 90, 30).map(dec=>({label:dec+\"°\", radec:[decLabelRa,dec]}));\n return raLabels.concat(decLabels);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":803,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":866,"value":"mapFeatures = [\"Grid\", \"GridLabels\", \"MilkyWay\", \"Ecliptic\", \"DESI_N\", \"DESI_S\", \"EBOSS_N\", \"EBOSS_S\", \"Shapes\"]","pinned":false,"hidden":false,"mode":"ojs"},{"id":858,"value":"projectionNames = [\"Eckert4\", \"Aitoff\", \"Hammer\", \"Mollweide\", \"CylindricalEqualArea\", \"Healpix\", \"Orthographic\"]","pinned":false,"hidden":false,"mode":"ojs"},{"id":274,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":192,"value":"DESInorth = await FileAttachment(\"desi-north.json\").json().then(F => createPath(F, {type:\"Polygon\"}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":185,"value":"DESIsouth = await FileAttachment(\"desi-south.json\").json().then(F => createPath(F, {type:\"Polygon\"}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":942,"value":"EBOSSnorth = await FileAttachment(\"eboss-north@1.json\").json().then(F => createPath(F, {type:\"Polygon\"}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":939,"value":"EBOSSsouth = await FileAttachment(\"eboss-south.json\").json().then(F => createPath(F, {type:\"Polygon\"}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":77,"value":"Ecliptic = await FileAttachment(\"Ecliptic@1.json\").json().then(F => createPath(F))","pinned":false,"hidden":false,"mode":"ojs"},{"id":108,"value":"MilkyWay = await FileAttachment(\"MilkyWay@2.json\").json().then(F => createPath(F))","pinned":false,"hidden":false,"mode":"ojs"},{"id":79,"value":"createPath = function(data, {type=\"LineString\"}={}) {\n const RA = data[\"dRA\"].reduce((all,delta) => (all.push(all[all.length-1]+delta/10), all), [data[\"RA0\"]]);\n const DEC = data[\"dDEC\"].reduce((all,delta) => (all.push(all[all.length-1]+delta/10), all), [data[\"DEC0\"]]);\n RA.push(data[\"RA0\"]);\n DEC.push(data[\"DEC0\"]);\n let coords = d3.transpose([RA.reverse(),DEC.reverse()]);\n if(type==\"Polygon\") coords = [coords];\n return {type:type, coordinates: coords};\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":17,"value":"grid = {\n const graticule = d3.geoGraticule();\n graticule.stepMajor([360, 360]);\n graticule.stepMinor([30, 15]);\n graticule.extentMinor([[-180,-75.001],[180,75.001]]);\n return graticule();\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":14,"value":"outline = ({type: \"Sphere\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":7,"value":"md`---`","pinned":false,"hidden":false,"mode":"ojs"},{"id":913,"value":"import {zoom} from \"@d3/versor-zooming\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":896,"value":"d3gp = require(\"d3-geo-projection@4\")","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W34","stdlib":"1","titling":"auto"}