% ***********************************************************************
% plasma_ri.m
% ***********************************************************************
%
% Data description:
%
% The refractive indices of plasma, assuming 4.6g/dL of albumin, 2.6g/dL 
% of globulin, and 0.38g/dL of fibrinogen for 250-1000nm with a 
% resolution of 5nm.
%
% Prepared by: Daniel Yim on December 2011
%
% ***********************************************************************
%
% Original data source:
%     Specific refraction increment of plasma proteins:
%     - G.S. Adair and M.E. Robinson. The specic refraction increments of 
%       serum-albumin and serum-globulin. Biochem. J., 24(4):993-1011, 1930.
%     - J. Voros. The density and refractive index of adsorbing protein
%       layers. Biophysical Journal, 87:553-561, 2004.
%
%     Index of refraction of water:
%     K.F. Palmer and D. Williams. Optical properties of water in the 
%     near infrared. J. Opt. Soc. Am., 64:1107-1110, August 1974.
% 
% ***********************************************************************
Mplasma = [ ...
  250.0, 1.379; ...
  255.0, 1.376; ...
  260.0, 1.374; ...
  265.0, 1.372; ...
  270.0, 1.37; ...
  275.0, 1.368; ...
  280.0, 1.366; ...
  285.0, 1.365; ...
  290.0, 1.363; ...
  295.0, 1.362; ...
  300.0, 1.36; ...
  305.0, 1.359; ...
  310.0, 1.358; ...
  315.0, 1.357; ...
  320.0, 1.356; ...
  325.0, 1.355; ...
  330.0, 1.354; ...
  335.0, 1.353; ...
  340.0, 1.352; ...
  345.0, 1.351; ...
  350.0, 1.351; ...
  355.0, 1.35; ...
  360.0, 1.349; ...
  365.0, 1.349; ...
  370.0, 1.348; ...
  375.0, 1.347; ...
  380.0, 1.347; ...
  385.0, 1.346; ...
  390.0, 1.346; ...
  395.0, 1.345; ...
  400.0, 1.346; ...
  405.0, 1.345; ...
  410.0, 1.344; ...
  415.0, 1.343; ...
  420.0, 1.342; ...
  425.0, 1.341; ...
  430.0, 1.34; ...
  435.0, 1.339; ...
  440.0, 1.338; ...
  445.0, 1.338; ...
  450.0, 1.337; ...
  455.0, 1.336; ...
  460.0, 1.336; ...
  465.0, 1.335; ...
  470.0, 1.335; ...
  475.0, 1.334; ...
  480.0, 1.334; ...
  485.0, 1.334; ...
  490.0, 1.334; ...
  495.0, 1.333; ...
  500.0, 1.333; ...
  505.0, 1.333; ...
  510.0, 1.333; ...
  515.0, 1.333; ...
  520.0, 1.333; ...
  525.0, 1.333; ...
  530.0, 1.333; ...
  535.0, 1.333; ...
  540.0, 1.333; ...
  545.0, 1.333; ...
  550.0, 1.333; ...
  555.0, 1.333; ...
  560.0, 1.333; ...
  565.0, 1.333; ...
  570.0, 1.333; ...
  575.0, 1.333; ...
  580.0, 1.333; ...
  585.0, 1.333; ...
  590.0, 1.333; ...
  595.0, 1.333; ...
  600.0, 1.333; ...
  605.0, 1.333; ...
  610.0, 1.333; ...
  615.0, 1.333; ...
  620.0, 1.333; ...
  625.0, 1.332; ...
  630.0, 1.332; ...
  635.0, 1.332; ...
  640.0, 1.332; ...
  645.0, 1.332; ...
  650.0, 1.332; ...
  655.0, 1.332; ...
  660.0, 1.332; ...
  665.0, 1.332; ...
  670.0, 1.332; ...
  675.0, 1.332; ...
  680.0, 1.332; ...
  685.0, 1.332; ...
  690.0, 1.332; ...
  695.0, 1.332; ...
  700.0, 1.332; ...
  705.0, 1.332; ...
  710.0, 1.332; ...
  715.0, 1.332; ...
  720.0, 1.332; ...
  725.0, 1.332; ...
  730.0, 1.332; ...
  735.0, 1.332; ...
  740.0, 1.332; ...
  745.0, 1.332; ...
  750.0, 1.332; ...
  755.0, 1.332; ...
  760.0, 1.332; ...
  765.0, 1.331; ...
  770.0, 1.331; ...
  775.0, 1.331; ...
  780.0, 1.331; ...
  785.0, 1.331; ...
  790.0, 1.331; ...
  795.0, 1.331; ...
  800.0, 1.331; ...
  805.0, 1.331; ...
  810.0, 1.331; ...
  815.0, 1.331; ...
  820.0, 1.33; ...
  825.0, 1.33; ...
  830.0, 1.33; ...
  835.0, 1.33; ...
  840.0, 1.33; ...
  845.0, 1.33; ...
  850.0, 1.33; ...
  855.0, 1.33; ...
  860.0, 1.33; ...
  865.0, 1.33; ...
  870.0, 1.33; ...
  875.0, 1.33; ...
  880.0, 1.33; ...
  885.0, 1.33; ...
  890.0, 1.329; ...
  895.0, 1.329; ...
  900.0, 1.329; ...
  905.0, 1.329; ...
  910.0, 1.329; ...
  915.0, 1.329; ...
  920.0, 1.329; ...
  925.0, 1.329; ...
  930.0, 1.329; ...
  935.0, 1.329; ...
  940.0, 1.329; ...
  945.0, 1.329; ...
  950.0, 1.328; ...
  955.0, 1.328; ...
  960.0, 1.328; ...
  965.0, 1.328; ...
  970.0, 1.328; ...
  975.0, 1.328; ...
  980.0, 1.328; ...
  985.0, 1.328; ...
  990.0, 1.328; ...
  995.0, 1.328; ...
  1000.0, 1.328];

plot(Mplasma(:,1), Mplasma(:,2), 'b', 'linewidth',1.5);
set(gca,'FontSize',14);
title('plasma index of refraction','FontSize',14);
ylabel('refractive index','fontsize',14);
xlabel('wavelength (nm)','fontsize',14);
axis([250 1000 1.25 1.45]);
set(gcf,'PaperUnits','inches','PaperPosition',[0 0 6 5]);
grid on;

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