wassname/scikit-image
1
0
Fork 0
mirror of https://github.com/wassname/scikit-image.git synced 2026-07-08 03:46:28 +08:00
Code Issues Packages Projects Releases Wiki Activity

transform.iradon: Refactoring for shorter functions.

Browse Source 
Will facilitate testing.
This commit is contained in:
Jostein Bø Fløystad
2013-06-19 00:12:14 +02:00
parent da423931b5
commit b61ff7513e
1 changed files with 9 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files
skimage/transform/radon_transform.py
+9 -6
View File
@@ -113,6 +113,14 @@ def radon(image, theta=None, circle=False):
return out
def _sinogram_circle_to_square(sinogram):
size = int(np.ceil(np.sqrt(2) * sinogram.shape[0]))
sinogram_padded = np.zeros((size, sinogram.shape[1]))
pad = int(np.ceil((size - sinogram.shape[0] - 1) / 2))
sinogram_padded[pad:pad + sinogram.shape[0], :] = sinogram
return sinogram_padded
def iradon(radon_image, theta=None, output_size=None,
filter="ramp", interpolation="linear", circle=False):
"""
@@ -177,12 +185,7 @@ def iradon(radon_image, theta=None, output_size=None,
output_size = int(np.floor(np.sqrt((radon_image.shape[0])**2
/ 2.0)))
if circle:
radon_size = int(np.ceil(np.sqrt(2) * radon_image.shape[0]))
radon_image_padded = np.zeros((radon_size, radon_image.shape[1]))
radon_pad = int(np.floor((radon_size - radon_image.shape[0] - 1) / 2.))
radon_image_padded[radon_pad:radon_pad + radon_image.shape[0], :] \
= radon_image
radon_image = radon_image_padded
radon_image = _sinogram_circle_to_square(radon_image)
n = radon_image.shape[0]