X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=SimuPBI_3WiresPhant_func.m;h=e1250e7d460887760000e54337617ba5acf0e82f;hb=HEAD;hp=65b550ff0084a76b1b299495d63b7af294b1e371;hpb=3fc60640f9d26ac5d192a16a6df766c3cc850073;p=CreaPhase.git diff --git a/SimuPBI_3WiresPhant_func.m b/SimuPBI_3WiresPhant_func.m index 65b550f..e1250e7 100644 --- a/SimuPBI_3WiresPhant_func.m +++ b/SimuPBI_3WiresPhant_func.m @@ -28,7 +28,8 @@ function [] = SimuPBI_3WiresPhant_func (vers, oversamp, basename_output, dist, e addpath('utilities_ESRF') addpath('utilities_LW') -addpath(genpath('octave_packages')) +addpath('octave_packages/image-1.0.15') +addpath('octave_packages/signal-1.1.3') ################################################ ############### INPUT parameters ############### @@ -166,13 +167,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