基于小波变换的边缘检测matlab程序

基于小波变换的边缘检测matlab程序，亲调可用......................................................................................................................................................

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