%$Release 4.0
% Copyright (c) Copyright 2004 - 2006 - The General Hospital Corporation and
% President and Fellows of Harvard University.
%
% All rights reserved.
%
% Redistribution and use in source and binary forms, with or without
% modification, are permitted provided that the following conditions are met:
%
% * Redistributions of source code must retain the above copyright
% notice, this list of conditions and the following disclaimer.
% * Redistributions in binary form must reproduce the above copyright
% notice, this list of conditions and the following disclaimer in the
% documentation and/or other materials provided with the distribution.
% * Neither the name of The General Hospital Corporation and Harvard
% University nor the names of its contributors may be used to endorse or
% promote products derived from this software without specific prior written
% permission.
%
% The Software has been designed for research purposes only and has not been
% reviewed or approved by the Food and Drug Administration or by any other
% agency. YOU ACKNOWLEDGE AND AGREE THAT CLINICAL APPLICATIONS ARE NEITHER
% RECOMMENDED NOR ADVISED.
%
% THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
% AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
% IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
% ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
% LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
% CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
% SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
% INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
% CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
% ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
% POSSIBILITY OF SUCH DAMAGE.
function [exs,varargout] = GetExtinctions( lambda,WhichSpectrum )
%
% GetExtinctions( lambda )
%
% Returns the specific absorption coefficients for
% [HbO Hb H2O lipid aa3]
% for the specified wavelengths. Note that the specific
% absorption coefficient (defined base e) is equal to the
% specific extinction coefficient (defined base 10) times 2.303.
%
% These values for the molar extinction coefficient e
% in [cm-1/(moles/liter)] were compiled by Scott Prahl
% (prahl@ece.ogi.edu) using data from
% W. B. Gratzer, Med. Res. Council Labs, Holly Hill, London
% N. Kollias, Wellman Laboratories, Harvard Medical School, Boston
% To convert this data to absorbance A, multiply by the
% molar concentration and the pathlength. For example, if x is the
% number of grams per liter and a 1 cm cuvette is being used,
% then the absorbance is given by
%
% (e) [(1/cm)/(moles/liter)] (x) [g/liter] (1) [cm]
% A = ---------------------------------------------------
% 66,500 [g/mole]
%
% using 66,500 as the gram molecular weight of hemoglobin.
% To convert this data to absorption coefficient in (cm-1), multiply
% by the molar concentration and 2.303,
%
% 琦 = (2.303) e (x g/liter)/(66,500 g Hb/mole)
% where x is the number of grams per liter. A typical value of x
% for whole blood is x=150 g Hb/liter.
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% $Author: dboasg $
%
% $Date: 2013-07-13 14:19:04 $
%
% $Revision: 1.2 $
%
% $Log: not supported by cvs2svn $
% Revision 1.1 2009/01/02 16:55:16 dboasg
% initial
%
% Revision 1.2 2006/12/11 18:40:58 dboasg
% updated copy right info
%
% Revision 1.1 2006/09/21 17:17:29 thuppert
% initial version
%
% Revision 1.2 2005/07/29 17:12:08 thuppert
% Added MGH copyright header to all *.m files!
%
% Revision 1.1.1.1 2005/07/01 14:25:11 dboas
% Initial import of homer version 4.0
%
% Revision 1.2 2003/08/29 18:49:59 bc
% Updated chimeric version (dboas + thuppert)
%
% Revision 1.4 2002/02/18 20:20:02 dboas
% Corrected the units of lipid absorption to be per cm rather than per mm.
% Also, now interpolate for wavelengths between those specified.
%
% Revision 1.3 2001/04/27 20:46:53 dboas
% Maria Angela Updated the extinction coefficients for Hb and HbO
%
% Revision 1.2 2000/09/06 12:10:12 dboas
% Added extinction coefficients for lipid and aa3.
%
% Revision 1.1.1.1 2000/05/25 13:14:47 dboas
% initial
%
% Revision 1.2 2000/01/10 00:14:14 dboas
% Storing the source and detector lists for use by other functions
%
% Revision 1.1 1999/11/18 14:21:00 tgaudett
% Initial Routines for Chromophores
% HB, HbO, H2O
%
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
num_lambda = length(lambda);
if ~exist('WhichSpectrum')
WhichSpectrum=1;
end
switch(WhichSpectrum)
%**************************************************************************
%**************************************************************************
%**************************************************************************
%**************************************************************************
case 1
Citation=sprintf('%s\n%s','W. B. Gratzer, Med. Res. Council Labs, Holly Hill,London',...
'N. Kollias, Wellman Laboratories, Harvard Medical School, Boston');
%These values for the molar extinction coefficient e in [cm-1/(moles/liter)] were compiled by Scott Prahl (prahl@ece.ogi.edu) using data from
%
%W. B. Gratzer, Med. Res. Council Labs, Holly Hill, London
%N. Kollias, Wellman Laboratories, Harvard Medical School, Boston
%To convert this data to absorbance A, multiply by the molar concentration and the pathlength. For example, if x is the number of grams per liter and a 1 cm cuvette is being used, then the absorbance is given by
%
% (e) [(1/cm)/(moles/liter)] (x) [g/liter] (1) [cm]
% A = ---------------------------------------------------
% 66,500 [g/mole]
%
%using 66,500 as the gram molecular weight of hemoglobin.
%To convert this data to absorption coefficient in (cm-1), multiply by the molar concentration and 2.303,
%
%琦 = (2.303) e (x g/liter)/(66,500 g Hb/mole)
%where x is the number of grams per liter. A typical value of x for whole blood is x=150 g Hb/liter.
vLambdaHbOHb = [
250 106112 112736;
252 105552 112736;
254 107660 112736;
256 109788 113824;
258 112944 115040;
260 116376 116296;
262 120188 117564;
264 124412 118876;
266 128696 120208;
268 133064 121544;
270 136068 122880;
272 137232 123096;
274 138408 121952;
276 137424 120808;
278 135820 119840;
280 131936 118872;
282 127720 117628;
284 122280 114820;
286 116508 112008;
288 108484 107140;
290 104752 98364;
292 98936 91636;
294 88136 85820;
296 79316 77100;
298 70884 69444;
300 65972 64440;
302 63208 61300;
304 61952 58828;
306 62352 56908;
308 62856 57620;
310 63352 59156;
312 65972 62248;
314 69016 65344;
316 72404 68312;
318 75536 71208;
320 78752 74508;
322 82256 78284;
324 85972 82060;
326 89796 85592;
328 93768 88516;
330 97512 90856;
332 100964 93192;
334 103504 95532;
336 104968 99792;
338 106452 104476;
340 107884 108472;
342 109060 110996;
344 110092 11352