Лабораторная работа по "Программированию"
Автор: guhf57 • Ноябрь 9, 2022 • Лабораторная работа • 492 Слов (2 Страниц) • 152 Просмотры
Лабораторная работа 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);
...