load('Salina.mat');

for i = 0:16
    % initialize matrices with length of 100000 to avoid reallocation
    eval(sprintf('class%02d = zeros(100000,204);', i));
    eval(sprintf('classIter%02d = 1;', i));
end
for row = 1:size(salinas_corrected, 1)
    for col = 1:size(salinas_corrected, 2)
        i = salinas_gt(row, col);
        % load data into corresponding classes
        eval(sprintf('class%02d(classIter%02d, :) = reshape(salinas_corrected(row, col, :), [1, size(salinas_corrected, 3)]);', i, i));
        eval(sprintf('classIter%02d = classIter%02d + 1;', i, i));
    end
end
for i = 0:16
    % cut off excessive chunk
    eval(sprintf('class%02d = class%02d(1:classIter%02d-1,:);', i, i, i));
end

population = [];
for i = 16:-1:1
    eval(sprintf('population = [class%02d;population];', i));
end

% clear up workspace
clear -regexp ^classIter\d{2}$ ^salinas_
clear row col i;