コード例 #1
0
ファイル: DataProcessing.cs プロジェクト: tcidy/TF_Tool
        //
        // Neuro Temp Data Processing:
        // Look for the time series of data, Starting after a line that starts with 'Scan,'
        // For each time series (in odd columns):
        //    start_val = mean(first 5 elements)
        // start_elem_ix = index of ANY channel that is 0.2 higher than the start_val
        // stop_elem_ix = start_elem_ix+interval
        // analysis_range = [stop_elem_ix-10:stop_elem_ix-1]         // This is 10 elements!
        // output_value = mean(analysis_range) - start_val

        public static void ProcessCRMTempData(string[] infiles, string outfile, double interval,
                                              bool doPlots = true)
        {
            using (StreamWriter tw = new StreamWriter(outfile))
            {
                string[] headers = new string[] {
                    "" /* line number */,
                    "Measurement Filename",
                    String.Join(",",
                                Enumerable.Range(1, 4).Select(x => "Probe " + x.ToString())
                                )
                };
                tw.WriteLine(String.Join(",", headers));
                int i = 0;

                foreach (string fn in infiles)
                {
                    StreamReader    sr       = new StreamReader(fn);
                    string          namepart = Path.GetFileName(fn);
                    string          line;
                    List <double>[] values = new List <double> [4];
                    for (int j = 0; j < 4; j++)
                    {
                        values[j] = new List <double>();
                    }
                    // Skip header line
                    while ((line = sr.ReadLine()) != null)
                    {
                        line = line.Trim();
                        if (line == "")
                        {
                            continue;
                        }
                        string[] parts = line.Split(new char[] { ',' }, StringSplitOptions.RemoveEmptyEntries);
                        if (parts.Length != 5)
                        {
                            throw new FormatException("Wrong number of element in data row of: " + fn);
                        }
                        for (int j = 1; j <= 4; j++) // skip first element
                        {
                            values[j - 1].Add(Double.Parse(parts[j]));
                        }
                    }
                    ;
                    if (values[0].Count == 0)
                    {
                        throw new FormatException("Data file is too short. No data samples" + fn);
                    }
                    var startvals = values.Select(
                        v => v.GetRange(0, 5).Mean()
                        ).ToArray();
                    int start_elem_ix = -1;
                    for (int j = 0; j < 4; j++)
                    {
                        int ix2 = values[j].FindIndex(x => x > startvals[j] + 0.2);

                        if (ix2 < values[j].Count)
                        {
                            if (ix2 >= 0 && (start_elem_ix == -1 || start_elem_ix > ix2) && values[j][ix2 + 1] > values[j][ix2] + 0.2)
                            {
                                start_elem_ix = ix2;
                            }
                        }
                    }
                    if (start_elem_ix < 6 || start_elem_ix == -1)
                    {
                        MessageBox.Show("Start element was found to be invalid [<=10 || not found], ix=" +
                                        start_elem_ix.ToString() + ", in file " +
                                        fn);
                        continue;
                    }
                    //throw new FormatException("Start element was found to be invalid [<=10 || not found], ix=" +
                    //    start_elem_ix.ToString() + ", in file " +
                    //    fn);
                    int end_elem_ix = start_elem_ix + ((int)interval);
                    var endvals     = values.Select(
                        v => v.GetRange(end_elem_ix - 10, 10).Mean()
                        ).ToArray();
                    var deltaT    = startvals.Zip(endvals, (x1, x2) => (x2 - x1));
                    var lineparts = new string[]
                    {
                        i.ToString(), // line number
                        namepart,     // filename
                        String.Join(",", deltaT.Select(x => x.ToString()))
                    };
                    tw.WriteLine(
                        String.Join(",", lineparts)
                        );
                    if (doPlots)
                    {
                        PlotTempSeriesForm ptsf = new PlotTempSeriesForm(
                            title: namepart, data: values, startvals: startvals,
                            startx2: start_elem_ix,
                            endx1: end_elem_ix - 10, endx2: end_elem_ix - 1
                            );
                        ptsf.Show();
                        i++;
                    }
                }
            }
        }
