CT 重建的 MATLAB 代码,挺适合你这种喜欢钻研图像算法的人。代码结构清晰,功能分块做得还不错:从数据读取、预、重建、后到图像显示都安排得妥妥的。用的是Radon 变换,配上常见的FBP或者ART、MLEM这类迭代方法,跑起来还蛮稳定的。
MATLAB 的fft2、ifft2这些函数在里面用得挺多,代码里也贴心配了几个滤波器的实现,比如Hamming
和Shepp-Logan
。如果你对图像质量要求高,稍微改改滤波器参数就行。响应也快,适合做算法原型。
需要注意一点,循环结构比较密集,尤其是用ART那块,for
和while
控制的更新迭代逻辑最好理清楚,不然容易卡逻辑。你可以从程序-改过
文件入手,一步步跑一下每段逻辑,比较容易上手。
哦对了,图像显示部分也没忽视,imshow
、volshow
这些都集成好了。跑完能直接看重建结果,调试起来方便。如果你有自己的 CT 投影数据,直接替换文件路径就能用了。
如果你对滤波反投影还不太熟,推荐去看看CT 重建算法探索滤波反投影与直接反投影对比这篇文章;另外MATLAB Radon 变换图像重建助手也挺适合新手上手的,图解讲得清楚。
,如果你熟 MATLAB,想快速搭一套 CT 图像重建流程,这份资源蛮值得试试的。