data1: collect manual data information
clear close all load filelist
% from lines.xlsx
zinches = [ 31 46.5 65 86 107 128 149 170 193 216 235 255]';
% cycle through images sz = size(files); N = sz(1); for m=1:N rgb = imread(files(m,:)); sz = size(rgb); small = imcrop(rgb,[sz(1)/4 sz(2)/4 512 512]); g = rgb2gray(im2double(small)); imshow(g); cdata = ginput; data(:,:,m) = cdata(1:2,:); fprintf('image %d threshold %g width %g\n',m,graythresh(g),data(2,1,m)-data(1,1,m)); end save imgdata files data zinches