diff options
| author | Picca Frédéric-Emmanuel <picca@synchrotron-soleil.fr> | 2016-05-02 17:18:20 +0200 |
|---|---|---|
| committer | Picca Frédéric-Emmanuel <picca@synchrotron-soleil.fr> | 2016-05-02 17:18:20 +0200 |
| commit | ddc148d462ee43384ea1f60f3e017211ace7dc5a (patch) | |
| tree | 53cb134fd56e89b4a3ac4c58de1b848ec017bf2d | |
| parent | 02ce8a2704110b759cd3a7f23483a9bee4dcf2fd (diff) | |
attenuation correction
| -rw-r--r-- | binoculars/backends/sixs.py | 15 |
1 files changed, 10 insertions, 5 deletions
diff --git a/binoculars/backends/sixs.py b/binoculars/backends/sixs.py index 218cba2..f573f57 100644 --- a/binoculars/backends/sixs.py +++ b/binoculars/backends/sixs.py @@ -399,6 +399,7 @@ class FlyScanUHV(SIXS): "omega": "UHV_OMEGA", "delta": "UHV_DELTA", "gamma": "UHV_GAMMA", + "attenuation": "attenuation", } def get_pointcount(self, scanno): @@ -412,9 +413,9 @@ class FlyScanUHV(SIXS): omega = h5_nodes['omega'][index] delta = h5_nodes['delta'][index] gamma = h5_nodes['gamma'][index] + attenuation = h5_nodes['attenuation'][index] - return (image, - (mu, omega, delta, gamma)) + return (image, attenuation, (mu, omega, delta, gamma)) def process_image(self, index, dataframe, pixels): util.status(str(index)) @@ -428,7 +429,9 @@ class FlyScanUHV(SIXS): # extract the data from the h5 nodes h5_nodes = dataframe.h5_nodes - intensity, values = self.get_values(index, h5_nodes) + intensity, attenuation, values = self.get_values(index, h5_nodes) + + intensity *= 3.2 ** attenuation weights = numpy.ones_like(intensity) weights *= ~mask @@ -477,6 +480,7 @@ class FlyScanUHV2(FlyScanUHV): "omega": "omega", "delta": "delta", "gamma": "gamma", + "attenuation": "attenuation", } @@ -487,6 +491,7 @@ class SBSMedH(FlyScanUHV): "mu": "data_18", "gamma": "data_20", "delta": "data_19", + "attenuation": "data_xx", } def get_pointcount(self, scanno): @@ -500,9 +505,9 @@ class SBSMedH(FlyScanUHV): mu = h5_nodes['mu'][index] gamma = h5_nodes['gamma'][index] delta = h5_nodes['delta'][index] + attenuation = h5_nodes['attenuation'][index] - return (image, - (pitch, mu, gamma, delta)) + return (image, attenuation, (pitch, mu, gamma, delta)) def load_matrix(filename): |