I=imread('heykel.png');
figure,imshow(I)
H=I(250,:);
I(250,:)=255;
[h,w,ch]=size(I);
figure,plot(1:800,H)

x=1:100;
y=x.^2;
figure,plot(x,y)

noiseI=imnoise(I,'salt & pepper', 0.02);
figure,imshow(noiseI)



F=fspecial('average',10);
F
size(F)
If=conv2(I,F);
If=conv2(double(I),F);
figure,imshow(If)
figure,imshow(uint8(If))

G=fspecial('gaussian',50,2);
figure,imagesc(G)
G=fspecial('gaussian',50,5);
figure,imagesc(G)
surf(G)
figure,surf(G)
close all


I=imread('heykel.png');
H1=I(250,:);
 I(250,:)=255;
 figure,plot(1:800,H1)
 F=fspecial('average',11);
If=conv2(double(I),F);
If=If(6:h-5,6:w-5);
H2=If(250,:);
 figure,plot(1:790,H2)
[h,w]=size(I);
nI=I;
for i=1:10000
    
    xh=randsample(h,1);
    yw=randsample(w,1);
    nI(xh,yw)=0;
end

figure,imshow(I)
% fI=zeros(h,w);
fI=nI;
radius=1;
for i=1:1:h
    [i h]
    for j=1:1:w
        
        if nI(i,j)==0
        window=nI(max(i-radius,1):min(i+radius,h),max(j-radius,1):min(j+radius,w));
        fI(i,j)=max(window(:));
        end
    end
end