コード例 #1
0
 public MZData[] PeakDetect(MZData[] Data)
 {
     PeakDetecting.PeakDetector pd    = new PeakDetecting.PeakDetector();
     PeakDetecting.peakinside[] Peaks = new PeakDetecting.peakinside[1];
     pd.PeaksDetecting(ref Data, ref Peaks);
     MZData[] OutData = new MZData[Peaks.GetLength(0)];
     for (int i = 0; i < Peaks.GetLength(0); i++)
     {
         OutData[i].Intensity = Peaks[i].Value;
         OutData[i].Mass      = Peaks[i].Center;
     }
     return(OutData);
 }
コード例 #2
0
        public MZData[] Centroid(MZData[] Data,int Len, bool StickMode /* former "in" */)
        {
            int total = 0, u;
            int o = 0, i = 0, count = Len;
            double sumIi, sumI, last = 0.0;
            double du = 0.0;
            bool goingdown = false;
            MZData[] OutData;

            if (StickMode) {
                //считаем пока не начнутся нули или пока следующий не станет меньше помассе
                for ( i = 1 ; i<count ; i++){
                    if (Data[i].Mass < Data[i-1].Mass || Data[i].Mass == 0){
                        break;
                    }
                }
                OutData = new MZData[i];
                count = i;
                for (i=0; i<count ; i++){
                    OutData[i].Intensity = Data[i].Intensity;
                    OutData[i].Mass = Data[i].Mass;
                }
                return OutData;
            }

            //пропуск начальных нулей
            while(i < count && Data[i].Intensity == 0.0) ++i;

            //считает области больше нуля
            while(i < count)
            {
                while(i < count && Data[i].Intensity != 0.0)
                {
                    if(last > Data[i].Intensity) {
                        goingdown = true;
                    }else{
                        if(goingdown) {
                            ++total;
                            goingdown = false;
                        }
                    }

                    last = Data[i].Intensity;
                    ++i;
                }

                last = 0.0;
                goingdown = false;

                while(i < count && Data[i].Intensity == 0.0)
                    i++;

                total++;
            }

            //запасает память на подсчитанные области
            OutData = new MZData[total];
            i = 0; o = 0; total = 0; last = 0.0; goingdown = false;

            while(i < count && Data[i].Intensity == 0.0) i++;

            while(i < count)
            {
                sumIi = sumI = 0.0;
                o = i -1;
                while(i < count && Data[i].Intensity != 0.0){

                    //если пошло на спад
                    if(last > Data[i].Intensity) {
                        goingdown = true;
                    }else{
                        if(goingdown) {
                            u = Convert.ToInt32((sumIi / sumI)/* + 0.5*/);
                            OutData[total].Intensity = sumI;
                            OutData[total].Mass = Data[o+u].Mass;
                            ++total;

                            sumIi = sumI = 0.0;
                            o = i -1;
                            goingdown = false;
                        }
                    }

                    sumIi += Data[i].Intensity*(i-o);
                    sumI += Data[i].Intensity;

                    last = Data[i].Intensity;
                    i++;
                }

                u = Convert.ToInt32((sumIi / sumI) /*+0.5*/ );
                du = sumIi / sumI - (double)u;
                //интенсивность по интегралу
                OutData[total].Intensity = sumI;
                //сентроид - по апексу
                //OutData[total].Mass = Data[o+u].Mass;
                //центроид по центру
                OutData[total].Mass = Data[o+u].Mass*(1-du) + Data[o+u+1].Mass*du;

                last = 0.0;
                goingdown = false;

                while(i < count && Data[i].Intensity == 0.0)
                    i++;
                total++;
            }
            return OutData;
        }
コード例 #3
0
 public MZData[] PeakDetect(MZData[] Data )
 {
     PeakDetecting.PeakDetector pd = new PeakDetecting.PeakDetector();
     PeakDetecting.peakinside[] Peaks = new PeakDetecting.peakinside[1];
     pd.PeaksDetecting(ref Data, ref Peaks);
     MZData[] OutData = new MZData[Peaks.GetLength(0)];
     for (int i = 0 ; i < Peaks.GetLength(0) ; i++){
         OutData[i].Intensity = Peaks[i].Value;
         OutData[i].Mass = Peaks[i].Center;
     }
     return OutData;
 }
コード例 #4
0
        public MZData[] Centroid(MZData[] Data, int Len, bool StickMode /* former "in" */)
        {
            int    total = 0, u;
            int    o = 0, i = 0, count = Len;
            double sumIi, sumI, last = 0.0;
            double du        = 0.0;
            bool   goingdown = false;

            MZData[] OutData;

            if (StickMode)
            {
                //считаем пока не начнутся нули или пока следующий не станет меньше помассе
                for (i = 1; i < count; i++)
                {
                    if (Data[i].Mass < Data[i - 1].Mass || Data[i].Mass == 0)
                    {
                        break;
                    }
                }
                OutData = new MZData[i];
                count   = i;
                for (i = 0; i < count; i++)
                {
                    OutData[i].Intensity = Data[i].Intensity;
                    OutData[i].Mass      = Data[i].Mass;
                }
                return(OutData);
            }


            //пропуск начальных нулей
            while (i < count && Data[i].Intensity == 0.0)
            {
                ++i;
            }

            //считает области больше нуля
            while (i < count)
            {
                while (i < count && Data[i].Intensity != 0.0)
                {
                    if (last > Data[i].Intensity)
                    {
                        goingdown = true;
                    }
                    else
                    {
                        if (goingdown)
                        {
                            ++total;
                            goingdown = false;
                        }
                    }

                    last = Data[i].Intensity;
                    ++i;
                }

                last      = 0.0;
                goingdown = false;

                while (i < count && Data[i].Intensity == 0.0)
                {
                    i++;
                }

                total++;
            }

            //запасает память на подсчитанные области
            OutData = new MZData[total];
            i       = 0; o = 0; total = 0; last = 0.0; goingdown = false;

            while (i < count && Data[i].Intensity == 0.0)
            {
                i++;
            }

            while (i < count)
            {
                sumIi = sumI = 0.0;
                o     = i - 1;
                while (i < count && Data[i].Intensity != 0.0)
                {
                    //если пошло на спад
                    if (last > Data[i].Intensity)
                    {
                        goingdown = true;
                    }
                    else
                    {
                        if (goingdown)
                        {
                            u = Convert.ToInt32((sumIi / sumI) /* + 0.5*/);
                            OutData[total].Intensity = sumI;
                            OutData[total].Mass      = Data[o + u].Mass;
                            ++total;

                            sumIi     = sumI = 0.0;
                            o         = i - 1;
                            goingdown = false;
                        }
                    }

                    sumIi += Data[i].Intensity * (i - o);
                    sumI  += Data[i].Intensity;

                    last = Data[i].Intensity;
                    i++;
                }

                u  = Convert.ToInt32((sumIi / sumI) /*+0.5*/);
                du = sumIi / sumI - (double)u;
                //интенсивность по интегралу
                OutData[total].Intensity = sumI;
                //сентроид - по апексу
                //OutData[total].Mass = Data[o+u].Mass;
                //центроид по центру
                OutData[total].Mass = Data[o + u].Mass * (1 - du) + Data[o + u + 1].Mass * du;

                last      = 0.0;
                goingdown = false;

                while (i < count && Data[i].Intensity == 0.0)
                {
                    i++;
                }
                total++;
            }
            return(OutData);
        }