bug in radon transform?
Emmanuelle Gouillart
emmanuelle.gouillart at nsup.org
Wed Apr 4 17:31:49 EDT 2012
Hello,
I have the impression that the center of rotation used for
computing the radon transform (skimage.transform.radon) is not correct:
>>> import numpy as np
>>> from skimage.transform import radon, iradon
>>>
>>> a = np.zeros((3, 3))
>>> a[1, 1] = 1
>>> projs = radon(a, [0, 45, 90, 135])
>>> projs
array([[ 0. , 0. , 0. , 0. ],
[ 0. , 0. , 0. , 0.13603897],
[ 0. , 0.46446609, 1. , 0.79289322],
[ 1. , 0.46446609, 0. , 0. ],
[ 0. , 0. , 0. , 0. ],
[ 0. , 0. , 0. , 0. ],
[ 0. , 0. , 0. , 0. ]])
(the projection at 90�C, that is projs[2], should have the non-zero
coefficient centered).
Before I try to modify this, I'd like to be sure that it is not the
expected behaviour, and a bug indeed...
Cheers,
Emmanuelle
More information about the scikit-image
mailing list