private static List <double> PreProcess(List <double[]> signalsList, double Q, int m, double r, int N, int iterations, int locality) { double[] f3 = new double[signalsList.Count]; bool first = true; for (int i = 0; i < f3.Length; i++) { if (first) { //f3[i] = BetaBandpass(signalsList[i][0], true); f3[i] = signalsList[i][0]; first = false; } else { //f3[i] = BetaBandpass(signalsList[i][0], false); f3[i] = signalsList[i][0]; } } string internalDirectory = _dataAccessFacade.GetGeneralSettings().GetModalDirectory("Emotion"); var file_F3 = File.AppendText(Path.Combine(internalDirectory, "F3.json")); file_F3.WriteLine(f3.ToJsonString()); file_F3.Flush(); file_F3.Close(); //f3 = FilterRLC.LCHP(f3, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 5, Q); //f3 = FilterRLC.LCLP(f3, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 50, Q); for (int i = 0; i < f3.Length; i++) { if (double.IsInfinity(f3[i]) || double.IsNaN(f3[i])) { Console.WriteLine(i + " f3 es " + f3[i]); f3[i] = f3[i - 1] * (1 + (_xrand.Next(-1, 1) / 10)); } } var file_F3Filtered = File.AppendText(Path.Combine(internalDirectory, "F3Filtered.json")); file_F3Filtered.WriteLine(f3.ToJsonString()); file_F3Filtered.Flush(); file_F3Filtered.Close(); double[] c4 = new double[signalsList.Count]; first = true; for (int i = 0; i < c4.Length; i++) { if (first) { //c4[i] = BetaBandpass(signalsList[i][1], true); c4[i] = signalsList[i][1]; first = false; } else { //c4[i] = BetaBandpass(signalsList[i][1], false); c4[i] = signalsList[i][1]; } } var file_C4 = File.AppendText(Path.Combine(internalDirectory, "C4.json")); file_C4.WriteLine(c4.ToJsonString()); file_C4.Flush(); file_C4.Close(); //c4 = FilterRLC.LCHP(c4, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 5, Q); //c4 = FilterRLC.LCLP(c4, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 50, Q); for (int i = 0; i < c4.Length; i++) { if (double.IsInfinity(c4[i]) || double.IsNaN(c4[i])) { Console.WriteLine(i + " c4 es " + c4[i]); c4[i] = c4[i - 1] * (1 + (_xrand.Next(-1, 1) / 10)); } } var file_C4Filtered = File.AppendText(Path.Combine(internalDirectory, "C4Filtered.json")); file_C4Filtered.WriteLine(c4.ToJsonString()); file_C4Filtered.Flush(); file_C4Filtered.Close(); var emdF3 = new Emd(); var imfsF3 = emdF3.GetImfs(f3, 4, iterations, locality); Console.WriteLine("obtenido imfs F3"); var emdC4 = new Emd(); var imfsC4 = emdC4.GetImfs(c4, 4, iterations, locality); Console.WriteLine("obtenido imfs C4"); var file_imfsF3 = File.AppendText(Path.Combine(internalDirectory, "imfsF3.json")); file_imfsF3.WriteLine(imfsF3.ToJsonString()); file_imfsF3.Flush(); file_imfsF3.Close(); var file_imfsC4 = File.AppendText(Path.Combine(internalDirectory, "imfsC4.json")); file_imfsC4.WriteLine(imfsC4.ToJsonString()); file_imfsC4.Flush(); file_imfsC4.Close(); //return MH(imfsF3, imfsC4); return(CalcEntropy(imfsF3, imfsC4, N, m, r)); }
private List <double> PreProcess(List <double[]> signalsList, double Q, int m, double r, int N, int iterations, int locality) { double[] f3 = new double[signalsList.Count]; bool first = true; for (int i = 0; i < f3.Length; i++) { if (first) { f3[i] = BetaBandpass(signalsList[i][0], true); //f3[i] = signalsList[i][0]; first = false; } else { f3[i] = BetaBandpass(signalsList[i][0], false); //f3[i] = signalsList[i][0]; } } //f3 = FilterRLC.LCHP(f3, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 5, Q); //f3 = FilterRLC.LCLP(f3, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 50, Q); for (int i = 0; i < f3.Length; i++) { if (double.IsInfinity(f3[i]) || double.IsNaN(f3[i])) { Console.WriteLine("F3 es NaN o Infinityyyyyy"); f3[i] = f3[i - 1] * (1 + (_xrand.Next(-1, 1) / 10)); } } double[] c4 = new double[signalsList.Count]; first = true; for (int i = 0; i < c4.Length; i++) { if (first) { c4[i] = BetaBandpass(signalsList[i][1], true); //c4[i] = signalsList[i][1]; first = false; } else { c4[i] = BetaBandpass(signalsList[i][1], false); //c4[i] = signalsList[i][1]; } } //c4 = FilterRLC.LCHP(c4, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 5, Q); //c4 = FilterRLC.LCLP(c4, EEGEmoProc2ChSettings.Instance.SamplingHz.Value, 50, Q); for (int i = 0; i < c4.Length; i++) { if (double.IsInfinity(c4[i]) || double.IsNaN(c4[i])) { Console.WriteLine("C4 es NaN o Infinityyyyyy"); c4[i] = c4[i - 1] * (1 + (_xrand.Next(-1, 1) / 10)); } } var emdF3 = new Emd(); var imfsF3 = emdF3.GetImfs(f3, 4, iterations, locality); var emdC4 = new Emd(); var imfsC4 = emdC4.GetImfs(c4, 4, iterations, locality); //return MH(imfsF3, imfsC4); return(CalcEntropy(imfsF3, imfsC4, N, m, r)); }