31:T4c2,// Initialize plot plot = Plotly.newPlot( plotDiv, [ { x: df.Date, y: detrended(d3.range(df.Date.length), df.Series0), type: "scatter", mode: markerstyle, marker: { size: "5" }, name: "CH0" }, { x: df.Date, y: detrended(d3.range(df.Date.length), df.Series1), mode: markerstyle, marker: { size: "5" }, name: "CH1" }, { x: df.Date, y: detrended(d3.range(df.Date.length), df.Series2), mode: markerstyle, marker: { size: "5" }, name: "CH2" }, { x: df.Date, y: detrended(d3.range(df.Date.length), df.Series3), mode: markerstyle, marker: { size: "5" }, name: "CH3" }, { x: [], y: [], mode: "markers", marker: { color: "red", size: 7 }, name: "Selected Points" } ], { title: "Time Series Plot with Interactive Selection", xaxis: { title: "Date" }, yaxis: { title: "Counts" }, showlegend: false, xaxis: { rangeslider: { visible: true, type: "date" } }, hovermode: "closest", // dragmode: "select", updatemenus: [ { type: "buttons", showactive: true } ] } )32:[] 2f:{"title":"BSM Time Series Tool","cells":[{"id":0,"value":"# BSM Time Series Tool ","pinned":false,"hidden":false,"mode":"md"},{"id":1013,"value":"viewof network = Inputs.select(['PB','IV','GF'],{label: 'Network', value: 'PB'})","pinned":false,"hidden":false,"mode":"ojs"},{"id":657,"value":"viewof station = Inputs.select(stationNames, {\n label: \"Station\",\n value: \"B026\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1008,"value":"viewof start = Inputs.text({ label: \"Start\", value: \"2023-03-01T00:00:00Z\" }) //viewof start = Inputs.datetime({\n// label: \"Start\",\n// value: \"2023-11-10T00:00:00Z\"\n//}) // note there is an issue with the datetime picker, where the display defaults to local time and I haven't figured out how to fix it","pinned":false,"hidden":false,"mode":"ojs"},{"id":1018,"value":"viewof end = Inputs.text({ label: \"End\", value: \"2023-04-01T00:00:00Z\" }) //viewof end = Inputs.datetime({ label: \"End\", value: \"2023-11-11T00:00:00Z\" })","pinned":false,"hidden":false,"mode":"ojs"},{"id":1021,"value":"viewof rate = Inputs.select(['LS','RS'],{ label: \"Sample Rate\", value: \"RS\" })","pinned":false,"hidden":false,"mode":"ojs"},{"id":1576,"value":"viewof markerstyle = Inputs.select([\"lines\", \"markers\", \"markers+lines\"], {\n label: \"Plot Style\",\n value: \"lines\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":992,"value":"show = html`${plotDiv}${textArea}${saveButton}${clearButton}`","pinned":true,"hidden":false,"mode":"ojs"},{"id":988,"value":"## Data Import and Preparation","pinned":false,"hidden":false,"mode":"md"},{"id":1591,"value":"stationNames = ['AVN2', 'B001', 'B003', 'B004', 'B005', 'B006', 'B007', 'B009',\n 'B010', 'B011', 'B012', 'B013', 'B014', 'B017', 'B018', 'B019',\n 'B020', 'B022', 'B023', 'B024', 'B026', 'B027', 'B028', 'B030',\n 'B031', 'B032', 'B033', 'B035', 'B036', 'B039', 'B040', 'B045',\n 'B054', 'B057', 'B058', 'B065', 'B066', 'B067', 'B072', 'B073',\n 'B075', 'B076', 'B078', 'B079', 'B081', 'B082', 'B084', 'B086',\n 'B087', 'B088', 'B089', 'B093', 'B201', 'B202', 'B203', 'B204',\n 'B205', 'B206', 'B207', 'B208', 'B900', 'B901', 'B916', 'B917',\n 'B918', 'B921', 'B926', 'B927', 'B928', 'B933', 'B934', 'B935',\n 'B941', 'B943', 'B944', 'B946', 'B950', 'BOZ1', 'BUY1', 'ESN1',\n 'HALK', 'TEPE', 'SIV1', 'TSM1', 'TSM2', 'TSM3', 'TSM4', 'TSM5',\n 'TSM6']","pinned":false,"hidden":false,"mode":"ojs"},{"id":503,"value":"function parseText(gtxt) {\n \n return d3.csvParse(\n gtxt\n .toString()\n .split(\"\\n\")\n .filter((row) => !row.includes(\"#\"))\n .join(\"\\n\")\n)\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":546,"value":"g1txt = d3.text(\n \"https://service.iris.edu/fdsnws/dataselect/1/query?net=\" +\n network +\n \"&sta=\" +\n station +\n \"&loc=T0&cha=\" +\n rate +\n \"1&starttime=\" +\n start +\n \"&endtime=\" +\n end +\n \"&format=geocsv&nodata=404\"\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":434,"value":"g2txt = d3.text(\n \"https://service.iris.edu/fdsnws/dataselect/1/query?net=\" +\n network +\n \"&sta=\" +\n station +\n \"&loc=T0&cha=\" +\n rate +\n \"2&starttime=\" +\n start +\n \"&endtime=\" +\n end+\n \"&format=geocsv&nodata=404\"\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":435,"value":"g3txt = d3.text(\n \"https://service.iris.edu/fdsnws/dataselect/1/query?net=\" +\n network +\n \"&sta=\" +\n station +\n \"&loc=T0&cha=\" +\n rate +\n \"3&starttime=\" +\n start +\n \"&endtime=\" +\n end +\n \"&format=geocsv&nodata=404\"\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":437,"value":"g4txt = d3.text(\n \"https://service.iris.edu/fdsnws/dataselect/1/query?net=\" +\n network +\n \"&sta=\" +\n station +\n \"&loc=T0&cha=\" +\n rate +\n \"4&starttime=\" +\n start+\n \"&endtime=\" +\n end+\n \"&format=geocsv&nodata=404\"\n)","pinned":false,"hidden":false,"mode":"ojs"},{"id":576,"value":"g0 = parseText(g1txt)","pinned":false,"hidden":false,"mode":"ojs"},{"id":574,"value":"g1 = parseText(g2txt)","pinned":false,"hidden":false,"mode":"ojs"},{"id":572,"value":"g2 = parseText(g3txt)","pinned":false,"hidden":false,"mode":"ojs"},{"id":559,"value":"g3 = parseText(g4txt)","pinned":false,"hidden":false,"mode":"ojs"},{"id":13,"value":"date_rng = g1.map((d) => d[\"Time\"])","pinned":false,"hidden":false,"mode":"ojs"},{"id":1464,"value":"data0_notZerod = linearInterpolate(g0.map((d) => parseFloat(d[\" Sample\"])))","pinned":true,"hidden":false,"mode":"ojs"},{"id":16,"value":"data0 = data0_notZerod.map((d) => d - data0_notZerod[0])","pinned":true,"hidden":false,"mode":"ojs"},{"id":1454,"value":"data1_notZerod = linearInterpolate(g1.map((d) => parseFloat(d[\" Sample\"]))) ","pinned":true,"hidden":false,"mode":"ojs"},{"id":19,"value":"data1 = data1_notZerod.map((d) => d - data1_notZerod[0])","pinned":true,"hidden":false,"mode":"ojs"},{"id":1448,"value":"data2_notZerod = linearInterpolate(g2.map((d) => parseFloat(d[\" Sample\"])))","pinned":true,"hidden":false,"mode":"ojs"},{"id":22,"value":"data2 = data2_notZerod.map((d) => d - data2_notZerod[0])","pinned":true,"hidden":false,"mode":"ojs"},{"id":629,"value":"data3_notZerod = linearInterpolate(g3.map((d) => parseFloat(d[\" Sample\"])))","pinned":true,"hidden":false,"mode":"ojs"},{"id":1434,"value":"data3 = data3_notZerod.map((d) => d - data3_notZerod[0])","pinned":true,"hidden":false,"mode":"ojs"},{"id":1367,"value":"// Linear interpolation function\nfunction linearInterpolate(arr) {\n let result = [];\n \n for (let i = 0; i < arr.length; i++) {\n if (arr[i] === 999999) {\n let startIndex = i;\n let endIndex = i;\n\n // Move forward to find the first non-999999 value\n while (endIndex < arr.length && arr[endIndex] === 999999) {\n endIndex++;\n }\n\n // If there are surrounding non-999999 values, set the starting value to the next non-999999 value\n if (endIndex < arr.length) {\n const nextValue = arr[endIndex];\n result.push(nextValue);\n } else {\n // If there are no subsequent non-999999 values, keep the original value\n result.push(arr[i]);\n }\n\n // Skip consecutive 999999 values\n i = endIndex - 1;\n } else {\n // If the value is not 999999, keep the original value\n result.push(arr[i]);\n }\n }\n\n return result;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":10,"value":"df = Object({\n Date: date_rng,\n Series0: data0,\n Series1: data1,\n Series2: data2,\n Series3: data3\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1036,"value":"function detrended(x,y){\n var lr = {};\n var n = y.length;\n var sum_x = 0;\n var sum_y = 0;\n var sum_xy = 0;\n var sum_xx = 0;\n var sum_yy = 0;\n\n for (var i = 0; i < y.length; i++) {\n\n sum_x += x[i];\n sum_y += y[i];\n sum_xy += (x[i]*y[i]);\n sum_xx += (x[i]*x[i]);\n sum_yy += (y[i]*y[i]);\n } \n\n lr['slope'] = (n * sum_xy - sum_x * sum_y) / (n*sum_xx - sum_x * sum_x);\n lr['intercept'] = (sum_y - lr.slope * sum_x)/n;\n lr['r2'] = Math.pow((n*sum_xy - sum_x*sum_y)/Math.sqrt((n*sum_xx-sum_x*sum_x)*(n*sum_yy-sum_y*sum_y)),2);\n\n return x.map(d => y[d] - d*lr.slope);\n}\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":98,"value":"## Plotting","pinned":false,"hidden":false,"mode":"md"},{"id":1029,"value":"plotDiv = html`

