m基于shepp-Logan模型和滤波反投影的医学图像多尺度全局重建和局部重建matlab仿真

本文涉及的产品
模型在线服务 PAI-EAS,A10/V100等 500元 1个月
模型训练 PAI-DLC,5000CU*H 3个月
交互式建模 PAI-DSW,5000CU*H 3个月
简介: m基于shepp-Logan模型和滤波反投影的医学图像多尺度全局重建和局部重建matlab仿真

1.算法描述

    从投影重建物体的截面图像是图像处理中非常重要的技术此技术在物体的无损伤性检测其内部缺陷的应用中能起很大作用从投影重建图像的技术早在20世纪中期就已经制成常规医疗诊断设备的商品1917年奥地利数学家J.Radon发表的论文证明了二维物体或三维物体可以从许多投影来重建其内部数据英国的EMI公司中央研究所1972年由G.N.Hounsfield设计出X射线断层分析仪此后于1974年5月在加拿大蒙特利尔第一次国际会议上正式命名为Computerized Tomography简称XCT或CT当时这台仪器可以获得人体的各个部位鲜明清晰的断层图像另外美国塔夫茨大学A.M.Cormack在1963年就提出了精确的非迭代的级数展开法用于图像重建的算法中并指出在诊断医学中图像重建的可能性以上两人由于对CT研制工作的开创性贡献获得了1979年诺贝尔医学奖。

   计算机层析成像技术(Computed Tomography,简称CT)是利用具有一定能量的射线源(X射线,γ射线)对物体进行断层扫描,并根据物体外部的探测器获得的物理量(指物质对射线的衰减系数)生成的一维投影数据,通过特定的重建算法,得到所扫描断层的二维图像。CT断层图像具有无影像重叠、空间和密度分辨率高、可直接进行数字化处理等优点,通过近十几年的发展已成为非接触无损检测的主流技术,是关键部件检测、机械仿型设计、安全检查等方面强有力的手段,并广泛应用于航空、航天、机械、汽车、船舶、公安等领域。此外,在医学领域,医用CT机已成为疾病诊断的重要辅助工具,对病灶的位置及病变程度的良好再现已使其成为医学影像领域里的主要医疗器械。CT机有三个组成部分。一是数据采集系统,二是图像重建(主要是图像重建算法),三是图像后处理及显示系统。CT图像重建算法主要分为两类,即变换法和级数展开法。变换法包括滤波反投影算法、傅立叶变换法以及ρ滤波法(rho-filtered layergrams),小波变换法;级数展开法主要是代数重建算法。滤波反投影是比较常用的变换算法,它具有速度快,空间和密度分辨率好的优点。

   滤波反投影法是目前CT图像重建领域使用得最为广泛的一种算法,其有效性和准确性也在临床上得到了验证。

9238b013cb4c9d3f3a6b1fcc0311a7f7_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png

 通常在计算机中,采用的离散反投影算法,因此离散情形下的卷积反投影算法可以如下描述:

2435175573d3d4c2587ffc511dee4c08_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png

   其中O(x,y)代表原始图像,R(x,y)代表重建图像。当然孤立地去看相对误差并不是一个好的衡量标准,它必须同人们的视觉评估结合起来进行考虑。 

shepp-Logan模型:

    用S-L头模型进行计算机仿真研究的主要好处之一是可以获得该模型投影数据的解析表达式。一个中心位置在原点且未经旋转的椭圆,其长轴与X轴重合,短轴与Y轴重合。假设椭圆内的密度为r、椭圆外密度为零,则该椭圆图像可用以下方程表示:

1996dba716d4379891155f9a90431d56_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png

式中A,B分别为椭圆长轴与短轴的长度。

6ab7983e2f589335998a1aed48dd32dc_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png

    在研究从投影重建图像的算法时,为了比较客观的评价各种重建算法的有效性,人们常选用公认的Sheep-Logan 头模型作为研究对象。该模型由10个位置、大小、方向、密度各异的椭圆组成,象征一个脑断层图像。

该模型的主要参数如下表所示:

de9b2f30a792ddb4ff5482dc09848a2d_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png

多尺度全局重建算法计算流程:

586a9fcc803ddc419f7d0ccb11c4a172_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png

2.仿真效果预览
matlab2022a仿真结果如下:

9dd9a1f0feb4ae8c86ce249a848c11c9_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png
4638f91a03108dd03379b42b8ba46eea_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png
5921fd914e56382002cf53d387d26ca5_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png
daacf522ca5713588d703ac58e8500ac_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png
17a8eb4c7de5b936a36ab0ce4aab1ec9_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png
a85571ada69627481929974571ca07ca_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png
39fe355fd7b79f63e285b749d00646d8_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png
1c20a1b8903661d174b0cf7e2a35ff03_watermark,size_14,text_QDUxQ1RP5Y2a5a6i,color_FFFFFF,t_100,g_se,x_10,y_10,shadow_20,type_ZmFuZ3poZW5naGVpdGk=.png

3.MATLAB核心程序

 
n = size(PR,1);
sideSize = n;
a = 1;
[Length, Count] = size(PR);
w = [-pi:(2*pi)/Length:pi-(2*pi)/Length];
 
rn1 = abs(2/a*sin(a.*w./2));
rn2 = sin(a.*w./2);
rd = (a*w)./2;
r = rn1*(rn2/rd)^2;
 
f = fftshift(r);
for i = 1:Count
        IMG = fft(PR(:,i));
        fimg = IMG.*f';
        g(:,i) = ifft(fimg);%利用FFT,IFFT进行滤波
