Muc数字图像处理-频域处理+复原-雷登变换实验代码

srcImage = imread('C:/Users/lenovo/Documents/MATLAB/正弦图/牙齿.png'); % 读入图像 srcImage = rgb2gray(srcImage); % 将图像转换为灰度图像 theta1 = 0:0.5:179.5; theta2 = 0:2:178; theta3 = 0:10:170; [R1, x1] = radon(srcImage, theta1); [R2, x2] = radon(srcImage, theta2); [R3, x3] = radon(srcImage, theta3); figure; imagesc(theta1, x1, R1); xlabel('Parallel Rotation Angle - \theta (degree)'); ylabel('Parallel Sensor Position - x'' (pixels)'); colormap(hot); colorbar; I1 = iradon(R1,theta1,'none');%defulat:Ram-Lak滤波器 I2 = iradon(R1,theta1,'Hamming'); I3 = iradon(R1,theta1); I4 = iradon(R2,theta2,'none');%defulat:Ram-Lak滤波器 I5 = iradon(R2,theta2,'Hamming'); I6 = iradon(R2,theta2); I7 = iradon(R3,theta3,'none');%defulat:Ram-Lak滤波器 I8 = iradon(R3,theta3,'Hamming'); I9 = iradon(R3,theta3); figure; subplot(3,3,1);imshow(I1,[]);title('theta1=0:0.5:179.5 无滤波'); subplot(3,3,2);imshow(I2,[]);title('theta1=0:0.5:179.5 Hamming滤波'); subplot(3,3,3);imshow(I3,[]);title('theta1=0:0.5:179.5 Ram-Lak滤波'); subplot(3,3,4);imshow(I4,[]);title('theta2=0:2:178 无滤波'); subplot(3,3,5);imshow(I5,[]);title('theta2=0:2:178 Hamming滤波'); subplot(3,3,6);imshow(I6,[]);title('theta2=0:2:178 Ram-Lak滤波'); subplot(3,3,7);imshow(I7,[]);title('theta3=0:10:170 无滤波'); subplot(3,3,8);imshow(I8,[]);title('theta3=0:10:170 Hamming滤波'); subplot(3,3,9);imshow(I9,[]);title('theta3=0:10:170 Ram-Lak滤波');