hanyoseob / python-art

Algebraic Reconstruction Technique (ART)

ART

Title

Algebraic Reconstruction Technique (ART)

Abstract

The algebraic reconstruction technique (ART) is an iterative reconstruction technique used in computed tomography. It reconstructs an image from a series of angular projections (a sinogram). Gordon, Bender and Herman first showed its use in image reconstruction; whereas the method is known as Kaczmarz method in numerical linear algebra.

An advantage of ART over other reconstruction methods (such as filtered backprojection (FBP)) is that it is relatively easy to incorporate prior knowledge into the reconstruction process.

Definition

Results

python-art's Issues

ART

hello dear;
when I ran the code, I got this problem .

TypeError Traceback (most recent call last)
/var/folders/04/6zwtc3x93jz3v9tkgfbfsk7c0000gn/T/ipykernel_35174/3701552822.py in
5 niter = 2e2
6 bpos = True
----> 7 x_art = ART(A, AT, y, x0, mu, niter, bpos)

/var/folders/04/6zwtc3x93jz3v9tkgfbfsk7c0000gn/T/ipykernel_35174/2609942379.py in ART(A, AT, b, x, mu, niter, bpos)
27
28 def ART(A, AT, b, x, mu=1e0, niter=1e2, bpos=True):
---> 29 ATA = AT(A(np.ones_like(x)))
30
31 for i in range(int(niter)):

/var/folders/04/6zwtc3x93jz3v9tkgfbfsk7c0000gn/T/ipykernel_35174/3663427198.py in (y)
7
8 A = lambda x: radon(x, THETA, circle=False).astype(np.float32)
----> 9 AT = lambda y: iradon(y, THETA, circle=False, filter=None, output_size=N).astype(np.float32)/(np.pi/(2*len(THETA)))
10 AINV = lambda y: iradon(y, THETA, circle=False, output_size=N).astype(np.float32)

TypeError: iradon() got an unexpected keyword argument 'filter'