private void btnProcessAccelerationTimeSeries_Click(object sender, EventArgs e)
{
if ((bgwCalculate != null) && (bgwCalculate.IsBusy))
{
bgwCalculate.CancelAsync();
return;
}
//simpleMultipleImagesShow imagesRepresentingForm = new simpleMultipleImagesShow();
//imagesRepresentingForm.Show();
DoWorkEventHandler bgwCalculate_DoWorkHandler = delegate(object currBGWsender, DoWorkEventArgs args)
{
BackgroundWorker selfWorker = currBGWsender as BackgroundWorker;
//simpleMultipleImagesShow multImagesRepresentingForm = (simpleMultipleImagesShow)((args.Argument as object[])[0]);
//Type theShowImagesType = multImagesRepresentingForm.GetType();
//MethodInfo thePicturePlacingMethodInfo = theShowImagesType.GetMethod("PlaceAPicture");
int imageRepresentingCounter = 0;
DateTime dbgDT = new DateTime(2014, 7, 9, 9, 0, 0);
dbgDT = dbgDT.AddMinutes(33);
DateTime dtSeriesStart = accSubseries[0].StartTime;
double tsOverallSeriesDurationMillisec = (accSubseries[accSubseries.Count - 1].EndTime - accSubseries[0].StartTime).TotalMilliseconds;
string strToWrite = " fileName ; lat ; lon ; date ; time ; time(s) since start ; period(s) ; spectrum amplitude";
ServiceTools.logToTextFile(strOutputDirectory + "\\100sData-spectra-maximums.dat", strToWrite + Environment.NewLine, true);
foreach (TimeSeries <double> accSubseria in accSubseries)
{
if (selfWorker.CancellationPending)
{
break;
}
int startindex = 0;
while (true)
{
int endIndex;
if (selfWorker.CancellationPending)
{
break;
}
TimeSeries <double> currTimeSeria = accSubseria.SubSeria(startindex, new TimeSpan(1000000000),
out endIndex); //100s
currTimeSeria = currTimeSeria.InterpolateSeria(new TimeSpan(500000));
currTimeSeria = currTimeSeria.ExtractDataDeviationValues();
//обработать и оценить наличие выраженных периодов
Complex[] sourceSignalArray = currTimeSeria.DataRealValuesComplexArray();
Fourier.Forward(sourceSignalArray);
//Transform.FourierForward(sourceSignalArray);
List <Complex> FourierTransformedSignal = new List <Complex>(sourceSignalArray);
List <double> FourierTransformedSignalAmplitudes =
FourierTransformedSignal.ConvertAll <double>(
cVal =>
((double.IsNaN(cVal.Magnitude)) || (double.IsInfinity(cVal.Magnitude)))
? (0.0d)
: (cVal.Magnitude));
List <double> FourierTransformedSignalPeriods = new List <double>();
for (int ind = 0; ind < FourierTransformedSignalAmplitudes.Count; ind++)
{
FourierTransformedSignalPeriods.Add(currTimeSeria.TotalSeriaDuration.TotalSeconds /
(double)ind);
}
FourierTransformedSignalAmplitudes =
new List <double>(
FourierTransformedSignalAmplitudes.Zip <double, double, double>(
FourierTransformedSignalPeriods,
(amp, periodSec) =>
((double.IsNaN(periodSec)) || (double.IsInfinity(periodSec))) ? (0.0d) : (amp)));
FourierTransformedSignalPeriods =
FourierTransformedSignalPeriods.ConvertAll <double>(
dval => ((double.IsNaN(dval)) || (double.IsInfinity(dval))) ? (0.0d) : (dval));
//проанализируем этот участок - есть ли выраженные пики по амплитуде конкретных частот
//найти максимум в спектре и выдать данные об этом максимуме в файл
// сначала отфильтруем периоды меньше 1с - для данных по динамике судна они несущественны
FourierTransformedSignalAmplitudes =
new List <double>(
FourierTransformedSignalAmplitudes.Zip <double, double, double>(
FourierTransformedSignalPeriods,
(amp, periodSec) => (periodSec <= 1.0d) ? (0.0d) : (amp)));
FourierTransformedSignalPeriods =
FourierTransformedSignalPeriods.ConvertAll <double>(dVal => (dVal <= 1.0d) ? (0.0d) : (dVal));
DescriptiveStatistics currAmpsStat =
new DescriptiveStatistics(FourierTransformedSignalAmplitudes);
List <double> lAmpsOutstanding = FourierTransformedSignalAmplitudes.ConvertAll <double>(dVal =>
{
if (dVal / currAmpsStat.Mean >= 100.0d)
{
return(dVal);
}
return(0.0d);
});
List <double> lPeriodsOutstanding =
new List <double>(FourierTransformedSignalPeriods.Zip <double, double, double>(lAmpsOutstanding,
(per, amp) => (amp == 0.0d) ? (0.0d) : (per)));
if (lAmpsOutstanding.Sum() > 0.0d)
{
MultipleScatterAndFunctionsRepresentation renderer =
new MultipleScatterAndFunctionsRepresentation(2048, 1536);
renderer.dvScatterXSpace.Add(currTimeSeria.TimeStampsValuesSeconds);
renderer.dvScatterFuncValues.Add(currTimeSeria.dvDoubleDataValues);
renderer.scatterLineColors.Add(new Bgr(255, 50, 50));
renderer.scatterDrawingVariants.Add(SequencesDrawingVariants.polyline);
renderer.dvScatterXSpace.Add(DenseVector.OfEnumerable(FourierTransformedSignalPeriods));
renderer.dvScatterFuncValues.Add(DenseVector.OfEnumerable(FourierTransformedSignalAmplitudes));
renderer.scatterLineColors.Add(new Bgr(50, 255, 50));
renderer.scatterDrawingVariants.Add(SequencesDrawingVariants.squares);
renderer.dvScatterXSpace.Add(DenseVector.OfEnumerable(lPeriodsOutstanding));
renderer.dvScatterFuncValues.Add(DenseVector.OfEnumerable(lAmpsOutstanding));
renderer.scatterLineColors.Add(new Bgr(50, 50, 255));
renderer.scatterDrawingVariants.Add(SequencesDrawingVariants.circles);
renderer.Represent();
if (strOutputDirectory != "")
{
double maxAmp = lAmpsOutstanding.Max();
int idx = lAmpsOutstanding.FindIndex(dval => dval == maxAmp);
double maxAmpPeriod = lPeriodsOutstanding[idx];
GPSdata gpsMark = gpsSeriaData.GetMostClose(currTimeSeria.StartTime).Item2;
string fName = currTimeSeria.StartTime.ToString("s").Replace(":", "-") +
"-100sData-spectrum.jpg";
renderer.SaveToImage(strOutputDirectory + "\\" + fName, true);
strToWrite = "" + fName + " ; ";
strToWrite += gpsMark.LatDec + " ; ";
strToWrite += gpsMark.LonDec + " ; ";
strToWrite += currTimeSeria.StartTime.Date.ToString("yyyy-MM-dd") + " ; ";
strToWrite += currTimeSeria.StartTime.ToString("HH-mm-ss") + " ; ";
strToWrite += (currTimeSeria.StartTime - dtSeriesStart).TotalSeconds + " ; ";
strToWrite += maxAmpPeriod.ToString() + " ; ";
strToWrite += maxAmp.ToString() + " ; ";
ServiceTools.logToTextFile(strOutputDirectory + "\\100sData-spectra-maximums.dat", strToWrite + Environment.NewLine, true);
ThreadSafeOperations.SetText(lblStatusString, "processing: " + currTimeSeria.StartTime.ToString("s"), false);
}
}
if ((currTimeSeria.StartTime >= dbgDT) || (currTimeSeria.EndTime >= dbgDT))
{
startindex++;
startindex--;
}
if (endIndex == accSubseria.Count - 1)
{
break;
}
Application.DoEvents();
selfWorker.ReportProgress(Convert.ToInt32(100.0d * (currTimeSeria.EndTime - dtSeriesStart).TotalMilliseconds / tsOverallSeriesDurationMillisec));
startindex += Convert.ToInt32((endIndex - startindex) / 2.0d);
}
}
};
RunWorkerCompletedEventHandler bgwCalculate_CompletedHandler = delegate(object currBGWCompletedSender, RunWorkerCompletedEventArgs args)
{
ThreadSafeOperations.ToggleButtonState(btnProcessAccelerationTimeSeries, true, "Process acceleration timeseries", false);
};
ProgressChangedEventHandler bgwCalculate_ProgressChanged = delegate(object bgwDataReaderSender, ProgressChangedEventArgs args)
{
ThreadSafeOperations.UpdateProgressBar(prbReadingProcessingData, args.ProgressPercentage);
};
ThreadSafeOperations.ToggleButtonState(btnProcessAccelerationTimeSeries, true, "STOP", true);
bgwCalculate = new BackgroundWorker();
bgwCalculate.WorkerSupportsCancellation = true;
bgwCalculate.WorkerReportsProgress = true;
bgwCalculate.DoWork += bgwCalculate_DoWorkHandler;
bgwCalculate.RunWorkerCompleted += bgwCalculate_CompletedHandler;
bgwCalculate.ProgressChanged += bgwCalculate_ProgressChanged;
//object[] bgwCalculateArgs = new object[] { imagesRepresentingForm };
object[] bgwCalculateArgs = new object[] { };
bgwCalculate.RunWorkerAsync(bgwCalculateArgs);
}