end
filtPR= real(g);
............................................................
I=phantom(256);  %生产头部模型图
figure(1);
imshow(I);       %显示图像
IMG=double(I);   %双精度显示
[cod_a,cod_h,cod_v,cod_d,map]=wtest(IMG);%一层小波变换
figure(2);
subplot(221)%显示分解系数1
imshow(cod_a,map);
subplot(222)%显示分解系数2
imshow(cod_h,map);
subplot(223)%显示分解系数3
imshow(cod_v,map);
subplot(224)%显示分解系数4
imshow(cod_d,map);
Z1=idwt2(cod_a,cod_h,cod_v,cod_d,'db1');%重建
figure(3);
imshow(Z1,map);
for i=1:256
    for j=1:256
        
        error(i,j)=(I(i,j)-Z1(i,j))^2/I(i,j)^2;%误差求解
    end
end
for i=1:256
    for j=1:256
        if I(i,j)==0
        error(i,j)=255;    
        end
       
    end
end
for i=1:256
    for j=1:256
        if error(i,j)>255;
        error(i,j)=255;    
        end
       
    end
end
相关文章
|
7天前
|
算法 5G 数据安全/隐私保护
大规模MIMO通信系统信道估计matlab性能仿真,对比LS,OMP,MOMP以及CoSaMP
本文介绍了大规模MIMO系统中的信道估计方法,包括最小二乘法(LS)、正交匹配追踪(OMP)、多正交匹配追踪(MOMP)和压缩感知算法CoSaMP。展示了MATLAB 2022a仿真的结果,验证了不同算法在信道估计中的表现。最小二乘法适用于非稀疏信道,而OMP、MOMP和CoSaMP更适合稀疏信道。MATLAB核心程序实现了这些算法并进行了性能对比。以下是部分
152 84
|
6天前
|
算法
基于GA遗传优化的TSP问题最优路线规划matlab仿真
本项目使用遗传算法(GA)解决旅行商问题(TSP),目标是在访问一系列城市后返回起点的最短路径。TSP属于NP-难问题,启发式方法尤其GA在此类问题上表现出色。项目在MATLAB 2022a中实现,通过编码、初始化种群、适应度评估、选择、交叉与变异等步骤,最终展示适应度收敛曲线及最优路径。
|
7天前
|
算法 BI Serverless
基于鱼群算法的散热片形状优化matlab仿真
本研究利用浴盆曲线模拟空隙外形,并通过鱼群算法(FSA)优化浴盆曲线参数,以获得最佳孔隙度值及对应的R值。FSA通过模拟鱼群的聚群、避障和觅食行为,实现高效全局搜索。具体步骤包括初始化鱼群、计算适应度值、更新位置及判断终止条件。最终确定散热片的最佳形状参数。仿真结果显示该方法能显著提高优化效率。相关代码使用MATLAB 2022a实现。
|
6天前
|
监控 算法 数据安全/隐私保护
基于视觉工具箱和背景差法的行人检测,行走轨迹跟踪,人员行走习惯统计matlab仿真
该算法基于Matlab 2022a,利用视觉工具箱和背景差法实现行人检测与轨迹跟踪,通过构建背景模型(如GMM),对比当前帧与模型差异,识别运动物体并统计行走习惯,包括轨迹、速度及停留时间等特征。演示三维图中幅度越大代表更常走的路线。完整代码含中文注释及操作视频。
|
1月前
|
安全
【2023高教社杯】D题 圈养湖羊的空间利用率 问题分析、数学模型及MATLAB代码
本文介绍了2023年高教社杯数学建模竞赛D题的圈养湖羊空间利用率问题,包括问题分析、数学模型建立和MATLAB代码实现,旨在优化养殖场的生产计划和空间利用效率。
102 6
【2023高教社杯】D题 圈养湖羊的空间利用率 问题分析、数学模型及MATLAB代码
|
1月前
|
存储 算法 搜索推荐
【2022年华为杯数学建模】B题 方形件组批优化问题 方案及MATLAB代码实现
本文提供了2022年华为杯数学建模竞赛B题的详细方案和MATLAB代码实现,包括方形件组批优化问题和排样优化问题,以及相关数学模型的建立和求解方法。
74 3
【2022年华为杯数学建模】B题 方形件组批优化问题 方案及MATLAB代码实现
|
1月前
|
数据采集 存储 移动开发
【2023五一杯数学建模】 B题 快递需求分析问题 建模方案及MATLAB实现代码
本文介绍了2023年五一杯数学建模竞赛B题的解题方法,详细阐述了如何通过数学建模和MATLAB编程来分析快递需求、预测运输数量、优化运输成本,并估计固定和非固定需求,提供了完整的建模方案和代码实现。
58 0
【2023五一杯数学建模】 B题 快递需求分析问题 建模方案及MATLAB实现代码
|
4月前
|
数据安全/隐私保护
耐震时程曲线,matlab代码,自定义反应谱与地震波,优化源代码,地震波耐震时程曲线
地震波格式转换、时程转换、峰值调整、规范反应谱、计算反应谱、计算持时、生成人工波、时频域转换、数据滤波、基线校正、Arias截波、傅里叶变换、耐震时程曲线、脉冲波合成与提取、三联反应谱、地震动参数、延性反应谱、地震波缩尺、功率谱密度
基于混合整数规划的微网储能电池容量规划(matlab代码)
基于混合整数规划的微网储能电池容量规划(matlab代码)
|
4月前
|
算法 调度
含多微网租赁共享储能的配电网博弈优化调度(含matlab代码)
含多微网租赁共享储能的配电网博弈优化调度(含matlab代码)

相关实验场景

更多