% ***********************************************************************
% plasma_extinction.m
% ***********************************************************************
%
% Data description:
% The effective absorption coefficient of plasma in [um^-1] 
% (350 - 1000nm) with a resolution of 5nm:
%
% Prepared by: Daniel Yim on December 2011
%
% ***********************************************************************
%
% Original data source:
%
%     M. Meinke, G. Mller, M. Friebel, and J. Helfmann. Optical 
%     properties of platelets and blood plasma and their inuence on the 
%     optical behavior of whole blood in the visible to near infrared 
%     wavelength range. J. Biomed. Opt., 12(1):014024-1-014024-9, 2007.
%
% ***********************************************************************
M = [ ...
  350.0, 7.661400000000001E-5; ...
  355.0, 7.064400018550002E-5; ...
  360.0, 6.279662541455936E-5; ...
  365.0, 5.503143812183899E-5; ...
  370.0, 5.0263750261853285E-5; ...
  375.0, 4.686231233627968E-5; ...
  380.0, 4.491137493457812E-5; ...
  385.0, 4.3651249989068745E-5; ...
  390.0, 4.3544E-5; ...
  395.0, 4.55038749345781E-5; ...
  400.0, 4.881374982543124E-5; ...
  405.0, 5.24948747053531E-5; ...
  410.0, 5.629150019643137E-5; ...
  415.0, 6.0482125294646845E-5; ...
  420.0, 6.404812515287187E-5; ...
  425.0, 6.669931254364218E-5; ...
  430.0, 6.8972E-5; ...
  435.0, 7.089918747822033E-5; ...
  440.0, 7.12944999781375E-5; ...
  445.0, 7.1476E-5; ...
  450.0, 7.308974989085306E-5; ...
  455.0, 7.51396876637203E-5; ...
  460.0, 7.695962506542187E-5; ...
  465.0, 7.821975001093125E-5; ...
  470.0, 7.72065E-5; ...
  475.0, 7.41766250872844E-5; ...
  480.0, 7.099775010898129E-5; ...
  485.0, 6.874143763092657E-5; ...
  490.0, 6.561249978170615E-5; ...
  495.0, 5.9519187107220306E-5; ...
  500.0, 5.321162467255942E-5; ...
  505.0, 4.49937498254312E-5; ...
  510.0, 3.6825E-5; ...
  515.0, 3.167012510914688E-5; ...
  520.0, 2.7305875218128155E-5; ...
  525.0, 2.3073312827192198E-5; ...
  530.0, 1.923575021812829E-5; ...
  535.0, 1.5283124738146945E-5; ...
  540.0, 1.1413249803568704E-5; ...
  545.0, 7.485624901784413E-6; ...
  550.0, 4.3475E-6; ...
  555.0, 2.775062532710925E-6; ...
  560.0, 1.893625021696871E-6; ...
  565.0, 1.779E-6; ...
  570.0, 1.779E-6; ...
  575.0, 1.779E-6; ...
  580.0, 1.779E-6; ...
  585.0, 1.647E-6; ...
  590.0, 1.647E-6; ...
  595.0, 1.647E-6; ...
  600.0, 1.647E-6; ...
  605.0, 1.647E-6; ...
  610.0, 1.647E-6; ...
  615.0, 1.647E-6; ...
  620.0, 1.647E-6; ...
  625.0, 1.647E-6; ...
  630.0, 1.647E-6; ...
  635.0, 1.647E-6; ...
  640.0, 1.647E-6; ...
  645.0, 1.647E-6; ...
  650.0, 1.647E-6; ...
  655.0, 1.6387499672062495E-6; ...
  660.0, 1.5150000000000001E-6; ...
  665.0, 1.5150000000000001E-6; ...
  670.0, 1.5150000000000001E-6; ...
  675.0, 1.383E-6; ...
  680.0, 1.383E-6; ...
  685.0, 1.383E-6; ...
  690.0, 1.383E-6; ...
  695.0, 1.383E-6; ...
  700.0, 1.383E-6; ...
  705.0, 1.383E-6; ...
  710.0, 1.383E-6; ...
  715.0, 1.383E-6; ...
  720.0, 1.252E-6; ...
  725.0, 1.1942500327937412E-6; ...
  730.0, 1.1199999999999999E-6; ...
  735.0, 1.1199999999999999E-6; ...
  740.0, 1.1199999999999999E-6; ...
  745.0, 1.1199999999999999E-6; ...
  750.0, 1.1199999999999999E-6; ...
  755.0, 1.1199999999999999E-6; ...
  760.0, 1.1199999999999999E-6; ...
  765.0, 9.88E-7; ...
  770.0, 9.88E-7; ...
  775.0, 8.56E-7; ...
  780.0, 8.56E-7; ...
  785.0, 8.56E-7; ...
  790.0, 8.56E-7; ...
  795.0, 8.56E-7; ...
  800.0, 8.56E-7; ...
  805.0, 8.56E-7; ...
  810.0, 8.56E-7; ...
  815.0, 8.56E-7; ...
  820.0, 8.56E-7; ...
  825.0, 7.659374891515529E-7; ...
  830.0, 7.249999999999999E-7; ...
  835.0, 7.249999999999999E-7; ...
  840.0, 7.249999999999999E-7; ...
  845.0, 7.249999999999999E-7; ...
  850.0, 7.249999999999999E-7; ...
  855.0, 7.249999999999999E-7; ...
  860.0, 7.249999999999999E-7; ...
  865.0, 7.249999999999999E-7; ...
  870.0, 7.249999999999999E-7; ...
  875.0, 7.249999999999999E-7; ...
  880.0, 7.249999999999999E-7; ...
  885.0, 7.249999999999999E-7; ...
  890.0, 8.232500216968844E-7; ...
  895.0, 8.56E-7; ...
  900.0, 8.56E-7; ...
  905.0, 8.56E-7; ...
  910.0, 8.56E-7; ...
  915.0, 9.88E-7; ...
  920.0, 1.1199999999999999E-6; ...
  925.0, 1.1199999999999999E-6; ...
  930.0, 1.1199999999999999E-6; ...
  935.0, 1.252E-6; ...
  940.0, 1.383E-6; ...
  945.0, 1.5397500109312458E-6; ...
  950.0, 1.647E-6; ...
  955.0, 2.1244999781375055E-6; ...
  960.0, 2.469999956274987E-6; ...
  965.0, 2.964E-6; ...
  970.0, 3.8210000437250515E-6; ...
  975.0, 5.1131876637203145E-6; ...
  980.0, 7.608625131009303E-6; ...
  985.0, 1.1026375098215737E-5; ...
  990.0, 1.6930000000000002E-5; ...
  995.0, 2.2389812380004615E-5; ...
  1000.0, 2.676224980356869E-5];

plot(M(:,1), M(:,2), 'b', 'linewidth',1.5);
title('plasma absorption coefficient','FontSize',14);
set(gca,'FontSize',14);

% Molar extinction coefficient
ylabel('effective absorption coefficient (um^{-1})','fontsize',14);
xlabel('wavelength (nm)','fontsize',14);
axis([350 1000 0 1*10^-4]);
set(gcf,'PaperUnits','inches','PaperPosition',[0 0 6 5]);
grid on;

pause;
print -dpng plasma_extinction.png -r100
close;
