vid=videoinput('winvideo')
for  i=1:400,
obr = getsnapshot(vid); obrazszary=rgb2gray(obr); figure(1); subplot(3,3,1);  imshow(obrazszary); title('Obraz kamery');
obraz=imresize(obrazszary,0.2);
obr_sp=imnoise(obraz,'salt & pepper');
obr_g=imnoise(obraz,'gaussian');
obr_p=pasy(obraz, 1, 30,'poziome');

%filtracja uśredniająca
h_usr3=fspecial('average',3), h_usr7=fspecial('average',7),
h_g3=fspecial('gaussian',3), h_usr7=fspecial('gaussian',7),
%uśredniajaca
obr_usr_sp=imfilter(obr_sp,h_usr3),
obr_usr_g=imfilter(obr_g,h_usr3),
obr_usr_p=imfilter(obr_p,h_usr3),
obr_usr_sp=imfilter(obr_sp,h_usr7),
obr_usr_g=imfilter(obr_g,h_usr7),
obr_usr_p=imfilter(obr_p,h_usr7),

%gaus
obr_g_sp=imfilter(obr_sp,h_usr3),
obr_g_g=imfilter(obr_g,h_usr3),
obr_g_p=imfilter(obr_p,h_usr3),
obr_g_sp=imfilter(obr_sp,h_usr7),
obr_g_g=imfilter(obr_g,h_usr7),
obr_g_p=imfilter(obr_p,h_usr7),

%medianowa

obr_med_sp=medfilt2(obr_sp,[3 3]),
obr_med_g=medfilt2(obr_g,[3 3]),
obr_med_p=medfilt2(obr_p,[3 3]),
obr_med_sp=medfilt2(obr_sp,[7 7]),
obr_med_g=medfilt2(obr_g,[7 7]),
obr_med_p=medfilt2(obr_p,[7 7]),
% detekcja krawędzi

h_sobel1=fspecial('sobel'), h_sobel2=-h_sobel1;
subplot(3,3,2);
kraw_sob1=filter2(h_sobel1,obraz), imshow(mat2gray(kraw_sob1)); title('Detekcja krawedzi ');
subplot(3,3,3);
kraw_sob2=filter2(h_sobel2,obraz), imshow(mat2gray(kraw_sob2)); title('Detekcja krawedzi 1');
subplot(3,3,4);
kraw_sob_e=sqrt(kraw_sob1.^2 + kraw_sob2.^2); imshow(mat2gray(kraw_sob_e)); title('Detekcja krawedzi 2');
subplot(3,3,5);
kraw_sob_m=abs(kraw_sob1) + abs(kraw_sob2); imshow(mat2gray(kraw_sob_m));
title('Detekcja krawedzi medianowa');

%laplasian
h_lap=fspecial('laplacian');
kraw_lap=filter2(h_lap,obraz);
subplot(3,3,6);
imshow(mat2gray(kraw_lap));
title('Filtr laplasian');
%filtr log

h_log=fspecial('log');
kraw_log=filter2(h_log,obraz);
subplot(3,3,7);
imshow(mat2gray(kraw_log));
title('filtr LOG');
end