PET图像重建OSEM算法

hxq_晴天

于 2024-09-18 15:18:36 发布

阅读量102

点赞数

分类专栏： PET图像重建文章标签： python numpy 机器学习

本文链接：https://blog.csdn.net/weixin_51870331/article/details/142334348

版权


from skimage.io import imread
from skimage import data_dir
from skimage.transform import radon, rescale, iradon
from skimage.draw import disk
# from skimage.draw import  circle
from scipy.ndimage import gaussian_filter
import numpy as np
import matplotlib.pyplot as plt



def osem(sinogram,nsub=1,niter=1,sfwhm=1.333,obj=True):


    shape = sinogram.shape
    assert shape[0] == shape[1],"osem sino must have same size"
    theta = np.linspace(0., 180., shape[0], endpoint=False)
    recon = np.zeros(shape)
    rr, cc = disk((shape[0] / 2, shape[1] / 2), shape[0] / 2 - 1)
    recon[rr, cc] = 1

    # normalization matrix
    nview = len(theta)
    norm = np.ones(shape)
    wgts = []
    for sub in range(nsub):
        views = range(sub, nview, nsub)
        wgt = iradon(norm[:, views], theta=theta[views], filter_name=None, circle=True)
        wgts.append(wgt)

    #