精华内容
下载资源
问答
  • 图像重建

    2019-11-19 18:56:06
    图像重建流程 Gate仿真 得到Root数据,用Analysis.C将数据转化成listmode数据 将listmode数据转化为3D的sinogram数据 将3Dsinogram用SSRB转化成2D的slice切片sinogram进行重建 ...

    图像重建流程
    Gate仿真
    得到Root数据,用Analysis.C将数据转化成listmode数据
    将listmode数据转化为3D的sinogram数据
    将3Dsinogram用SSRB转化成2D的slice切片sinogram进行重建

    展开全文
  • CT图像重建:使用MATLAB的计算机断层扫描图像重建项目
  • 图像重建程序

    2018-12-27 14:35:45
    医学图像重建算法很好的代码,自用运行没有问题,包括FDK、MLEM等
  • 相对于MR,CT在参数及扫描方面并不是太难,但是CT的图像重建及各种算法则是非常难的,也是比较抽象的。这篇文章介绍CT图像重建算法等相关内容。CT技术是CT诊断的基础,帮助医务工作者充分掌握CT技术是我们的责任和...

    本文作者:飞利浦资深CT专家,Colin Han,韩太林。感谢CT Clinical Evidence授权转载。

    相对于MR,CT在参数及扫描方面并不是太难,但是CT的图像重建及各种算法则是非常难的,也是比较抽象的。这篇文章介绍CT图像重建算法等相关内容。

    6b1b2d395c07bc374bdc0c8ee0c9784b.png

    CT技术是CT诊断的基础,帮助医务工作者充分掌握CT技术是我们的责任和义务!b1b43beca5e8d9b4a44631fa8fa29dea.png

    CT基础知识

    Image Reconstruction Method

    图像重建算法

    在上一期内容CT原理1中,我们主要介绍了X线与物质的作用,也就是X线的衰减是如何发生的,正是由于这种衰减的存在,X线才可以被用于CT成像,那么探测器获得的衰减信号最终是如何被转换成CT图像的呢?这就涉及到图像的重建算法了,所以今天我们接着X线的衰减,继续介绍CT图像重建的原理。

    1

    重建算法的分类

    8d13966336fb828f301c2e6a7ebb7248.png

    CT重建算法共有3类,如下图

    (1)反投影法

    (2)迭代重建算法

    (3)解析法:包括滤波反投影法和傅里叶变换法

    58bce0fc90d5c7da3dd5d157f8571cb4.png

    其中(3)解析法中的A滤波反投影法是在(1)反投影法的基础上发展起来的,通过加入滤波函数解决了图像锐利度的问题,如下图,没有加入滤波函数(Filter Function)时重建的图像是模糊的(红色),而加入滤波函数(Filter Function)后重建的图像是清晰的(黄色),由于这种方法(FBP)重建速度快,且图像质量比较高,成为了最常用的CT图像重建方法;而(2)迭代重建算法虽然很早就提出了(早于FBP法),但是由于其计算量非常大,重建速度慢,有赖于计算机性能的突破,所以这种方法的普及反而在FBP算法之后。

    64750c81a7d1837b0635716a51b97859.png

    4995fffa37e62910c68586c1451f1254.png

    2

    反投影法原理

    8d13966336fb828f301c2e6a7ebb7248.png

    反投影法是最简单的一种CT图像重建方法

    CT图像重建的基础:相同的X射线强度,穿过不同物质,发生的衰减不同,利用这一规律可以区别人体内的不同物质。

    CT图像重建的过程: 简单来说,CT扫描中X线穿透人体的每个层面内的结构可以被分成多个小立方体(被称为体素:Voxel),每个小立方体都对应一个单独的衰减信号,把这个信号输入到图像平面矩阵中相应的小格子(被称为像素:Pixel)中,把每一个体素的衰减信号都输入到相应的像素内,然后以不同的灰度反映出来,这就是CT图像重建的过程。

    744ed2cdf2ff19abaa611e216549d177.png

    下面以4个体素为例对反投影法的原理进行简单说明:

    假设四个体素的衰减系数分别为1、2、3、4,那么直接把相应体素的衰减系数直接输入到对应的像素内就可以生成图像了。是不是觉得太简单了呢?

    d53c185fcde7cbc3547b4937f376ae33.png

    但实际上,每个体素具体的衰减系数我们是不知道,通过扫描获得的是不同方向上,各体素衰减系数的总和,如下图,投影1-6(方向为从球管到探测器)获得的总衰减系数分别为:7、3、5、4、6、5,而每个体素具体的衰减系数是多少呢?87c42d8303bcde7e67949f105926c6a7.png接下来,我们利用已知的数据进行反投影(方向为从探测器到球管):

    根据探测器获得的数据反过来对该方向上的每个体素进行相同的赋值,称之为反投影。

    方向1反投影结果为:

    99f048dfdc1e93f45e2cf85923b35964.png

    方向2反投影结果为:

    99fb251ebff9b2f3b9a7dcc192c48107.png

    方向3反投影结果为:

    76025575972dbdb1ef8387638c4d3700.png

    方向4反投影结果为:

    56246c165fb16e7260d9fb0214988850.png

    方向5反投影结果为:

    2851a906400cca7747496bf8b6ebd28d.png

    方向6反投影结果为:

    5a8d49642ad420abb370e13e00103651.png

    把这6次反投影的结果叠加到一起为:

    5b73c552d334891005ed40f4bc0c2fa9.png

    叠加的过程导致每个像素的所得数据为3次自身衰减系数+另外三个体素的衰减系数(三个方向的三次反投影),所以需要减去一个相同的基数(四个体素的衰减系数之和),并除以2,才是确切的衰减系数.

    由投影1+投影2或者投影4+投影5可得出4个体素总的衰减系数为10,那么减去基数后结果为:

    1c2391ec1733690ad74dba4c2fe72f53.png

    然后再除以2结果为:

    8c6596118cd4ca83a1ba5eb9a2925d2d.png

    最终获得的这个像素矩阵就是对上述四个体素的准确还原.而这种方法就是CT图像重建中反投影法的基本原理,由于这个过程不需要对每一个体素的具体衰减系数进行求解,而是直接进行反投影获得叠加数据,然后再对获得叠加数据进行线性的调整就能实现对每个体素的还原,因此重建速度非常快,但是这种方法要求各个方向的投影数据必须是完备的,而且噪声与辐射剂量呈线性关系。

    64750c81a7d1837b0635716a51b97859.png

    3

    迭代重建算法原理

    8d13966336fb828f301c2e6a7ebb7248.png

    迭代重建算法是计算量非常大的一种CT图像重建方法

    简单来说就是利用求解线性方程组的方式来重建图像,实际上是从—幅假设的初始图像出发,采用逐步逼近的方法,将理论投影值同实际测量投影值不断进行比较和迭代更新,直到最终获得最优解,最早被用在PET成像中,如下图[1]。

    509f9c9745e30520221b4aca992a03c8.png下面还是以4体素为例对迭代重建算法进行简单说明:

    假设扫描获得的投影值结果为:76e1f58fd242084e5aea93d2a1880e6a.png

    先对图像进行猜测,如每个像素都为10:

    ef91f7c4d9add7bfb6765d7c70183609.png

    然后对这个猜测结果获得的投影值与实际扫描获得的投影值进行比较,获得校准系数(为实际测量的投影值与猜测所得的投影值之比),如下图:

    3095d65da34cc5a55e4ea62dab7028fd.png

    然后根据本次校准系数进行一次迭代:

    b46513592a3577c9db41e9837ca77ee5.png

    然后又获得新的校准系数,并进行二次迭代:

    93df289c5872f5f9c3cc339a24f495ca.png

    然后又获得了新的校准系数,并进行三次迭代:

    02db85aa11705b9c9a83af8d06c16c45.png

    ......

    以此类推,不断迭代,直到最终校准后的各体素衰减值对应的投影值与实际扫描获得的投影值差距为最小(Error=min),则满足了迭代终止条件,然后才根据最后校准获得的体素衰减值来重建图像.这种方法就是CT图像重建中迭代重建法的基本原理,由于这种方法需要对每一个体素所对应的具体衰减系数进行求解,而医学图像所用的矩阵至少为512x512,因此每层图像需要求解的未知数为262144个,所以这种方法的计算量非常之大,故重建速度慢,但是这种方法的噪声不再受到辐射剂量的限制,能够在辐射剂量较低的条件下,获得高质量图像。

    64750c81a7d1837b0635716a51b97859.png

    重点回顾

    CT图像重建的过程就是把每一个体素的衰减系数对应到每一个像素中去,反投影法不需要求解也能进行还原,重建速度快,而迭代重建需要先求得最优解,然后再进行图像重建,速度慢但在低剂量时也可获得高质量图像

    注意事项!

     相同的X射线强度,穿过不同物质,发生的衰减不同,利用这一规律可以区别人体内的不同物质。

    ✔ 反投影法不需要对每一个体素的具体衰减系数进行求解,因此重建速度非常快,但是这种方法要求各个方向的投影数据必须是完备的,而且噪声与辐射剂量呈线性关系.

     迭代重建需要对每一个体素所对应的具体衰减系数进行求解,因此重建速度慢,但能够在辐射剂量较低的条件下,获得高质量图像。

    [1].Knesaurek K, Machac J, Vallabhajosula S, Buchsbaum MS. A new iterative reconstruction technique for attenuation correction in high-resolution positron emission tomography. Eur J Nucl Med. 1996;23(6):656‐661. doi:10.1007/BF00834527.

    展开全文
  • ART图像重建

    2014-09-18 20:24:11
    ART应用于图像重建工作 代数重建 matlab代码
  • matlab ART重建算法,该算法可用于CT图像重建或EST图像重建-Matlab ART reconstruction algorithm, which can be used for CT image reconstruction or image reconstruction EST
  • ART迭代算法图像重建

    2021-02-11 09:00:18
    采用代数迭代算法进行图像重建的MATLAB程序,具有借鉴的价值,欢迎大家分享,谢谢! 采用代数迭代算法进行图像重建的MATLAB程序,具有借鉴的价值,欢迎大家分享,谢谢!
  • CT图像重建代码

    2018-06-02 15:53:24
    CT图像重建的MATLAB代码,根据自己设计的问题更改之后可用
  • MATLAB图像重建代码

    2018-03-29 21:35:54
    MATLAB的图像重建代码,图像重建UI设计是本课题最后呈现出来的系统,需要具有美观,舒适,简约且功能全面等特点。图(a)展示的是GUI设计时的代码编辑面,里面的内容刚好就是读取图片的代码。图(b)显示的是系统的...
  • matlab编写的图像重建迭代算法,ART重建,给出了重建结果。
  • DCT图像重建

    2013-06-09 15:13:28
    图像分为8*8的块 DCT图像重建 自己编写 代码
  • 医学图像重建

    2013-06-01 18:06:47
    医学图像重建,曾更生著,高等教育出版社。
  • ctpet医学图像重建.rar

    2021-02-07 15:09:25
    主要针对从事医学图像重建研究和开发的人员
  • ct图像重建

    2012-12-27 14:27:27
    用于医学ct图像重建,所使用的开发语言是MATLAB
  • 图像重建的统计自适应子集算法
  • 图像重建的分块迭代算法,赵双任,杨新铁,图像重建和重新投影(Reconstruction and Re-projection)的迭代算法已在投影视野受限(Limited Field of View, LFOV)的情下使用。该情况下投影数据被�
  • 英文名称为Algebraic reconstruction ...Gordon、Bender和Herman首次证明了它在图像重建中的应用。而这种方法被称为数值线性代数中的Kaczmarz方法。相对于其他重建方法(如滤波反投影),ART的一个优点是,将先验...

    英文名称为Algebraic reconstruction technique:即代数重建算法

    代数重建技术(ART)是一种用于计算机层析成像的迭代重建技术。它从一系列的角度投影(正弦图)中重建一幅图像。Gordon、Bender和Herman首次证明了它在图像重建中的应用。而这种方法被称为数值线性代数中的Kaczmarz方法。相对于其他重建方法(如滤波反投影),ART的一个优点是,将先验知识(已知的约束条件)纳入重建过程是相对容易的。该算法的实质是用迭代法求解线性方程组的解。

    clc;clear all;close all;N = 180; %图像大小N2 = N^2;I = phantom(N);% 产生头模型图像theta = linspace(0,180,181);theta = theta(1:180);%% = = = = = = 产生投影数据 = = = = = = %%P_num = 260; %探测器通道个数P = ParallelBeam(theta ,N ,P_num); %产生投影数据%P = radon(I,theta);%% = = = = = = 获取投影矩阵 = = = = = =%%delta = 1;% 网格大小,角度增量[W_ind,W_dat] = SystemMatrix(theta,N,P_num,delta);%% = = = = = = 进行ART迭代 = = = = = = %%F = zeros(N2,1);  %初始图像向量lambda = 0.25;  %松弛因子c = 0;  %迭代计数器irt_num = 5;  while(c    for j = 1:length(theta)        for i = 1:1:P_num            % 取得一条射线所穿过的网格编号和长度            u = W_ind((j-1)*P_num + i,:);  % 编号            v = W_dat((j-1)*P_num + i,:);  % 长度                        if any(u) == 0                continue;            end            %恢复投影矩阵中与这一条射线对应的行向量 w            w = zeros(1,N2);            ind = u > 0;            w(u(ind))=v(ind);            % 图像进行一次迭代            PP = w * F;   %前向投影            C = (P(i,j)-PP)/sum(w.^2) * w';  % 修正项                       F = F + lambda * C;        end    end    F(F<0) = 0;   % 小于0的像素值置为0    c = c+1;endF = reshape(F,N,N)'; % 转换成N x N的图像矩阵%% = = = = = = 仿真结果显示 = = = = = = %%figure(1);imshow(I);xlabel('(a)180x180头模型图像');figure(2);%A = imadjust(F);imshow(F,[]);xlabel('(b)ART算法重建的图像');     

    26c63d337c114115d101cb11dc2b2b36.png

    c90a8672f3d3febe785767e61fb4b224.png

    往期回顾>>>>>>

    【模式识别】Matlab指纹识别【优化求解】A*算法解决三维路径规划问题  matlab自动识别银行卡号【优化问题】基于栅格地图——遗传算法的机器人最优路径规划【优化求解】模拟退火遗传实现带时间窗的车辆路径规划问题【数学建模】Matlab实现SEIR模型分享到朋友圈获取6个赞即可获取完整代码

    25a92fa72661737e0b8313fc2e1ea0fa.png

    展开全文
  • CT图像重建例子

    2014-03-07 09:10:27
    基于C++编写的CT图像重建例子,包含CT原始数据
  • 提出一种正则化广义逆ERT 图像重建算法, 利用ERT 仿真软件得到的数据进行图像重建。与 常用的ERT 图像重建算法进行比较, 重建出的图像经过统一的门限滤波后, 反投影算法、 灵敏度系数算 法和正则化广义逆...
  • CT图像重建 matlab 程序

    2013-12-30 22:20:23
    MATLAB图像重建,CT图像重建;matlab代码及实验报告模板
  • 具有感知压缩感的彩色图像重建
  • 空间图像重建的快速碰撞检测算法
  • 代数图像重建算法的充要条件
  • 图像重建技术matlab源代码,原文件为一个老鼠的图片,经重建后达到比较好的效果,入口为test
  • matlab_radon图像重建

    2016-11-08 22:59:59
    matlab代码 利用拉东(radon)变换求解图像重建
  • 这是一个医学图像重建作业matlab源码,通过滤波反投影实现的,可运行,是医学图像计算课程作业
  • 磁共振图像重建

    2014-05-01 15:17:09
    有关磁共振图像重建的代码,希望对大家有用
  • 为了提升稀疏采样环境下的图像重建质量,针对广义全变分模型重建图像时不能充分利用图像本身结构自相似性信息的不足,建立了一个非局部约束下的改进广义全变分图像重建模型。该模型引入了变化域非局部自相似性作为...
  • 基于opencv+ c++ 的图像重建系统,可以参考一下,有数据,有源码。

空空如也

空空如也

1 2 3 4 5 ... 20
收藏数 5,357
精华内容 2,142
关键字:

图像重建