Main Page | Data Structures | File List | Data Fields | Globals

peaks.c File Reference

contains function for finding guassian-like peaks in a vector of intensity More...

#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "peaks.h"

Functions

int specFindPeak (INT32 *spec, int length, float *pcentroid, float *psize, int loch, int hich, int ifwhm, float sigma, int maxpk)
 find all the peaks in an intensity profile.


Detailed Description

contains function for finding guassian-like peaks in a vector of intensity


Function Documentation

int specFindPeak INT32 *  spec,
int  length,
float *  pcentroid,
float *  psize,
int  loch,
int  hich,
int  ifwhm,
float  sigma,
int  maxpk
 

find all the peaks in an intensity profile.

the algorithm looks for gaussian-like peaks by looking for local maxima in a fit score. the fit score is the difference of the integrated intensity within 1 fwhm and the integrated intensity between 1 and 2 fwhm away from the point, normalizded by the integrated intensity within 2 fwhm of the point.

Parameters:
spec intensity trace in which peaks are found
length length of intensity trace
pcentroid array where centroid of located peaks are written
psize array where the area of the peaks are written
loch lower limit of the index for the search
hich upper limit of the index for the search
ifwhm integer full width at half maximum of expected peak
sigma minimum "sigma" (fit score) value a point must yield for the algorithm to look for a peak
maxpk maximum number of peak to be located
Returns:
the number of peaks found


Generated on Fri Apr 20 16:43:07 2007 for LaneRuler by  doxygen 1.3.9.1