- 博客(6)
- 收藏
- 关注
RSS订阅原创 迭代重建算法中投影矩阵的计算
在前面学习的重建算法都是属于解析法，它是以Radon变换为理论基础，首先对投影数据在连续域进行一些处理，然后将其离散化进行重建。接下来学习到的是图像重建的迭代算法，该算法的主要思路是：从一幅假设的初始图象出发，采用逐步逼近的方法，将理论投影值与实际测量投影值进行比较，在某种最优准则下寻找最优解。 该算法本质上类似于解方程组，但是在实际的图像重建过程中，由于运算量巨大、方程的欠定性以及测量误差、噪声的影响等原因，直接通过求逆矩阵来解方程组是比较困难的，很难将其应用在实际生活中...
2021-08-10 15:31:36
1826
1
原创 扇束图像重建存在的一些问题
在扇束图像重建的仿真实验中（包括等角和等距扇束），刚开始都遇见了同样的问题，就是重建出来的图像如下图所示： 在该重建结果中，我们能够很清楚的看到该图像出现了重叠现象，与原始图像相差过大，严重影响到医生的诊断。经过反复的实验，终于略微了解了一些出现该问题的原因。 我认为主要原因是扇束只能覆盖一小部分的物体，因此重建之后会出现图中的结果。然后我将焦距和探测器的长度增大之后，在很大程度上解决了该问题。下图展示了增大相关参数之后的效果： ...
2021-08-08 16:56:21
705
1
原创 matlab中使用分段函数对矩阵中的元素进行操作
最近编写程序的时候遇见一个问题，如何对一个矩阵中属于不同范围的元素进行不同的操作？刚开始使用if--elseif--else来对自变量进行划分并对其进行不同的运算，但是发现结果是错误的，程序没有按照自己预想的那样运行。查阅一些资料后发现了一种方法，我觉得还是比较简单的，接下来演示一下。分段函数表达式如下： 我们使用分段函数，对矩阵中大于0的元素加1，对矩阵中小于等于0的元...
2021-08-06 18:05:50
1693
原创 图像重建中常用的滤波器的设计
在医学图像重建的过程中，如果不使用滤波器，重建出来的图像与使用滤波器重建出来的图像的质量相差甚远，效果如下图所示： 由上图我们可以看出，使用滤波器能极大的提高重建图像的质量，其中最为常见就是RL滤波器与SL滤波器。首先介绍RL滤波器，RL滤波器是使用窗函数对斜坡滤波器进行截断产生的。如下图所示： 在使用该滤波器的时候，需要将其离散化，离散化之后的滤波器的函数...
2021-08-05 15:41:36
3305
5
原创 平行束反投影重建
在医学图像重建的过程中，平行束直接反投影重建是最基本的重建方法，掌握直接反投影以后，可以帮助我们很轻松的学会滤波反投影重建的原理。 对于一幅256*256的shepp-logan头模型，通过radon函数可以轻松的得到该模型的平行束投影数据。假设探测器的通道数为256，在进行反投影之前，我们需要对进行简单的修正，使其在任何一个角度下，数据的个数与探测器的通道数相同。代码如下：P=radon(I,theta);[mm,nn] = size(P);%计算投影矩阵的行和列...
2021-08-05 09:01:10
1756
原创 在Matlab中将一幅图片的中心设置为x-y轴坐标原点
我们知道，一幅图像由许许多多的像素组成，这些像素可以理解为一个矩阵（假设该矩阵为img）中的元素。但是在matlab中索引是从1开始，并且坐标原点位于左上角，这里我们将此时的坐标轴记为i-j坐标轴。使用img(i,j)表示其中的任何一个像素。 也许我们会觉得这样也挺好呀，对图像操作没有什么问题。但是我们要想遍历整个图片进行一些操作的时候，我们需要两个for循环。对于matlab而言，运行for循环会大大降低程序的运行效率。。同时，我们在做图像重建的时候，经常要用到坐标原点位...
2021-08-05 07:42:20
4445
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人