matlab包络检波用椭圆滤波器的频谱代码.pdf

matlab包络检波用椭圆滤波器的频谱代码

clearall;clc;

�T变换

I=imread(cameraman.tif);

A1=fft2(I);

A2=ifft2(A1);

A3=uint8(A2);

figure,subplot(221),imshow(I);

subplot(222);imshow(A3);

�T变换

I=imread(cameraman.tif);

figure,subplot(223);imshow(I);

I=im2double(I);

T=dctmtx(8);

B=blkproc(I,[88],P1*x*P2,T,T);

Mask=

;

B2=blkproc(B,[8

8],P1.*x,Mask);%此处为点乘(.*)

I2=blkproc(B2,[88],P1*x*P2,T,T);

subplot(224);imshow(I2);

%直方图均衡

I=imread(pout.tif);

figure,subplot(321);imshow(I);

subplot(322);imhist(I);[J,T]=histeq(I,64);%图像灰度扩展

到0~255，但是只有64个灰度级

subplot(323);imshow(J);

subplot(324);imhist(J);

J=histeq(I,32);

subplot(325);imshow(J);%

图像灰度扩展到0~255，但是只有32个灰度级

subplot(326);imhist(J);

figure,plot((0:255)/255,T);%转移函数的变换曲线

%灰度拉伸

I=imread(pout.tif);

figure,subplot(221);imshow(I);

subplot(222);imhist(I);

[X,map]=imread(pout.tif);

Y=imadjust(X,[0.10.99],[01]);

subplot(223);imshow(Y,map);

subplot(224);imhist(Y);

%图像空间域处理

[I,map]=imread(pout.tif);

J1=imnoise(I,gaussian,0,0.02);

%叠加均值为0，方差为0.02的高斯噪声，可以用localvar代

替

%J2=imnoise(I,saltpepper,0.04);%

叠加密度为0.04的椒盐噪声

M4=[010;101;010];

M4=M4/4;%4邻域平均滤波

I_filter1=filter2(M4,J1);

figure,subplot(321);imshow(I);

subplot(322);imshow(J1);

subplot(323);imshow(I_filter1,map);

M8=[111;101;11

1];%8邻域平均滤波

M8=M8/8;

I_filter2=filter2(M8,J1);

subplot(324);imshow(I_filter2,map);

K1=medfilt2(J1,[3,3]);

subplot(325);imshow(K1);

K2=medfilt2(J2,[77]);

subplot(326);imshow(K2);

%频域处理

I=imread(pout.tif);

figure,subplot(221);imshow(I);

J1=imnoise(I,salt

pepper);%叠加椒盐噪声

subplot(222);imshow(J1);

f=double