Tag peaks that are touching the detector mask
jonwright opened this issue · comments
For old data, take a detector mask, grow it and cast to boolean, and just do:
cf.filter( mask[ cf.s_raw.astype(int), cf.f_raw.astype(int) ] )
Perhaps this could be added as a columnfile method, something like columnfile.detectormask( mask, growth_cycles=3 )
?
For new data, we could look at the pixel intensity which is adjacent to a masked pixel. This means counting up the integrated intensity twice. Places to look at would include:
- sparse properties props(): ImageD11/sinograms/properties.py
- labelimage.measurepeaks. Also peaksearch.py would need to read a mask from the command line.
- new frelon frelon_peaksearch