我试图在4pi成像仪的滤波反投影中实现Tikhonov正则化,但我编写的代码根本没有改变最终图像。我使用lambda和k的任意值来尝试查看过滤器的工作方式,但到目前为止还没有骰子
在python中,我尝试过:
theta = np.linspace(latra[0], latra[1], min(image.shape)-1)#, endpoint=True)
print len(theta)
sinogram = radon(image, theta=theta, circle=False)
im_fft = fft2(sinogram)
L = 4.5e-2
k = L/5
nX = im_fft.shape[0]
kxMax = k*nX/2
nY = im_fft.shape[1]
kyMax = k*nY/2
x = np.arange(-kxMax,kxMax,nX)
y = np.arange(-kyMax,kyMax,nY)
xx, yy = np.meshgrid(x,y)
tkF = np.sqrt(xx**2 + yy**2)/(1 + L**4*np.sqrt(xx**2+yy**2)**4)
Fk = np.conv2d(im_fft,tkF)
result = ifft2(Fk)
其中image
是我使用4pi成像仪从辐射数据进行的简单反投影。问题是过滤器根本没有改变简单的反投影图像
目前没有回答
相关问题 更多 >
编程相关推荐