Wizja lab3

info=imaqhwinfo
vid=videoinput('winvideo')
%preview(vid)
get(vid)
get(getselectedsource(vid))
set(vid, 'TriggerRepeat', Inf)
set(gcf, 'doublebuffer', 'on')
start(vid)
figure(1);
while(vid.FramesAcquired <=10)
    data=getsnapshot(vid);
    obraz=rgb2gray(data);
    subplot(2,2,1);
    imshow(obraz);

end
 subplot(2,2,2);
 imhist(obraz);
 subplot(2,2,3);
 histeq(obraz);
 %level = graythresh(obraz);
 %bin = im2bw(obraz,level);
 %subplot(2,2,4);
 %imshow(bin);
 figure(2);
 obr_g=imnoise(obraz,'gaussian');
 subplot(3,3,1);
 imshow(obr_g);
 hg3=fspecial('gaussian',[3 3]);
 hg4=fspecial('gaussian',[7 7]);
 obr_g2=imfilter(obr_g,hg3);
 obr_g3=imfilter(obr_g,hg4);
 subplot(3,3,2);
 imshow(obr_g2);
 subplot(3,3,3);
 imshow(obr_g3);
 hg5=fspecial('average',[3 3]);
 obr_p2=imfilter(obr_g,hg4);
 hg6=fspecial('average',[7 7]);
 obr_p3=imfilter(obr_g,hg5);
 subplot(3,3,5);
 imshow(obr_p2);
 subplot(3,3,6);
 imshow(obr_p3);
 obr_sp2=medfilt2(obr_g,[3 3]);
 obr_sp3=medfilt2(obr_g,[7 7]);
 subplot(3,3,8);
 imshow(obr_sp2);
 subplot(3,3,9);
 imshow(obr_sp3);

 figure(3);
 obr_p=pasy(obraz,1,20,'poziome');
 subplot(3,3,1);
 imshow(obr_p);
 hg3=fspecial('gaussian',[3 3]);
 hg4=fspecial('gaussian',[7 7]);
 obr_p2=imfilter(obr_p,hg3);
 obr_p3=imfilter(obr_p,hg4);
 subplot(3,3,2);
 imshow(obr_p2);
 subplot(3,3,3);
 imshow(obr_p3);
 hg5=fspecial('average',[3 3]);
 obr_p2=imfilter(obr_p,hg4);
 hg6=fspecial('average',[7 7]);
 obr_p3=imfilter(obr_p,hg5);
 subplot(3,3,5);
 imshow(obr_p2);
 subplot(3,3,6);
 imshow(obr_p3);
 obr_sp2=medfilt2(obr_p,[3 3]);
 obr_sp3=medfilt2(obr_p,[7 7]);
 subplot(3,3,8);
 imshow(obr_sp2);
 subplot(3,3,9);
 imshow(obr_sp3);

 figure(4);
 obr_sp=imnoise(obraz,'salt & pepper');
 subplot(3,3,1);
 imshow(obr_sp);
 hg3=fspecial('gaussian',[3 3]);
 hg4=fspecial('gaussian',[7 7]);
 obr_p2=imfilter(obr_sp,hg3);
 obr_p3=imfilter(obr_sp,hg4);
 subplot(3,3,2);
 imshow(obr_p2);
 subplot(3,3,3);
 imshow(obr_p3);
 hg5=fspecial('average',[3 3]);
 obr_p2=imfilter(obr_sp,hg4);
 hg6=fspecial('average',[7 7]);
 obr_p3=imfilter(obr_sp,hg5);
 subplot(3,3,5);
 imshow(obr_p2);
 subplot(3,3,6);
 imshow(obr_p3);
 obr_sp2=medfilt2(obr_sp,[3 3]);
 obr_sp3=medfilt2(obr_sp,[7 7]);
 subplot(3,3,8);
 imshow(obr_sp2);
 subplot(3,3,9);
 imshow(obr_sp3);

 %operator sobela
 figure(5);
 h_sob1=fspecial('sobel')
 h_sob2=-h_sob1'
 kraw_sob1=filter2(h_sob1,obraz);
 subplot(2,2,1);
 imshow(mat2gray(kraw_sob1))
 kraw_sob2=filter2(h_sob2,obraz);
 subplot(2,2,2);
 imshow(mat2gray(kraw_sob2))
 kraw_sob_c=sqrt(kraw_sob1.^2+kraw_sob2.^2);
 subplot(2,2,3);
 imshow(mat2gray(kraw_sob_c))
 kraw_sob_m=abs(kraw_sob1)+abs(kraw_sob2);
 subplot(2,2,4);
 imshow(mat2gray(kraw_sob_m))

 %operator laplace'a(są jednokierunkowe, nie potrzeba drugiego wektora)
 h_lap=fspecial('laplacian')

 %operator LOG
 h_log=fspecial('log')

stop(vid)
delete(vid)
clear

close(gfc)