コード例 #2
0
ファイル: DataProcessing.cs プロジェクト: tcidy/TF_Tool
        // FIXME: The 0.2 constant threshold for RF on seems perhaps too low
        // FIXME: Assumption that the sample rate is exactly 1 Hz
        //
        // Neuro Temp Data Processing:
        // Look for the time series of data, Starting after a line that starts with 'Scan,'
        // For each time series (in odd columns):
        //    start_val = mean(first 30 elements)           [ or maybe it should be 30???]
        // start_elem_ix = index of first element of probe 1 that is its start_val+0.2
        // stop_elem_ix = start_elem_ix+interval
        // analysis_range = [stop_elem_ix-10:stop_elem_ix-1]         // This is 10 elements!
        // output_value = mean(analysis_range) - start_val

        public static void ProcessNeuroTempData(string[] infiles, string outfile, double interval,
                                                bool doPlots = true)
        {
            using (StreamWriter tw = new StreamWriter(outfile))
            {
                int i = 0;
                foreach (string fn in infiles)
                {
                    StreamReader sr       = new StreamReader(fn);
                    string       namepart = Path.GetFileName(fn);
                    string       line;
                    int          number_channels;
                    do
                    {
                        line = sr.ReadLine();
                    } while (line != null && !line.StartsWith("Total Channels", StringComparison.InvariantCultureIgnoreCase));
                    string[] line_break = line.Split(new char[] { ',' }, StringSplitOptions.RemoveEmptyEntries);
                    number_channels = Int32.Parse(line_break[1]);
                    List <double>[] values = new List <double> [number_channels];
                    for (int j = 0; j < number_channels; j++)
                    {
                        values[j] = new List <double>();
                    }
                    do
                    {
                        line = sr.ReadLine();
                    } while (line != null && !line.StartsWith("scan,", StringComparison.InvariantCultureIgnoreCase));

                    // Skip header line
                    if (line == null)
                    {
                        throw new FormatException("Data file is too short. No data samples" + fn);
                    }

                    while ((line = sr.ReadLine()) != null)
                    {
                        line = line.Trim();
                        if (line == "")
                        {
                            continue;
                        }
                        string[] parts = line.Split(new char[] { ',' }, StringSplitOptions.RemoveEmptyEntries);
                        if (parts.Length != 2 * (number_channels + 1))
                        {
                            throw new FormatException("Wrong number of element in data row of: " + fn);
                        }
                        for (int j = 1; j <= number_channels; j++) // skip first element
                        {
                            values[j - 1].Add(Double.Parse(parts[j * 2]));
                        }
                    }
                    ;
                    if (values[0].Count == 0)
                    {
                        throw new FormatException("Data file is too short. No data samples" + fn);
                    }
                    var startvals = values.Select(
                        v => v.GetRange(0, 30).Mean()
                        ).ToArray();

                    //int start_elem_ix = values[0].Skip(30).ToList().FindIndex(x => x > startvals[0] + 0.2) + 30; // Ignore the first 30 elements
                    int start_elem_ix = -1;

                    for (int j = 0; j < number_channels; j++)
                    {
                        int ix2 = values[j].FindIndex(x => x > startvals[j] + 0.3);
                        if (ix2 >= 30 && (start_elem_ix == -1 || start_elem_ix > ix2))
                        {
                            start_elem_ix = ix2;
                        }
                    }
                    if (start_elem_ix == -1)
                    {
                        MessageBox.Show("Start element was found to be invalid [<=10 || not found], ix=" +
                                        start_elem_ix.ToString() + ", in file " +
                                        fn);
                        continue;
                    }
                    //throw new FormatException("Start element was found to be invalid, ix=" +
                    //    start_elem_ix.ToString() + ", in file " +
                    //    fn);
                    int end_elem_ix = start_elem_ix + ((int)interval);
                    if (end_elem_ix >= values[0].Count)
                    {
                        MessageBox.Show("In file " + fn + ", the RF pulse start + interval position was past the end of the data");
                        continue;
                    }
                    //throw new FormatException("In file " + fn + ", the RF pulse start + interval position was past the end of the data");
                    var endvals = values.Select(
                        v => v.GetRange(end_elem_ix - 10, 10).Mean()
                        ).ToArray();
                    var deltaT    = startvals.Zip(endvals, (x1, x2) => (x2 - x1));
                    var lineparts = new string[]
                    {
                        i.ToString(), // line number
                        namepart,     // filename
                        String.Join(",", deltaT.Select(x => x.ToString()))
                    };
                    string[] headers = new string[] {
                        "" /* line number */,
                        "Filename",
                        String.Join(",",
                                    Enumerable.Range(1, number_channels).Select(x => "Probe " + x.ToString())
                                    )
                    };
                    if (i == 0)
                    {
                        tw.WriteLine(String.Join(",", headers));
                    }
                    tw.WriteLine(
                        String.Join(",", lineparts)
                        );
                    if (doPlots)
                    {
                        PlotTempSeriesForm ptsf = new PlotTempSeriesForm(
                            title: namepart, data: values, startvals: startvals,
                            startx2: start_elem_ix,
                            endx1: end_elem_ix - 10, endx2: end_elem_ix - 1
                            );
                        ptsf.Show();
                    }
                    i++;
                }
            }
        }