如何为4pi反投影实现tikhonov滤波器?

2024-03-29 10:26:33 发布

您现在位置:Python中文网/ 问答频道 /正文

我试图在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成像仪从辐射数据进行的简单反投影。问题是过滤器根本没有改变简单的反投影图像


Tags: 图像imagefft过滤器npnxshapexx