打开医学图像,使用Radon生成雷登变换图(正弦图),并使用iradon进行重建
时间: 2023-06-14 14:05:09
浏览: 288
好的,这里提供一个Python的示例代码,需要安装NumPy和matplotlib库:
```python
import numpy as np
import matplotlib.pyplot as plt
from skimage.transform import radon, iradon
# 生成一个圆形模板
N = 128
R = 40
x, y = np.indices((N, N))
mask = ((x-N//2)**2 + (y-N//2)**2) < R**2
# 生成雷登变换图
theta = np.linspace(0, 180, 180, endpoint=False)
sinogram = radon(mask.astype(float), theta=theta, circle=True)
# 显示正弦图
fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(8, 4.5))
ax1.imshow(mask, cmap='gray')
ax1.set_title('Original')
ax2.imshow(sinogram, cmap='gray', extent=(-90, 90, 0, sinogram.shape[0]), aspect='auto')
ax2.set_title('Sinogram')
plt.show()
# 重建图像
reconstruction = iradon(sinogram, theta=theta, circle=True)
reconstruction *=
```
最低
0.47元/天
开通会员,查看完整答案

成为会员后, 你将解锁


相关推荐
















