addpath('utilities_ESRF')
addpath('utilities_LW')
-addpath('octave_packages')
+addpath('octave_packages/image-1.0.15')
+addpath('octave_packages/signal-1.1.3')
################################################
############### INPUT parameters ###############
%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
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