Changeset 1077 for trunk/src/uvmat.m

Timestamp:

Mar 20, 2020, 11:29:01 PM (4 years ago)

Author:

sommeria

Message:

various bug repair

File:

• trunk/src/uvmat.m (modified) (3 diffs)

trunk/src/uvmat.m

@@ -1415 +1415 @@
 end
 
-prompt = {'offset luminosity value in the absence of illumination';'smoothing width for the reference line (in pixels)'};
+prompt = {'offset luminosity value in the absence of illumination';'smoothing width for the reference line (in pixels)';...
+    'LIF Temperature Coeff1';'LIF Temperature Coeff2';'LIF Temperature Coeff3'};
     dlg_title = 'set the parameters for LIF';
-    num_lines= 2;
-    def     = { num2str(BlackOffset);num2str(RefLineWidth)};
+    num_lines= 5;
+    def     = { num2str(BlackOffset);num2str(RefLineWidth);'0';'0';'0'};
     answer = inputdlg(prompt,dlg_title,num_lines,def);
 if isempty(answer) 
@@ -1425 +1426 @@
     XmlData.LIFCalib.BlackOffset=str2double(answer{1}) ;% image value for black background, to be determined by taking images with a cover on the objective lens
     XmlData.LIFCalib.RefLineWidth=str2double(answer{2}) ;% smoothing width used for the reference line
+    XmlData.LIFCalib.TLIFCoeff=[str2double(answer{3}) str2double(answer{4}) str2double(answer{5})];
 end
 
@@ -1443 +1445 @@
 XmlData.LIFCalib.Ray2Coord=LineData{2}.Coord;
 XmlData.LIFCalib.RefLineCoord=LineData{3}.Coord;
+
+%% rescale the image
+[nby,nbx]=size(UvData.Field.A);
+x=linspace(UvData.Field.Coord_x(1),UvData.Field.Coord_x(2),nbx)-nbx/2;
+y=linspace(UvData.Field.Coord_y(1),UvData.Field.Coord_y(2),nby)-nby/2;
+[X,Y]=meshgrid(x,y);
+coeff_quad=0.15*4/(nbx*nbx);% image luminosity reduced by 10% at the edge
+UvData.Field.A=double(UvData.Field.A).*(1+coeff_quad*(X.*X+Y.*Y));
 
 %% display the current image in polar axes with origin at the  illumination source