private Image <Bgr, byte> FillGraphImage(Size imgSize)
{
string curDirPath = Directory.GetCurrentDirectory() + "\\logs\\";
DirectoryInfo dirInfo = new DirectoryInfo(curDirPath);
List <Dictionary <string, object> > lReadData = new List <Dictionary <string, object> >();
if (defaultGraphsTimeSpan)
{
graphsTimeSpan = new Tuple <DateTime, DateTime>(DateTime.UtcNow.AddDays(-1), DateTime.UtcNow);
}
GraphVariablesTypes currVarType = GraphVariablesTypes.none;
if (rbtnPressureGraph.Checked)
{
currVarType = GraphVariablesTypes.Pressure;
}
if (rbtnAirTempGraph.Checked)
{
currVarType = GraphVariablesTypes.AirTemp;
}
if (rbtnWaterTempGraph.Checked)
{
currVarType = GraphVariablesTypes.WaterTemp;
}
if (rbtnWindSpeedGraph.Checked)
{
currVarType = GraphVariablesTypes.WindSpeed;
}
if ((fRenderer == null) || (!prevGraphsTimeSpan.Equals(graphsTimeSpan) || currVarType != prevGraphVariable))
{
fRenderer = new MultipleScatterAndFunctionsRepresentation(imgSize);
switch (currVarType)
{
case GraphVariablesTypes.Pressure:
{
currValueColor = new Bgr(Color.Blue);
break;
}
case GraphVariablesTypes.AirTemp:
{
currValueColor = new Bgr(Color.Red);
break;
}
case GraphVariablesTypes.WaterTemp:
{
currValueColor = new Bgr(Color.RoyalBlue);
break;
}
case GraphVariablesTypes.WindSpeed:
{
currValueColor = new Bgr(Color.Gray);
break;
}
default:
{
currValueColor = new Bgr(Color.Blue);
break;
}
}
}
else if (fRenderer != null)
{
if (fRenderer.TheImage.Size != imgSize)
{
fRenderer.ResizeCanvas(imgSize);
}
}
if (!prevGraphsTimeSpan.Equals(graphsTimeSpan))
{
IEnumerable <string> ncFileNames = Directory.EnumerateFiles(curDirPath,
"IoffeVesselInfoStream-MeteoDataLog-*.nc",
SearchOption.TopDirectoryOnly);
foreach (string ncFileName in ncFileNames)
{
Tuple <DateTime, DateTime> currFileDateTimeRange = null;
try
{
currFileDateTimeRange = ServiceTools.GetNetCDFfileTimeStampsRange(ncFileName);
}
catch (Exception ex)
{
#region report
#if DEBUG
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
theLogWindow = ServiceTools.LogAText(theLogWindow,
"an exception has been thrown during file reading: " + Environment.NewLine + ncFileName +
Environment.NewLine + "message: " + ex.Message + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription());
});
#else
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
ServiceTools.logToTextFile(errorLogFilename,
"an exception has been thrown during file reading: " + Environment.NewLine + ncFileName +
Environment.NewLine + "message: " + ex.Message + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription(), true, true);
});
#endif
#endregion report
}
if (currFileDateTimeRange == null)
{
continue;
}
if ((currFileDateTimeRange.Item1 >= graphsTimeSpan.Item1) &&
(currFileDateTimeRange.Item1 <= graphsTimeSpan.Item2) ||
(currFileDateTimeRange.Item2 >= graphsTimeSpan.Item1) &&
(currFileDateTimeRange.Item2 <= graphsTimeSpan.Item2))
{
Dictionary <string, object> dictFileData = null;
try
{
dictFileData = NetCDFoperations.ReadDataFromFile(ncFileName);
}
catch (Exception ex)
{
#region report
#if DEBUG
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
theLogWindow = ServiceTools.LogAText(theLogWindow,
"an exception has been thrown during file reading: " + Environment.NewLine +
ncFileName +
Environment.NewLine + "message: " + ex.Message + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription());
});
#else
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
ServiceTools.logToTextFile(errorLogFilename,
"an exception has been thrown during file reading: " + Environment.NewLine + ncFileName +
Environment.NewLine + "message: " + ex.Message + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription(), true, true);
});
#endif
#endregion report
}
if (dictFileData != null)
{
lReadData.Add(dictFileData);
}
}
}
foreach (Dictionary <string, object> currFileDataDict in lReadData)
{
if (currFileDataDict == null)
{
continue;
}
string varNameDateTime = "DateTime";
List <long> currFileDateTimeLongTicksList =
new List <long>((currFileDataDict[varNameDateTime] as long[]));
List <DateTime> currFileDateTimeList =
currFileDateTimeLongTicksList.ConvertAll(longVal => new DateTime(longVal));
string varNameMeteoData = "MeteoData";
List <MeteoData> currFileMeteoDataList =
MeteoData.OfDenseMatrix(currFileDataDict[varNameMeteoData] as DenseMatrix);
if (tsMeteoDataForGraphs == null)
{
try
{
tsMeteoDataForGraphs = new TimeSeries <MeteoData>(currFileMeteoDataList, currFileDateTimeList,
true);
}
catch (Exception ex)
{
#region report
#if DEBUG
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
theLogWindow = ServiceTools.LogAText(theLogWindow,
"couldn`t create timeseries: exception has been thrown" + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription() + Environment.NewLine + "message: " +
ex.Message);
});
#else
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
ServiceTools.logToTextFile(errorLogFilename,
"couldn`t create timeseries: exception has been thrown" + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription() + Environment.NewLine + "message: " +
ex.Message, true, true);
});
#endif
#endregion report
}
}
else
{
try
{
tsMeteoDataForGraphs.AddSubseriaData(currFileMeteoDataList, currFileDateTimeList, true);
}
catch (Exception ex)
{
#region report
#if DEBUG
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
theLogWindow = ServiceTools.LogAText(theLogWindow,
"couldn`t create timeseries: exception has been thrown" + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription() + Environment.NewLine + "message: " +
ex.Message);
});
#else
ServiceTools.ExecMethodInSeparateThread(this, () =>
{
ServiceTools.logToTextFile(errorLogFilename,
"couldn`t create timeseries: exception has been thrown" + Environment.NewLine +
ServiceTools.CurrentCodeLineDescription() + Environment.NewLine + "message: " +
ex.Message, true, true);
});
#endif
#endregion report
}
}
}
if (tsMeteoDataForGraphs == null)
{
return(null);
}
tsMeteoDataForGraphs.SortByTimeStamps();
tsMeteoDataForGraphs.RemoveDuplicatedTimeStamps();
DateTime utcNow = DateTime.UtcNow;
if (defaultGraphsTimeSpan)
{
tsMeteoDataForGraphs.RemoveValues(dt => (utcNow - dt).TotalSeconds > 86400);
}
else
{
tsMeteoDataForGraphs.RemoveValues(
dt => !((dt >= graphsTimeSpan.Item1) && (dt <= graphsTimeSpan.Item2)));
}
List <TimeSeries <MeteoData> > subSeriesBy1Minute =
tsMeteoDataForGraphs.SplitByTimeSpan(new TimeSpan(0, 1, 0));
List <double> lSubSeriesEntriesCount = subSeriesBy1Minute.ConvertAll(subs => (double)subs.Count);
DescriptiveStatistics statsCounts = new DescriptiveStatistics(lSubSeriesEntriesCount);
aveMinuteEntriesCount = Convert.ToInt32(statsCounts.Mean);
// = tsMeteoData.TimeStamps.ConvertAll(dt => (dt - maxDateTime).TotalSeconds);
}
List <MeteoData> meteoDataList = tsMeteoDataForGraphs.DataValues;
DateTime maxDateTime = tsMeteoDataForGraphs.TimeStamps.Max();
if ((currVarType != prevGraphVariable) || !prevGraphsTimeSpan.Equals(graphsTimeSpan))
{
double minVarValue = 0.0d;
double maxVarValue = 1.0d;
List <double> currVarToShowValues = new List <double>();
switch (currVarType)
{
case GraphVariablesTypes.Pressure:
{
currVarToShowValues = meteoDataList.ConvertAll(mdt => mdt.pressure);
TimeSeries <double> currVarTS = new TimeSeries <double>(currVarToShowValues,
tsMeteoDataForGraphs.TimeStamps);
currVarTS.RemoveValues(dVal => dVal <= 900.0d);
currVarToShowValues = new List <double>(currVarTS.DataValues);
currFileSecondsList = currVarTS.TimeStamps.ConvertAll(dt => (dt - maxDateTime).TotalSeconds);
fRenderer.yAxisValuesConversionToRepresentTicksValues =
new Func <double, string>(dVal => dVal.ToString("F1"));
break;
}
case GraphVariablesTypes.AirTemp:
{
currVarToShowValues = meteoDataList.ConvertAll(mdt => mdt.airTemperature);
TimeSeries <double> currVarTS = new TimeSeries <double>(currVarToShowValues,
tsMeteoDataForGraphs.TimeStamps);
currVarTS.RemoveValues(dVal => ((dVal < -20.0d) || (dVal > 50.0d)));
currVarToShowValues = new List <double>(currVarTS.DataValues);
currFileSecondsList = currVarTS.TimeStamps.ConvertAll(dt => (dt - maxDateTime).TotalSeconds);
fRenderer.yAxisValuesConversionToRepresentTicksValues =
new Func <double, string>(dVal => dVal.ToString("F2"));
break;
}
case GraphVariablesTypes.WaterTemp:
{
currVarToShowValues = meteoDataList.ConvertAll(mdt => mdt.waterTemperature);
TimeSeries <double> currVarTS = new TimeSeries <double>(currVarToShowValues,
tsMeteoDataForGraphs.TimeStamps);
currVarTS.RemoveValues(dVal => ((dVal < -20.0d) || (dVal > 50.0d)));
currVarToShowValues = new List <double>(currVarTS.DataValues);
currFileSecondsList = currVarTS.TimeStamps.ConvertAll(dt => (dt - maxDateTime).TotalSeconds);
fRenderer.yAxisValuesConversionToRepresentTicksValues =
new Func <double, string>(dVal => dVal.ToString("F2"));
break;
}
case GraphVariablesTypes.WindSpeed:
{
currVarToShowValues = meteoDataList.ConvertAll(mdt => mdt.windSpeed);
TimeSeries <double> currVarTS = new TimeSeries <double>(currVarToShowValues,
tsMeteoDataForGraphs.TimeStamps);
currVarTS.RemoveValues(dVal => ((dVal < 0.0d) || (dVal > 50.0d)));
currVarToShowValues = new List <double>(currVarTS.DataValues);
currFileSecondsList = currVarTS.TimeStamps.ConvertAll(dt => (dt - maxDateTime).TotalSeconds);
fRenderer.yAxisValuesConversionToRepresentTicksValues =
new Func <double, string>(dVal => dVal.ToString("F1"));
break;
}
default:
return(null);
}
dvVarValues = DenseVector.OfEnumerable(currVarToShowValues);
dvVarValues = dvVarValues.Conv(Extensions.StandardConvolutionKernels.gauss, aveMinuteEntriesCount * 10);
}
fRenderer.dvScatterFuncValues.Add(dvVarValues);
fRenderer.dvScatterXSpace.Add(DenseVector.OfEnumerable(currFileSecondsList));
fRenderer.xAxisValuesConversionToRepresentTicksValues = (dValSec) =>
{
DateTime currDT = tsMeteoDataForGraphs.TimeStamps.Max().AddSeconds(dValSec);
return(currDT.ToString("yyyy-MM-dd" + Environment.NewLine + "HH:mm"));
};
fRenderer.scatterLineColors.Add(currValueColor);
fRenderer.scatterDrawingVariants.Add(SequencesDrawingVariants.polyline);
fRenderer.xSpaceMin = currFileSecondsList.Min();
fRenderer.xSpaceMax = currFileSecondsList.Max();
fRenderer.overallFuncMin = dvVarValues.Min();
fRenderer.overallFuncMax = dvVarValues.Max();
fRenderer.fixSpecifiedMargins = true;
fRenderer.Represent();
Image <Bgr, byte> retImg = fRenderer.TheImage;
// расположим надпись
string strSign = "current value: " + dvVarValues.Last().ToString("F2");
List <TextBarImage> textBarsCases = new List <TextBarImage>();
TextBarImage tbimTopLeftSign = new TextBarImage(strSign, retImg);
tbimTopLeftSign.PtSurroundingBarStart =
new Point(fRenderer.LeftServiceSpaceGapX + tbimTopLeftSign.textHalfHeight,
fRenderer.TopServiceSpaceGapY + tbimTopLeftSign.textHalfHeight);
textBarsCases.Add(tbimTopLeftSign);
TextBarImage tbimBtmLeftSign = new TextBarImage(strSign, retImg);
tbimBtmLeftSign.PtSurroundingBarStart = new Point(fRenderer.LeftServiceSpaceGapX + tbimBtmLeftSign.textHalfHeight,
retImg.Height - fRenderer.BtmServiceSpaceGapY - tbimBtmLeftSign.textHalfHeight - tbimBtmLeftSign.textHeight * 2);
textBarsCases.Add(tbimBtmLeftSign);
TextBarImage tbimTopRightSign = new TextBarImage(strSign, retImg);
tbimTopRightSign.PtSurroundingBarStart =
new Point(
retImg.Width - fRenderer.RightServiceSpaceGapX - tbimTopRightSign.textHalfHeight -
tbimTopRightSign.textBarSize.Width, fRenderer.TopServiceSpaceGapY + tbimTopLeftSign.textHalfHeight);
textBarsCases.Add(tbimTopRightSign);
TextBarImage tbimBtmRightSign = new TextBarImage(strSign, retImg);
tbimBtmRightSign.PtSurroundingBarStart =
new Point(
retImg.Width - fRenderer.RightServiceSpaceGapX - tbimBtmRightSign.textHalfHeight -
tbimBtmRightSign.textBarSize.Width,
retImg.Height - fRenderer.BtmServiceSpaceGapY - tbimBtmRightSign.textHalfHeight -
tbimBtmRightSign.textHeight * 2);
textBarsCases.Add(tbimBtmRightSign);
textBarsCases.Sort((case1, case2) => (case1.SubImageInTextRect.CountNonzero().Sum() > case2.SubImageInTextRect.CountNonzero().Sum()) ? 1 : -1);
MCvFont theFont = new MCvFont(Emgu.CV.CvEnum.FONT.CV_FONT_HERSHEY_PLAIN, 2.0d, 2.0d)
{
thickness = 2,
};
// retImg.Draw(strSign, textBarsCases[0].ptTextBaselineStart, Emgu.CV.CvEnum.FontFace.HersheyPlain, 2.0d, new Bgr(Color.Green), 2);
retImg.Draw(strSign, ref theFont, textBarsCases[0].ptTextBaselineStart, new Bgr(Color.Green));
retImg.Draw(textBarsCases[0].rectSurroundingBar, new Bgr(Color.Green), 2);
prevGraphsTimeSpan = graphsTimeSpan;
prevGraphVariable = currVarType;
return(retImg);
}
public static DenseVector Conv(this DenseVector dvSource, StandardConvolutionKernels kernelType, int kernelHalfWidth = 10)
{
DenseVector dvKernel = ConvKernel(kernelType, kernelHalfWidth);
return(dvSource.Conv(dvKernel));
}
/// <summary>
/// Cartesians to polar.
/// </summary>
/// <param name="dmData">input densematrix data.</param>
/// <param name="centerPoint">The center point of new polar system - in old cartesian coordinate system.</param>
/// <param name="dmMask">The mask matrix.</param>
/// <param name="angleGridNodesCount">The angle grid nodes count.</param>
/// <returns>MathNet.Numerics.LinearAlgebra.Double.DenseMatrix.</returns>
/// todo: регулировать размер сетки для свертки по гауссу
public static DenseMatrix CartesianToPolar(DenseMatrix dmData, PointD centerPoint, DenseMatrix dmMask, int angleGridNodesCount = 144)
{
double angleValueDiff = 2 * Math.PI / (angleGridNodesCount - 1);
DenseMatrix dmDistances = DenseMatrix.Create(dmData.RowCount, dmData.ColumnCount,
(r, c) => centerPoint.Distance(new PointD(c, r)));
dmDistances = (DenseMatrix)dmDistances.PointwiseMultiply(dmMask);
double maxDistance = dmDistances.Values.Maximum();
int distanceNodesCount = Convert.ToInt32(maxDistance) + 1;
DenseMatrix dmCountOfElementsSummedToThePolarPointValue = DenseMatrix.Create(angleGridNodesCount,
distanceNodesCount, 0.0d);
DenseMatrix dmAngles = DenseMatrix.Create(dmData.RowCount, dmData.ColumnCount,
(row, column) =>
{
double dx = (double)column - centerPoint.X;
double dy = (double)row - centerPoint.Y;
double r = dmDistances[row, column];
if (r == 0.0d)
{
return(0.0d);
}
double cosPhi = dx / r;
double phi = Math.Acos(cosPhi);
if (dy > 0)
{
phi = 2.0d * Math.PI - phi;
}
return(phi);
});
dmAngles = (DenseMatrix)dmAngles.PointwiseMultiply(dmMask);
DenseMatrix dmOutData = DenseMatrix.Create(angleGridNodesCount, distanceNodesCount, 0.0d);
//Vector2D vCenterPoint = new Vector2D(centerPoint);
//DenseMatrix dmSmoothed = dmData.Conv2(StandardConvolutionKernels.gauss, 7);
//int cartesianRowCount = dmSmoothed.RowCount;
//int cartesianColCount = dmSmoothed.ColumnCount;
//DenseMatrix dmOutData = DenseMatrix.Create(angleGridNodesCount, distanceNodesCount, (row, col) =>
//{
// PointPolar ptPolCurrPoint = new PointPolar((double) col, (double) row*angleValueDiff);
// PointD ptdCurrPoint = (ptPolCurrPoint.PointD() + vCenterPoint).ToPointD();
// int cartesianRow = Convert.ToInt32(ptdCurrPoint.Y);
// int cartesianCol = Convert.ToInt32(ptdCurrPoint.X);
// if ((cartesianRow < 0) || (cartesianRow >= cartesianRowCount) || (cartesianCol < 0) ||
// (cartesianCol >= cartesianColCount))
// {
// return 0.0d;
// }
// else
// return dmSmoothed[cartesianRow, cartesianCol];
//});
for (int cartesianRow = 0; cartesianRow < dmData.RowCount; cartesianRow++)
{
for (int cartesianCol = 0; cartesianCol < dmData.ColumnCount; cartesianCol++)
{
double currAngle = dmAngles[cartesianRow, cartesianCol];
double currDistance = dmDistances[cartesianRow, cartesianCol];
if (dmMask[cartesianRow, cartesianCol] == 0.0d)
{
continue;
}
int angleRow = Convert.ToInt32(currAngle / angleValueDiff);
int distanceCol = Convert.ToInt32(currDistance);
dmOutData[angleRow, distanceCol] += dmData[cartesianRow, cartesianCol];
if (dmMask[cartesianRow, cartesianCol] > 0.0d)
{
dmCountOfElementsSummedToThePolarPointValue[angleRow, distanceCol] += 1.0d;
}
}
}
dmCountOfElementsSummedToThePolarPointValue.MapInplace(x => (x == 0.0d) ? (0.0d) : (1.0d / x));
dmOutData = (DenseMatrix)dmOutData.PointwiseMultiply(dmCountOfElementsSummedToThePolarPointValue);
dmOutData.MapInplace(x => (x == 0.0d) ? (1.0d) : (x));
//сгладим
foreach (Tuple <int, MathNet.Numerics.LinearAlgebra.Vector <double> > tplRowVectorAndIndex in dmOutData.EnumerateRowsIndexed())
{
DenseVector currVector = (DenseVector)tplRowVectorAndIndex.Item2;
DenseVector currVectorSmoothed = currVector.Conv(StandardConvolutionKernels.gauss, 5);
// и восстановим краевые значения
for (int i = 0; i < 5; i++)
{
currVectorSmoothed[i] = currVector[i];
currVectorSmoothed[currVectorSmoothed.Count - i - 1] = currVector[currVector.Count - i - 1];
}
dmOutData.SetRow(tplRowVectorAndIndex.Item1, currVectorSmoothed.ToArray());
}
return(dmOutData);
}
