Лабораторная работа по "Программированию"

Лабораторная работа 2

Задание 1

Image_2D=imread('lena512.bmp');

Bloc_Size_Y=270;

Bloc_Size_X=200;

Start_Y=150;

Start_X=190;

close all;

%Image_2D_=(Image_2D);

Image_2D_=double(Image_2D);

[siz_y, siz_x]=size(Image_2D_);

Bloc_=zeros(Bloc_Size_Y, Bloc_Size_X);

y=1;

while (y<=siz_y)&&(y<=Bloc_Size_Y+Start_Y)

    x=1;

    while (x<=siz_x)&&(x<=Bloc_Size_X+Start_X)

        if y>=Start_Y && x>=Start_X

            Bloc_(y-Start_Y+1, x-Start_X+1)=Image_2D_(y, x);

        end

        x=x+1;

    end

    y=y+1;

end

%figure, imshow(Image_2D,256);

figure, imshow(Image_2D);

%figure, hist(Image_2D); % Для GNU Octave (CLI)

figure, imhist(Image_2D,256); % Для MATLAB R2021a

saveas(gcf,'lena512_Histogram.jpg','jpg')

Bloc=uint8(Bloc_); imwrite(Bloc,'Bloc.bmp','bmp');

%figure, imshow(Bloc,256);

figure, imshow(Bloc);

%figure, hist(Bloc);   % Для GNU Octave (CLI)

%saveas(gcf,'Bloc_Histogram.jpg','jpg')   % Для GNU Octave (CLI)

figure, imhist(Bloc,256);  % Для MATLAB R2021a

saveas(gcf,'Bloc_Histogram.jpg','jpg')    % Для MATLAB R2021a

[pic 1][pic 2]

[pic 3][pic 4]

Рисунок 1 – Выполнение задания 1

Задание 2

close all;

N=300;

x=zeros(N,N);

y1=N/6;

y2=N-y1;

x1=N/3;

x2=N-x1;

x(y1:y2,x1:x2) = 1;

%x(5:24,13:17) = 1;

%figure, imshow(x,'notruesize')

figure, imshow(x)

%

F = fft2(x);

F2 = log(abs(F));

figure,

%imshow(F2,[-1 5],'notruesize');

imshow(F2,[-1 5]);

%imshow(F2,'DisplayRange',[-1 5]);

colormap(jet); colorbar

%

F = fft2(x,256,256);

figure, imshow(log(abs(F)),[-1 5]); colormap(jet); colorbar

%

F = fft2(x,256,256);

F2 = fftshift(F);

figure, imshow(log(abs(F2)),[-1 5]); colormap(jet); colorbar

[pic 5][pic 6]

[pic 7][pic 8]

Рисунок 2 – Выполнение задания 2

Задание 3

TH=500;

%TH=1000;

N=64;

%

close all

[f, map] = imread('frame_64.bmp'); aaa(:,:)=f(:,:);

figure, imshow(f,map);

%

r=double(f);

f2=fft2(r,N,N); J_(:,:)=f2(:,:);

f3=fftshift(f2);

f4=abs(f3);

f5 = log(1+f3.*conj(f3)/N); % Power spectrum

%

umin=(min(min(f5))); umax=(max(max(f5)));

as=(f5-umin)/(umax-umin);

f1_=uint8(255*as);

figure,

imshow(f1_);

%imshow(f1_,256);