% Roy test

path(path,'/home/S00/schestr0/NRR_Smith');

%% Section 1

[a,bump_images,c,d,e]=make_1d_images(2,50,0.1);

figure(1);
hold on;
title('Changing bump width');
for i=1:2,
  plot(bump_images(:,i),'R');
end

[a,bump_images,c,d,e]=make_1d_images(2,50,0.3);

for i=1:2,
  plot(bump_images(:,i),'B');
end

[a,bump_images,c,d,e]=make_1d_images(2,50,0.5);

for i=1:2,
  plot(bump_images(:,i),'G');
end
hold off;


%% Section 2

[a,bump_images,c,d,e]=make_1d_images(2,10,0.2);

figure(2);
hold on;
title('Changing image width');
for i=1:2,
  plot(bump_images(:,i),'R');
end

[a,bump_images,c,d,e]=make_1d_images(2,50,0.2);

for i=1:2,
  plot(bump_images(:,i),'B');
end

[a,bump_images,c,d,e]=make_1d_images(2,200,0.2);

for i=1:2,
  plot(bump_images(:,i),'G');
end
hold off;

%% Section 3

[a,bump_images,c,d,e]=make_1d_images(4,100,0.3);

average_smoothed_bump_images = average_smooth(bump_images, 9);

figure(3);
hold on;
title('Average smooth effect');
for i=1:4,
  plot(average_smoothed_bump_images(:,i),'R');
end
hold off;
average_smoothed_bump_images = average_smooth(average_smoothed_bump_images, 9);

figure(4);
hold on;
title('Average smooth of order 2 effect');
for i=1:4,
  plot(average_smoothed_bump_images(:,i),'R');
end
hold off;

%% Section 4

gaussian_smoothed_bump_images = gaussian_smooth(bump_images, 40);

figure(5);
hold on;
title('Gaussian smooth effect');
for i=1:4,
  plot(gaussian_smoothed_bump_images(:,i),'R');
end
hold off;
