FBP.zip_CBP FBP_teaj4m_卷积反投影_图像处理

clear;clc;close all; %% 图像输入部分 N=input('选择输入圆大小:\n'); %大圆小圆 [X,Y]=meshgrid(linspace(-N/2+0.5,N/2-0.5,N)); image=zeros(N); image(X.^2+Y.^2<(0.4*N)^2)=0.3; image((X-(0.1*N)).^2+(Y-(0.1*N)).^2<(0.1*N)^2)=1; figure(); imshow(image,[]); theta = 0:0.5:179.5; %% 正投影 projectionImg = RadonImrotate(theta,image); figure(),imshow(projectionImg,[]) %% %flag='RL','SL'0.16,'Hamming' reconImg = CBP(projectionImg,theta,'SL'); figure(),imshow(reconImg); figure,imshow(reconImg,[]); %% 反投影 %FBP:s=0.16;%0:1:179.5; %FBP:s=0.32;0:0.5:179.5;l %CBP:s=0.643;0:0.5:179.5; figure,plot(reconImg(:,round(N/2+0.1*N)),'r') hold on,plot(image(:,round(N/2+0.1*N)),'b') legend('重建图像','原始图像'); %% 结果 rss=reconImg-image; figure,imshow(rss) figure,imshow(rss,[]) rss=round(rss,2);