call_phansalkar.m

function call_phansalkar()

choroidIm = imread('norm.tif');
newIm = im2double(choroidIm);

%% *uncomment below for image resizing option*

newIm = imresize(newIm, [304 304]);
BW = phansalkar(newIm, [15 15]); 

%% *automated window size adjustment based on input image size*

% windowSize = round(size(newIm ,1)/(304/15));
% BW = phansalkar(newIm, [windowSize windowSize]);
BW2 = imcomplement(BW);

figure; imshow([newIm BW2],[])
imwrite(BW2, 'norm_BW_MATLAB.tif')