close all;
clear;
%%______________________________FANTOM ZABURZENIA ARCHITEKTURY  z tłem _______________________________________%%
%% bez tła
img = zeros(512:512);
% / \
Idist1 = drawLine(img, [130 130], [176 176], 1,3);Idist1 = Idist1*0.4056;
Idist2 = drawLine(img, [275 130], [210 190], 1, 3);Idist2 = Idist2*0.4056;
Idist9 = drawLine(img, [310 250], [220 212], 1, 3);Idist9 = Idist9*0.4056;
Idist10 = drawLine(img, [213 195], [290 170], 1, 3);Idist10 = Idist10*0.4056;
Idist11 = drawLine(img, [150 110], [200 180], 1, 3);Idist11 = Idist11*0.4056;
Idist12 = drawLine(img, [120 175], [170 190], 1, 3);Idist12 = Idist12*0.4056;
Idist13 = drawLine(img, [120 245], [187 213], 1, 3);Idist13 = Idist13*0.4056;
Idist14 = drawLine(img, [190 245], [167 302], 1, 3);Idist14 = Idist14*0.4056;
Idist3 = drawLine(img, [217 222], [260 265], 1, 3);Idist3 = Idist3*0.4056;
Idist4 = drawLine(img, [195 215], [120 310], 1, 3);Idist4 = Idist4*0.4056;
% |
Idist5 = drawLine(img, [210 200], [305 200], 1, 3);Idist5 = Idist5*0.4056;
Idist6 = drawLine(img, [120 200], [175 200], 1, 3);Idist6 = Idist6*0.4056;
% -
Idist7 = drawLine(img, [200 120], [200 180], 1, 3);Idist7 = Idist7*0.4056;
Idist8 = drawLine(img, [200 235], [200 300], 1, 3);Idist8 = Idist8*0.4056;

Idist = Idist1+Idist2+Idist3+Idist4+Idist5+Idist6+Idist7+Idist8+Idist9+Idist10+Idist11+Idist12+Idist13+Idist14;
% figure; imagesc(Idist); title('fantom zaburzenia bez tła');colormap(gray);
% zaszumienie fantomu 
% gDist10 = imnoise(Idist,'gaussian',0.10); figure; imagesc(gDist10); colormap(gray); title('zaszumiony gaussian 10%');
% gDist25 = imnoise(Idist,'gaussian',0.25); figure; imagesc(gDist25); colormap(gray); title('zaszumiony gaussian 25%');
%  gDist45 = imnoise(Idist,'gaussian',0.45); figure; imagesc(gDist45); colormap(gray); title('zaszumiony gaussian 45%');
% % 
% % spDist10 = imnoise(Idist,'salt & pepper',0.10);   figure; imagesc(spDist10); colormap(gray); title('zaszumiony salt & paper 10%');
% % spDist25 = imnoise(Idist,'salt & pepper',0.25); figure; imagesc(spDist25); colormap(gray); title('zaszumiony salt & paper 25%');
% spDist45 = imnoise(Idist,'salt & pepper',0.45); figure; imagesc(spDist45); colormap(gray); title('zaszumiony salt & paper 45%');
% 

%% z tłem
Idist111 = drawLine(img, [400 1], [1 1], 1, 3);Idist111 = Idist111*0.4056;
Idist112 = drawLine(img, [390 12], [1 50], 1, 3);Idist112 = Idist112*0.4056;
Idist113 = drawLine(img, [380 3], [1 100], 1, 3);Idist113 = Idist113*0.4056;
Idist114 = drawLine(img, [372 14], [1 150], 1, 3);Idist114 = Idist114*0.4056;
Idist115 = drawLine(img, [389 5], [1 200], 1, 3);Idist115 = Idist115*0.4056;
Idist116 = drawLine(img, [363 16], [1 250], 1, 3);Idist116 = Idist116*0.4056;
Idist117 = drawLine(img, [378 17], [1 300], 1, 3);Idist117 = Idist117*0.4056;
Idist118 = drawLine(img, [391 8], [1 350], 1, 3);Idist118 = Idist118*0.4056;

Idist119 = drawLine(img, [398 7], [1 400], 1, 3);Idist119 = Idist119*0.4056;
Idist120 = drawLine(img, [393 12], [50 400], 1, 3);Idist120 = Idist120*0.4056;
Idist121 = drawLine(img, [394 16], [100 400], 1, 3);Idist121 = Idist121*0.4056;
Idist122 = drawLine(img, [395 22], [150 400], 1, 3);Idist122 = Idist122*0.4056;
Idist123 = drawLine(img, [396 27], [200 400], 1, 3);Idist123 = Idist123*0.4056;
Idist124 = drawLine(img, [397 31], [250 400], 1, 3);Idist124 = Idist124*0.4056;
Idist125 = drawLine(img, [398 38], [300 400], 1, 3);Idist125 = Idist125*0.4056;
Idist126 = drawLine(img, [399 45], [350 400], 1, 3);Idist126 = Idist126*0.4056;
Idist127 = drawLine(img, [400 50], [400 400], 1, 3);Idist127 = Idist127*0.4056;

Itlo = Idist111+Idist112+Idist113+Idist114+Idist115+Idist116+Idist117+Idist118+Idist118+Idist120+Idist121+Idist122+Idist123+Idist124+Idist125+Idist126+Idist127;
% I =  Itlo;
I = Idist + Itlo;
% figure; imagesc(I); title('fantom zaburzenia z tłem');colormap(gray);
% zaszumiony 
% gDist10 = imnoise(I,'gaussian',0.10); %figure; imagesc(gDist10); colormap(gray); title('zaszumiony gaussian 10% z tłem');
% gDist25 = imnoise(I,'gaussian',0.25); %figure; imagesc(gDist25); colormap(gray); title('zaszumiony gaussian 25% z tłem');
gDist45 = imnoise(I,'gaussian',0.45); %figure; imagesc(gDist45); colormap(gray); title('zaszumiony gaussian 45% z tłem');
% gDist65 = imnoise(I,'gaussian',0.65);
% 
% % 
% spDist10 = imnoise(I,'salt & pepper',0.10);%   figure; imagesc(spDist10); colormap(gray); title('zaszumiony salt & paper 10% z tłem');
% spDist25 = imnoise(I,'salt & pepper',0.25); %figure; imagesc(spDist25); colormap(gray); title('zaszumiony salt & paper 25% z tłem');
% spDist45 = imnoise(I,'salt & pepper',0.45); %figure; imagesc(spDist45); colormap(gray); title('zaszumiony salt & paper 45% z tłem');