Page History
Table of Contents |
---|
Note |
---|
Images are shown for DIFFERENT EVENTS for demonstration purpose only |
...
Algorithm
Chuck suggested merging algorithm with sequence of procedures for photons split between pixels as follows.
Original array
This algorithm consumes (float) calibrated data expressed in terms of number of photons. For test purpose use image shape=(50,50) random standard normal distribution with sigma=10
...
Input array with float number of photons per pixel is split for two
- floor of input (
...
- uint) and
- difference between input and previuos uint array - (float) fractional array [0,1).
Pixels with intensity below zero and masked are set to 0 (photon) in both arrays.
Find local maximums
Fractional float array is used to search for local maximums in the rank=1 cross region (among 4 pixels closest to the central):
...
- C++ algorithm is implemented in ImgAlgos::AlgImgPros::mapOfPhotonNumbersV1 for 2-d arrays only
- Boost-python wrapper is in ImgAlgos.pyImgAlgos.cpp
- Pure python wrapper ImgAlgos.PyAlgos.py - supports 2-d and n-d arrays.
Presentation in the Detector package
...