function [Output] = LU_pix_by_pix(Image,Library)

% Linear Unmixing
% Silas J. Leavesley, PhD
% Сòòò½APP
% Last Updated: 2/7/2017

% To visualize output data:
% figure(1)
% for i = 1:length(Output.data(1,1,:))
%   subplot(1,length(Output.data(1,1,:)),i)
%   imagesc(Output.data(:,:,i))
% end
% colormap(gray) % or any other colormap, as desired

tic

Image=double(Image);

[L M N]=size(Image);
[O P]=size(Library);

h=waitbar(0,'Please wait...');
for i=1:L
    waitbar(i/L,h,sprintf('Unmixing row %5.0f out of %5.0f',i,L))
    for j=1:M
        [x(:,i,j),resnorm,residual]=lsqnonneg(Library,permute(double(Image(i,j,:)),[3 1 2]));
        Output.Error.RMS(i,j)=(sum(residual.^2)/O)^0.5;
        RMS_signal(i,j)=(sum((Image(i,j,:).^2))/O)^0.5;
    end
end

Output.data=permute(x,[2 3 1]);
Output.Error.Percent_RMS=Output.Error.RMS./RMS_signal;
delete(h)
toc