![](https://img-blog.csdnimg.cn/20201014180756925.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
CT重建
文章平均质量分 56
xtphysics
莱顿大学LUMC在读博士
展开
-
基于Operator Discretization Library (ODL)的CT图像正反投
持续更新…1. FBPimport odlimport matplotlib.pyplot as pltreco_space = odl.uniform_discr(min_pt=[-20, -20], max_pt=[20, 20], shape=[512, 512])geometry = odl.tomo.parallel_beam_geometry(reco_space)projector = odl.tomo.RayTransform(reco_space,geometry, impl=原创 2021-07-15 18:00:13 · 2890 阅读 · 10 评论 -
快速批量Radon变换的Python实现——使用PyTorch函数affine_grid和grid_sample
Radon变换的旋转过程Radon变换的旋转过程可以使用ndimage.rotate()对图像进行特定角度的旋转,但这个过程比较耗时,不利于投影角度较多或者需要对批量图像进行Radon变换的情况。而通过利用PyTorch中的仿射变换相关函数affine_grid()和grid_sample()也可以实现旋转的过程,可以在GPU上进行运算,Radon变换所需的时间大幅减少。快速批量Radon变换的Python实现from scipy import ndimageimport numpy as npi原创 2021-05-12 14:37:17 · 1895 阅读 · 2 评论 -
扇形束CT重建快速生成系统矩阵(system matrix)的Python实现
import numpy as npimport mathSDD = 962.9SOD = 570ReconCentre = [0,0]pixel_spacing = 0.6836channel_spacing = 0.5channel_num = int(1024)view_num = 720image_size = 512'''每个探测器单元的中心(每条Ray在探测器上的landpoint)到探测器中心的距离'''Ray2DetCentre = []for i in ran原创 2021-05-05 15:52:03 · 3039 阅读 · 13 评论 -
Python实现平行束滤波反投影——Inverse Radon Transformation
参考博文进行了平行束滤波反投影的修改,将时域滤波修改为频域滤波,重建后消除原博文中图像的竖条状伪影。https://blog.csdn.net/hsyxxyg/article/details/106433940频域平行束滤波反投影(反radon变换)产生频域滤波信号:nextpow2函数的实现可参见github:https://github.com/freenowill/Denoise-/blob/ba99babb685f6c80f07a23275c8c2127de601815/Denois原创 2021-04-21 17:39:49 · 2236 阅读 · 7 评论