反Rad变换

Question

谁能帮助我实现不带过滤器的iradon变换（反投影）？？？ 我有点知道这个概念，但是很难实现。

这是我的ra变换（正投影）代码

def sinogram(img):    
    radon = np.zeros((len(img), len(img)), dtype='float64')

    for i in range(len(img)):
        rot = misc.imrotate(img, -i * 180 / len(img)).astype('float64')
        radon[:, i] = sum(rot) # Sum is the integral in Radon transform   

    return radon