图像的空间域转换为频率域

load woman; figure, imshow(X,map); title('原图'); nbcol = size(map,1); [cA1,cH1,cV1,cD1] = dwt2(X,'db1'); cod_X = wcodemat(X,nbcol); cod_cA1 = wcodemat(cA1,nbcol); cod_cH1 = wcodemat(cH1,nbcol); cod_cV1 = wcodemat(cV1,nbcol); cod_cD1 = wcodemat(cD1,nbcol); dec2d = [... cod_cA1, cod_cH1; ... cod_cV1, cod_cD1 ... ]; figure; imshow(dec2d,map); title('频率分量分布图'); for i=1:128 for j=1:128 cod_cH1(i,j)=0; end end % for i=1:128 %test % for j=1:128 % cod_cV1(i,j)=0; % end % end % for i=1:128 % for j=1:128 % cod_cD1(i,j)=0; % end % end % for i=1:128 % for j=1:128 % cod_cA1(i,j)=0; % end % end x=idwt2(cod_cA1, cod_cH1,cod_cV1, cod_cD1 ,'db1'); figure, imshow(x,map); title('去掉其中的一个高频');