-
Notifications
You must be signed in to change notification settings - Fork 6
/
GetExtinctions.m
1600 lines (1549 loc) · 52.2 KB
/
GetExtinctions.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
%$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 extinction coefficients for
% [HbO Hb H2O lipid aa3]
% for the specified wavelengths.
%
% These values for the molar extinction coefficient e
% in [cm-1/(moles/liter)] were compiled by Scott Prahl
% ([email protected]) 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,
%
% µa = (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: 2009/01/02 16:55:16 $
%
% $Revision: 1.1 $
%
% $Log: GetExtinctions.m,v $
% 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 ([email protected]) 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,
%
%µa = (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 113524;
346 109032 116052;
348 107984 118752;
350 106576 122092;
352 105040 125436;
354 103696 128776;
356 101568 132120;
358 97828 133632;
360 94744 134940;
362 92248 136044;
364 89836 136972;
366 88484 137900;
368 87512 138856;
370 88176 139968;
372 91592 141084;
374 95140 142196;
376 98936 143312;
378 103432 144424;
380 109564 145232;
382 116968 145232;
384 125420 148668;
386 135132 153908;
388 148100 159544;
390 167748 167780;
392 189740 180004;
394 212060 191540;
396 231612 202124;
398 248404 212712;
400 266232 223296;
402 284224 236188;
404 308716 253368;
406 354208 270548;
408 422320 287356;
410 466840 303956;
412 500200 321344;
414 524280 342596;
416 521880 363848;
418 515520 385680;
420 480360 407560;
422 431880 429880;
424 376236 461200;
426 326032 481840;
428 283112 500840;
430 246072 528600;
432 214120 552160;
434 165332 552160;
436 132820 547040;
438 119140 501560;
440 102580 413280;
442 92780 363240;
444 81444 282724;
446 76324 237224;
448 67044 173320;
450 62816 103292;
452 58864 62640;
454 53552 36170;
456 49496 30698.8;
458 47496 25886.4;
460 44480 23388.8;
462 41320 20891.2;
464 39807.2 19260.8;
466 37073.2 18142.4;
468 34870.8 17025.6;
470 33209.2 16156.4;
472 31620 15310;
474 30113.6 15048.4;
476 28850.8 14792.8;
478 27718 14657.2;
480 26629.2 14550;
482 25701.6 14881.2;
484 25180.4 15212.4;
486 24669.6 15543.6;
488 24174.8 15898;
490 23684.4 16684;
492 23086.8 17469.6;
494 22457.6 18255.6;
496 21850.4 19041.2;
498 21260 19891.2;
500 20932.8 20862;
502 20596.4 21832.8;
504 20418 22803.6;
506 19946 23774.4;
508 19996 24745.2;
510 20035.2 25773.6;
512 20150.4 26936.8;
514 20429.2 28100;
516 21001.6 29263.2;
518 22509.6 30426.4;
520 24202.4 31589.6;
522 26450.4 32851.2;
524 29269.2 34397.6;
526 32496.4 35944;
528 35990 37490;
530 39956.8 39036.4;
532 43876 40584;
534 46924 42088;
536 49752 43592;
538 51712 45092;
540 53236 46592;
542 53292 48148;
544 52096 49708;
546 49868 51268;
548 46660 52496;
550 43016 53412;
552 39675.2 54080;
554 36815.2 54520;
556 34476.8 54540;
558 33456 54164;
560 32613.2 53788;
562 32620 52276;
564 33915.6 50572;
566 36495.2 48828;
568 40172 46948;
570 44496 45072;
572 49172 43340;
574 53308 41716;
576 55540 40092;
578 54728 38467.6;
580 50104 37020;
582 43304 35676.4;
584 34639.6 34332.8;
586 26600.4 32851.6;
588 19763.2 31075.2;
590 14400.8 28324.4;
592 10468.4 25470;
594 7678.8 22574.8;
596 5683.6 19800;
598 4504.4 17058.4;
600 3200 14677.2;
602 2664 13622.4;
604 2128 12567.6;
606 1789.2 11513.2;
608 1647.6 10477.6;
610 1506 9443.6;
612 1364.4 8591.2;
614 1222.8 7762;
616 1110 7344.8;
618 1026 6927.2;
620 942 6509.6;
622 858 6193.2;
624 774 5906.8;
626 707.6 5620;
628 658.8 5366.8;
630 610 5148.8;
632 561.2 4930.8;
634 512.4 4730.8;
636 478.8 4602.4;
638 460.4 4473.6;
640 442 4345.2;
642 423.6 4216.8;
644 405.2 4088.4;
646 390.4 3965.08;
648 379.2 3857.6;
650.0 506.0 3743.0;
652.0 488.0 3677.0;
654.0 474.0 3612.0;
656.0 464.0 3548.0;
658.0 454.3 3491.3;
660.0 445.0 3442.0;
662.0 438.3 3364.7;
664.0 433.8 3292.8;
666.0 431.3 3226.3;
668.0 429.0 3133.0;
670.0 427.0 3013.0;
672.0 426.5 2946.0;
674.0 426.0 2879.0;
676.0 424.0 2821.7;
678.0 423.0 2732.3;
680.0 423.0 2610.8;
682.0 422.0 2497.3;
684.0 420.0 2392.0;
686.0 418.0 2292.7;
688.0 416.5 2209.3;
690.0 415.5 2141.8;
692.0 415.0 2068.7;
694.0 415.0 1990.0;
696.0 415.5 1938.5;
698.0 416.0 1887.0;
700.0 419.3 1827.7;
702.0 422.5 1778.5;
704.0 425.5 1739.5;
706.0 429.7 1695.7;
708.0 435.0 1647.0;
710.0 441.0 1601.7;
712.0 446.5 1562.5;
714.0 451.5 1529.5;
716.0 458.0 1492.0;
718.0 466.0 1450.0;
720.0 472.7 1411.3;
722.0 479.5 1380.0;
724.0 486.5 1356.0;
726.0 494.3 1331.7;
728.0 503.0 1307.0;
730.0 510.0 1296.5;
732.0 517.0 1286.0;
734.0 521.0 1286.0;
736.0 530.7 1293.0;
738.0 546.0 1307.0;
740.0 553.5 1328.0;
742.0 561.0 1349.0;
744.0 571.0 1384.3;
746.0 581.3 1431.3;
748.0 592.0 1490.0;
750.0 600.0 1532.0;
752.0 608.0 1574.0;
754.0 618.7 1620.7;
756.0 629.7 1655.3;
758.0 641.0 1678.0;
760.0 645.5 1669.0;
762.0 650.0 1660.0;
764.0 666.7 1613.3;
766.0 681.0 1555.0;
768.0 693.0 1485.0;
770.0 701.5 1425.0;
772.0 710.0 1365.0;
774.0 722.0 1288.3;
776.0 733.7 1216.3;
778.0 745.0 1149.0;
780.0 754.0 1107.5;
782.0 763.0 1066.0;
784.0 775.0 1021.3;
786.0 787.0 972.0;
788.0 799.0 918.0;
790.0 808.0 913.0;
792.0 817.0 908.0;
794.0 829.0 887.3;
796.0 840.7 868.7;
798.0 852.0 852.0;
800.0 863.3 838.7;
802.0 873.3 828.0;
804.0 881.8 820.0;
806.0 891.7 812.0;
808.0 903.0 804.0;
810.0 914.3 798.7;
812.0 924.7 793.7;
814.0 934.0 789.0;
816.0 943.0 787.0;
818.0 952.0 785.0;
820.0 962.7 783.0;
822.0 973.0 781.0;
824.0 983.0 779.0;
826.0 990.5 778.5;
828.0 998.0 778.0;
830.0 1008.0 778.0;
832.0 1018.0 777.7;
834.0 1028.0 777.0;
836.0 1038.0 777.0;
838.0 1047.7 777.0;
840.0 1057.0 777.0;
842.0 1063.5 777.5;
844.0 1070.0 778.0;
846.0 1079.3 779.3;
848.0 1088.3 780.3;
850.0 1097.0 781.0;
852.0 1105.7 783.0;
854.0 1113.0 785.0;
856.0 1119.0 787.0;
858.0 1126.0 789.3;
860.0 1134.0 792.0;
862.0 1142.0 795.3;
864.0 1149.7 799.0;
866.0 1157.0 803.0;
868.0 1163.7 807.7;
870.0 1170.3 812.3;
872.0 1177.0 817.0;
874.0 1182.0 820.5;
876.0 1187.0 824.0;
878.0 1193.0 830.0;
880.0 1198.7 835.7;
882.0 1204.0 841.0;
884.0 1209.3 847.0;
886.0 1214.3 852.7;
888.0 1219.0 858.0;
890.0 1223.7 863.3;
892.0 1227.5 867.8;
894.0 1230.5 871.3;
896.0 1234.0 875.3;
898.0 1238.0 880.0;
900.0 1241.3 883.3;
902 1202 765.04;
904 1206 767.44;
906 1209.2 769.8;
908 1211.6 772.16;
910 1214 774.56;
912 1216.4 776.92;
914 1218.8 778.4;
916 1220.8 778.04;
918 1222.4 777.72;
920 1224 777.36;
922 1225.6 777.04;
924 1227.2 776.64;
926 1226.8 772.36;
928 1224.4 768.08;
930 1222 763.84;
932 1219.6 752.28;
934 1217.2 737.56;
936 1215.6 722.88;
938 1214.8 708.16;
940 1214 693.44;
942 1213.2 678.72;
944 1212.4 660.52;
946 1210.4 641.08;
948 1207.2 621.64;
950 1204 602.24;
952 1200.8 583.4;
954 1197.6 568.92;
956 1194 554.48;
958 1190 540.04;
960 1186 525.56;
962 1182 511.12;
964 1178 495.36;
966 1173.2 473.32;
968 1167.6 451.32;
970 1162 429.32;
972 1156.4 415.28;
974 1150.8 402.28;
976 1144 389.288;
978 1136 374.944;
980 1128 359.656;
982 1120 344.372;
984 1112 329.084;
986 1102.4 313.796;
988 1091.2 298.508;
990 1080 283.22;
992 1068.8 267.932;
994 1057.6 252.648;
996 1046.4 237.36;
998 1035.2 222.072;
1000 1024 206.784 ];
vLambdaHbOHb(:,2) = vLambdaHbOHb(:,2) * 2.303;
vLambdaHbOHb(:,3) = vLambdaHbOHb(:,3) * 2.303;
%**************************************************************************
%**************************************************************************
%**************************************************************************
%**************************************************************************
case 2
Citation=sprintf('%s\n%s','J.M. Schmitt, "Optical Measurement of Blood Oxygenation by Implantable Telemetry," Technical Report G558-15, Stanford."',...
'M.K. Moaveni, "A Multiple Scattering Field Theory Applied to Whole Blood," Ph.D. dissertation, Dept. of Electrical Engineering, University of Washington, 1970');
% These values for the molar extinction coefficient e in [cm-1/(moles/liter)] were compiled by Scott Prahl ([email protected]) using data from
%
% J.M. Schmitt, "Optical Measurement of Blood Oxygenation by Implantable Telemetry," Technical Report G558-15, Stanford."
% M.K. Moaveni, "A Multiple Scattering Field Theory Applied to Whole Blood," Ph.D. dissertation, Dept. of Electrical Engineering, University of Washington, 1970.
% 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,
%
% µa = (2.303) e (x g/liter)/(66,500 g Hb/mole)
vLambdaHbOHb = [ 630 680 4280;
640 440 3640;
650 380 3420;
660 320 3200;
670 320 3080;
680 320 2960;
690 280 2560;
700 320 2160;
710 340 1840;
720 360 1520;
730 400 1500;
740 440 1520;
750 520 1620;
760 600 1720;
770 660 1420;
780 720 1120;
790 760 1020;
800 800 920;
810 860 880;
820 920 840;
830 980 840;
840 1040 840;
850 1060 800;
860 1080 840;
870 1120 840;
880 1160 840;
890 1180 860;
900 1200 880;
910 1220 920;
920 1240 880;
930 1240 800;
940 1200 800;
950 1200 720 ];
vLambdaHbOHb(:,2) = vLambdaHbOHb(:,2) * 2.303;
vLambdaHbOHb(:,3) = vLambdaHbOHb(:,3) * 2.303;
%**************************************************************************
%**************************************************************************
%**************************************************************************
%******************************************************************
case 3
Citation=sprintf('%s','S. Takatani and M. D. Graham, "Theoretical analysis of diffuse reflectance from a two-layer tissue model," IEEE Trans. Biomed. Eng., BME-26, 656--664, (1987). ');
% These values for the molar extinction coefficient e in [cm-1/(moles/liter)] were compiled by Scott Prahl ([email protected]) using data from
%
% S. Takatani and M. D. Graham, "Theoretical analysis of diffuse reflectance from a two-layer tissue model," IEEE Trans. Biomed. Eng., BME-26, 656--664, (1987).
% 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,
%
% µa = (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 = [450 68000 58000;
460 45040 20600;
480 27360 13360;
500 20200 16360;
507 19240 19240;
510 19040 20000;
520 23520 25080;
522 25680 25680;
540 57080 41120;
542 57480 44000;
549 49840 49840;
555 36000 52160;
560 33880 50160;
569 45080 45080;
577 61480 36800;
579 54920 35440;
586 28920 28920;
600 3200 14600;
605 1860 9496;
615 1152 5776;
625 732 4400;
635 488 3796;
645 396 3436;
655 340 3244;
665 292 3156;
675 288 3028;
685 272 2796;
695 280 2424;
705 300 1988;
715 328 1628;
725 368 1464;
735 412 1464;
745 480 1616;
755 556 1756;
765 616 1640;
775 684 1340;
785 736 1040;
795 776 964;
805 880 896;
815 880 880;
825 952 832;
835 996 820;
845 1048 820;
855 1068 820;
865 1116 820;
875 1140 848;
885 1168 832;
895 1188 884;
905 1208 896;
915 1220 924;
925 1228 860;
935 1216 848;
945 1212 756;
955 1196 704;
965 1176 616;
975 1148 552;
985 1108 424;
995 1052 372];
vLambdaHbOHb(:,2) = vLambdaHbOHb(:,2) * 2.303;
vLambdaHbOHb(:,3) = vLambdaHbOHb(:,3) * 2.303;
end
%
% ABSORPTION SPECTRUMOF H20
% FROM G. M. Hale and M. R. Querry, "Optical constants of water in the 200nm to
% 200µm wavelength region," Appl. Opt., 12, 555--563, (1973).
%
% ON THE WEB AT
% http://omlc.ogi.edu/spectra/water/abs/index.html
%
vLambdaH2O = [
200.00 0.069000;
225.00 0.027400;
250.00 0.016800;
275.00 0.010700;
300.00 0.0067000;
325.00 0.0041800;
350.00 0.0023300;
375.00 0.0011700;
400.00 0.00058000;
425.00 0.00038000;
450.00 0.00028000;
475.00 0.00024700;
500.00 0.00025000;
525.00 0.00032000;
550.00 0.00045000;
575.00 0.00079000;
600.00 0.0023000;
625.00 0.0028000;
650.00 0.0032000;
675.00 0.0041500;
700.00 0.0060000;
725.00 0.015900;
750.00 0.026000;
775.00 0.024000;
800.00 0.020000;
810.00 0.019858;
820.00 0.023907;
825.00 0.028000;
830.00 0.029069;
840.00 0.034707;
850.00 0.043000;
860.00 0.046759;
870.00 0.051999;
875.00 0.056000;
880.00 0.055978;
890.00 0.060432;
900.00 0.068000;
910.00 0.072913;
920.00 0.10927;
925.00 0.14400;
930.00 0.17296;
940.00 0.26737;
950.00 0.39000;
960.00 0.42000;
970.00 0.45000;
975.00 0.45000;
980.00 0.43000;
990.00 0.41000;
1000.0 0.36000];
%
% Extinction coefficient for lipid.
% I got this from Brian Pogue who got this from Matcher and Cope (DAB)
% In units of per mm and convert to per cm.
%
vLambdaLipid = [
650 0.000080;
652 0.000080;
654 0.000080;
656 0.000080;
658 0.000080;
660 0.000080;
662 0.000080;
664 0.000080;
666 0.000080;
668 0.000080;
670 0.000080;
672 0.000080;
674 0.000080;
676 0.000080;
678 0.000080;
680 0.000080;
682 0.000080;
684 0.000080;
686 0.000080;
688 0.000080;
690 0.000080;
692 0.000080;
694 0.000080;
696 0.000080;
698 0.000080;
700 0.000080;
702 0.000080;
704 0.000080;
706 0.000080;
708 0.000080;
710 0.000080;
712 0.000080;
714 0.000080;
716 0.000080;
718 0.000080;
720 0.000096;
722 0.000101;
724 0.000096;
726 0.000100;
728 0.000090;
730 0.000089;
732 0.000093;
734 0.000105;
736 0.000123;
738 0.000148;
740 0.000179;
742 0.000214;
744 0.000254;
746 0.000296;
748 0.000341;
750 0.000385;
752 0.000426;
754 0.000462;
756 0.000491;
758 0.000510;
760 0.000515;
762 0.000498;
764 0.000458;
766 0.000399;
768 0.000333;
770 0.000267;
772 0.000206;
774 0.000155;
776 0.000106;
778 0.000063;
780 0.000033;
782 0.000021;
784 0.000023;
786 0.000029;
788 0.000027;
790 0.000017;
792 0.000006;
794 0.000000;
796 0.000002;
798 0.000009;
800 0.000021;
802 0.000037;
804 0.000055;
806 0.000073;
808 0.000091;
810 0.000109;
812 0.000128;
814 0.000146;
816 0.000163;
818 0.000178;
820 0.000193;
822 0.000205;
824 0.000214;
826 0.000219;
828 0.000221;
830 0.000219;
832 0.000213;
834 0.000205;
836 0.000193;
838 0.000180;
840 0.000167;
842 0.000155;
844 0.000145;
846 0.000138;
848 0.000135;
852 0.000130;
854 0.000130;
856 0.000137;
858 0.000153;
860 0.000179;
862 0.000216;
864 0.000265;
866 0.000329;
868 0.000408;
870 0.000505;
872 0.000618;
874 0.000758;
876 0.000919;
878 0.001103;
880 0.001314;
882 0.001552;
884 0.001809;
886 0.002082;
888 0.002365;
890 0.002653;
892 0.002939;
894 0.003224;
896 0.003521;
898 0.003833;
900 0.004168;
902 0.004545;
904 0.004976;
906 0.005451;
908 0.005969;
910 0.006522;
912 0.007106;
914 0.007686;
916 0.008235;
918 0.008744;
920 0.009179;
922 0.009484;
924 0.009644;
926 0.009602;
928 0.009363;
930 0.008895;
932 0.008275;
934 0.007497;
936 0.006648;
938 0.005792;
940 0.004947;
942 0.004166;
944 0.003467;
946 0.002859;
948 0.002337;
950 0.001898;
952 0.001544;
954 0.001258;
956 0.001022;
958 0.000830;
960 0.000678;
962 0.000553;
964 0.000451;
966 0.000365;
968 0.000298;
970 0.000244;
972 0.000204;
974 0.000176;
976 0.000163;
978 0.000165;
980 0.000185;
982 0.000221;
984 0.000278;
986 0.000351;
988 0.000448;
990 0.000564;
992 0.000698;
994 0.000853;
996 0.001027;
998 0.001215;
1000 0.001420;
1002 0.001643;
1004 0.001882;
1006 0.002127;
1008 0.002371;
1010 0.002619;
1012 0.002858;
1014 0.003077;
1016 0.003279;
1018 0.003463;
1020 0.003635;
1022 0.003793;
1024 0.003941;
1026 0.004084;
1028 0.004220;
1030 0.004338;
1032 0.004438;
1034 0.004505;
1036 0.004537;
1038 0.004524;
1040 0.004468;
1042 0.004365;
1044 0.004229;
1046 0.004065;
1048 0.003880;
1050 0.003677;
1052 0.003469;
1054 0.003259;
1056 0.003051;
1058 0.002843] * 10;
vLambdaAA3 = [
6.5000000e+002 1.1361272e+000;
6.5050000e+002 1.1318779e+000;
6.5100000e+002 1.1276285e+000;
6.5150000e+002 1.1233792e+000;
6.5200000e+002 1.1191298e+000;
6.5250000e+002 1.1135125e+000;
6.5300000e+002 1.1078952e+000;
6.5350000e+002 1.1022779e+000;
6.5400000e+002 1.0966605e+000;
6.5450000e+002 1.0903251e+000;
6.5500000e+002 1.0839896e+000;
6.5550000e+002 1.0776542e+000;
6.5600000e+002 1.0713187e+000;
6.5650000e+002 1.0643820e+000;
6.5700000e+002 1.0574453e+000;
6.5750000e+002 1.0505086e+000;
6.5800000e+002 1.0435719e+000;
6.5850000e+002 1.0362062e+000;
6.5900000e+002 1.0288404e+000;
6.5950000e+002 1.0214747e+000;
6.6000000e+002 1.0141090e+000;
6.6050000e+002 1.0062588e+000;
6.6100000e+002 9.9840857e-001;
6.6150000e+002 9.9055837e-001;
6.6200000e+002 9.8270816e-001;
6.6250000e+002 9.7439221e-001;
6.6300000e+002 9.6607626e-001;
6.6350000e+002 9.5776031e-001;
6.6400000e+002 9.4944435e-001;
6.6450000e+002 9.4080158e-001;
6.6500000e+002 9.3215880e-001;
6.6550000e+002 9.2351603e-001;