Halftomo: change CoR rounding
Currently, when using half-tomography with a given CoR:
- Reconstructed slices have shape
2*c
wherec = 2*int(round(CoR))
- CoR is set to
int(round(CoR))
(ascor_halftomo
).
Both points are problematic. Point (2) entails a loss in accuracy when a subpixel CoR is needed (which is often the case).
Point (1) prevents from picking a predictable output size, see _get_end_xy()
in dataset_validator.py
.
Proposed:
- Always use
2*ceil(CoR)
for slice size - Reconstruct with a possibly-floating-point
cor_halftomo
(as should be!)