反Radon变换

Question

有谁能帮我实现无过滤器的iradon变换（反投影）？？？ 我有点了解这个概念，但很难实现。在

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

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