Add parameters for fitting edge spread function with sum of erfc at 52 keV

Signed-off-by: Wolfgang Ludwig <wolfgang.ludwig@esrf.fr>

zUtil_ID19/mncf.m

@@ -21,7 +21,7 @@ function mt = mncf(opt, varargin)
 %   dimmensions: 1 | {2} create 1-dim mtf plot or 2-dim mtf map
 
     conf = struct( ...
-        'size', 11, ...
+        'size', 31, ...
         'oversampling', 1, ...
         'dimensions', 2, ...
         'p', []);
@@ -75,10 +75,9 @@ function mt = mncf(opt, varargin)
             ma = 1;
             p = [0.3978, 1.0493, 0.4608, 3.4978, 0.0674, 17.8457, 679.3300];
         case 8
-            % for testing
-            mi = 0;
-            ma = 1;
-            p = [0.3246, 1.1299, 0.4810, 3.4915, 0.1187, 15.4062, 729.8940];
+            mi=0.035;
+            ma=0.9875;
+            p = [0.7288    4.3464    0.1172   50.8488    0.4826  461.0864  493.5530];
         case 9
             % for fits
             mi = 0;

zUtil_Imaging/gtOTFGetMeasured.m

@@ -23,7 +23,7 @@ function [otf, psf] = gtOTFGetMeasured(opt, varargin)
 % Wolfgang/Peter's matlab functions
 
     conf = struct( ...
-        'size', 11, ...
+        'size', 21, ...
         'oversampling', 1, ...
         'p', []);
     conf = parse_pv_pairs(conf, varargin);
@@ -81,10 +81,11 @@ function [otf, psf] = gtOTFGetMeasured(opt, varargin)
             ma = 1;
             p = [0.3978, 1.0493, 0.4608, 3.4978, 0.0674, 17.8457, 679.3300];
         case 8
-            % for testing
-            mi = 0;
-            ma = 1;
-            p = [0.3246, 1.1299, 0.4810, 3.4915, 0.1187, 15.4062, 729.8940];
+            % 0.77 micron pixels, 10X NA=0.23 WBID11, 2x eyepiece, 25
+            % microns freestanding GGG, 52keV (2017) p = [0.7288    4.3464    0.1172   50.8488    0.4826  461.0864  493.5530]
+            mi=0.035;
+            ma=0.9875;
+            p = [0.7288    2.3464    0.1172   50.8488    0.4826  461.0864  493.5530];
         case 9
             % for fits
             mi = 0;