精华内容
下载资源
问答
  • 利用迭代重建法SIRT,用户可以通过扫描得到的数据或者读取有关的数据内容,输入到工具箱,然后通过工具箱的计算得到该事物的雏形。 图象重建工具箱的优势在于应用范围很广,无论医疗,建筑,工艺制造等等,都能用...
  • 提出了一种基于代数迭代的Mojette变换层析重建算法。在确定最佳投影角度的基础上, 结合传统层析技术中的乘性代数迭代算法进行重建。利用该算法对轴对称火焰进行了层析重建。数值模拟结果表明, 与基于角的重建算法...
  • 基于线积分模型的代数迭代算法程序,整个程序是基于matlab开发。压缩包内有两个*.m文件,ite.m是执行程序,medfuncMatirx.m是计算正投影矩阵的程序。正投影矩阵的计算是基于线积分模型,结合西北大学张顺利老师的...
  • 针对传统重建算法对火焰重建精度低、重建速度慢的问题,提出了基于正则先验的全变差代数迭代(ARTTV)算法,以提高对称与非对称火焰的重建精度。同时,为了提高重建速度,建立了基于“ARTTV-粒子群算法(PSO)内核”的极限...
  • 采用傅里叶变换相位展开技术处理莫尔条纹,提取轴向任意截面的投影信息,使用一种新的莫尔层析迭代算法重建温度分布,该迭代算法与已有的代数重建算法不同,它由莫尔偏折基本公式和网格重建技术推导得出,直接使用...
  • 构造了平滑约束矩阵作为先验信息引入到重建迭代过程, 建立了一种平滑约束OSEM(SC-OSEM)迭代重建算法。分别将中值滤波、全变差最小(TVM)方法作为平滑约束条件, 通过数值模拟, 针对不完备理想投影数据、含金属不完备...
  • 为研究层析成像算法对温度场二维重建质量的影响,实现了两种典型重建算法:代数迭代重建算法(ART)和模拟退火(SA)算法。在不同的射线分布和吸收谱线数目情况下,使用两种算法对给定单峰温度场和双峰温度场分别...
  • ART迭代算法图像重建

    2021-02-11 09:00:18
    采用代数迭代算法进行图像重建的MATLAB程序,具有借鉴的价值,欢迎大家分享,谢谢! 采用代数迭代算法进行图像重建的MATLAB程序,具有借鉴的价值,欢迎大家分享,谢谢!
  • 软件工程论文基于图像旋转的迭代式CT重建算法软件工程研究 本文是一篇研究软件工程的论文本文采用可以稀疏重建的迭代重建算法提出了一种基于图像旋转的全变分迭代重建算法主要研究内容如下1. 深入研究并实现了迭代...
  • 为此提出了一种新的光学层析技术的代数迭代重建算法,在算法中引入了包含先验知识的属性矩阵,并摒弃了通常所采用的对超松弛系数人为的确定取法,采用了变超松弛系数。实验计算结果表明,引入属性矩阵和变超松弛系数的...
  • 广义代数重建算法(E-ART)可以重建得到高质量的双基材料图像,但是该算法收敛速度较慢。针对该问题,提出一种加速收敛算法(AE-ART),其基本思想是通过为基函数添加权重来增大高能和低能投影曲线的夹角或减小高能和低能...
  • 英文名称为Algebraic reconstruction technique:即代数重建算法代数重建技术(ART)是一种用于计算机层析成像的迭代重建技术。它从一系列的角度投影(正弦图)中重建一幅图像。Gordon、Bender和Herman首次证明了它在...

    英文名称为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算法重建的图像');     

    547d5d9d4dcd257d83571c76c8cfce95.png

    b32ad22b15caf4b82c7d6c7448d9ee76.png

    往期回顾>>>>>>

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

    f4c728b3d434fef2f9781b0d4602cd1c.png

    展开全文
  • 详细分析了投影噪声、投影方向数、场分布性质对重建精度的影响,并与代数迭代重建算法结果进行对比.结果表明,该算法以两个正交方向投影数据重建单峰余弦模拟场平均误差仅为0.03%,而代数迭代重建算法为3.81%;该算法以...
  • 针对运用最小二乘法求解DV-Hop定位算法带来的节点定位误差较大的问题,提出基于代数重建法的DV-Hop定位算法,运用一种由图象重建问题而引入的逐次迭代算法——代数重建法。仿真结果表明,改进算法能降低无线...
  • 在分析了相位衬度CT的特点之后,将压缩感知理论引入相位衬度CT重建中,并在该理论框架下将L1约束融入代数迭代重建(ART)算法中,提出了一种微分相位衬度CT重建算法。数值模拟和实际实验表明,该方法可以根据少量投影...
  • 分析了不同分布先验概率参数和不同噪声信噪比等影响因子对算法重建结果的影响,并与高斯牛顿迭代、贝叶斯迭代重建算法相比较。结果表明,所提算法重建速度有很大提高,并能较好反映被测介质的相对位置。
  • 为了解决这一问题,针对三维锥束CT情况,研究了一种增加惩罚项的联合代数迭代算法,同时研究了该算法与有序子集结合的方法。计算机仿真试验表明:在锥束CT图像重建中,该方法能够在抑制噪声的同时提高重建图像的收敛...
  • 采用代数迭代(ART)算法和最大似然期望最大化(MLEM)算法,利用开路傅里叶变换红外(OP-FTIR)光谱仪的测量结果,通过仿真模拟了高斯空间分布模型下的气体二维浓度场重建,并利用重建评价指标——逼近度和相关系数,分析了这...
  • 建立了基于参考温度的多光谱测温法数学模型,提出了一种基于模拟退火理论的变松弛因子的光学层析技术新算法,通过计算机数值模拟,详细考察了该算法对非对称温度场分布的重建效果并与传统的代数迭代重建算法及滤波反...
  • 采用数值模拟方法,利用代数迭代重建算法(ART)与最小二乘QR分解算法进行重建与对比,并考察了气体浓度边界条件对重建结果的影响。实验室内利用体积分数为4 %的NH3混合气体作为测量介质,重建了测量区域内圆形烟道出口...
  • 迭代算法 Iterative methods KA:Kaczmarz算法 KA: Kaczmarz’s Algorithm ART:代数重建技术 ...SIRT:同步迭代重建技术 SIRT: Simultaneous Iterative Reconstruction Technique CG:共轭梯度算法 CG: Conju...

    迭代算法

    Iterative methods

    KA:Kaczmarz算法

    KA: Kaczmarz’s Algorithm

    ART:代数重建技术

    ART: Algebraic Reconstruction Technique

    SIRT:同步迭代重建技术

    SIRT: Simultaneous Iterative Reconstruction Technique

    CG:共轭梯度算法

    CG: Conjugate Gradient method

    LM:Levenberg-Marquardt算法

    LM: Levenberg-Marquardt method

    全局优化算法

    Global optimization methods

    GS:网格搜索方法

    GS: Grid-Search method

    MC:蒙特卡罗法

    MC: Monte-Carlo method

    SA:模拟退火法,包括Metropolis和HeatBath

    SA: Simulated Annealing method, including Metropolis and HeatBath

    GA:遗传算法

    GA: Genetic Algorithm

    参考文献:

    Aster R C, Borchers B, Thurber C H. Parameter estimation and inverse problems [M]. Academic Press, 2011.

    Sen M K, Stoffa P L. Global optimization methods in geophysical inversion [M]. Cambridge University Press, 2013.

    完整源码下载地址:

    http://page2.dfpan.com/fs/4lcje2219291061f358/

    更多精彩文章请关注微信号:在这里插入图片描述

    展开全文
  • 英文名称为Algebraic reconstruction technique:即代数重建算法代数重建技术(ART)是一种用于计算机层析成像的迭代重建技术。它从一系列的角度投影(正弦图)中重建一幅图像。Gordon、Bender和Herman首次证明了它在...

    英文名称为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

    展开全文
  • 在有无噪声情况下比较代数重建算法、联合代数重建算法及同时迭代重建算法的重建图像,证明了SIRT迭代重建具有较高的图像重建精度与较强的抗噪声性能,因此将SIRT作为高噪声中子投影图像CT迭代重建算法的保真项。...
  • MATLAB-ART算法

    千次阅读 2019-04-14 21:03:24
    ART(Algebra Reconstruction Technique, ART),即代数重建法。 在图像重建方法中,迭代重建法的经典方法是Gorden R....提出的联合迭代重建算法(Simultaneous Iterative Reconstruction Technique...

    ART(Algebra Reconstruction Technique, ART),即代数重建法。

    在图像重建方法中,迭代重建法的经典方法是Gorden R.等提出的代数重建法(Algebra Reconstruction Technique, ART),及Gilbert P.提出的联合迭代重建算法(Simultaneous Iterative Reconstruction Technique, SIRT)。

    (1) 联合代数重建方法(SART) 代数重建算法在迭代过程中,每次投影计算的修正值并不是完全相同,穿过同一像素网格时,图像的模糊误差修正将会引起重建区域的严重噪声,且算法需要较多的迭代次数才能得到较好的重建结果,重建效率不高。针对这些问题Anderson和Kak于1984年提出了联合代数重建算法。该算法对于每个像素是同一投影角度内通过该像素的所有射线误差值之累加,其实质就是对ART中的噪声进行了平滑,因此可以获得较为理想的重建结果。

    (2) 乘型代数重建方法(Multiplicative ATR ,MATR) 上面所介绍的ART算法在求解过程中,每个像素的校正过程都是附加上一个校正值来实现的,因此通常称为“加型”ATR算法。如果当每一个像素的校正是通过乘上一个校正值来实现的,便得到所谓的“乘型”代数重建方法,MATR算法的特点是:初始估计值产中的每个分量值必须大于零。并且在迭代过程中像素的值变为零的点,它将始终保持为零。

    (3) 改进的联合代数重建方法(Modified SATR, MSATR) 联合代数重建算法在针对简单的中心对称图像进行重建时,若投影角度较少,会出现严重的边缘效应。虽然边缘噪声分布区域并非图像感兴趣区域,而且也不是影响图象质量的主要矛盾,但其使得重建图像和原始图像之间的误差趋近于零,导致迭代过程无法进行,最终得到的重建结果中间区域失真,无法达到重建要求。改进的联合代数重建算法就是为了解决边缘效应问题而出现的,是由华中科技大学的李春芳等人提出,由于投影和重建数据计算过程均来自于加权因子和估值的乘积,误差来源于迭代过程中估值的前后不一致,而权因子只是起到了放大误差的作用,因此不应采取联合代数重建中以加权因子作为误差分配的唯一准则的方法,而应对图像边缘部分和中间部分采取相同量级的修正。MSART的初始值不能采用零,而应设置为非零的相等值。

    其本质意义很艾略特理论相似,涨完了跌,跌完了涨,没什么使用价值 ,不过就是多了一层数字转换,算是障眼法

    展开全文
  • 首先基于滤波反投影(FBP)重建算法和联合代数迭代重建技术(SART)算法分别重建图像,得到两种重建算法的重建数值区间范围,然后将FBP的数值区间映射到SART的数值区间,最后将经过映射之后的重建图像结合内容自适应图像...
  • ART算法——matlab实现

    千次阅读 2019-11-30 17:50:46
    英文名称为Algebraic reconstruction technique:即代数重建算法 代数重建技术(ART)是一种用于计算机层析成像的迭代重建技术。它从一系列的角度投影(正弦图)中重建一幅图像。Gordon、Bender和Herman首次证明了它在...
  • 分析了DPC-CT的特点,在凸集投影(POCS)的理论框架下,将L1范数、曲波系数约束和经典的代数迭代算法(ART)相结合提出了一种适合DPC-CT的重建算法。数值模拟和实验的结果表明,该方法可以根据少量投影数据获得较好的...

空空如也

空空如也

1 2 3
收藏数 43
精华内容 17
关键字:

代数迭代重建算法