小波多尺度边缘检测主程序

%小波多尺度边缘检测主程序
close all
clc
I = imread（‘cameraman.tif‘）;
[rowcol]=size（I）;

%使用的小波函数
wname=‘bior4.4‘;

% 显示原图
subplot（231）;
imshow（I）;
title（‘原图‘）;

%canny算法
bw= edge（I‘canny‘）;
subplot（232）;
imshow（bw）;
title（‘Canny法‘）;

I = double（I）;

% 对图像进行三层小波变换
[ca1ch1cv1cd1] = dwt2（Iwname）;
[ca2ch2cv2cd2] = dwt2（ca1wname）;
[ca3ch3cv3cd3] = dwt2（ca2wname）;

% 各层小波系数求极值点及极值点的梯度值（角度）；
thr=4;a0=1;n=7;
[edge_mf1grads1] = local_max_mode（cv1ch1thra0n）;

thr=15;a0=1;n=5;
[edge_mf2grads2] = local_max_mode（cv2ch2thra0n）;

thr=30;a0=1;n=3;
[edge_mf3grads3] = local_max_mode（cv3ch3thra0n）;

% 对极值点矩阵归一化
edge_mf1=guiyi（edge_mf1）;
edge_mf2=guiyi（edge_mf2）;