matlab小波变换图像去噪

小波变换可以用于图像去噪。以下是使用MATLAB进行小波去噪的简单步骤：

1. 读取图像并将其转换为灰度图像。

img = imread('image.jpg'); % 读取图像
grayImg = rgb2gray(img); % 将图像转换为灰度图像

2. 对图像进行小波变换。

[cA,cH,cV,cD] = dwt2(grayImg,'haar'); % 对图像进行2D小波变换

其中，`cA`是近似系数，`cH`是水平细节系数，`cV`是垂直细节系数，`cD`是对角细节系数。

3. 对小波系数进行阈值处理。

% 设定阈值
sigma = median(abs(cD(:))) / 0.6745;
threshold = sigma * sqrt(2*log10(numel(cD)));

% 阈值处理
cDthresh = cD .* (abs(cD) > threshold);
cHthresh = cH .* (abs(cH) > threshold);
cVthresh = cV .* (abs(cV) > threshold);

这里使用的是基于中值绝对偏差（MAD）的阈值处理方法。

4. 对处理后的小波系数进行反变换并显示图像。

% 反变换
denoisedImg = idwt2(cA,cHthresh,cVthresh,cDthresh,'haar');

% 显示图像
imshow(denoisedImg);

这里使用的是Haar小波基。

完整代码如下：

% 读取图像并将其转换为灰度图像
img = imread('image.jpg');
grayImg = rgb2gray(img);

% 对图像进行小波变换
[cA,cH,cV,cD] = dwt2(grayImg,'haar');

% 对小波系数进行阈值处理
sigma = median(abs(cD(:))) / 0.6745;
threshold = sigma * sqrt(2*log10(numel(cD)));
cDthresh = cD .* (abs(cD) > threshold);
cHthresh = cH .* (abs(cH) > threshold);
cVthresh = cV .* (abs(cV) > threshold);

% 对处理后的小波系数进行反变换并显示图像
denoisedImg = idwt2(cA,cHthresh,cVthresh,cDthresh,'haar');
imshow(denoisedImg);