图像去噪

分数阶各向异性扩散图像去噪，实现了文献Fractional-Order Anisotropic Diffusion for Image Denoising， Jian Bai， Xiang Chu， IEEE T. Image Process.，2007， 16(10): 2492-3502

img = double（imread（‘standard_lena.bmp‘））;
n=512;
uo=zeros（nn）;
uo=img（1:n1:n）;
init = 2055615866;
randn（‘seed‘ init）;
sigma = 1;
figure;imagesc（uo）;colormap（gray）;
un =sigma * randn（size（uo））;
u = uo +25*un;
SNR1（uou）
figure;imagesc（u）;colormap（gray）;
u1=fft2（u）;
a=1.5;
t=4^（-a）;
j=sqrt（-1）;
g1=zeros（nn）;
g2=zeros（nn）;
g3=zeros（nn）;
g4=zeros（nn）;
for i=1:n
    p（i）=（1-exp（-2*3.1415926535*j*（i-1）/n））^a*exp（2*3.1415926535*j*（i-1）/n）^（a/2）;
    q（i）=conj（p（i））;
end
for k=1:3300
    k
    ty=0;
    for i=1:n
        g1（:i）=p（i）*u1（:i）;
        g2（i:）=p（i）*u1（i:）;
    end
    g1=real（ifft2（g1））;
    g2=real（ifft2（g2））;
    g3=g1./abs（1+g1.^2+g2.^2）;
    g4=g2./abs（1+g1.^2+g2.^2）;
    g3=fft2（g3）;
    g4=fft2（g4）;
    for i=1:n
        g1（:i）=q（i）*g3（:i）;
        g2（i:）=q（i）*g4（i:）;
    end
    u1=u1-（g1+g2）*t;
    ut=ifft2（u1）;
    SNR1（uoreal（ut））
end
ut=ifft2（u1）;
figure;imagesc（real（ut））;colormap（gray）;
SNR1（uout）
SNR1（uereal（ut））

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-11-22 22:56  frac\
     文件      263222  2009-06-11 13:13  frac\BARB.BMP
     文件         249  2009-06-11 13:13  frac\SNR1.m
     文件        7168  2011-07-18 12:45  frac\Thumbs.db
     文件        1012  2009-11-27 15:29  frac\dqzb.m
     文件        1114  2009-06-11 13:13  frac\fdqz.m
     文件         164  2009-06-11 13:13  frac\fractiond2.m
     文件         172  2009-06-11 13:13  frac\fractiondc2.m
     文件      263222  2009-06-11 13:13  frac\standard_lena.bmp