Changeset 1148 for trunk/src/geometry_calib.m

Timestamp:

Jun 21, 2024, 4:51:59 PM (7 days ago)

Author:

sommeria

Message:

civ_3D corrected

File:

• trunk/src/geometry_calib.m (modified) (4 diffs)

trunk/src/geometry_calib.m

@@ -982 +982 @@
 end
 
+
+
+%% read the current image, displayed in the GUI uvmat
+huvmat=findobj(allchild(0),'Name','uvmat');
+UvData=get(huvmat,'UserData');
+A=UvData.Field.A;%currently displayed image
+npxy=size(A);
+
 %% initiate the grid in phys coordinates
 CalibData=get(handles.geometry_calib,'UserData');%get information stored on the GUI geometry_calib
@@ -987 +995 @@
 if isfield(CalibData,'grid')
     grid_input=CalibData.grid;%retrieve the previously used grid
+else
+   %S=skewness(double(reshape(A,1,[])));
+   A=double(A);
+   A=A-mean(mean(A));
+   S=mean(mean(A.*A.*A))/(mean(mean(A.*A)))^1.5
+   grid_input.CheckWhite=sign(S);%propose white markers if image skewness>0, black markers otherwise
 end
 [T,CalibData.grid,CalibData.grid.CheckWhite,CalibData.grid.FilterWindow]=create_grid(grid_input,'detect_grid');%display the GUI create_grid, read the set of phys coordinates T
@@ -992 +1006 @@
 X=[CalibData.grid.x_0 CalibData.grid.x_1 CalibData.grid.x_0 CalibData.grid.x_1]';%corner absissa in the phys coordinates (cm)
 Y=[CalibData.grid.y_0 CalibData.grid.y_0 CalibData.grid.y_1 CalibData.grid.y_1]';%corner ordinates in the phys coordinates (cm)
-
-%% read the current image, displayed in the GUI uvmat
-huvmat=findobj(allchild(0),'Name','uvmat');
-UvData=get(huvmat,'UserData');
-A=UvData.Field.A;%currently displayed image
-npxy=size(A);
 
 %% calculate transform matrices for plane projection: rectangle assumed to be viewed in perspective
@@ -1034 +1042 @@
 Rangx=DataOut.Coord_x;% x coordinates of first and last pixel centres in phys
 Rangy=DataOut.Coord_y;% y coordinates of first and last pixel centres in phys
+
+
+%% inverse the image in the case of black lines
 if CalibData.grid.CheckWhite
     Amod=double(Amod);%case of white grid markers: will look for image maxima