matlab求给代码注释，会多少给多少就行

Question

clear all;
load wbarb; %小波变换边缘提取程序
I = ind2gray(X,map);%检索图转成灰度图
imshow(I);
I1 = imadjust(I,stretchlim(I),[0,1]);%调整图像的像素值，可以改变对比度和颜色
figure;
imshow(I1);
[N,M] = size(I); %m：图像的宽度 n：图像的高度m
h = [0.125,0.375,0.375,0.125];
g = [0.5,-0.5];
delta = [1,0,0]; %差值
J = 3;
a(1:N,1:M,1,1:J+1) = 0;
dx(1:N,1:M,1,1:J+1) = 0;
dy(1:N,1:M,1,1:J+1) = 0;
d(1:N,1:M,1,1:J+1) = 0;
a(:,:,1,1) = conv2(h,h,I,'same'); %二维卷积
dx(:,:,1,1) = conv2(delta,g,I,'same');
dy(:,:,1,1) = conv2(g,delta,I,'same');
x = dx(:,:,1,1);
y = dy(:,:,1,1);
d(:,:,1,1) = sqrt(x.^2+y.^2);
I1 = imadjust(d(:,:,1,1),stretchlim(d(:,:,1,1)),[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(:,:,1,j+1) = conv2(hj,hj,a(:,:,1,j),'same');
dx(:,:,1,j+1) = conv2(delta,gj,a(:,:,1,j),'same');
dy(:,:,1,j+1) = conv2(gj,delta,a(:,:,1,j),'same');
x = dx(:,:,1,j+1);
y = dy(:,:,1,j+1);
dj(:,:,1,j+1) = sqrt(x.^2+y.^2);
I1 = imadjust(dj(:,:,1,j+1),stretchlim(dj(:,:,1,j+1)),[0 1]);figure;imshow(I1);

Answer 1

clear all; //清除以前历史的数据
load wbarb;%波变换边缘提取程序
I = ind2gray(X,map);%检索图转成灰度图
imshow(I); %显示转换后的灰度图
I1 = imadjust(I,stretchlim(I),[0,1]);%调整图像的像素值，可以改变对比度和颜色
figure;
imshow(I1); %显示调整后的图像
[N,M] = size(I); %m：图像的宽度 n：图像的高度m
h = [0.125,0.375,0.375,0.125]; %h和g可能是小波变换的系数设置，自己查询是不是
g = [0.5,-0.5];
delta = [1,0,0]; %差值
J = 3;
a(1:N,1:M,1,1:J+1) = 0; %a,d,dx,dy设置出四个矩阵（MATLAB里面图像就是矩阵形式）
dx(1:N,1:M,1,1:J+1) = 0;
dy(1:N,1:M,1,1:J+1) = 0;
d(1:N,1:M,1,1:J+1) = 0;
a(:,:,1,1) = conv2(h,h,I,'same'); %二维卷积，下同（图像是2维数据）。注意括号里面的I相同，其他参数不同，看前面设置
dx(:,:,1,1) = conv2(delta,g,I,'same');
dy(:,:,1,1) = conv2(g,delta,I,'same');
x = dx(:,:,1,1); %二维卷积之后，赋给X,Y
y = dy(:,:,1,1);
d(:,:,1,1) = sqrt(x.^2+y.^2); %X,Y的平方和，再开根号，赋给d.
I1 = imadjust(d(:,:,1,1),stretchlim(d(:,:,1,1)),[0 1]);figure;imshow(I1); %调整图像后显示
lh = length(h); %获取图像的宽度和高度
lg = length(g);
for j = 1:J+1 %for循环，可能用于求像素值？不太确定
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(:,:,1,j+1) = conv2(hj,hj,a(:,:,1,j),'same'); %2维卷积
dx(:,:,1,j+1) = conv2(delta,gj,a(:,:,1,j),'same');
dy(:,:,1,j+1) = conv2(gj,delta,a(:,:,1,j),'same');
x = dx(:,:,1,j+1);
y = dy(:,:,1,j+1);
dj(:,:,1,j+1) = sqrt(x.^2+y.^2);
I1 = imadjust(dj(:,:,1,j+1),stretchlim(dj(:,:,1,j+1)),[0 1]);figure;imshow(I1);

追问

您好，我还想问一下，这个程序的小波基函数是什么呀，结果七张图片都是什么意思

追答

小波变换你可以去搜下，具体的我也不是知道太多。