[bug] fix the potential typo in Custo et al for CSF mua, fix #232

we noticed that the absorption coefficient in Custo et al. 2006, mua=0.004/mm is 10x larger than the cited reference, Strangman 2015, which shows a value of 0.0004/mm. we are now moving to the upstream value to avoid further propagating a possible typo
fangq · Sep 13, 2024 · 6fbf3db · 6fbf3db
1 parent c9456a7
commit 6fbf3db
Show file tree

Hide file tree

Showing 7 changed files with 7 additions and 7 deletions.
diff --git a/example/usc19-5/USC_19-5_atlas.json b/example/usc19-5/USC_19-5_atlas.json
@@ -46,7 +46,7 @@
 				"n": 1.37
 			},
 			{
-				"mua": 0.004,
+				"mua": 0.0004,
 				"mus": 0.009,
 				"g": 0.89,
 				"n": 1.37

diff --git a/mcxlab/examples/demo_4layer_head.m b/mcxlab/examples/demo_4layer_head.m
@@ -35,7 +35,7 @@
 
 cfg.prop = [0 0 1 1            % medium 0: the environment
             0.019 7.8   0.89 1.37     % medium 1: skin & skull
-            0.004 0.009 0.89 1.37     % medium 2: CSF
+            0.0004 0.009 0.89 1.37    % medium 2: CSF
             0.02  9.0   0.89 1.37     % medium 3: gray matter
             0.08 40.9   0.84 1.37];   % medium 4: white matter
 

diff --git a/mcxlab/examples/demo_colin27_atlas.m b/mcxlab/examples/demo_colin27_atlas.m
@@ -48,7 +48,7 @@
 cfg.prop = [0         0    1.0000    1.0000 % background/air
             0.0190    7.8182    0.8900    1.3700 % scalp
             0.0190    7.8182    0.8900    1.3700 % skull
-            0.0040    0.0090    0.8900    1.3700 % csf
+            0.0004    0.0090    0.8900    1.3700 % csf
             0.0200    9.0000    0.8900    1.3700 % gray matters
             0.0800   40.9000    0.8400    1.3700 % white matters
             0         0    1.0000    1.0000]; % air pockets

diff --git a/mcxlab/examples/demo_fullhead_atlas.m b/mcxlab/examples/demo_fullhead_atlas.m
@@ -25,7 +25,7 @@
 
 % tissue labels:0-ambient air,1-scalp,2-skull,3-csf,4-gray matter,5-white matter,6-air cavities
 cfg.vol = USC_atlas;
-cfg.prop = [0, 0, 1, 1; 0.019 7.8 0.89 1.37; 0.02 9.0 0.89 1.37; 0.004 0.009 0.89 1.37; 0.019 7.8 0.89 1.37; 0.08 40.9 0.84 1.37; 0, 0, 1, 1];
+cfg.prop = [0, 0, 1, 1; 0.019 7.8 0.89 1.37; 0.02 9.0 0.89 1.37; 0.0004 0.009 0.89 1.37; 0.019 7.8 0.89 1.37; 0.08 40.9 0.84 1.37; 0, 0, 1, 1];
 
 % light source
 cfg.srcnum = 1;

diff --git a/mcxlab/examples/demo_svmc_brain19_5.m b/mcxlab/examples/demo_svmc_brain19_5.m
@@ -37,7 +37,7 @@
 cfg.prop = [0.0   0.0   1.0  1.0    % background medium
             0.019 7.8   0.89 1.37   % scalp
             0.019 7.8   0.89 1.37   % skull
-            0.004 0.009 0.89 1.37   % CSF
+            0.0004 0.009 0.89 1.37  % CSF
             0.02  9.0   0.89 1.37   % gray matter
             0.08  40.9  0.84 1.37   % white matter
             0.0   0.0   1.0  1.0];  % air cavity

diff --git a/mcxlab/examples/demo_test_csf.m b/mcxlab/examples/demo_test_csf.m
@@ -30,7 +30,7 @@
 cfg.prop = [0 0 1 1            % medium 0: the environment
             0.019 7.8   0.89 1.37     % medium 1: scalp
             0.019 7.8   0.89 1.37     % medium 2: skull
-            0.004 0.009 0.89 1.37     % medium 3: CSF
+            0.0004 0.009 0.89 1.37    % medium 3: CSF
             0.02  9.0   0.89 1.37     % medium 4: gray matter
             0.08 40.9   0.84 1.37];   % medium 5: white matter (not used)
 

diff --git a/src/mcx_bench.c b/src/mcx_bench.c
@@ -490,7 +490,7 @@ MSTR(
              {"mua": 0.00, "mus": 0.0, "g": 1.00, "n": 1.0},
              {"mua": 0.019, "mus": 7.8182, "g": 0.89, "n": 1.37},
              {"mua": 0.019, "mus": 7.8182, "g": 0.89, "n": 1.37},
-             {"mua": 0.004, "mus": 0.009, "g": 0.89, "n": 1.37},
+             {"mua": 0.0004, "mus": 0.009, "g": 0.89, "n": 1.37},
              {"mua": 0.02, "mus": 9.0, "g": 0.89, "n": 1.37},
              {"mua": 0.08, "mus": 40.9, "g": 0.84, "n": 1.37},
              {"mua": 0.0,"mus": 0.0, "g": 1.0, "n": 1.0}