X-Git-Url: https://git.creatis.insa-lyon.fr/pubgit/?a=blobdiff_plain;f=SimuPBI_3WiresPhant_func.m;h=e1250e7d460887760000e54337617ba5acf0e82f;hb=05be8ac4091470e0513b4eaec95fe80969b0ce1c;hp=46b7eaf0bf53002f0cf4678277017481b7e3ca69;hpb=e1e85a737673e93b6b4b01bb0b072cd35dfa3c72;p=CreaPhase.git

diff --git a/SimuPBI_3WiresPhant_func.m b/SimuPBI_3WiresPhant_func.m
index 46b7eaf..e1250e7 100644
--- a/SimuPBI_3WiresPhant_func.m
+++ b/SimuPBI_3WiresPhant_func.m
@@ -26,8 +26,10 @@
 
 function [] = SimuPBI_3WiresPhant_func (vers, oversamp, basename_output, dist, energy, ps, nbProj, range_angle, model_ctf, model_Fresnel, height, dir_out, noise_type, noise_amount)
 
-addpath('/users/lweber/Matlab/SimulationsPBI/CreaPhase/utilities_ESRF')
-addpath('/users/lweber/Matlab/SimulationsPBI/CreaPhase/utilities_LW')
+addpath('utilities_ESRF')
+addpath('utilities_LW')
+addpath('octave_packages/image-1.0.15')
+addpath('octave_packages/signal-1.1.3')
 
 ################################################
 ############### INPUT parameters ############### 
@@ -66,8 +68,15 @@ addpath('/users/lweber/Matlab/SimulationsPBI/CreaPhase/utilities_LW')
 %range_angle=180; 
 
 ## Which model do you want to use for the propagation? use 1 or 0
-%model_ctf=0;
-%model_Fresnel=1;
+## model_ctf refers to the Contrast Transfer Function propagation model. use 1 if you want to simulate propagation with the CTF model, or 0 otherwise.
+## should be equal to 0 or 1 
+% model_ctf = 1
+% model_ctf = 0
+
+## model_Fresnel refers to the Fresnel propagation model. Use 1 if you want to simulate propagation with the fresnel model, or 0 otherwise.
+## should be equal to 0 or 1 
+% model_Fresnel = 1
+% model_Fresnel = 0
 
 ################################
 ############# Parameters related to the object
@@ -158,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