Changeset 17 for trunk/src/geometry_calib.m

Timestamp:

Feb 18, 2010, 7:37:48 PM (15 years ago)

Author:

sommeria

Message:

possibilite de créer une grille physique dans geometry_calib (le faire apres avoir introduit les points sur l'image avec la souris)
petites corrections mineures sur les autres fichiers

File:

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

trunk/src/geometry_calib.m

@@ -109 +109 @@
 %set(htable,'ColumnNames',{'x','y','z','X(pixels)','Y(pixels)'})
 
+%------------------------------------------------------------------------
 % --- Outputs from this function are returned to the command line.
 function varargout = geometry_calib_OutputFcn(hObject, eventdata, handles)
-% varargout  cell array for returning output args (see VARARGOUT);
-% hObject    handle to figure
-% eventdata  reserved - to be defined in a future version of MATLAB
-% handles    structure with handles and user data (see GUIDATA)
-
+%------------------------------------------------------------------------
 % Get default command line output from handles structure
 varargout{1} = handles.output;
@@ -123 +120 @@
 function Phi_Callback(hObject, eventdata, handles)
 
-
-
-
-%--------------------------------------------------
+%------------------------------------------------------------------------
 %read input xml file and update the edit boxes
 function loadfile(handles,fileinput)
-
+%------------------------------------------------------------------------
 %read the input xml file
 t=xmltree(fileinput);
@@ -267 +261 @@
 unit=unitlist{get(handles.CoordUnit,'value')};
 GeometryCalib.CoordUnit=unit;
-
+GeometryCalib.SourceCalib.PointCoord=Object.Coord;
 huvmat=findobj(allchild(0),'Name','uvmat');
 hhuvmat=guidata(huvmat);%handles of elements in the GUI uvmat
@@ -285 +279 @@
     outputfile=answer{1};
 end
-testappend=0;
-if exist(outputfile,'file');%=1 if the output file already exists, 0 else 
-    t=xmltree(outputfile); %read the file
-    backupfile=outputfile;
-    testexist=2;
-    while testexist==2
-        backupfile=[backupfile '~'];% make a backup name by adding  ~ to the xml file name
-        testexist=exist(backupfile,'file');
-    end
-    [success,message]=copyfile(outputfile,backupfile);%make backup   
-    t=xmltree(outputfile); %read the file
-    uid=find(t,'ImaDoc');
-    if ~isequal(uid,1)%if the xml file is not ImaDoc, delete it (after backup)
-        if isequal(success,1)
-            delete(outputfile)
-        else
-            msgbox_uvmat('ERROR',['error in the backup of the existing xml file: ' message])
-            return
-        end
-    else
-        uid_calib=find(t,'ImaDoc/GeometryCalib');
-        testappend=1;
-        if isempty(uid_calib)
-            [t,uid_calib]=add(t,1,'element','GeometryCalib');
-        else %if GeometryCalib already exists, delete its content
-            uid_child=children(t,uid_calib);
-            t=delete(t,uid_child);
-%             testappend=1;
-        end
-    end
-end
-if ~testappend %create a new xml file for calibration data
-    t=xmltree;
-    t=set(t,1,'name','ImaDoc');
-    [t,uid_calib]=add(t,1,'element','GeometryCalib');
-end
-% hgrid=get(handles.REPLICATE,'parent');%read the calibration image source on the interface userdata
-% imagename=get(hgrid,'UserData');
-% if exist(imagename,'file')
-%     GeometryCalib.SourceCalib.ImageCalib=imagename;
+update_imadoc(GeometryCalib,outputfile)
+% testappend=0;
+% if exist(outputfile,'file');%=1 if the output file already exists, 0 else 
+%     t=xmltree(outputfile); %read the file
+%     backupfile=outputfile;
+%     testexist=2;
+%     while testexist==2
+%         backupfile=[backupfile '~'];% make a backup name by adding  ~ to the xml file name
+%         testexist=exist(backupfile,'file');
+%     end
+%     [success,message]=copyfile(outputfile,backupfile);%make backup   
+%     t=xmltree(outputfile); %read the file
+%     uid=find(t,'ImaDoc');
+%     if ~isequal(uid,1)%if the xml file is not ImaDoc, delete it (after backup)
+%         if isequal(success,1)
+%             delete(outputfile)
+%         else
+%             msgbox_uvmat('ERROR',['error in the backup of the existing xml file: ' message])
+%             return
+%         end
+%     else
+%         uid_calib=find(t,'ImaDoc/GeometryCalib');
+%         testappend=1;
+%         if isempty(uid_calib)
+%             [t,uid_calib]=add(t,1,'element','GeometryCalib');
+%         else %if GeometryCalib already exists, delete its content
+%             uid_child=children(t,uid_calib);
+%             t=delete(t,uid_child);
+% %             testappend=1;
+%         end
+%     end
 % end
-GeometryCalib.SourceCalib.PointCoord=Object.Coord;
-t=struct2xml(GeometryCalib,t,uid_calib); 
-save(t,outputfile);
+% if ~testappend %create a new xml file for calibration data
+%     t=xmltree;
+%     t=set(t,1,'name','ImaDoc');
+%     [t,uid_calib]=add(t,1,'element','GeometryCalib');
+% end
+% % hgrid=get(handles.REPLICATE,'parent');%read the calibration image source on the interface userdata
+% % imagename=get(hgrid,'UserData');
+% % if exist(imagename,'file')
+% %     GeometryCalib.SourceCalib.ImageCalib=imagename;
+% % end
+% GeometryCalib.SourceCalib.PointCoord=Object.Coord;
+% t=struct2xml(GeometryCalib,t,uid_calib); 
+% save(t,outputfile);
 msgbox_uvmat('CONFIRMATION',{[outputfile ' updated with calibration data'];...
     ['Error rms (along x,y)=' num2str(GeometryCalib.ErrorRms) ' pixels'];...
@@ -364 +359 @@
 end
 % %record image source
-% hgrid=get(handles.REPLICATE,'parent');%read the calibration image source on the interface userdata
-% imagename=get(hgrid,'UserData');
-% if exist(imagename,'file')
-%     GeometryCalib.SourceCalib.ImageCalib=imagename;
-% end
 GeometryCalib.SourceCalib.PointCoord=Object.Coord;
-
-
-%root PROJETS
 
 %open and read the dataview GUI
@@ -380 +367 @@
 end
 CalibData=get(handles.figure1,'UserData');%read the calibration image source on the interface userdata
-% filename='PROJETS';%default
-% if isfield(CalibData,'XmlInput')
-%      [pp,filename]=fileparts(CalibData.XmlInput);
-% end
-% while ~isequal(filename,'PROJETS') && numel(filename)>1
-%     filename_1=filename;
-%     pp_1=pp;
-%     [pp,filename]=fileparts(pp)
-% end
-% projinput=fullfile(pp_1,filename_1)
-% dd=dataview(projinput)
-
-% 
-% Device=[];%default
-% 
-% h_dataview=dataview;
-% hhdataview=guidata(h_dataview);
-% drawnow
 
 if isfield(CalibData,'XmlInput')
@@ -443 +412 @@
 end
 if testinput
-    outcome=dataview(XmlInput,SubCampaignTest,GeometryCalib)%,SubCampaignTest)
+    outcome=dataview(XmlInput,SubCampaignTest,GeometryCalib)
 end
 %     %A COMPLETER
@@ -955 +924 @@
     set(handles.ListCoord,'String',Tabchar)  
     ListCoord_Callback(hObject, eventdata, handles)
-    PLOT_Callback(hObject,eventdata,handles)
+    MenuPlot_Callback(hObject,eventdata,handles)
 end
 
@@ -1051 +1020 @@
 end
 
-
-
-% --------------------------------------------------------------------
+%------------------------------------------------------------------------
 function MenuCreateGrid_Callback(hObject, eventdata, handles)
+%------------------------------------------------------------------------
 hcalib=get(handles.calib_type,'parent');%handles of the GUI geometry_calib
 CalibData=get(hcalib,'UserData');
@@ -1061 +1029 @@
     Tinput=CalibData.grid;
 end
-T=create_grid(Tinput);%display translate_points GUI and get shift parameters 
+T=create_grid(Tinput)%display the GUI create_grid 
 CalibData.grid=T;
 set(hcalib,'UserData',CalibData)
@@ -1067 +1035 @@
 %grid in phys space
 Coord_cell=get(handles.ListCoord,'String');
-data=read_geometry_calib(Coord_cell);
-data.Coord(:,1)=T(1)+data.Coord(:,1);
-data.Coord(:,2)=T(2)+data.Coord(:,2);
-data.Coord(:,3)=T(3)+data.Coord(:,3);
-data.Coord(:,[4 5])=data.Coord(:,[4 5]);
-for i=1:size(data.Coord,1)
-    for j=1:5
-          Coord{i,j}=num2str(data.Coord(i,j),4);%phys x,y,z
+data=read_geometry_calib(Coord_cell)
+nbpoints=size(data.Coord,1); %nbre of calibration points
+data.Coord(1:size(T,1),1:3)=T;
+for i=1:nbpoints
+   for j=1:5
+          Coord{i,j}=num2str(data.Coord(i,j),4);%display coordiantes with 4 digits
    end
 end
+for i=nbpoints+1:size(data.Coord,1)
+    for j=1:3
+          Coord{i,j}=num2str(data.Coord(i,j),4);%display coordiantes with 4 digits
+    end
+    for j=4:5
+          Coord{i,j}='';%display coordiantes with 4 digi
+    end
+end
+
+
+%size(data.Coord,1)
 Tabchar=cell2tab(Coord,'    |    ');
 set(handles.ListCoord,'Value',1)
 set(handles.ListCoord,'String',Tabchar)
 
-
-
-% --------------------------------------------------------------------
+%-----------------------------------------------------------------------
 function MenuTranslatePoints_Callback(hObject, eventdata, handles)
+%-----------------------------------------------------------------------
 hcalib=get(handles.calib_type,'parent');%handles of the GUI geometry_calib
 CalibData=get(hcalib,'UserData')