clear all
close all
clc

% I00 = imread("0mm.png");
% I01 = imread("1mm.png");
% I05 = imread("5mm.png");
% I10 = imread("5mm.png");

[I00p BW00 D00 IDX00] =Improcess_V1("0mm.tif");
[I01p BW01 D01 IDX01]=Improcess_V1("1mm.tif");
[I05p BW05 D05 IDX05]=Improcess_V1("5mm.tif");
[I10p BW10 D10 IDX10]=Improcess_V1("10mm.tif");


N = 8; %number of contours
% figure(1)
% imshow((I00_8b))
% 
% figure(2)
% imshow(I00_8b_mean)
% 
% figure(3)
% imshow(I00_8b_med )
% 
% figure(4)
% imshow(I00_8b_med_mean )

figure(1)
subplot(1,4,1)
imshow(I00p)
subplot(1,4,2)
imshow(I01p)
subplot(1,4,3)
imshow(I05p)
subplot(1,4,4)
imshow(I10p)
% colormap(hot)


% figure(2)
% subplot(1,4,1)
% contourf(flip(I00p),N, 'square' )
% subplot(1,4,2)
% % imcontour(I01p,N )
% contourf(flip(I01p),N, 'square' )
% subplot(2,2,3)
% contourf(flip(I05p),N, 'square' )
% % imcontour(I05p,N )
% subplot(2,2,4)
% contourf(flip(I10p),N, 'square' )
% % imcontour(I10p,N )
% 

figure(3);
subplot(1,4,1)
imshow(BW00)
subplot(1,4,2)
imshow(BW01)
subplot(1,4,3)
imshow(BW05)
subplot(1,4,4)
imshow(BW10)

figure(4);
subplot(1,4,1)
imshow(D00, [0,40])
colorbar

subplot(1,4,2)
imshow(D01, [0,40])
colorbar

subplot(1,4,3)
imshow(D05, [0,40])
colorbar

subplot(1,4,4)
imshow(D10, [0,40])
colorbar

colormap(turbo)


% [D01, IDX01] = bwdist(BW01);
% % D01(D01>=10)=NaN;
% % sdm00 = signedDistanceMap(1-BW00,InterpolationMethod="none");
% figure()
% imshow(D01,[0,10]);
% colormap("thermal-2")
% show(sdm00,BoundaryColor=[0 0 0],Colorbar="on");

% I00_8b_med_mean(I00_8b_med_mean<110)=NaN;