diff options
Diffstat (limited to 'libcrystfel')
-rw-r--r-- | libcrystfel/src/filters.c | 55 |
1 files changed, 0 insertions, 55 deletions
diff --git a/libcrystfel/src/filters.c b/libcrystfel/src/filters.c index fc8738c8..2ce782d9 100644 --- a/libcrystfel/src/filters.c +++ b/libcrystfel/src/filters.c @@ -42,61 +42,6 @@ #include "image.h" -static int compare_vals(const void *ap, const void *bp) -{ - const signed int a = *(signed int *)ap; - const signed int b = *(signed int *)bp; - - if ( a > b ) return 1; - if ( a < b ) return -1; - return 0; -} - - -static void clean_panel(struct image *image, int sx, int sy) -{ - int x, y; - const int s = sizeof(signed int); - - for ( x=0; x<512; x++ ) { - - signed int vals[128]; - double m; - - for ( y=0; y<128; y++ ) { - vals[y] = image->data[(x+sx)+(y+sy)*image->width]; - } - - qsort(&vals[0], 128, s, compare_vals); - - m = gsl_stats_int_median_from_sorted_data(vals, 1, 128); - - for ( y=0; y<128; y++ ) { - image->data[(x+sx)+(y+sy)*image->width] -= m; - } - - } -} - - -/* Pre-processing to make life easier */ -void filter_cm(struct image *image) -{ - int px, py; - - if ( (image->width != 1024) || (image->height != 1024) ) return; - - for ( px=0; px<2; px++ ) { - for ( py=0; py<8; py++ ) { - - clean_panel(image, 512*px, 128*py); - - } - } - -} - - void filter_noise(struct image *image) { int x, y; |