public void Plot(StylisedPeak stylisedPeak)
{
// Get observations
Vector vector;
if (stylisedPeak.ForceObservations != null)
{
vector = stylisedPeak.ForceObservations;
}
else
{
IntensityMatrix matrix = this._core.Options.SelectedMatrix;
if (matrix == null)
{
vector = null;
}
else
{
vector = matrix.Find(stylisedPeak.Peak);
}
}
Debug.WriteLine("PeakPlot: " + stylisedPeak);
Dictionary <string, MCharting.Series> seriesNames = new Dictionary <string, MCharting.Series>();
Peak peak = stylisedPeak?.Peak;
// Clear plot
MCharting.Plot plot = this.PrepareNewPlot(stylisedPeak != null && !stylisedPeak.IsPreview, peak, vector?.Source);
try // <- CompletNewPlot
{
// Get selection
this.SelectedPeak = peak;
this.SetCaption("Plot of {0}.", peak);
if (peak == null)
{
return;
}
// Get options
StylisedPeakOptions opts = stylisedPeak.OverrideDefaultOptions ?? new StylisedPeakOptions(this._core);
// Get order data
ObservationInfo[] obsOrder = this._core.Observations.ToArray();
// Group legends
IEnumerable <GroupInfoBase> order = opts.ShowAcqisition ? (IEnumerable <GroupInfoBase>)opts.ViewBatches : (IEnumerable <GroupInfoBase>)opts.ViewGroups;
Dictionary <GroupInfoBase, MCharting.Series> groupLegends = this.DrawLegend(plot, order);
if (vector == null)
{
return;
}
// Show acquisition and batches?
if (opts.ShowAcqisition)
{
// --- RAW DATA (points) ---
MCharting.Series legendEntry = new MCharting.Series();
legendEntry.Name = "Observations";
legendEntry.Style.DrawPoints = new SolidBrush(Color.Black);
plot.LegendEntries.Add(legendEntry);
this.AddToPlot(plot, peak, seriesNames, vector, "Raw data", opts, EPlot.ByBatch, groupLegends, legendEntry);
// --- TREND (thick line) ---
if (stylisedPeak.ForceTrend != null)
{
MCharting.Series legendEntry2 = new MCharting.Series();
legendEntry2.Name = "Trend";
legendEntry2.Style.DrawLines = new Pen(Color.Black, this._core.Options.LineWidth);
legendEntry2.Style.DrawLines.Width = 4;
plot.LegendEntries.Add(legendEntry2);
this.AddToPlot(plot, peak, seriesNames, stylisedPeak.ForceTrend, "Trend data", opts, EPlot.ByBatch | EPlot.DrawLine | EPlot.DrawBold, groupLegends, legendEntry2);
}
this.DrawLabels(plot, opts.ConditionsSideBySide, order, opts.DrawExperimentalGroupAxisLabels);
return;
}
// Sort data
ConfigurationTrend trend = opts.SelectedTrend;
Vector avg = stylisedPeak.ForceTrend ?? trend.CreateTrend(this._core, vector);
Vector min = MinSmoother.CreateTrend(this._core, vector);
Vector max = MaxSmoother.CreateTrend(this._core, vector);
// --- PLOT MEAN & SD (lines across)
if (opts.ShowVariableMean)
{
this.AddMeanAndSdLines(plot, opts, vector.Values, peak, groupLegends);
}
// --- RANGE (shaded area) ---
if (opts.ShowRanges)
{
this.AddUpperAndLowerShade(plot, opts, seriesNames, peak, min, max, groupLegends);
}
// --- RAW DATA (points) ---
if (opts.ShowPoints && !stylisedPeak.IsPreview)
{
MCharting.Series legendEntry = new MCharting.Series();
legendEntry.Name = "Observations";
legendEntry.Style.DrawPoints = new SolidBrush(Color.Black);
plot.LegendEntries.Add(legendEntry);
this.AddToPlot(plot, peak, seriesNames, vector, "Raw data", opts, EPlot.None, groupLegends, legendEntry);
}
// --- RANGE (lines) ---
if (opts.ShowMinMax)
{
MCharting.Series legendEntry = new MCharting.Series();
legendEntry.Name = "Range min/max";
legendEntry.Style.DrawLines = new Pen(Color.Gray, this._core.Options.LineWidth);
plot.LegendEntries.Add(legendEntry);
this.AddToPlot(plot, peak, seriesNames, min, "Min value", opts, EPlot.DrawLine, groupLegends, legendEntry);
this.AddToPlot(plot, peak, seriesNames, max, "Max value", opts, EPlot.DrawLine, groupLegends, legendEntry);
}
// --- TREND (thick line) ---
if (opts.ShowTrend)
{
MCharting.Series legendEntry = new MCharting.Series();
legendEntry.Name = "Trend";
legendEntry.Style.DrawLines = new Pen(Color.Black, this._core.Options.LineWidth);
legendEntry.Style.DrawLines.Width = this._core.Options.LineWidth * 4;
plot.LegendEntries.Add(legendEntry);
this.AddToPlot(plot, peak, seriesNames, avg, "Trend data", opts, EPlot.DrawLine | EPlot.DrawBold, groupLegends, legendEntry);
}
// --- LABELS ---
this.DrawLabels(plot, opts.ConditionsSideBySide, order, opts.DrawExperimentalGroupAxisLabels);
}
finally
{
this.CompleteNewPlot(plot);
}
}
/// <summary>
/// Adds the mean and standard deviation lines to the plot.
/// </summary>
/// <param name="plot">Target plot</param>
/// <param name="o">Options</param>
/// <param name="observations">Set of observations</param>
/// <param name="peak">Peak being plotted</param>
/// <param name="groupLegends">Dictionary to get group legends from</param>
private void AddMeanAndSdLines(
MCharting.Plot plot,
StylisedPeakOptions o,
IReadOnlyList <double> observations,
Peak peak,
Dictionary <GroupInfoBase, MCharting.Series> groupLegends)
{
MCharting.Series legendEntry = new MCharting.Series();
legendEntry.Name = "Std. Dev. Min/Max";
legendEntry.Style.DrawLines = new Pen(Color.Gray, this._core.Options.LineWidth);
legendEntry.Style.DrawLines.DashStyle = DashStyle.Dot;
plot.LegendEntries.Add(legendEntry);
MCharting.Series legendEntry2 = new MCharting.Series();
legendEntry2.Name = "Mean";
legendEntry2.Style.DrawLines = new Pen(Color.Black, this._core.Options.LineWidth);
plot.LegendEntries.Add(legendEntry2);
// Iterate the types
foreach (GroupInfo group in o.ViewGroups)
{
// Get the Y values
double yMean = Maths.Mean(observations);
double sd = Maths.StdDev(observations, yMean);
double yMin = yMean - sd / 2;
double yMax = yMean + sd / 2;
// Get the X values
int xLeft = @group.Range.Min;
int xRight = @group.Range.Max;
if (o.ConditionsSideBySide)
{
int typeOffset = this.GetTypeOffset(o.ViewGroups, @group);
xLeft += typeOffset;
xRight += typeOffset;
}
// Create the series
MCharting.Series sMean = plot.Series.Add(@group.DisplayName + ": Mean");
MCharting.Series sMin = plot.Series.Add(@group.DisplayName + ": StdDevMin");
MCharting.Series sMax = plot.Series.Add(@group.DisplayName + ": StdDevMax");
sMean.Style.StrictOrder = SeriesStyle.EOrder.X;
sMin.Style.StrictOrder = SeriesStyle.EOrder.X;
sMax.Style.StrictOrder = SeriesStyle.EOrder.X;
sMean.ApplicableLegends.Add(groupLegends[group]);
sMin.ApplicableLegends.Add(groupLegends[group]);
sMax.ApplicableLegends.Add(groupLegends[group]);
sMean.ApplicableLegends.Add(legendEntry2);
sMin.ApplicableLegends.Add(legendEntry);
sMax.ApplicableLegends.Add(legendEntry);
sMean.Tag = peak;
sMin.Tag = peak;
sMax.Tag = peak;
Color c = @group.ColourLight;
sMean.Style.DrawLines = new Pen(c, this._core.Options.LineWidth);
sMin.Style.DrawLines = new Pen(c, this._core.Options.LineWidth);
sMax.Style.DrawLines = new Pen(c, this._core.Options.LineWidth);
sMean.Style.DrawLines.DashStyle = DashStyle.Dash;
sMin.Style.DrawLines.DashStyle = DashStyle.Dash;
sMax.Style.DrawLines.DashStyle = DashStyle.Dash;
// Add the points
AddDataPoint(sMean, xLeft, yMean, @group);
AddDataPoint(sMean, xRight, yMean, @group);
AddDataPoint(sMin, xLeft, yMin, @group);
AddDataPoint(sMin, xRight, yMin, @group);
AddDataPoint(sMax, xLeft, yMax, @group);
AddDataPoint(sMax, xRight, yMax, @group);
}
}
/// <summary>
/// Plots a line of conditions on the graph
/// </summary>
/// <param name="peak">Peak to plot</param>
/// <param name="seriesNames">Names of existing series</param>
/// <param name="intensities">Intensity values (y)</param>
/// <param name="seriesName">Name of this series</param>
/// <param name="xInfo">Info for x position of intentises (in same order as intensities)</param>
/// <param name="line">Line or dots/</param>
/// <param name="bold">Bold line?</param>
/// <param name="isConditions">Order by conditions or obervations</param>
private void AddToPlot(
MCharting.Plot plot,
Peak peak,
Dictionary <string, MCharting.Series> seriesNames,
Vector intensities,
string seriesName,
StylisedPeakOptions o,
EPlot draw,
Dictionary <GroupInfoBase, MCharting.Series> groupLegends,
MCharting.Series legend)
{
// Iterate whatever it is we're iterating
for (int i = 0; i < intensities.Observations.Length; ++i)
{
// Get the values
ObservationInfo obs = intensities.Observations[i];
// Name the series
if (draw.HasFlag(EPlot.ByBatch))
{
if (!o.ViewBatches.Contains(obs.Batch))
{
continue;
}
}
else
{
if (!o.ViewGroups.Contains(obs.Group))
{
continue;
}
}
bool colorByBatch = draw.HasFlag(EPlot.ByBatch);
GroupInfoBase seriesUsing = colorByBatch ? (GroupInfoBase)obs.Batch : (GroupInfoBase)obs.Group;
string name = seriesName + " for " + seriesUsing.DisplayName;
MCharting.Series series;
// Create the series (if required)
if (!seriesNames.ContainsKey(name))
{
series = plot.Series.Add(name);
series.Style.StrictOrder = SeriesStyle.EOrder.X;
series.ApplicableLegends.Add(groupLegends[seriesUsing]);
series.ApplicableLegends.Add(legend);
seriesNames.Add(name, series);
series.Tag = peak;
Color colour = (draw.HasFlag(EPlot.DrawBold) || draw.HasFlag(EPlot.ByBatch) || !o.ShowTrend) ? seriesUsing.Colour : seriesUsing.ColourLight;
if (draw.HasFlag(EPlot.DrawLine))
{
series.Style.DrawLines = new Pen(colour);
series.Style.DrawLines.Width = draw.HasFlag(EPlot.DrawBold) ? this._core.Options.LineWidth * 4 : this._core.Options.LineWidth;
}
else
{
UiControls.CreateIcon(series, seriesUsing);
}
}
else
{
series = seriesNames[name];
}
// Get the X position
double xPos;
if (draw.HasFlag(EPlot.ByBatch))
{
xPos = obs.Order;
if (o.ConditionsSideBySide)
{
xPos += this.GetBatchOffset(o.ViewBatches, obs.Batch);
}
}
else
{
xPos = obs.Time;
if (o.ConditionsSideBySide)
{
xPos += this.GetTypeOffset(o.ViewGroups, obs.Group);
}
}
// Get the Y position
double yPos = intensities.Values[i];
if (double.IsNaN(yPos) || double.IsInfinity(yPos))
{
yPos = 0;
}
// Create the point
MCharting.DataPoint cdp = new MCharting.DataPoint(xPos, yPos);
IntensityInfo tag = new IntensityInfo(obs.Time, obs.Rep, obs.Group, yPos);
cdp.Tag = tag;
series.Points.Add(cdp);
}
}
private void AddUpperAndLowerShade(
MCharting.Plot plot,
StylisedPeakOptions o,
Dictionary <string, MCharting.Series> seriesNames,
Peak peak,
Vector min,
Vector max,
Dictionary <GroupInfoBase, MCharting.Series> groupLegends)
{
MCharting.Series legendEntry = new MCharting.Series();
legendEntry.Name = "Range";
legendEntry.Style.DrawVBands = new SolidBrush(Color.Gray);
plot.LegendEntries.Add(legendEntry);
// Iterate the conditions
for (int i = 0; i < min.Observations.Length; i++)
{
ObservationInfo obs = min.Observations[i];
Debug.Assert(max.Observations[i] == min.Observations[i], "Expected max and min trends to match sequences.");
if (o.ViewGroups.Contains(obs.Group))
{
// Name the series
string name = "Range for " + obs.Group.DisplayName;
MCharting.Series series;
// Create the series (if required)
if (!seriesNames.ContainsKey(name))
{
series = plot.Series.Add(name);
series.Style.StrictOrder = SeriesStyle.EOrder.X;
Color c = obs.Group.ColourLight;
c = Color.FromArgb(0x80, c.R, c.G, c.B);
series.Tag = peak;
series.Style.DrawVBands = obs.Group.CreateBrush(c);
series.ApplicableLegends.Add(groupLegends[obs.Group]);
series.ApplicableLegends.Add(legendEntry);
seriesNames.Add(name, series);
}
else
{
series = seriesNames[name];
}
// Get the X coordinate
int typeOffset = this.GetTypeOffset(o.ViewGroups, obs.Group);
double xVal = obs.Time;
if (o.ConditionsSideBySide)
{
xVal += typeOffset;
}
// Get the Y coordinates
double yMin = min.Values[i];
double yMax = max.Values[i];
if (double.IsNaN(yMin) || double.IsInfinity(yMin))
{
yMin = 0;
}
if (double.IsNaN(yMax) || double.IsInfinity(yMax))
{
yMax = 0;
}
// Create the point
IntensityInfo info1 = new IntensityInfo(obs.Time, null, obs.Group, yMin);
IntensityInfo info2 = new IntensityInfo(obs.Time, null, obs.Group, yMax);
MCharting.DataPoint cdp = new MCharting.DataPoint(xVal, new[] { yMin, yMax });
cdp.Tag = new[] { info1, info2 };
series.Points.Add(cdp);
}
}
}
