Changeset 84 for trunk/src/set_grid.m


Ignore:
Timestamp:
Apr 24, 2010, 2:07:22 PM (15 years ago)
Author:
sommeria
Message:

-update_imadoc: copy the timing information of the movie file avi in the new xml file when a geometry_calibration is performed
-set_grid: modified to produce grids in px coordiantes and to display the produced grid
-geometry_calib: display the point with max error
-civ: bug corrected for copying avi movies to png files
-uvmat: small bug fixes
-px_XYZ: introduce default values for Calib parameters

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/set_grid.m

    r73 r84  
    2222function varargout = set_grid(varargin)
    2323
    24 % Last Modified by GUIDE v2.5 04-Feb-2008 16:05:02
    25 
    26 % Begin initialization code - DO NOT EDIT
     24% Last Modified by GUIDE v2.5 23-Apr-2010 15:44:47
     25
     26% Begin initialization code - DO NOT PLOT
    2727gui_Singleton = 1;
    2828gui_State = struct('gui_Name',       mfilename, ...
    … …  
    4141    gui_mainfcn(gui_State, varargin{:});
    4242end
    43 % End initialization code - DO NOT EDIT
     43% End initialization code - DO NOT PLOT
    4444
    4545%-------------------------------------------------------------------
    … …  
    5858%  if =[] or absent, no plot (mask mode in uvmat)
    5959% parameters on the uvmat interface (obtained by 'get_plot_handle.m')
    60 function set_grid_OpeningFcn(hObject, eventdata, handles,inputfile)
     60function set_grid_OpeningFcn(hObject, eventdata, handles,inputfile,CoordType)
    6161
    6262% Choose default command line output for set_grid
    … …  
    7070% set(hObject,'Position',[0.7 0.1 0.25 0.5])%set the position of the set_grid interface
    7171set(hObject,'DeleteFcn',@closefcn)
    72 set(handles.TITLE,'Value',1)
    73 set(handles.ObjectStyle,'Value',1)
    74 set(handles.ProjMode,'Value',1)
     72% set(handles.TITLE,'Value',1)
     73%set(handles.ObjectStyle,'Value',1)
     74%set(handles.ProjMode,'Value',1)
    7575set(handles.MenuCoord,'ListboxTop',1)
    7676set(handles.MenuCoord,'Value',1);
    … …  
    8080   set(handles.image_2,'String',inputfile)
    8181end
    82 
     82if exist('CoordType','var')
     83    if strcmp(CoordType,'px')
     84        set(handles.MenuCoord,'Value',2)
     85    end
     86end
    8387
    8488% --- Outputs from this function are returned to the command line.
    8589function varargout = set_grid_OutputFcn(hObject, eventdata, handles)
    86 % varargout  cell array for returning output args (see VARARGOUT);
    87 % hObject    handle to figure
    88 % eventdata  reserved - to be defined in a future version of MATLAB
    89 % handles    structure with handles and user data (see GUIDATA)
    90 
    9190% Get default command line output from handles structure
    9291varargout{1} = handles.output;
    9392varargout{2}=handles;
    94 
    95 % --- Executes on selection change in ObjectStyle.
    96 function ObjectStyle_Callback(hObject, eventdata, handles)
    97 
    98 ProjMode_Callback(hObject, eventdata, handles)
    99 
    100 %----------------------------------------------
    101 function xObject_Callback(hObject, eventdata, handles)
    102 
    103 
    104 function yObject_Callback(hObject, eventdata, handles)
    105 
    106 
    107 % --- Executes on selection change in zObject.
    108 function zObject_Callback(hObject, eventdata, handles)
    109 
    110 
    111 %---------------------------------------------------
    112 % --- Executes on selection change in ProjMode.
    113 function ProjMode_Callback(hObject, eventdata, handles)
    114 menu=get(handles.ProjMode,'String');
    115 value=get(handles.ProjMode,'Value');
    116 ProjMode=menu{value};
    117 menu=get(handles.ObjectStyle,'String');
    118 value=get(handles.ObjectStyle,'Value');
    119 ObjectStyle=menu{value};
    120 test3D=isequal(get(handles.ZObject,'Visible'),'on');%3D case
    121 if isequal(ObjectStyle,'plane')||isequal(ObjectStyle,'volume')
    122     set(handles.Phi,'Visible','on')
    123     if test3D%3D case
    124         set(handles.Theta,'Visible','on')
    125         set(handles.Psi,'Visible','on')
    126     end
    127     set(handles.XMin,'Visible','on')
    128     set(handles.XMax,'Visible','on')
    129     set(handles.YMin,'Visible','on')
    130     set(handles.YMax,'Visible','on')
    131     if test3D
    132         set(handles.Theta,'Visible','on')
    133         set(handles.Psi,'Visible','on')
    134         set(handles.ZMin,'Visible','on')
    135         set(handles.ZMax,'Visible','on')
    136     end
    137 else
    138     set(handles.Phi,'Visible','off')
    139     set(handles.Theta,'Visible','off')
    140     set(handles.Psi,'Visible','off')
    141     set(handles.XMin,'Visible','off')
    142     set(handles.XMax,'Visible','off')
    143     set(handles.YMin,'Visible','off')
    144     if isequal(ProjMode,'interp')
    145         set(handles.YMax,'Visible','off')
    146     else
    147         set(handles.YMax,'Visible','on')
    148     end
    149     if isequal(ObjectStyle,'rectangle')|isequal(ObjectStyle,'ellipse')
    150         set(handles.XMax,'Visible','on')
    151     else
    152        set(handles.XMax,'Visible','off')
    153     end
    154     set(handles.ZMin,'Visible','off')
    155     set(handles.ZMax,'Visible','off')
    156 end
    157 if isequal(ProjMode,'projection')|isequal(ProjMode,'inside')|isequal(ProjMode,'outside')|isequal(ObjectStyle,'points')
    158     set(handles.DX,'Visible','off')
    159     set(handles.DY,'Visible','off')
    160     set(handles.DZ,'Visible','off')   
    161 else
    162     set(handles.DX,'Visible','on')
    163     set(handles.DY,'Visible','on')
    164     if test3D%3D case
    165         set(handles.DZ,'Visible','on')
    166     end
    167 end
    168 
    169 %---------------------------------------------
    170 % --- Executes on selection change in TITLE.
    171 function TITLE_Callback(hObject, eventdata, handles)
    172 hsetobject=get(handles.TITLE,'parent');
    173 SetData=get(hsetobject,'UserData');%get the hidden interface data
    174 %      function named CALLBACK in UNTITLED.M with the given input arguments.
    175 menu=get(handles.TITLE,'String');
    176 value=get(handles.TITLE,'Value');
    177 titl=menu{value};
    178 if isequal(titl,'POINTS')
    179      menu_style={'points'};
    180      menu_proj={'projection';'interp';'filter';'none'};
    181 elseif isequal(titl,'LINE')
    182      menu_style={'line';'polyline';'rectangle';'polygon';'ellipse'};%'line' =default
    183      menu_proj={'projection';'interp';'filter';'none'};
    184 elseif isequal(titl,'PATCH')
    185      menu_style={'rectangle';'polygon';'ellipse'};%'line' =default
    186      menu_proj={'inside';'outside';'none'};
    187  elseif isequal(titl,'PLANE')
    188      menu_style={'plane'};
    189      menu_proj={'projection';'interp'};
    190 elseif isequal(titl,'VOLUME')
    191      menu_style={'volume'};
    192      menu_proj={'none'};
    193  
    194 end
    195 set(handles.ObjectStyle,'String',menu_style)
    196 set(handles.ObjectStyle,'Value',1)
    197 set(handles.ProjMode,'String',menu_proj)
    198 set(handles.ProjMode,'Value',1)
    199 if isfield(SetData,'ParentButton')
    200     update_parentbutton(SetData.ParentButton,titl)
    201 end
    202 ObjectStyle_Callback(hObject, eventdata, handles) 
    203 
    204 %-----------
    205 function update_parentbutton(ParentButton,titl)
    206 
    207 if isstruct(ParentButton)
    208     parentfields=fields(ParentButton);
    209     for ibutton=1:length(parentfields)
    210         buttonhandle=eval(['ParentButton.' parentfields{ibutton}]);
    211         if ishandle(buttonhandle)
    212             set(buttonhandle,'Value',0)
    213             set(buttonhandle,'BackgroundColor',[0 1 0])%put unactivated buttons to green
    214         end
    215     end
    216     if isfield(ParentButton,titl)
    217        buttonhandle=eval(['ParentButton.' titl]);
    218        if ishandle(buttonhandle)
    219             set(buttonhandle,'Value',1)
    220             set(buttonhandle,'BackgroundColor',[1 1 0])%put activated button to yellow
    221        end
    222     end
    223 end
    224 %------------
    225 function Phi_Callback(hObject, eventdata, handles)
    226 update_slider(hObject, eventdata,handles)
    227 
    228 function Theta_Callback(hObject, eventdata, handles)
    229 update_slider(hObject, eventdata,handles)
    230 
    231 function update_slider(hObject, eventdata,handles)
    232 %rotation angles
    233 Phi=(pi/180)*str2num(get(handles.Phi,'String'));%first Euler angle in radian
    234 Theta=(pi/180)*str2num(get(handles.Theta,'String'));%second Euler angle in radian
    235 
    236 %components of the unitiy vector normal to the projection plane
    237 NormVec_X=-sin(Phi)*sin(Theta);
    238 NormVec_Y=cos(Phi)*sin(Theta);
    239 NormVec_Z=cos(Theta);
    240 huvmat=findobj('Tag','uvmat');%find the current uvmat interface handle
    241 UvData=get(huvmat,'UserData');%Data associated to the current uvmat interface
    242 Z=NormVec_X *(UvData.X)+NormVec_Y *(UvData.Y)+NormVec_Z *(UvData.Z);
    243 set(handles.z_slider,'Min',min(Z))
    244 set(handles.z_slider,'Max',max(Z))
    245 ZMax_Callback(hObject, eventdata, handles)
    246 
    247 function DX_Callback(hObject, eventdata, handles)
    248 
    249 
    250 function DY_Callback(hObject, eventdata, handles)
    251 
    252 
    253 function DZ_Callback(hObject, eventdata, handles)
    254 
    25593
    25694
    … …  
    362200
    363201%-----------------------------------------------------------------------
    364 % --- Executes on button press in edit: PLOT the defined object and its projected field
    365 function edit_Callback(hObject, eventdata, handles)
    366 hsetobject=get(hObject,'parent');
    367 SetData=get(hsetobject,'UserData');%get the hidden interface data
    368 %IndexObj=SetData.IndexObj%index of the current projection object in the list of projection objects (UvData.ProjObject)
    369 huvmat=findobj('Tag','uvmat');%find the current uvmat interface handle
    370 UvData=get(huvmat,'UserData');%Data associated to the current uvmat interface
    371 if isfield(UvData,'CuurentObjectIndex')
    372     IndexObj=UvData.CurrentObjectIndex;
    373 else
    374     IndexObj=[];
    375 end
    376 ObjectData=read_set_grid(handles);%read the interface input parameters defining the object
    377 [UvData,IndexObj]=update_obj(UvData,IndexObj,ObjectData,SetData.PlotHandles);
    378 uvmat('write_plot_param',PlotHandles,UvData.Object{IndexObj}.PlotParam); %update the display of plotting parameters for the current object
    379 SetData.IndexObj=IndexObj;
    380 set(gcbf,'UserData',SetData)%update object index in the set_grid interface
    381 set(huvmat,'UserData',UvData)%update the data in the uvmat interface
     202% --- Executes on button press in plot: PLOT the defined object and its projected field
     203function plot_Callback(hObject, eventdata, handles)
     204grid_pix_A=get_grid(handles);
     205huvmat=uvmat(get(handles.image_1,'String'));
     206hhuvmat=guidata(huvmat);
     207set(hhuvmat.transform_fct,'Value',1)
     208uvmat('run0_Callback',hObject,eventdata,hhuvmat); %file input with xml reading  in uvmat
     209axes(hhuvmat.axes3);
     210hold on
     211plot(grid_pix_A(:,1),grid_pix_A(:,2),'.')
     212
     213% --- Executes on button press in plot_2.
     214function plot_2_Callback(hObject, eventdata, handles)
     215[grid_pix_A,grid_pix_B]=get_grid(handles);
     216huvmat=uvmat(get(handles.image_2,'String'));
     217hhuvmat=guidata(huvmat);
     218set(hhuvmat.transform_fct,'Value',1)
     219uvmat('run0_Callback',hObject,eventdata,hhuvmat); %file input with xml reading  in uvmat
     220axes(hhuvmat.axes3);
     221hold on
     222plot(grid_pix_B(:,1),grid_pix_B(:,2),'.')
     223
    382224
    383225
    … …  
    401243delete_object(IndexObj);
    402244
    403 %----------------------------------------------------
    404 function YMin_Callback(hObject, eventdata, handles)
    405 % hObject    handle to YMin (see GCBO)
    406 % eventdata  reserved - to be defined in a future version of MATLAB
    407 % handles    structure with handles and user data (see GUIDATA)
    408 
    409 % Hints: get(hObject,'String') returns contents of YMin as text
    410 %        str2double(get(hObject,'String')) returns contents of YMin as a double
    411 
    412 
    413 function ZMin_Callback(hObject, eventdata, handles)
    414 
    415 
    416 function ZMax_Callback(hObject, eventdata, handles)
    417 DZ=str2num(get(handles.ZMax,'String'));
    418 ZMin=get(handles.z_slider,'Min');
    419 ZMax=get(handles.z_slider,'Max');
    420 rel_step(1)=DZ/(ZMax-ZMin);
    421 rel_step(2)=0.2;
    422 set(handles.z_slider,'SliderStep',rel_step)
    423 
    424 function YMax_Callback(hObject, eventdata, handles)
    425 
    426 
    427 function XMin_Callback(hObject, eventdata, handles)
    428 
    429 
    430 function XMax_Callback(hObject, eventdata, handles)
    431 
    432245
    433246% ------------------------------------------------------
    434247function save_Callback(hObject, eventdata, handles)
    435248% ------------------------------------------------------
    436 Object=read_set_object(handles);%read the set_grid interface;
    437 DX=Object.DX;
    438 DY=Object.DY;
    439 RangeX=Object.RangeX;
    440 RangeY=Object.RangeY;
    441  array_realx=[RangeX(2):DX:RangeX(1)];
    442  array_realy=[RangeY(2):DY:RangeY(1)];
    443  nx_patch=length(array_realx);
    444  ny_patch=length(array_realy);
    445  [grid_realx,grid_realy]=meshgrid(array_realx,array_realy);
    446  grid_real(:,1)=reshape(grid_realx,nx_patch*ny_patch,1);
    447  grid_real(:,2)=reshape(grid_realy,nx_patch*ny_patch,1);
    448  grid_real(:,3)=zeros(nx_patch*ny_patch,1);
     249[grid_pix_A,grid_pix_B]=get_grid(handles);
     250
     251 %ECRIRE FICHIERS
     252nbpointsA=size(grid_pix_A);
     253XA=grid_pix_A(:,1);
     254YA=grid_pix_A(:,2);
     255unitcolumn=32*ones(size(XA));
     256Xchar=num2str(XA);
     257blanc=char(unitcolumn);
     258Ychar=num2str(YA);
     259tete=['1 ' num2str(nbpointsA(1))];
     260txt=[Xchar blanc Ychar];
     261textgrid={tete;txt};
     262textout=char(textgrid);
     263imageA=get(handles.image_1,'String');
     264[Pathsub]=name2display(imageA);
     265Answer = msgbox_uvmat('INPUT_TXT','grid file name (*.grid)',fullfile(Pathsub,'gridA.grid'));
     266% Answer = inputdlg('grid file name (*.grid)',' ',1,{fullfile(Pathsub,'gridA.grid')},'on');
     267dlmwrite(Answer,textout,'');
     268msgbox_uvmat('CONFIRMATION',[Answer ' written as ASCII text file']);
     269if ~isempty(grid_pix_B)
     270    nbpointsB=size(grid_pix_B);
     271    XB=grid_pix_B(:,1);
     272    YB=grid_pix_B(:,2);
     273    unitcolumn=32*ones(size(XB));
     274    Xchar=num2str(XB);
     275    blanc=char(unitcolumn);
     276    Ychar=num2str(YB);
     277    tete=['1 ' num2str(nbpointsB(1))];
     278    txt=[Xchar blanc Ychar];
     279    textgrid={tete;txt};
     280    textout=char(textgrid);
     281    Answer = msgbox_uvmat('INPUT_TXT','grid file name (*.grid)',fullfile(Pathsub,'gridB.grid'));
     282    dlmwrite(Answer,textout,'');
     283    msgbox_uvmat('CONFIRMATION',[Answer ' written as ASCII text file']);
     284end
     285
     286%-------------------------
     287function [grid_pix_A,grid_pix_B]=get_grid(handles);
     288%Object=read_set_object(handles);%read the set_grid interface;
     289grid_pix_B=[];%default
     290DX=str2num(get(handles.DX,'String'));
     291DY=str2num(get(handles.DY,'String'));
     292XMin=str2num(get(handles.XMin,'String'));
     293XMax=str2num(get(handles.XMax,'String'));
     294YMin=str2num(get(handles.YMin,'String'));
     295YMax=str2num(get(handles.YMax,'String'));
     296array_realx=[XMin:DX:XMax];
     297array_realy=[YMin:DY:YMax];
     298nx_patch=length(array_realx);
     299ny_patch=length(array_realy);
     300[grid_realx,grid_realy]=meshgrid(array_realx,array_realy);
     301grid_real(:,1)=reshape(grid_realx,nx_patch*ny_patch,1);
     302grid_real(:,2)=reshape(grid_realy,nx_patch*ny_patch,1);
     303grid_real(:,3)=zeros(nx_patch*ny_patch,1);
    449304 
    450305imageA=get(handles.image_1,'String');
    … …  
    482337     tsaiA=[];
    483338end
    484 size(grid_real)
    485 tsaiA
    486 if isempty(tsaiA)
     339MenuCoord=get(handles.MenuCoord,'String');
     340val=get(handles.MenuCoord,'Value');
     341if isempty(tsaiA)||strcmp(MenuCoord{val},'px')
    487342    grid_imaA(:,1)=grid_real(:,1);
    488343    grid_imaA(:,2)=grid_real(:,2);
    … …  
    544399grid_real2(:,2)=grid_real_y;
    545400grid_real2(:,3)=zeros(nx_patch_new,1);
    546 if isempty(tsaiA)
     401if isempty(tsaiA)||strcmp(MenuCoord{val},'px')
    547402    grid_pix_A(:,1)=grid_real2(:,1);
    548403   grid_pix_A(:,2)= grid_real2(:,2);
    … …  
    554409end
    555410
    556  %ECRIRE FICHIERS
    557 nbpointsA=size(grid_pix_A);
    558 XA=grid_pix_A(:,1);
    559 YA=grid_pix_A(:,2);
    560 unitcolumn=32*ones(size(XA));
    561 Xchar=num2str(XA);
    562 blanc=char(unitcolumn);
    563 Ychar=num2str(YA);
    564 tete=['1 ' num2str(nbpointsA(1))];
    565 txt=[Xchar blanc Ychar];
    566 textgrid={tete;txt};
    567 textout=char(textgrid);
    568 Answer = msgbox_uvmat('INPUT_TXT','grid file name (*.grid)',fullfile(Pathsub,'gridA.grid'));
    569 % Answer = inputdlg('grid file name (*.grid)',' ',1,{fullfile(Pathsub,'gridA.grid')},'on');
    570 dlmwrite(Answer,textout,'');
    571 msgbox_uvmat('CONFIRMATION',[Answer ' written as ASCII text file']);
    572 if testB
    573     nbpointsB=size(grid_pix_B);
    574     XB=grid_pix_B(:,1);
    575     YB=grid_pix_B(:,2);
    576     unitcolumn=32*ones(size(XB));
    577     Xchar=num2str(XB);
    578     blanc=char(unitcolumn);
    579     Ychar=num2str(YB);
    580     tete=['1 ' num2str(nbpointsB(1))];
    581     txt=[Xchar blanc Ychar];
    582     textgrid={tete;txt};
    583     textout=char(textgrid);
    584     Answer = msgbox_uvmat('INPUT_TXT','grid file name (*.grid)',fullfile(Pathsub,'gridB.grid'));
    585     dlmwrite(Answer,textout,'');
    586     msgbox_uvmat('CONFIRMATION',[Answer ' written as ASCII text file']);
    587 end
    588 
    589 
    590 % --- Executes on slider movement.
    591 function z_slider_Callback(hObject, eventdata, handles)
    592 %A ADAPTER
    593 Z_value=get(handles.z_slider,'Value');
    594 
    595 %rotation angles
    596 Phi=(pi/180)*str2num(get(handles.Phi,'String'));%first Euler angle in radian
    597 Theta=(pi/180)*str2num(get(handles.Theta,'String'));%second Euler angle in radian
    598 
    599 %components of the unity vector normal to the projection plane
    600 NormVec_X=-sin(Phi)*sin(Theta);
    601 NormVec_Y=cos(Phi)*sin(Theta);
    602 NormVec_Z=cos(Theta);
    603 
    604 %set new plane position and update graph
    605 set(handles.XObject,'String',num2str(NormVec_X*Z_value))
    606 set(handles.YObject,'String',num2str(NormVec_Y*Z_value))
    607 set(handles.ZObject,'String',num2str(NormVec_Z*Z_value))
    608 edit_Callback(hObject, eventdata, handles)
    609 
    610 
    611 
    612 function XObject_Callback(hObject, eventdata, handles)
    613 
    614 
    615 function YObject_Callback(hObject, eventdata, handles)
    616 
    617 
    618 
    619 
    620 function ZObject_Callback(hObject, eventdata, handles)
    621 
    622 
    623 function image_2_Callback(hObject, eventdata, handles)
    624 % hObject    handle to image_2 (see GCBO)
    625 % eventdata  reserved - to be defined in a future version of MATLAB
    626 % handles    structure with handles and user data (see GUIDATA)
    627 
    628 % Hints: get(hObject,'String') returns contents of image_2 as text
    629 %        str2double(get(hObject,'String')) returns contents of image_2 as a double
    630 
    631 
    632 
    633 function image_1_Callback(hObject, eventdata, handles)
    634 % hObject    handle to image_1 (see GCBO)
    635 % eventdata  reserved - to be defined in a future version of MATLAB
    636 % handles    structure with handles and user data (see GUIDATA)
    637 
    638 % Hints: get(hObject,'String') returns contents of image_1 as text
    639 %        str2double(get(hObject,'String')) returns contents of image_1 as a double
    640 
    641 
     411
     412
     413%------------------------------------------------------------------------
    642414% --- Executes on button press in HELP.
    643415function HELP_Callback(hObject, eventdata, handles)
     416%------------------------------------------------------------------------
    644417path_to_uvmat=which ('uvmat');% check the path of uvmat
    645418pathelp=fileparts(path_to_uvmat);
    646 helpfile=fullfile(pathelp,'UVMAT_DOC','uvmat_doc.html');
     419helpfile=fullfile(pathelp,'uvmat_doc','uvmat_doc.html');
    647420if isempty(dir(helpfile)), errordlg('Please put the help file uvmat_doc.html in the directory UVMAT/UVMAT_DOC')
    648421else
    … …  
    650423end
    651424
     425
     426
Note: See TracChangeset for help on using the changeset viewer.