function?[ESMmaxEDM]=gauss_gamma（I）
%made?by?白國(guó)樂(lè)，2017/6/20
%雙窗函數(shù)濾波，計(jì)算sar圖像的梯度強(qiáng)度和梯度方向
[mn]=size（I）;
a=2;????????????????????????????????????????%a是雙窗間距
c=2;????????????????????????????????????????%c是gauss向窗口長(zhǎng)度參數(shù)??x向，橫向
b=2;????????????????????????????????????????%b是gamma向窗口寬度參數(shù)??y向，豎向
global?p?q;?????????????????????????????
p=5;q=5;????????????????????????????????????%pq分別是單窗的行數(shù)和列數(shù)
%%%%%%%%%%%%%%%%%%%%%%%%%%%%????W01W02
W0=[];
??for?x=-p:1:p
????for?y=-q:1:q
?????W=abs（y）^（a-1）/（sqrt（2*pi）*c*gamma（a）*b^a）*exp（-（x^2/（2*c^2）+abs（y）/b））;??
?????W0=[W0W];?????????????????????????????????????????????????????????????
????end
???end
W0=reshape（W02*p+12*q+1）;?W01=W0;W02=W0;?
??for?i=p+1:2*p+1
??????for?j=1:2*q+1
????????W01（ij）=0;
??????end
??end

??for?i=1:p+1
?????for?j=1:2*q+1
?????????W02（ij）=0;
?????end
??end

%%%%%%%%%%%%%%%%%%%%%%%%%%%???W11W12
W1=[];
??for?x=-p:1:p
????for?y=-q:1:q
??????x1=x*cos（pi/4）-y*sin（pi/4）;y1=x*sin（pi/4）+y*cos（pi/4）;
?????W=abs（y1）^（a-1）/（sqrt（2*pi）*c*gamma（a）*b^a）*exp（-（x1^2/（2*c^2）+abs（y1）/b））;??
?????W1=[W1W];?????????????????????????????????????????????????????????????
????end
???end
W1=reshape（W12*p+12*q+1）;?W11=W1;W12=W1;?
??for?i=1:2*p+1
????for?j=1:2*q+1
??????if?i+j>=p+q+2
????????W11（ij）=0;
??????end
????end
??end

??for?i=1:2*p+1
?????for?j=1:2*q+1
????????if?i+j<=p+q+2
?????????W12（ij）=0;
????????end
?????end
??end

%%%%%%%%%%%%%%%%%%%%%%%%%%%????????W21W22
W2=[];
???for?x=-p:1:p
????for?y=-q:1:q
??????x2=x*cos（pi/2）-y*sin（pi/2）;y2=x*sin（pi/2）+y*cos（pi/2）;
?????W=abs（y2）^（a-1）/（sqrt（2*pi）*c*gamma（a）*b^a）*exp（-（x2^2/（2*c^2）+abs（y2）/b））;??
?????W2=[W2W];?????????????????????????????????????????????????????????????
????end
???end
W2=reshape（W22*p+12*q+1）;?W21=W2;W22=W2;?
??for?i=1:2*p+1
??????for?j=q+1:2*q+1
????????W21（ij）=0;
??????end
??end

??for?i=1:2*p+1
?????for?j=1:q+1
?????????W22（ij）=0;
?????end
??end

%%%%%%%%%%%%%%%%%%%%%%%%%????????W31W32
W3=[];
???for?x=-p:1:p
??????for?y=-q:1:q
????????x3=x*cos（3*pi/4）-y*sin（3*pi/4）;y3=x*sin（3*pi/4）+y*cos（3*pi/4）;
???????W=abs（y3）^（a-1）/（sqrt（2*pi）*c*gamma（a）*b^a）*exp（-（x3^2/（2*c^2）+abs（y3）/b））;??
???????W3=[W3W];?????????????????????????????????????????????????????????????
??????end
???end
W3=reshape（W32*p+12*q+1）;?W31=W3;W32=W3;?
??for?i=1:2*p+1
??????for?j=1:2*q+1
???????if?i<=j
????????W31（ij）=0;
???????end
??????end
???end
?
??for?i=1:2*p+1
?????for?j=1:2*q+1
???????if?i>=j
?????????W32（ij）=0;
???????end
?????end
??end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%???W41W42
W4=[];
??for?x=-p:1:p
????for?y=-q:1:q
??????x4=x*cos（pi）-y*sin（pi）;y4=x*sin（pi）+y*cos（pi）;
?????W=abs（y4）^（a-1）/（sqrt（2*pi）*c*gamma（a）*b^a）*exp（-（x4^2/（2*c^2）+abs（y4）/b））;??
?????W4=[W4W];?????????????????????????????????????????????????????????????
????end
???end
W4=reshape（W42*p+12*q+1）;?W41=W4;W42=W4;?
??for?i=p+1:2*p+1
??????for?j=1:2*q+1
????????W41（ij）=

?屬性????????????大小?????日期????時(shí)間???名稱(chēng)
-----------?---------??----------?-----??----
?????文件????????6904??2018-03-13?16:25??gauss_gamma.m
?????文件?????????221??2017-06-20?14:31??main.m
?????文件??????172632??2017-06-07?15:50??邊緣檢測(cè)test.png