file_list = {
    'C:\tyo\20130624_sla\111.jpg',
    'C:\tyo\20130624_sla\112.jpg',
    'C:\tyo\20130624_sla\113.jpg',
    'C:\tyo\20130624_sla\121.jpg',
    'C:\tyo\20130624_sla\122.jpg',
    'C:\tyo\20130624_sla\123.jpg',
    'C:\tyo\20130624_sla\211.jpg',
    'C:\tyo\20130624_sla\212.jpg',
    'C:\tyo\20130624_sla\213.jpg',
    'C:\tyo\20130624_sla\221.jpg',
    'C:\tyo\20130624_sla\222.jpg',
    'C:\tyo\20130624_sla\223.jpg',
    'C:\tyo\20130624_sla\311.jpg',
    'C:\tyo\20130624_sla\312.jpg',
    'C:\tyo\20130624_sla\313.jpg',
    'C:\tyo\20130624_sla\321.jpg',
    'C:\tyo\20130624_sla\322.jpg',
    'C:\tyo\20130624_sla\323.jpg',
    'C:\tyo\20130624_sla\411.jpg',
    'C:\tyo\20130624_sla\412.jpg',
    'C:\tyo\20130624_sla\413.jpg',
    'C:\tyo\20130624_sla\421.jpg',
    'C:\tyo\20130624_sla\422.jpg',
    'C:\tyo\20130624_sla\423.jpg',
    'C:\tyo\20130624_sla\511.jpg',
    'C:\tyo\20130624_sla\512.jpg',
    'C:\tyo\20130624_sla\513.jpg',
    'C:\tyo\20130624_sla\521.jpg',
    'C:\tyo\20130624_sla\522.jpg',
    'C:\tyo\20130624_sla\523.jpg'
    };

file_list = {
    'C:\tyo\20130627_sla\1_1.jpg',
    'C:\tyo\20130627_sla\1_2.jpg',
    'C:\tyo\20130627_sla\1_3.jpg',
    'C:\tyo\20130627_sla\1_4.jpg',
    'C:\tyo\20130627_sla\1_5.jpg',
    'C:\tyo\20130627_sla\1_6.jpg',
    'C:\tyo\20130627_sla\1_7.jpg',
    'C:\tyo\20130627_sla\1_8.jpg',
    'C:\tyo\20130627_sla\1_9.jpg',
    'C:\tyo\20130627_sla\1_10.jpg',
    'C:\tyo\20130627_sla\2_1.jpg',
    'C:\tyo\20130627_sla\2_2.jpg',
    'C:\tyo\20130627_sla\2_3.jpg',
    'C:\tyo\20130627_sla\2_4.jpg',
    'C:\tyo\20130627_sla\2_5.jpg',
    'C:\tyo\20130627_sla\2_6.jpg',
    'C:\tyo\20130627_sla\2_7.jpg',
    'C:\tyo\20130627_sla\2_8.jpg',
    'C:\tyo\20130627_sla\2_9.jpg',
    'C:\tyo\20130627_sla\2_10.jpg'
    };

% Tables to store graylevel threshold and leaf area
area = zeros(20,1);
threshold = zeros(20,1);

% Image width and height in millimeters
width = 215.9;
height = 297.2;

for i = 1:20
    A = imread(file_list{i}, 'jpg');
    imshow(A(:,:,1));
    level = graythresh(A(:,:,1));
    BW = im2bw(A(:,:,1), level);
    imshow(BW);
    area(i) = (1 - sum(sum(BW))/(length(BW(1,:))*length(BW(:,1)))) * width * height;
    threshold(i) = level;
    pause(2);
end





