/*
 * peakdetect.h
 */

#include <gsl/gsl_matrix.h>

extern gsl_matrix* createImageMatrix(ControlContext *ctx, int16_t *image);

extern void sigma_filter(gsl_matrix *m, double k);

extern void local_sigma_filter(gsl_matrix *m, int r, double k);

extern void apply_kernel(gsl_matrix *m, gsl_matrix *kernel);

extern gsl_matrix* generate_flat_kernel(int half);

extern gsl_matrix* floodfill(gsl_matrix *m, double threshold, int* count);

extern gsl_matrix* iterate(gsl_matrix *m, int* count);