Update fitGauss_samerho.m

Author: quantitativeTEM

functions/FitProgram/fitGauss_samerho.m

@@ -51,26 +51,9 @@
 
 %% Recalculate size of boxes used to cut out parts of the original image
 Box = round(FP.dist/input.dx)/2*ones(max(type),1); %calcBoxSize(FP.obs_bs,input.coordinates,input.dx);
-% % ISTEM
-% Box = [30;30;15];
-% % ABF Dimi
-% Box = [8;5;8];
-% % ABF DF4 Nico Armand
-% Box = [12;12;6];
-% % ABF HAADF Nico Armand
-% Box = [23.5;23.5;16];
-% % LAADF Nico Armand
-% Box = [17.5;17.5;9];
-% % IDPC
-% Box(3) = Box(1)/2;
-
-%% Zeta columnwise
-% obs_back = smoothBack(Box,FP.K,FP.L,FP.X,FP.Y,input.obs,input.dx);
-obs_back = 0;
 
 %% startmodel
-input_obs_bs = input.obs-obs_back;
-reshape_obs_bs = reshape(input_obs_bs,FP.K*FP.L,1);
+reshape_obs_bs = reshape(input.obs,FP.K*FP.L,1);
 Ga = sparse(FP.K*FP.L,FP.n_c+FP.fitzeta);
 job = cell(FP.numWorkers,1);
 max_n = 70-6*FP.numWorkers;
@@ -146,7 +129,7 @@
     Estimated_eta = thetalin;
     Estimated_zeta = FP.zeta;
 end
-obs_bs = input_obs_bs - Estimated_zeta;
+obs_bs = input.obs - Estimated_zeta;
 eta_estimatedbackground = Estimated_eta;
 
 %% options for lsqnonlin
@@ -311,7 +294,7 @@
     if FP.fitzeta
         Estimated_eta = thetalin(1:end-1);
         Estimated_zeta = thetalin(end);
-        obs_bs = input_obs_bs - Estimated_zeta;
+        obs_bs = input.obs - Estimated_zeta;
     else
         Estimated_eta = thetalin;
     end
@@ -374,4 +357,4 @@
     % fit
     StartParametersnonlin_s = [betaX_estimatedbackground(i) betaY_estimatedbackground(i)];
     EstimatedParametersnonlin(:,n) = lsqnonlin('criterionGauss_samerho',StartParametersnonlin_s',[],[],options,reshapeX_s,reshapeY_s,K_s,L_s,obs_s_back,rho_estimatedbackground(i));
-end
+end