From: Loriane Weber Date: Tue, 13 Sep 2016 09:47:09 +0000 (+0200) Subject: change the path of the images contained in "images" directory (attenuation and phase... X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?p=CreaPhase.git;a=commitdiff_plain;h=f71159f19dfee5fa7a4ce6aa09f41ce2fb5e125b change the path of the images contained in "images" directory (attenuation and phase maps). --- diff --git a/README_SimuPBI_3WiresPhant_func.txt b/README_SimuPBI_3WiresPhant_func.txt index b669791..51c0b39 100644 --- 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 65b550f..98b11ed 100644 --- 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