数字图像处理实验：椭圆绘制与矩阵操作.pdfVIP

%画椭圆清

除命令窗口

全部;clc;

n=256;I

ones(n,n);

xrepmat((1:n),n,1);y

=repmat((1:n),1,n);

z=((x‑n/2)./80).^2+((y‑n/2)./100).^2;

I((z1.1)(z0.9))=0;

显示图像(I);

%实验一(1)clearallI=ones(256,256);I(78:178,58:198)=0.6;对于m=78:178对于

n=58:198如果(((m-128)/50)^2+((n-128)/70)^2)=1.1(((m-128)/50)^2+

((n-128)/70)^2)=0.9I(m,n)=0.2;结束结束

结束显示图像(I)N=100;F=fft2(I,N,N);F2

fftshift(abs(F));创建图形窗口2显示对数

变换后的F2图像创建图形窗口3x=1:N;y=1:

N;绘制三维网格图;设置灰度颜色映射;显示

颜色条

%实验一(2)

clearall

I=imread(LENA128.bmp);

N=100;

%画椭圆

clearall;

clc;

n=256;

I=ones(n,n);

x=repmat((1:n),n,1);

y=repmat((1:n),1,n);

z((x-n/2)./80).^2+((y-n/2)./100).^2;

I((z1.1)(z0.9))=0;

imshow(I);

%实验一(1)

clearall

I=ones(256,256);

I(78:178,58:198)=0.6;

form=78:178

forn=58:198

if(((m-128)/50)^2+((n-128)/70)^2)=1.1(((m-128)/50)^2+((n-

128)/70)^2)=0.9

I(m,n)=0.2;

end

end

end

imshow(I)

N=100;

F=fft2(I,N,N);

F2=fftshift(abs(F));

figure(2)

imshow(log(F2))

figure(3)

x=1:N;y=1:N;

mesh(x,y,F2(x,y));colormap(gray);colorbar

%实验一(2)

clearall

I=imread(LENA128.bmp);

N=100;

I1=imrotate(I,45);F1=fft2(I1,N,N)./(length(I1))

^2;M1fftshift(abs(F1));subplot(2,2,1),显

示图像