1.算法描述

直接由正弦图得到反投影图像，会存在严重的模糊，这是早期 CT 系统所存在的问题。傅立叶中心切片定理表明，投影的一维傅立叶变换是得到投影区域的二维傅立叶变换的切片。滤波反投影重建算法在反投影前将每一个采集投影角度下的投影进行卷积处理，从而改善点扩散函数引起的形状伪影，有效地改善了重建的图像质量。

括号 [] 内部是一个一维傅里叶反变换，可以认为这是一个一维滤波器的传递函数。由于 ∣ ω ∣ |\omega|∣ω∣ 是一个不可积的斜坡函数（Slope function），可以通过对斜坡加窗进行限制，典型地如汉明窗（Hamming window）、韩窗（Hann window）。该式也可以使用空间卷积来实现：

这表明，将对应的投影 g ( ρ , θ ) g(\rho, \theta)g(ρ,θ) 与斜坡滤波器传递函数 s ( ρ ) s(\rho)s(ρ) 的傅里叶反变换进行卷积，可以得到角度 θ \thetaθ 的各个反投影，整个反投影图像可以通过对所有反投影图像积分得到。

滤波反投影，运用滤波反投影的方法对图像进行重建matlab仿真+程序操作视频对比包括 Ram-Lak (Ramp),Shepp-Logan,Ram-Lak Cosine,Ram-Lak Hamming,Ram-Lak Hann。

Ram-Lak

斜坡滤波器（甚至在被加窗时）在频率域的直流项为零，故每一幅反投影图像的均值将为零。这将意味着，每一幅反投影图像都将有正像素和负像素值，当所有的反投影图像相加形成最终的重建图像时，一些负像素值位置可能变成正像素，而平均值可能不为零，但是，典型地，最终的图像将还是有负像素值。Ram-Lak滤波函数如下所示：

Shepp-Logan

Shepp-Logan头模型是由L. A. Shepp和B. F. Logan于1974年首次提出的，它由10 个位置、大小、方向、密度各异的椭圆叠加而成，模拟一个脑部的断层。Shepp-Logan滤波函数如下所示：

2.仿真效果预览

matlab2022a仿真结果如下：

3.MATLAB部分代码预览

  function PushButton2Selected(h, eventdata)switch list1_valcase 'None'filter_val = 'None';case 'Ram-Lak (Ramp)'filter_val = 'Ram-Lak';case 'Shepp-Logan'filter_val = 'Shepp-Logan';case 'Ram-Lak Cosine'filter_val = 'Cosine';case 'Ram-Lak Hamming'filter_val = 'Hamming';case 'Ram-Lak Hann'filter_val = 'Hann';case 'Special'filter_val = 'Special';endswitch list2_valcase 'Linear'interpol_val = 'linear';case 'Nearest'interpol_val = 'nearest';case 'Spline'interpol_val = 'spline';case 'Cubic'interpol_val = 'pchip';endtheta = round(0 : 180/antal_projektioner : 179);del_sinogram = data.sinogram(:, theta+1);if strcmp(filter_val, 'Special') == 1if strcmp(file_name(end-4:end), '3.mat') == 1data.SB = data.SB + 0.3 * rand(size(data.SB));endif strcmp(file_name(end-4:end), 'f.mat') == 1filter_val = 'Hamming';rec_img = iradon(del_sinogram, theta, interpol_val, filter_val);elseFFT_sinogram = fftshift(fft(del_sinogram),1);FFT_sinogram = FFT_sinogram .* data.SB(:, theta+1);del_sinogram = ifft(ifftshift(FFT_sinogram,1));rec_img = iradon(del_sinogram, theta, interpol_val);endelserec_img = iradon(del_sinogram, theta, interpol_val, filter_val);endfigure, imagesc(rec_img), colormap(pink)title([data.txt ': ' num2str(antal_projektioner) ' projections, ' list1_val ' filter, ' list2_val ' interpolation'])end
A_023

4.完整MATLAB程序

matlab源码说明_我爱C编程的博客-CSDN博客

运用滤波反投影的方法对图像进行重建matlab仿真相关推荐

图像滤波边界处理方法对图像滤波的影响
%边缘效应与滤波 %实现空间滤波领域处理时的一个重要考虑因素,当滤波掩模的中心靠近图像轮廓时发生的情况 % 当一个n*n的方形掩模,其中心距离图像边缘为(n-1)/2个像素时,该掩模至少有一条边与图像 ...
基于基于全局差错能量函数的双目图像立体匹配算法matlab仿真,并提取图像的深度信息
目录 1.算法概述 2.仿真效果预览 3.核心MATLAB代码预览 4.完整MATLAB程序 1.算法概述全局的能量函数公式如下: E(f)=Edata(f)+Esmooth(f) 其中,Edata ...
m基于POCS算法的空域序列图像超分辨率重建matlab仿真
目录 1.算法仿真效果 2.MATLAB核心程序 3.算法涉及理论知识概要 4.完整MATLAB 1.算法仿真效果 matlab2022a仿真结果如下: 2.MATLAB核心程序 .......... ...
基于Luca-Kanade光流算法的图像运动场提取matlab仿真
目录 1.算法概述 2.仿真效果 3.MATLAB仿真源码 1.算法概述光流是一种运动模式,这种运动模式指的是一个物体.表面.边缘在一个视角下由一个观察者(比如眼睛.摄像头等)和背景之间形成的明显移 ...
怎么把matlab仿真数据压缩,JPEG图像压缩编码及其MATLAB仿真实现(1)
信息技术! ! 电! 子! 测! 量! 技! 术 ! ! "#"$%&'()$ *"+,-&"*"(% %"$.('#'/0 ...
基于sift特征提取的图像配准算法matlab仿真
目录 1.算法描述 2.仿真效果预览 3.MATLAB核心程序 4.完整MATLAB 1.算法描述 SIFT 是一种从图像中提取独特不变特征的方法,其特点为基于图像的一些局部特征,而与图像整体的大小和 ...
【MATLAB教程案例47】基于双目相机拍摄图像的三维重建matlab仿真
欢迎订阅<FPGA学习入门100例教程>.<MATLAB学习入门100例教程> 本课程学习成果预览: 目录 1.软件版本 2.基于双目相机拍摄图像的三维重建原理概述
【IHS+稀疏表示】基于IHS小波变换和稀疏表示相结合的遥感图像融合算法matlab仿真
1.软件版本 matlab2013b 2.本算法理论知识小波变换具有多尺度特性,其时域和频域均具有表征信号局部特性的能力和多分辨率分析的特点.在低频时具有高的频谱分辨率,在高频时具有低的频谱分 ...
基于affine+sift特征提取的图像配准算法matlab仿真
目录一.理论基础二.核心程序三.仿真测试结果一.理论基础由于相机正面拍摄物体时,相机的光轴方向可能发生变化,带来扭曲,所以第一步是对每张图片进行变换,模拟所有可能的仿射扭曲.这些扭曲由两个参 ...

运用滤波反投影的方法对图像进行重建matlab仿真

1.算法描述

2.仿真效果预览

3.MATLAB部分代码预览

4.完整MATLAB程序

运用滤波反投影的方法对图像进行重建matlab仿真相关推荐

最新文章

热门文章