30:T1494,stampView={ const twopi=2*Math.PI; const size=540; const c=DOM.context2d(size,size,1); c.font="24px sansserif"; c.lineWidth=1; c.lineJoin="round"; c.capStyle="round"; c.imageSmoothingEnabled=false; c.strokeStyle ="#fff"; c.fillStyle="#fff"; return Object.assign(c.canvas, { update(pid,zoomLevel,drot,dx,dy,dr1,dr2) { const zoom=zoomLevel; const npix=size/zoom; const cal=idIndex.get(pid); const T0 = (cal.petal_loc - 3) * 36; const loc=1000*cal.petal_loc+cal.device_loc; const C=coords[loc]; // Theta axis in FVC coords. let [x0,y0]=C.FVC; x0 += centerdx; y0 += centerdy; // Theta axis in FP coords. const [xaxis,yaxis] = C.FP; const FP2FVC=([x,y])=>{ const [dx,dy] = [x-xaxis,y-yaxis]; return [x0 + C.dFVCdx[0]*dx + C.dFVCdy[0]*dy, y0 + C.dFVCdx[1]*dx + C.dFVCdy[1]*dy]; }; const [sx,sy]=[(x0-npix/2)/downsampling,(y0-npix/2)/downsampling] .map(Math.round) .map(xy=>xy*downsampling); const to=(x,y)=>[(x-sx)*zoom,(y-sy)*zoom]; // Arm lengths. const r1=cal.length_r1 + dr1; const r2=cal.length_r2 + dr2; // Draw the image with x and y flipped. c.save(); c.translate(size/2,size/2); c.scale(-1,-1); c.translate(-size/2,-size/2); c.clearRect(0,0,size,size); const args=[ sx/downsampling,sy/downsampling,npix/downsampling,npix/downsampling, 0,0,size,size]; const useColor=displayOptions.includes("color"); if(displayOptions.includes("front")) c.drawImage((useColor?colorizedImage:grayImage).src,...args); if(displayOptions.includes("back")) c.drawImage((useColor?tintedSpots:spotImage).src,...args); // Draw the path of the theta arm. const drawEllipse = function(x,y,r) { const [rx1,ry1] = C.dFVCdx.map(u=>r*u); const [rx2,ry2] = C.dFVCdy.map(u=>r*u); const [rX,rY] = [[rx1,ry1],[rx2,ry2]].map(([x,y])=>Math.sqrt(x**2+y**2)*zoom); const rot = Math.atan2(-ry1,rx1); const [xc,yc] = to(x,y); c.beginPath(), c.setLineDash([2,2]), c.ellipse(xc,yc,rX,rY,rot,0,twopi), c.stroke(); c.beginPath(), c.setLineDash([]), c.arc(xc,yc,4,0,twopi), c.stroke(); } if(displayOptions.includes("circles")) drawEllipse(x0,y0,r1); let rotdeg = (C?.T && C?.P) ? (T0 + C.T + C.P) : 0; rotdeg += drot; const rotrad = rotdeg*Math.PI/180; const [cos,sin] = [Math.cos(rotrad),Math.sin(rotrad)]; const drawPath = function(x0,y0,rot,pts) { rot *= Math.PI/180; const [cos,sin]=[Math.cos(rot),Math.sin(rot)]; const path=pts.map(([x,y])=>{ const [xr,yr] = [cos*x - sin*y, sin*x + cos*y]; return [x0 + C.dFVCdx[0]*xr + C.dFVCdy[0]*yr, y0 + C.dFVCdx[1]*xr + C.dFVCdy[1]*yr]; }); c.beginPath(); c.moveTo(...to(...path[0])); path.slice(1).forEach(([x,y])=>c.lineTo(...to(x,y))); c.stroke(); } let [T,P] = [0,0]; // fiber position in FP coords let [xfiber,yfiber] = (C?.FIBER_FP) ? C.FIBER_FP : [xaxis+2, yaxis]; xfiber += dx; yfiber += dy; // fiber position in FVC coords const [xf,yf] = FP2FVC([xfiber,yfiber]); if(displayOptions.includes("circles")) drawEllipse(xf,yf,r2); if(displayOptions.includes("phi arm")) { drawPath(xf,yf,rotdeg,ferruleShape.outer); drawPath(xf,yf,rotdeg,ferruleShape.trap); } // Elbow location in FP coords. const [xelbow,yelbow] = [xfiber-cos*r2, yfiber-sin*r2]; // Calculate external theta. T=Math.atan2(yelbow-yaxis,xelbow-xaxis)*180/Math.PI; P=Math.atan2(yfiber-yelbow,xfiber-xelbow)*180/Math.PI - T; T -= T0; while(T < 0) T+=360; while(P < 0) P+=360; // Draw visible parts of the theta arm if requested. if(displayOptions.includes("theta base")) drawPath(x0,y0,T+T0+180,upperHousingShape.lower); if(displayOptions.includes("theta upper")) drawPath(x0,y0,T+T0+180,upperHousingShape.upper); if(displayOptions.includes("phi stop")) drawPath(x0,y0,T+T0+180,upperHousingShape.stop); // Elbow location in FVC coords. const [xfvc,yfvc] = [ x0 + C.dFVCdx[0]*(xelbow-xaxis) + C.dFVCdy[0]*(yelbow-yaxis), y0 + C.dFVCdx[1]*(xelbow-xaxis) + C.dFVCdy[1]*(yelbow-yaxis)]; if(displayOptions.includes("lines")) { c.beginPath(); c.moveTo(...to(x0,y0)); c.lineTo(...to(xfvc,yfvc)); c.lineTo(...to(xf,yf)); c.save(); c.lineWidth=2; c.strokeStyle=tintColor; c.stroke(); c.restore(); } // Draw text labels c.restore(); c.fillStyle = "#fff"; c.fillText(`${pid} @ ${loc}`, 4, size-8); const S=snapshotIndex.get(pid); let [locT,locP] = [S?.locT,S?.locP]; while(locT<0) locT += 360; while(locP<0) locP += 360; c.fillStyle="#fff"; c.textAlign="right"; c.fillText(`DB: T=${Number(locT).toFixed(1)}° P=${Number(locP).toFixed(1)}°`, size-4, size-8); c.fillStyle=tintColor; c.fillText(`T=${Number(T).toFixed(1)}° P=${Number(P).toFixed(1)}°`, size-4, size-34); c.textAlign="left"; } }); }31:[{"id":"3adf2707a54c2a42e3a7d37ef276689709f2847b66d4bda40b5089729cd299a5d3f61de0a96d87b936677028f3b57533efc96cb59d66ebbe9783287cca3660e6","href":"https://static.observableusercontent.com/files/3adf2707a54c2a42e3a7d37ef276689709f2847b66d4bda40b5089729cd299a5d3f61de0a96d87b936677028f3b57533efc96cb59d66ebbe9783287cca3660e6","name":"calib-20210621-00095513.csv","type":"text/csv","lastModified":1624368605885.499,"size":1254379},{"id":"76acdd55decded1ff5d0e21d7b2080c3c1f1555180ffbaad527a620c783bfced4794442d457e64e61bf9481d0c75ca66bc35630617bb2e87ac1676588bb9c11f","href":"https://static.observableusercontent.com/files/76acdd55decded1ff5d0e21d7b2080c3c1f1555180ffbaad527a620c783bfced4794442d457e64e61bf9481d0c75ca66bc35630617bb2e87ac1676588bb9c11f","name":"back-95747.png","type":"image/png","lastModified":1624561171988.049,"size":9285509},{"id":"789454e2c50ab6a201799107798fbc8ac9d79a46cafbeb07e4b2eb44dd6a8a90e6e920ee26d5be4aa5ef2c917e6881e5987960e9d3ff8aff46a58d18120f215b","href":"https://static.observableusercontent.com/files/789454e2c50ab6a201799107798fbc8ac9d79a46cafbeb07e4b2eb44dd6a8a90e6e920ee26d5be4aa5ef2c917e6881e5987960e9d3ff8aff46a58d18120f215b","name":"front-95748-95757.png","type":"image/png","lastModified":1624561229119.779,"size":23472748},{"id":"b04cc98415550a3d1047a45d252b8043ab57aee1f2999955a2193f171c0e7bc141682f46211da370dab20ece80ebc8e9fa7b6046d788fd5dab933ff2d78c2fed","href":"https://static.observableusercontent.com/files/b04cc98415550a3d1047a45d252b8043ab57aee1f2999955a2193f171c0e7bc141682f46211da370dab20ece80ebc8e9fa7b6046d788fd5dab933ff2d78c2fed","name":"coords-9574.json","type":"application/json","lastModified":1624561261955.875,"size":1546854},{"id":"9be2876da1967312d1577de1d0d98fd76abc9e1dfd9bc6d2ae750ee2efcf72541953f651f9c177093d9c1d633cde7ea4efd898fd1f2fce818dbe55e7007d1124","href":"https://static.observableusercontent.com/files/9be2876da1967312d1577de1d0d98fd76abc9e1dfd9bc6d2ae750ee2efcf72541953f651f9c177093d9c1d633cde7ea4efd898fd1f2fce818dbe55e7007d1124","name":"snapshot-95747.csv","type":"text/csv","lastModified":1624562753659.017,"size":203702}] 2e:{"title":"Front-Illuminated Image Analysis","cells":[{"id":0,"value":"# Front-Illuminated Image Analysis","pinned":false,"hidden":false,"mode":"md"},{"id":926,"value":"{\n const C=idIndex.get(posID);\n return html`r1=${Number(C.length_r1+dr1).toFixed(2)} mm, r2=${Number(C.length_r2+dr2).toFixed(2)} mm`\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":24,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":883,"value":"DOM.download(await new Promise(resolve => stampView.toBlob(resolve)), `${posID}.png`, \"Save PNG\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":32,"value":"viewof posID = Inputs.text({\n label:\"Positioner ID\",\n placeholder: \"Mnnnnn\",\n width:\"10em\",\n value:\"M05273\",\n submit:\"OK\",\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":641,"value":"viewof angAdjust=Inputs.range([-180,180],{step:0.2,label:'angle adjust [deg]:'})","pinned":false,"hidden":false,"mode":"ojs"},{"id":581,"value":"viewof spotdx=Inputs.range([-10,10],{step:0.05,label:'spot dx [mm]:'});","pinned":false,"hidden":false,"mode":"ojs"},{"id":583,"value":"viewof spotdy=Inputs.range([-10,10],{step:0.05,label:'spot dy [mm]:'});","pinned":false,"hidden":false,"mode":"ojs"},{"id":964,"value":"viewof centerdx=Inputs.range([-10,10],{step:0.05,label:'center dx [mm]:'});","pinned":false,"hidden":false,"mode":"ojs"},{"id":966,"value":"viewof centerdy=Inputs.range([-10,10],{step:0.05,label:'center dy [mm]:'});","pinned":false,"hidden":false,"mode":"ojs"},{"id":932,"value":"viewof dr1=Inputs.range([-1,1],{step:0.01,label:'T arm dr1 [mm]:'});","pinned":false,"hidden":false,"mode":"ojs"},{"id":934,"value":"viewof dr2=Inputs.range([-1,1],{step:0.01,label:'P arm dr2 [mm]:'});","pinned":false,"hidden":false,"mode":"ojs"},{"id":508,"value":"viewof displayOptions={\n const all=[\"front\",\"back\",\"circles\",\"lines\",\"phi arm\",\"phi stop\",\"theta upper\",\"theta base\",\"color\"];\n return Inputs.checkbox(all,{label:\"Display:\", value:all});\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":192,"value":"viewof zoomLevel=Inputs.radio([1,2,3,4,5,6], {value:4, label:\"Zoom level:\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":316,"value":"viewof colorMap=Inputs.radio([\"Cividis\",\"Viridis\",\"Inferno\",\"Magma\",\"Plasma\",\"Turbo\"],\n {value:\"Cividis\",label:\"Color map:\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":550,"value":"viewof spotHue=Inputs.range([0,360],{step:5,value:45,label:\"Spot hue:\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":574,"value":"viewof spotWhite=Inputs.range([0,1],{step:0.1,value:0.5,label:\"Spot lightness:\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":683,"value":"Plot.plot({\n height:600,\n width:600,\n x:{ticks:false, domain:[-425,425]},\n y:{ticks:false, domain:[-425,425]},\n color:{domain:diffRange[whichDiff], scheme:'BrBG'},\n marks:[\n Plot.dot(locs, {\n x:loc=>locIndex.get(loc).offset_x_cs5, y:loc=>locIndex.get(loc).offset_y_cs5, r:3,\n r:loc=>loc==theLoc?10:(hasFiber.get(loc)==\"Y\"?3:1),\n stroke:\"#aaa\", strokeWidth:1, //loc=>hasFiber.get(loc)==\"Y\"?0:1,\n fill:loc=>!(diffs[loc]?.[whichDiff])?0:diffs[loc][whichDiff],\n title:loc=>`${locIndex.get(loc).pos_id} ${whichDiff}=${Number(diffs[loc][whichDiff]).toFixed(2)}`\n }),\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":696,"value":"viewof whichDiff=Inputs.radio([\"dr\",\"dT\",\"dP\"],{value:\"dr\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":778,"value":"Plot.plot({\n x:{domain:[-15,15]},\n //y:{domain:[-15,15]},\n grid:true,\n marks:[\n Plot.frame(),\n //Plot.dot(locs,{x:loc=>diffs[loc].dT, y:loc=>diffs[loc].dP, r:1})\n //Plot.dot(locs,{x:loc=>locIndex.get(loc).offset_t, y:loc=>diffs[loc].dT, r:1}),\n Plot.rectY(locs.filter(loc=>Math.abs(diffs[loc].dP)<15), Plot.binX({y: \"count\"}, {x: {value:loc=>diffs[loc].dP, thresholds:d3.range(-15,16,0.5)}, fill:'#f008'})),\n Plot.rectY(locs.filter(loc=>Math.abs(diffs[loc].dT)<15), Plot.binX({y: \"count\"}, {x: {value:loc=>diffs[loc].dT, thresholds:d3.range(-15,16,0.5)}, fill:'#00f8'})),\n Plot.ruleX([d3.median(locs, loc=>diffs[loc].dP), d3.median(locs, loc=>diffs[loc].dT)])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":895,"value":"wrap=A=>{\n while(A<0) A+= 360;\n while(A>=360) A-= 360;\n return Number(Number(A).toFixed(3));\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":890,"value":"ids=[\"M01686\",\"M01719\",\"M01865\",\"M02060\",\"M02065\",\"M02351\",\"M02433\",\"M02480\",\"M02547\",\"M02809\",\"M02905\",\"M03415\",\"M03454\",\"M03456\",\"M03601\",\"M03675\",\"M03802\",\"M03809\",\"M03972\",\"M04060\",\"M04138\",\"M04235\",\"M04316\",\"M04588\",\"M04596\",\"M04744\",\"M04763\",\"M04792\",\"M05018\",\"M05300\",\"M05439\",\"M05472\",\"M05478\",\"M06489\",\"M07119\",\"M07369\",\"M07373\",\"M07416\",\"M07480\",\"M07500\",\"M07763\",\"M06783\",\"M03120\",\"M03294\",\"M03757\",\"M03966\",\"M04243\",\"M05152\",\"M03599\",\"M03596\",\"M03493\",\"M07163\",\"M05638\",\"M06847\",\"M04217\",\"M03837\",\"M03554\",\"M02582\",\"M04556\",\"M05038\",\"M02956\",\"M04529\",\"M04133\",\"M03156\",\"M04283\",\"M05607\",\"M03990\",\"M04478\",\"M02792\",\"M01190\",\"M01727\",\"M06214\",\"M06009\",\"M03767\",\"M03159\",\"M03943\",\"M02732\",\"M04583\",\"M05070\",\"M07601\"]","pinned":false,"hidden":false,"mode":"ojs"},{"id":808,"value":"diffMu=[d3.median(locs.filter(loc=>Math.abs(diffs[loc].dP)<15), loc=>diffs[loc].dP),\n d3.median(locs.filter(loc=>Math.abs(diffs[loc].dT)<15), loc=>diffs[loc].dT)]","pinned":false,"hidden":false,"mode":"ojs"},{"id":817,"value":"diffRms=[d3.variance(locs.filter(loc=>Math.abs(diffs[loc].dP)<15), loc=>diffs[loc].dP),\n d3.variance(locs.filter(loc=>Math.abs(diffs[loc].dT)<15), loc=>diffs[loc].dT)].map(Math.sqrt)","pinned":false,"hidden":false,"mode":"ojs"},{"id":762,"value":"theLoc = {\n const C=idIndex.get(posID);\n return C.device_loc+1000*C.petal_loc;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":700,"value":"diffRange={\n const nsig=4;\n const [muP,winP] = [diffMu[0], nsig*diffRms[0]];\n const [muT,winT] = [diffMu[1], nsig*diffRms[1]];\n return {\n dr:[-0.05,0.05],\n dT:[muT-winT,muT+winT],\n dP:[muP-winP,muP+winP],\n };\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":38,"value":"live = {\n stampView.update(posID,zoomLevel,angAdjust,spotdx,spotdy,dr1,dr2);\n return true;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":295,"value":"interpolator = d3[`interpolate${colorMap}`]","pinned":false,"hidden":false,"mode":"ojs"},{"id":862,"value":"newLUT = {\n const c = DOM.context2d(256, 1, 1);\n for(let i = 0; i < 256; i++) c.fillStyle = interpolator(i / 255), c.fillRect(i, 0, 1, 1);\n return new Uint32Array(c.getImageData(0,0,256,1).data.buffer);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":346,"value":"colorizedImage={\n const src=grayImage.src;\n const c=DOM.context2d(src.width,src.height,1);\n c.drawImage(src,0,0);\n const imgData = c.getImageData(0,0,src.width,src.height);\n /*\n for(let i=0; i[loc,coords[loc]?.FIBER_FVC?'Y':'N']))","pinned":false,"hidden":false,"mode":"ojs"},{"id":659,"value":"diffs = Object.fromEntries(locs.map(loc => {\n const C=coords[loc];\n const S=snapshotIndex.get(C.POSID);\n const result={};\n const wrap = dA=>{\n while(dA>180) dA-=360;\n while(dA<=-180) dA+=360;\n return dA;\n }\n if(C?.T!=null && S?.locT!=null) result['dT'] = wrap(C.T - S.locT);\n if(C?.P!=null && S?.locP!=null) result['dP'] = wrap(C.P - S.locP);\n if(C?.FIBER_FP!=null && S?.fpX!=null && S?.fpY!=null) {\n const [x,y] = C.FIBER_FP;\n result['dr'] = Math.sqrt((x-S.fpX)**2 + (y-S.fpY)**2);\n }\n return [loc, result];\n}));","pinned":false,"hidden":false,"mode":"ojs"},{"id":8,"value":"coords = FileAttachment(\"coords-9574.json\").json()","pinned":false,"hidden":false,"mode":"ojs"},{"id":134,"value":"locs = Array.from(locIndex.keys())","pinned":false,"hidden":false,"mode":"ojs"},{"id":230,"value":"idIndex = d3.index(calib, d=>d.pos_id)","pinned":false,"hidden":false,"mode":"ojs"},{"id":222,"value":"locIndex = d3.index(calib, d=>1000*d.petal_loc+d.device_loc)","pinned":false,"hidden":false,"mode":"ojs"},{"id":216,"value":"calib = d3.csvParse(await FileAttachment(\"calib-20210621-00095513.csv\").text(), d3.autoType)","pinned":false,"hidden":false,"mode":"ojs"},{"id":647,"value":"snapshotIndex = d3.index(snapshot, d=>d.pos_id)","pinned":false,"hidden":false,"mode":"ojs"},{"id":634,"value":"snapshot = d3.csvParse(await FileAttachment(\"snapshot-95747.csv\").text(), d3.autoType)","pinned":false,"hidden":false,"mode":"ojs"},{"id":395,"value":"import {upperHousingShape,ferruleShape} from \"@dkirkby/desi-pos-geom\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":297,"value":"d3 = require(\"d3-dsv@3\", \"d3-array@2\", \"d3-color@2\", \"d3-scale-chromatic@2\")","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W31","stdlib":"1","titling":"auto"}