public static string Text(LineScanFolder ls)
{
// name, unit, comment, data...
int dataPoints = ls.imgG.height;
string[] csvLines = new string[dataPoints + 3];
// times (ms)
csvLines[0] = "Time, ";
csvLines[1] = "ms, ";
csvLines[2] = ls.folderName + ", ";
for (int i = 0; i < dataPoints; i++)
csvLines[i + 3] = Math.Round(ls.timesMsec[i], 3).ToString() + ", ";
// raw PMT values (R)
if (ls.curveR != null)
{
csvLines[0] += "R, ";
csvLines[1] += "AFU, ";
csvLines[2] += ", ";
for (int i = 0; i < dataPoints; i++)
csvLines[i + 3] += Math.Round(ls.curveR[i], 3).ToString() + ", ";
}
// raw PMT values (G)
if (ls.curveG != null)
{
csvLines[0] += "G, ";
csvLines[1] += "AFU, ";
csvLines[2] += ", ";
for (int i = 0; i < dataPoints; i++)
csvLines[i + 3] += Math.Round(ls.curveG[i], 3).ToString() + ", ";
}
// delta raw PMT values (G)
if (ls.curveDeltaG != null)
{
csvLines[0] += "dG, ";
csvLines[1] += "AFU, ";
csvLines[2] += ", ";
for (int i = 0; i < dataPoints; i++)
csvLines[i + 3] += Math.Round(ls.curveDeltaG[i], 3).ToString() + ", ";
csvLines[0] += "f(dG), ";
csvLines[1] += "AFU, ";
csvLines[2] += "filtered, ";
double[] filteredChopped = ls.GetFilteredYs(ls.curveDeltaG);
double[] filtered = new double[dataPoints];
for (int i = 0; i < dataPoints; i++)
filtered[i] = 0;
Array.Copy(filteredChopped, 0, filtered, ls.filterPx * 2, filteredChopped.Length);
for (int i = 0; i < dataPoints; i++)
if (i < ls.filterPx * 2 || i > (dataPoints - ls.filterPx * 2 * 2))
csvLines[i + 3] += ", ";
else
csvLines[i + 3] += Math.Round(filtered[i], 3).ToString() + ", ";
}
// Green over Red
if (ls.curveGoR != null)
{
csvLines[0] += "G/R, ";
csvLines[1] += "%, ";
csvLines[2] += ", ";
for (int i = 0; i < dataPoints; i++)
csvLines[i + 3] += Math.Round(ls.curveGoR[i], 3).ToString() + ", ";
}
// Delta Green over Red
if (ls.curveDeltaGoR != null)
{
csvLines[0] += "dG/R, ";
csvLines[1] += "%, ";
csvLines[2] += ", ";
for (int i = 0; i < dataPoints; i++)
csvLines[i + 3] += Math.Round(ls.curveDeltaGoR[i], 3).ToString() + ", ";
csvLines[0] += "f(dG/R), ";
csvLines[1] += "AFU, ";
csvLines[2] += "filtered, ";
double[] filteredChopped = ls.GetFilteredYs(ls.curveDeltaGoR);
double[] filtered = new double[dataPoints];
for (int i = 0; i < dataPoints; i++)
filtered[i] = 0;
Array.Copy(filteredChopped, 0, filtered, ls.filterPx * 2, filteredChopped.Length);
for (int i = 0; i < dataPoints; i++)
if (i < ls.filterPx * 2 || i > (dataPoints - ls.filterPx * 2 * 2))
csvLines[i + 3] += ", ";
else
csvLines[i + 3] += Math.Round(filtered[i], 3).ToString() + ", ";
}
// convert to CSV
string csv = "";
foreach (string line in csvLines)
csv += line + "\n";
return csv;
}
public void UpdateGraphs(bool skipUpdatesIfBusy = true)
{
// use pretty colors
Color blue = System.Drawing.ColorTranslator.FromHtml("#1f77b4");
Color lightBlue = System.Drawing.ColorTranslator.FromHtml("#a2d1f2");
Color red = System.Drawing.ColorTranslator.FromHtml("#d62728");
Color lightRed = System.Drawing.ColorTranslator.FromHtml("#e63738");
Color green = System.Drawing.ColorTranslator.FromHtml("#2ca02c");
Color lightGreen = System.Drawing.ColorTranslator.FromHtml("#98df8a");
Color colorBaselineMarks = System.Drawing.ColorTranslator.FromHtml("#666666");
Color colorPeak = System.Drawing.ColorTranslator.FromHtml("#d62728");
Color colorZero = System.Drawing.ColorTranslator.FromHtml("#000000");
if (!lsFolder.isValid)
{
return;
}
if (skipUpdatesIfBusy && busyUpdating)
{
return;
}
else
{
busyUpdating = true;
}
// do the analysis
lsFolder.GenerateAnalysisCurves();
// update the structure profile plot
double[] columnPixelIndices = new double[lsFolder.imgG.width];
for (int i = 0; i < columnPixelIndices.Length; i++)
{
columnPixelIndices[i] = i;
}
formsPlot2.plt.Clear();
if ((cbR.Enabled && cbR.Checked))
{
double[] columnBrightnessR = ImageDataTools.GetAverageLeftright(lsFolder.imgR);
formsPlot2.plt.PlotScatter(columnPixelIndices, columnBrightnessR, markerSize: 0, color: red);
}
if ((cbG.Enabled && cbG.Checked))
{
double[] columnBrightnessG = ImageDataTools.GetAverageLeftright(lsFolder.imgG);
formsPlot2.plt.PlotScatter(columnPixelIndices, columnBrightnessG, markerSize: 0, color: green);
}
formsPlot2.plt.AxisAuto(0, .1);
formsPlot2.plt.PlotVLine(lsFolder.structure1, color: Color.Black);
formsPlot2.plt.PlotVLine(lsFolder.structure2, color: Color.Black);
// update styling before the render
formsPlot2.plt.YLabel("Intensity");
formsPlot2.plt.XLabel("Position (pixel number)");
formsPlot2.Render();
// update the time series plot
curveToCopy = null;
formsPlot1.plt.Clear();
if (cbRatio.Checked && cbRatio.Enabled)
{
if (cbDelta.Enabled && cbDelta.Checked)
{
curveToCopy = lsFolder.GetFilteredYs(lsFolder.curveDeltaGoR);
formsPlot1.plt.PlotScatter(lsFolder.timesMsec, lsFolder.curveDeltaGoR, color: lightBlue, lineWidth: 0, markerSize: 2);
formsPlot1.plt.PlotScatter(lsFolder.GetFilteredXs(), curveToCopy, markerSize: 0, color: blue, lineWidth: 2);
formsPlot1.plt.PlotVLine(lsFolder.baseline1 * lsFolder.scanLinePeriod, color: colorBaselineMarks, lineWidth: 2);
formsPlot1.plt.PlotVLine(lsFolder.baseline2 * lsFolder.scanLinePeriod, color: colorBaselineMarks, lineWidth: 2);
formsPlot1.plt.PlotHLine(0, color: colorZero, lineWidth: 2);
formsPlot1.plt.PlotHLine(curveToCopy.Max(), color: colorPeak, lineWidth: 2);
formsPlot1.plt.YLabel("Delta G/R (%)");
}
else
{
curveToCopy = lsFolder.GetFilteredYs(lsFolder.curveGoR);
formsPlot1.plt.PlotScatter(lsFolder.timesMsec, lsFolder.curveGoR, color: lightBlue, lineWidth: 0, markerSize: 2);
formsPlot1.plt.PlotScatter(lsFolder.GetFilteredXs(), curveToCopy, markerSize: 0, color: blue, lineWidth: 2);
formsPlot1.plt.YLabel("G/R (%)");
formsPlot1.plt.PlotHLine(curveToCopy.Max(), color: colorPeak, lineWidth: 2);
}
}
else if ((cbR.Checked && cbR.Enabled) && (cbG.Checked && cbG.Enabled))
{
formsPlot1.plt.PlotScatter(lsFolder.timesMsec, lsFolder.curveG, markerSize: 0, color: green);
formsPlot1.plt.PlotScatter(lsFolder.timesMsec, lsFolder.curveR, markerSize: 0, color: red);
formsPlot1.plt.YLabel("G and R (AFU)");
}
else if (cbR.Checked && cbR.Enabled)
{
formsPlot1.plt.PlotScatter(lsFolder.timesMsec, lsFolder.curveR, markerSize: 0, color: red);
formsPlot1.plt.YLabel("R (AFU)");
}
else if (cbG.Checked && cbG.Enabled)
{
if (cbDelta.Enabled && cbDelta.Checked)
{
curveToCopy = lsFolder.GetFilteredYs(lsFolder.curveDeltaG);
formsPlot1.plt.PlotScatter(lsFolder.timesMsec, lsFolder.curveDeltaG, lineWidth: 0, color: lightGreen, markerSize: 2);
formsPlot1.plt.PlotScatter(lsFolder.GetFilteredXs(), curveToCopy, markerSize: 0, color: green, lineWidth: 2);
formsPlot1.plt.PlotVLine(lsFolder.baseline1 * lsFolder.scanLinePeriod, color: colorBaselineMarks, lineWidth: 2);
formsPlot1.plt.PlotVLine(lsFolder.baseline2 * lsFolder.scanLinePeriod, color: colorBaselineMarks, lineWidth: 2);
formsPlot1.plt.PlotHLine(0, color: colorZero, lineWidth: 2);
formsPlot1.plt.PlotHLine(curveToCopy.Max(), color: colorPeak, lineWidth: 2);
formsPlot1.plt.YLabel("Delta G (%)");
}
else
{
curveToCopy = lsFolder.GetFilteredYs(lsFolder.curveG);
formsPlot1.plt.PlotScatter(lsFolder.timesMsec, lsFolder.curveG, lineWidth: 0, color: lightGreen, markerSize: 2);
formsPlot1.plt.PlotScatter(lsFolder.GetFilteredXs(), curveToCopy, markerSize: 0, color: green, lineWidth: 2);
formsPlot1.plt.PlotHLine(curveToCopy.Max(), color: colorPeak, lineWidth: 2);
formsPlot1.plt.YLabel("G (AFU)");
}
}
formsPlot1.plt.AxisAuto(.05, .1);
// make ticks smaller than typical
var tickFont = new Font(FontFamily.GenericMonospace, (float)10.0);
//formsPlot1.plt.Ticks(font: tickFont);
//formsPlot2.plt.Ticks(font: tickFont);
// update styling before the render
formsPlot1.plt.XLabel("Time (milliseconds)");
formsPlot1.Render();
// update peak label
if (curveToCopy == null)
{
lblPeak.Text = $"";
}
else
{
lblPeak.Text = $"{Math.Round(curveToCopy.Max(), 2)}%";
}
Application.DoEvents();
busyUpdating = false;
}
