Beispiel #1
0
        public bool Test(double start_wl, double end_wl, double max_slope, double min_slope, double max_error,
                         out bool noisyFit, out double slope)
        {
            noisyFit = false;

            //normalize
            var curve = _spectrum.SpectrumData;
            var x     = curve.Select(p => p.X).ToArray();
            var y     = curve.Select(p => p.Y).ToArray();
            var minY  = y.Min();
            var maxY  = y.Max();

            y     = y.Select(p => (p - minY) / (maxY - minY)).ToArray();
            curve = x.Zip(y, (xp, yp) => new System.Windows.Point(xp, yp)).ToList();

            var roi = curve.Where(p => p.X >= start_wl && p.X <= end_wl).ToList();

            x = roi.Select(p => p.X).ToArray();
            y = roi.Select(p => p.Y).ToArray();

            var    linear        = Fit.Line(x, y);
            double intercept     = linear.Item1;
            double s             = linear.Item2;
            var    standardError = GoodnessOfFit.PopulationStandardError(x.Select(d => s * d + intercept), y);

            if (Math.Round(standardError, 3) > max_error)
            {
                noisyFit = true;
            }

            slope = s;
            bool res = (!noisyFit) && (slope >= min_slope) && (slope <= max_slope);

            return(res);
        }