clear?all;
load?wbarb;
%?I?=?ind2gray（Xmap）;
RGB=imread（‘lena.jpg‘）;
I=rgb2gray（RGB）;
imshow（I）;
I1?=?imadjust（Istretchlim（I）[01]）;
figure;
imshow（I1）;
[NM]?=?size（I）;
%??h?=?[0.1250.3750.3750.125];
%??g?=?[0.5-0.5];
h?=?[0.050.10.10.05];
g?=?[0.25-0.25];
delta?=?[100];?
J?=?3;
a（1:N1:M11:J+1）?=?0;
dx（1:N1:M11:J+1）?=?0;
dy（1:N1:M11:J+1）?=?0;
d（1:N1:M11:J+1）?=?0;?
a（::11）?=?conv2（hhI‘same‘）;
dx（::11）?=?conv2（deltagI‘same‘）;
dy（::11）?=?conv2（gdeltaI‘same‘）;
x?=?dx（::11）;
y?=?dy（::11）;
d（::11）?=?sqrt（x.^2+y.^2）;
I1?=?imadjust（d（::11）stretchlim（d（::11））[0?1]）;
figure;
imshow（I1）;
lh?=?length（h）;
lg?=?length（g）;
for?j?=?1:J+1
??lhj?=?2^j*（lh-1）+1;
??lgj?=?2^j*（lg-1）+1;
??hj（1:lhj）=0;
??gj（1:lgj）=0;
??for?n?=?1:lh
????hj（2^j*（n-1）+1）=h（n）;
??end
??for?n?=?1:lg
????gj（2^j*（n-1）+1）=g（n）;
??end
??
??a（::1j+1）?=?conv2（hjhja（::1j）‘same‘）;
??dx（::1j+1）?=?conv2（deltagja（::1j）‘same‘）;
??dy（::1j+1）?=?conv2（gjdeltaa（::1j）‘same‘）;
?
??x?=?dx（::1j+1）;
??y?=?dy（::1j+1）;
??dj（::1j+1）?=?sqrt（x.^2+y.^2）;
?
??I1?=?imadjust（dj（::1j+1）stretchlim（dj（::1j+1））[0?1]）;
??figure;
??imshow（I1）;
end??

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????1247??2017-09-06?09:49??XBBH.m

-----------?---------??----------?-----??----

?????????????????1247????????????????????1