极化雷达图像Lee滤波源代码MATLAB版

极化雷达图像Lee滤波源代码（MATLAB版）

clear
fid=fopen（‘E:\CAUC\毕业设计--TEACHERS\毕业设计--数据\旧金山1\C3\C11.bin‘‘rb‘）;C11=fread（fid‘single‘）;
fid1=fopen（‘E:\CAUC\毕业设计--TEACHERS\毕业设计--数据\旧金山1\C3\C22.bin‘‘rb‘）;C22=fread（fid1‘single‘）;
fid2=fopen（‘E:\CAUC\毕业设计--TEACHERS\毕业设计--数据\旧金山1\C3\C33.bin‘‘rb‘）;C33=fread（fid2‘single‘）;
C1=C11;C2=C22;C3=C33;
span1=C11+C22+C33;
fid=fopen（‘E:\CAUC\毕业设计--TEACHERS\毕业设计--数据\旧金山1\Lee\C3\C11.bin‘‘rb‘）;C11=fread（fid‘single‘）;
fid1=fopen（‘E:\CAUC\毕业设计--TEACHERS\毕业设计--数据\旧金山1\Lee\C3\C22.bin‘‘rb‘）;C22=fread（fid1‘single‘）;
fid2=fopen（‘E:\CAUC\毕业设计--TEACHERS\毕业设计--数据\旧金山1\Lee\C3\C33.bin‘‘rb‘）;C33=fread（fid2‘single‘）;
C4=C11;C5=C22;C6=C33;
span2=C11+C22+C33;
C1=reshape（C1300300）;C2=reshape（C2300300）;C3=reshape（C3300300）;
C4=reshape（C4300300）;C5=reshape（C5300300）;C6=reshape（C6300300）;
span1=reshape（span1300300）;span2=reshape（span2300300）;
m=10;n=10;
b1=C1（m:m+9n:n+9）;b2=C2（m:m+9n:n+9）;b3=C3（m:m+9n:n+9）;
b4=C4（m:m+9n:n+9）;b5=C5（m:m+9n:n+9）;b6=C6（m:m+9n:n+9）;
b7=span1（m:m+9n:n+9）;b8=span2（m:m+9n:n+9）;
epi1=0;epi2=0;epi3=0;epi4=0;epi5=0;epi6=0;epi7=0;epi8=0;
e1=reshape（b11100）;e2=reshape（b21100）;e3=reshape（b31100）;
e4=reshape（b41100）;e5=reshape（b51100）;e6=reshape（b61100）;
e7=reshape（b71100）;e8=reshape（b81100）;
enl1=mean（e1）^2/std（e1）^2;enl2=mean（e2）^2/std（e2）^2;enl3=mean（e3）^2/std（e3）^2;
enl4=mean（e4）^2/std（e4）^2;enl5=mean（e5）^2/std（e5）^2;enl6=mean（e6）^2/std（e6）^2;
enl7=mean（e7）^2/std（e7）^2;enl8=mean（e8）^2/std（e8）^2;
for i=1:109;
    for j=100:109;
        epi1=epi1+sqrt（（C1（ij）-C1（i+1j））^2+（C1（ij）-C1（ij+1））^2）;
        epi2=epi2+sqrt（（C4（ij）-C4（i+1j））^2+（C4（ij）-C4（ij+1））^2）;
    end
end
for i=20:209;
    for j=200:209;
        epi3=epi1+sqrt（（C2（ij）-C2（i+1j））^2+（C2（ij）-C2（ij+1））^2）;
        epi4=epi2+sqrt（（C5（ij）-C5（i+1j））^2+（C5（ij）-C5（ij+1））^2）;
    end
end
for i=50:59;
    for j=50:59;
        epi5=epi1+sqrt（（C3（ij）-C3（i+1j））^2+（C3（ij）-C3（ij+1））^2）;
        epi6=epi2+sqrt（（C6（ij）-C6（i+1j））^2+（C6（ij）-C6（ij+1））^2）;
    end
end
for i=50:53;
    for j=50:53;
        epi7=epi1+sqrt（（span1（ij）-span1（i+1j））^2+（span1（ij）-span1（ij+1））^2）;
        epi8=epi2+sqrt（（span2（ij）-span2（i+1j））^2+（span2（ij）-span2（ij+1））^2）;
    end
end
epi11=epi2/epi1;epi22=epi4/epi3;epi33=epi6/epi5;epispan=epi8/epi7;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2013-04-29 18:37  lee滤波\
     文件        2406  2012-05-06 15:54  lee滤波\biaozhun.asv
     文件        2406  2012-05-06 15:55  lee滤波\biaozhun.m
     文件          92  2012-06-03 15:41  lee滤波\change.m
     文件       20807  2012-05-04 10:29  lee滤波\dblclass.asv
     文件       20891  2012-06-03 15:32  lee滤波\dblclass.m
     文件        2360  2012-06-03 15:48  lee滤波\dbl_class.mat
     文件       15745  2012-06-03 15:48  lee滤波\dbl_coordinate.mat
     文件       25069  2012-06-03 15:49  lee滤波\image_class.mat
     文件         260  2012-06-03 15:49  lee滤波\lei_distance.mat
     文件       18212  2012-05-04 10:29  lee滤波\oddclass.asv
     文件       18296  2012-06-03 15:36  lee滤波\oddclass.m
     文件        6309  2012-06-03 15:48  lee滤波\odd_class.mat
     文件       65303  2012-06-03 15:48  lee滤波\odd_coordinate.mat
     文件        3553  2012-01-08 14:42  lee滤波\scatterfilter.asv
     文件        3973  2012-06-03 15:40  lee滤波\scatterfilter.m
     文件       20746  2012-05-04 10:39  lee滤波\volclass.asv
     文件       20830  2012-06-03 15:37  lee滤波\volclass.m
     文件        6427  2012-06-03 15:49  lee滤波\vol_class.mat
     文件       47631  2012-06-03 15:49  lee滤波\vol_coordinate.mat
     文件        4431  2012-05-19 09:17  lee滤波\yanfilter.asv
     文件        4439  2012-06-03 15:49  lee滤波\yanfilter.m
     文件         134  2012-05-16 20:29  lee滤波\yanscatter.asv
     文件        3194  2010-11-24 19:03  lee滤波\yanscatter.m
     文件        1341  2012-05-03 21:56  lee滤波\yanscatter1.asv
     文件        1341  2011-12-30 19:59  lee滤波\yanscatter1.m