Biomedical Optics. Chapter 1

clear all; close all; clc;
%
% filename = 'blood.txt';
% delimiterIn = ' ';
% headerlinesIn = 2;
% A = importdata(filename,delimiterIn,headerlinesIn);

filename = 'blood.xlsx';
sheet = 1;
xlRange = 'C3:G378';

subsetA = xlsread(filename,sheet,xlRange);

filename2 = 'water_absorption.xlsx';
sheet = 1;
xlRange = 'E1:F232';

subsetB = xlsread(filename2,sheet,xlRange);

figure()
semilogy(subsetA(:,1),subsetA(:,2));
title('Molar Extinction Coefficient:')
xlabel('Wavelength(nm)');
ylabel('Molar extinction coefficient (cm-1/M');
grid on
hold on
semilogy(subsetA(:,1),subsetA(:,3));
legend('HbO2','Hb')
hold off


figure()
loglog(subsetB(:,1),subsetB(:,2));
title('Absorption spectrum')
xlabel('Wavelength(nm)');
ylabel('Absorption (cm-1');
grid on
hold on
loglog(subsetA(:,1),subsetA(:,4));
loglog(subsetA(:,1),subsetA(:,5));
legend('Water','HbO2','Hb')
hold off

% hold on
% semilogy(subsetA(:,1),subsetA(:,3));
%
% legend('HbO2','Hb')
% hold off

% figure()
% figure()
% semilogy(subsetA(:,1),subsetA(:,2));
% title('Molar Extinction Coefficient:')
% xlabel('Wavelength(nm)');
% ylabel('Molar extinction coefficient (cm-1/M');
% grid on
%
% hold on
% semilogy(subsetA(:,1),subsetA(:,3));
%
% legend('HbO2','Hb')
% hold off

% for k = [3, 5]
%    disp(A.colheaders{1, k})
%    disp(A.data(:, k))
%    disp(' ')
% end