change the path of the images contained in "images" directory (attenuation and phase...

diff --git a/README_SimuPBI_3WiresPhant_func.txt b/README_SimuPBI_3WiresPhant_func.txt
index b66979180c0e144240ab6130284275218b5f7122..51c0b39224a4a9c04989621c59c572fa8684b884 100644 (file)
--- a/README_SimuPBI_3WiresPhant_func.txt
+++ b/README_SimuPBI_3WiresPhant_func.txt
@@ -7,6 +7,7 @@
 ## SimuPBI_3WirePhant_func('Analytical', 2, 'test', [0 0.1 0.5], 19, 2.7, 300, 180, 1, 0, '', 'gaussian', 35)
 
 ## Please note that dir_out, noise_type, noise_amount are optional arguments.
+## This method requires the files 'attenuation.edf', and 'delta_beta_map_1200.edf' in the "images"
 
 ################################################
 ############### INPUT parameters ############### 
@@ -68,7 +69,7 @@
 
 
 ############# Parameters related to the object #############
-## in this case, the obejct is a three wire phantom, described in [Frachon, Weber et al., PMB, 2016]. 
+## in this case, the object is a three wire phantom, described in [Frachon, Weber et al., PMB, 2016]. 
 ## height of the object. If height == 1, the simulation is computed in 1D. If height >=2 , the simulation turns into 2D.
 %height=10;
 

diff --git a/SimuPBI_3WiresPhant_func.m b/SimuPBI_3WiresPhant_func.m
index 65b550ff0084a76b1b299495d63b7af294b1e371..98b11ed2c7b9551eff3039fa012a026ad1392563 100644 (file)
--- a/SimuPBI_3WiresPhant_func.m
+++ b/SimuPBI_3WiresPhant_func.m
@@ -166,13 +166,13 @@ betash = lambda*D; % in meter^2
 pad_method='extend';
 
 % Absorption image ('mu', in cm-1)
-mu=edfread('attenuation.edf'); % map of mu in cm-1
+mu=edfread('images/attenuation.edf'); % map of mu in cm-1
 [m n]=size(mu)
 
 absorption=mu*100; % in m-1 , to be consistant with lambda!! /10  to have less attenuation
 beta=(lambda/(4*pi))*absorption; % one slice of beta /!\ units
 
-delta_beta=edfread('delta_beta_map_1200.edf');
+delta_beta=edfread('images/delta_beta_map_1200.edf');
 delta=delta_beta.*beta;
 
 ## Noise addition