`","pinned":true,"hidden":false,"mode":"ojs"},{"id":968,"value":"$31","pinned":true,"hidden":false,"mode":"ojs"},{"id":222,"value":"selectedTimestamps = []","pinned":false,"hidden":false,"mode":"ojs"},{"id":1099,"value":"datax = []","pinned":true,"hidden":false,"mode":"ojs"},{"id":1531,"value":"gname = []","pinned":true,"hidden":false,"mode":"ojs"},{"id":1528,"value":"datay = []","pinned":true,"hidden":false,"mode":"ojs"},{"id":971,"value":"plot.on(\"plotly_click\", (event) => {\n const selected = new Date(event.points[0].x);\n datax.push(event.points[0].x);\n datay.push(event.points[0].y);\n gname.push(event.points[0].data.name);\n selectedTimestamps.push(selected);\n // Create a copy of the existing selected points trace\n const selectedPointsTrace = {\n x: selectedTimestamps,\n y: datay,\n mode: \"markers\",\n marker: { size: 10, color: \"red\" },\n type: \"scatter\",\n name: \"Selected Points\"\n };\n\n // Add the new trace for selected points without affecting the existing traces\n Plotly.addTraces(plot, selectedPointsTrace);\n\n const timestampsText =\n datax[datax.length - 1] + \" \" + gname[gname.length - 1] + \"\\n\";\n\n textArea.value += timestampsText;\n})","pinned":true,"hidden":false,"mode":"ojs"},{"id":1005,"value":"function clearPlot() {\n // Clear the last selected timestamp from text area, but preserve annotations\n selectedTimestamps.pop();\n datax.pop();\n datay.pop();\n gname.pop();\n // Delete all traces from the plot\n Plotly.deleteTraces(plot, [-1]); // Delete the last trace, which is the \"Selected Points\" trace\n // Redraw the plot\n Plotly.redraw(plot);\n const timestampsText = textArea.value\n .split(\"\\n\")\n .splice(0, selectedTimestamps.length)\n .join(\"\\n\");\n textArea.value = timestampsText + \"\\n\";\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":1001,"value":"clearButton = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":995,"value":"// Function to clear the plot\nclearButton.onclick = () => clearPlot()","pinned":false,"hidden":false,"mode":"ojs"},{"id":94,"value":"## Text box","pinned":false,"hidden":false,"mode":"md"},{"id":978,"value":"textArea = html``","pinned":false,"hidden":false,"mode":"ojs"},{"id":33,"value":"textStyle = {\n textArea.style.border = \"3px solid #ccc\";\n textArea.style.padding = \"10px\";\n textArea.style.width = \"90%\";\n textArea.style.height = \"100px\";\n textArea.placeholder = \"Selected Timestamps and Custom Text\";\n}","pinned":true,"hidden":false,"mode":"ojs"},{"id":947,"value":"// Function to save selected timestamps and custom text to CSV\nfunction saveToCSV() {\n const headerText = \"# Created \" + Date() + \"\\n# Time\";\n const customText = textArea.value;\n const csvContent = `${headerText}\\nCustom Text:\\n${customText}`;\n\n const blob = new Blob([csvContent], { type: \"text/csv\" });\n const url = URL.createObjectURL(blob);\n\n const a = document.createElement(\"a\");\n a.href = url;\n a.download =\n network +\n \".\" +\n station +\n \".\" +\n rate +\n \".\" +\n start +\n \".\" +\n end +\n \".offsetTimes.txt\";\n document.body.appendChild(a);\n a.click();\n document.body.removeChild(a);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":80,"value":"saveButton = html``;","pinned":true,"hidden":false,"mode":"ojs"},{"id":82,"value":"saveButton.onclick = () => saveToCSV();\n","pinned":true,"hidden":false,"mode":"ojs"},{"id":940,"value":"## Import","pinned":false,"hidden":false,"mode":"md"},{"id":938,"value":"Plotly = require(\"https://cdn.plot.ly/plotly-latest.min.js\")","pinned":true,"hidden":false,"mode":"ojs"}],"files":"$W32","stdlib":"1","titling":"auto"}