数字图像处理radon matlab变换算法代码

本次试验只要是对radon变换的理解和应用，本试验用的是Matlab软件编写的处理程序。

x=imread（‘1.bmp‘）;
subplot（131）;
imshow（x）;
title（‘原图‘）;
[mn]=size（x）;
M=128;

      %       投影        %
      
proj1=zeros（128128）;     %     定义一个128X128维的零矩阵     %   
for N=1:128               %     设置128个投影角度             %
    t=N*pi/128;           %     定义每个投影角度大小          %
    for i=1:64
        for j=1:64
            y1=（i-n/2）*cos（t）+（j-m/2）*sin（t）;  %   radon 变换    %
            z1=floor（y1）;                       %   取整          %
            proj1（Nz1+64）=radon（Nz1+64）+double（x（ij））*（z1+1-y1）;   %     插值     %
            proj1（Nz1+1+64）=radon（Nz1+1+64）+double（x（ij））*（y1-z1）;
        end
    end
end
subplot（132）;
imshow（proj1/10000）;
title（‘投影图‘）;

              %     反投影       %
        
proj2=fftshift（fft2（proj1））;
proj3=zeros（1M）;
proj3=double（proj3）;

             %     构建Ramp滤波器       %
for i=1:M                        
    proj3（i）=abs（i-M/2.0）/（M/2.0）;
end
proj3=proj3（ones（1281）:）;
proj2=proj2.*proj3;   

             %     构建高斯低通滤波器   %
d=zeros（MM）;
for i=1:M
    for j=1:M
        d（ij）=（（i-（2*M+1）/4）^2+（j-（2*M+1）/4）^2）^.5;
    end
end  
           %       指定高斯滤波器截止频率 Do  %
Do=60;          
Gauss_L=exp（-（（d.^2）./（2*Do^2）））;
proj2=proj2.*Gauss_L;       %    Ramp滤波器与高斯低通滤波器相乘    %
proj4=abs（ifft2（proj2））;
A=max（max（proj4））-min（min（proj4））;
proj4=（proj4-min（min（proj4）））./A.*255;     %   滤波后投影proj4    %

proj5=zeros（128128）;
for i=1:128
    for j=1:128
        x=double（j）-（M+1）/2.0;
        y=（M+1）/2.0-double（i）;
        t=pi/128.0;
        for n=1:128
            s1=x*cos（（n-1）*t）+y*sin（（n-1）*t）+M/2.0;
            if s1<=M&s1>=1
               proj5（ij）=proj5（ij）+（proj4（nfloor（s1））+（s1-floor（s1））*（proj4（nfloor（s1）+1）））;
            end
        end
    end
end
B=max（max（proj5））-min（min（proj5））;    %灰度拉伸
subplot（133）;
proj5=（proj5-min（min（proj5）））.*（255/B）;
imshow（uint8（proj5））;
title（‘滤波反投影图D。=60.bmp‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       5176  2009-05-14 10:45  实验五\1.bmp

     文件       2040  2009-06-07 20:38  实验五\project.m

     文件     123392  2009-06-08 21:51  实验五\实验报告.doc

     文件        160  2009-06-08 21:48  实验五\实验说明.txt

     目录          0  2009-06-10 12:43  实验五

----------- ---------  ---------- -----  ----

               130768                    5