source: trunk/src/set_object.m @ 864

Last change on this file since 864 was 864, checked in by sommeria, 9 years ago

bug_repaired

File size: 34.7 KB
RevLine 
[204]1%'set_object': GUI to edit a projection object
2%------------------------------------------------------------------------
3% function hset_object= set_object(data, PlotHandles,ZBounds)
4% associated with the GUI set_object.fig
5%
6% OUTPUT:
7% hset_object: handle of the GUI figure
8%
[445]9% INPUT:u
[204]10% data: structure describing the object properties
11%    .Style=...
12%    .ProjMode
13%    .CoordType: 'phys' or 'px'
[379]14%    .num_DX,.num_DY,.num_DZ : mesh along each dirction
[204]15%    .RangeX, RangeY
16%    .Coord(j,i), i=1, 2, 3,  components x, y, z of j=1...n position(s) characterizing the object components
17% PlotHandles: handles for projection plots NO MORE USED
18% Zbounds: bounds on Z ( 3D case)
[809]19
20%=======================================================================
21% Copyright 2008-2014, LEGI UMR 5519 / CNRS UJF G-INP, Grenoble, France
22%   http://www.legi.grenoble-inp.fr
23%   Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
[204]24%
25%     This file is part of the toolbox UVMAT.
[809]26%
[204]27%     UVMAT is free software; you can redistribute it and/or modify
[809]28%     it under the terms of the GNU General Public License as published
29%     by the Free Software Foundation; either version 2 of the license,
30%     or (at your option) any later version.
31%
[204]32%     UVMAT is distributed in the hope that it will be useful,
33%     but WITHOUT ANY WARRANTY; without even the implied warranty of
34%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
[809]35%     GNU General Public License (see LICENSE.txt) for more details.
36%=======================================================================
[204]37
38function varargout = set_object(varargin)
39
[850]40% Last Modified by GUIDE v2.5 16-Jan-2015 11:03:00
[204]41
[680]42% Begin initialization code - DO NOT REFRESH
[204]43gui_Singleton = 1;
44gui_State = struct('gui_Name',       mfilename, ...
45                   'gui_Singleton',  gui_Singleton, ...
46                   'gui_OpeningFcn', @set_object_OpeningFcn, ...
47                   'gui_OutputFcn',  @set_object_OutputFcn, ...
48                   'gui_LayoutFcn',  [] , ...
49                   'gui_Callback',   []);
50if nargin & ischar(varargin{1})
51    gui_State.gui_Callback = str2func(varargin{1});
52end
53
54if nargout
55    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
56else
57    gui_mainfcn(gui_State, varargin{:});
58end
[680]59% End initialization code - DO NOT REFRESH
[204]60%------------------------------------------------------------------------
61%------------------------------------------------------------------------
62% --- Executes just before set_object is made visible.
63%INPUT:
64% handles: handles of the set_object interface elements
65%'IndexObj': NON USED ANYMORE (To suppress) index of the object (on the UvData list) that set_object will modify
66%        if =[] or absent: index still undefined (create mode in uvmat)
[680]67%        if=0; no associated object (used for series), the button 'REFRESH' is  then unvisible
[204]68%'data': read from an existing object selected in the interface
[379]69%      .Name : class of object ('POINTS','LINE',....)
70%      .num_DX,num_DY,num_DZ; meshes for regular grids
[204]71%      .Coord: object position coordinates
72%      .ParentButton: handle of the uicontrol object calling the interface
73% PlotHandles: set of handles of the elements contolling the plotting of the projected field:
[680]74%  if =[] or absent, no refresh (mask mode in uvmat)
[204]75% parameters on the uvmat interface (obtained by 'get_plot_handle.m')
76function set_object_OpeningFcn(hObject, eventdata, handles, data, PlotHandles,ZBounds)
77%-------------------------------------------------------------------
78% Choose default command line output for set_object
79handles.output = hObject;
80% Update handles structure
81guidata(hObject, handles);
82
[613]83%% position
[657]84set(0,'Unit','pixels')
[613]85ScreenSize=get(0,'ScreenSize');% get the size of the screen, to put the fig on the upper right
[705]86PosGUI=get(handles.set_object,'Position');% fig width in pixels
87Width=PosGUI(3);%width of the gui set_object in pixels
88Height=PosGUI(4);
[613]89Left=ScreenSize(3)- Width-40; %right edge close to the right, with margin=40
90Bottom=ScreenSize(4)-Height-40; %put fig at top right
[657]91set(handles.set_object,'Unit','pixels')
[613]92set(handles.set_object,'Position',[Left Bottom Width Height])
93
[204]94%default
[206]95if ~exist('ZBounds','var')
96    ZBounds=0; %default
[204]97end
[309]98set(hObject,'WindowButtonDownFcn',{'mouse_down'})%set mouse click action function
[402]99set(hObject,'DeleteFcn',{@closefcn})
[204]100
101% fill the interface as set in the input data:
102if exist('data','var')
[379]103    if isfield(data,'Coord') &&size(data.Coord,2)==3
104        set(handles.z_slider,'Visible','on')
105    else
[204]106        set(handles.z_slider,'Visible','off')
107    end
[387]108    if isfield(data,'TypeMenu')
109        set(handles.Type,'String',data.TypeMenu)
110    end
111    if isfield(data,'ProjModeMenu')
112        set(handles.ProjMode,'UserData',data.ProjModeMenu)% data.ProjModeMenu as default menu (used in Type_Callback)
113    end
[591]114    errormsg=fill_GUI(data,handles.set_object);
[410]115    if ~isempty(errormsg)
116        msgbox_uvmat('ERROR','bad data input in set_object')
117        return
118    end
[387]119    Type_Callback(hObject, eventdata, handles)% update the GUI set_object depending on the object type   
[681]120    set(handles.REFRESH,'BackgroundColor',[1 0 0])
[204]121    if isfield(data,'RangeZ') && length(ZBounds) >= 2
[379]122        set(handles.num_RangeZ_2,'String',num2str(max(data.RangeZ),3))
[204]123        DZ=max(data.RangeZ);%slider step
[206]124        if ~isnan(ZBounds(1)) && ZBounds(2)~=ZBounds(1)
[204]125            rel_step(1)=min(DZ/(ZBounds(2)-ZBounds(1)),0.2);%must be smaller than 1
126            rel_step(2)=0.1;
127            set(handles.z_slider,'Visible','on')
128            set(handles.z_slider,'Min',ZBounds(1))
129            set(handles.z_slider,'Max',ZBounds(2))
130            set(handles.z_slider,'SliderStep',rel_step)
131            set(handles.z_slider,'Value',(ZBounds(1)+ZBounds(2))/2)
132        end
133    end
134    if isfield(data,'RangeX')
135        if ischar(data.RangeX)
136            data.RangeX=str2num(data.RangeX);
137        end
[379]138        set(handles.num_RangeX_2,'String',num2str(max(data.RangeX),3))
139        set(handles.num_RangeX_1,'String',num2str(min(data.RangeX),3))
[204]140    end
141    if isfield(data,'RangeY')
142        if ischar(data.RangeY)
143            data.RangeY=str2num(data.RangeY);
144        end
[379]145        set(handles.num_RangeY_2,'String',num2str(max(data.RangeY),3))
146        set(handles.num_RangeY_1,'String',num2str(min(data.RangeY),3))
[204]147    end
148    if isfield(data,'RangeZ')
149        if ischar(data.RangeZ)
150            data.RangeZ=str2num(data.RangeZ);
151        end
[379]152        set(handles.num_RangeZ_2,'String',num2str(max(data.RangeZ),3))
[204]153        if numel(data.RangeZ)>=2
[379]154            set(handles.num_RangeZ_1,'String',num2str(min(data.RangeZ),3))
[204]155        end
156    end 
[625]157    if ~isfield(data,'Angle')
158        data.Angle=[0 0 0];
159    end
160%     if isfield(data,'Angle') && isequal(numel(data.Angle),3)
[379]161         set(handles.num_Angle_1,'String',num2str(data.Angle(1)))
162         set(handles.num_Angle_2,'String',num2str(data.Angle(2)))
163         set(handles.num_Angle_3,'String',num2str(data.Angle(3)))
[625]164%     end
[204]165end
[625]166set(get(handles.set_object,'children'),'enable','off')
167set(handles.SAVE,'enable','on')
[680]168% set(handles.REFRESH,'enable','off')
[204]169
[625]170
[204]171%------------------------------------------------------------------------
172% --- Outputs from this function are returned to the command line.
173function varargout = set_object_OutputFcn(hObject, eventdata, handles)
174%------------------------------------------------------------------------
175% Get default command line output from handles structure
176varargout{1} = handles.output;
177varargout{2}=handles;
178
179%------------------------------------------------------------------------
[402]180% executed when closing the GUI set_object
181function closefcn(gcbo,eventdata)
182%------------------------------------------------------------------------
183huvmat=findobj(allchild(0),'Tag','uvmat');%find the current uvmat interface handle
184if ~isempty(huvmat)
185    hhuvmat=guidata(huvmat);
[650]186    set(hhuvmat.CheckViewObject,'value',0)%
187    set(hhuvmat.CheckEditObject,'Value',0)% desactivate the edit option
[402]188    % deselect the object in ListObject when view_field is closed
189    if isempty(findobj(allchild(0),'Tag','view_field'))
190        ObjIndex=get(hhuvmat.ListObject,'Value');
191        ObjIndex=ObjIndex(1);%keep only the first object selected
192        set(hhuvmat.ListObject,'Value',ObjIndex)
193        % draw all object colors in blue (unselected) in uvmat
[511]194        hother=[findobj(hhuvmat.PlotAxes,'Tag','proj_object');findobj(hhuvmat.PlotAxes,'Tag','DeformPoint')];%find all the proj object and deform point representations
[402]195        for iobj=1:length(hother)
196            if isequal(get(hother(iobj),'Type'),'rectangle')||isequal(get(hother(iobj),'Type'),'patch')
197                set(hother(iobj),'EdgeColor','b')
198                if isequal(get(hother(iobj),'FaceColor'),'m')
199                    set(hother(iobj),'FaceColor','b')
200                end
201            elseif isequal(get(hother(iobj),'Type'),'image')
202                Acolor=get(hother(iobj),'CData');
203                Acolor(:,:,1)=zeros(size(Acolor,1),size(Acolor,2));
204                set(hother(iobj),'CData',Acolor);
205            else
206                set(hother(iobj),'Color','b')
207            end
208            set(hother(iobj),'Selected','off')
209        end
210    end
211end
212hseries=findobj(allchild(0),'Name','series');%find the current series interface handle
213if ~isempty(hseries)
214    hhseries=guidata(hseries);
[667]215    set(hhseries.EditObject,'Value',0)
[402]216end
217
218
219%------------------------------------------------------------------------
[379]220% --- Executes on selection change in Type.
221function Type_Callback(hObject, eventdata, handles)
[204]222%------------------------------------------------------------------------
[681]223
[379]224ListType=get(handles.Type,'String');
225Type=ListType{get(handles.Type,'Value')};
[204]226% make correspondance between different object styles
[379]227Coord=get(handles.Coord,'Data');
228
229%% set the number of lines in the Coord table depending on object type
230switch Type
231    case{'line'}
232        if size(Coord,1)<2
233            if isequal(size(Coord,2),3)
234                Coord=[Coord; 0 0 0];%add a line for edition (3D case)
235            else
236                Coord=[Coord; 0 0]; %add a line for edition (2D case)
237            end
238        else
239            Coord=Coord(1:2,:);
240        end
241    case{'rectangle','ellipse','plane','volume'}
242        Coord=Coord(1,:);
[204]243end
[379]244set(handles.Coord,'Data',Coord)
[204]245
[379]246%% set the projection menu and the corresponding options
[387]247if isempty(get(handles.ProjMode,'UserData'))
248    switch Type
[648]249        case {'points','line','plane'}
[581]250            menu_proj={'projection';'interp_lin';'interp_tps';'none'};
[648]251        case 'polyline'
252            menu_proj={'interp_lin';'interp_tps';'none'};
[387]253        case {'polygon','rectangle','ellipse'}
[648]254            menu_proj={'inside';'outside';'mask_inside';'mask_outside';'interp_lin';'interp_tps';'none'};
[387]255        case 'volume'
[581]256            menu_proj={'interp_lin';'none'};
[387]257        otherwise
[581]258            menu_proj={'projection';'interp_lin';'interp_tps';'none'};%default
[387]259    end
260else
261    menu_proj=get(handles.ProjMode,'UserData');
262end
[410]263ProjModeList=get(handles.ProjMode,'String');
264menu_index=find(strcmp(ProjModeList{get(handles.ProjMode,'Value')},menu_proj));
265if isempty(menu_index)
266    menu_index=1;%
[204]267end
[410]268set(handles.ProjMode,'Value',menu_index);% value index must not exceed the menu length
[204]269set(handles.ProjMode,'String',menu_proj)
270ProjMode_Callback(hObject, eventdata, handles)
[379]271
[204]272%------------------------------------------------------------------------
273% --- Executes on selection change in ProjMode.
[675]274%------------------------------------------------------------------------
[681]275function ProjMode_Callback(hObject, eventdata, handles)
[675]276
[681]277set(handles.REFRESH,'BackgroundColor',[1 0 1])
[204]278menu=get(handles.ProjMode,'String');
279value=get(handles.ProjMode,'Value');
280ProjMode=menu{value};
[379]281menu=get(handles.Type,'String');
282value=get(handles.Type,'Value');
[204]283ObjectStyle=menu{value};
[379]284%%%%%%%%% TODO
285test3D=0; %TODO: update  test3D=isequal(get(handles.ZObject,'Visible'),'on');%3D case
286%%%%%%%%%
[204]287%default setting
[379]288set(handles.num_Angle_1,'Visible','off')
289set(handles.num_Angle_2,'Visible','off')
290set(handles.num_Angle_3,'Visible','off')
291set(handles.num_RangeX_1,'Visible','off')
292set(handles.num_RangeX_2,'Visible','off')
293set(handles.num_RangeY_1,'Visible','off')
[589]294if isequal(ProjMode,'interp_lin')|| isequal(ProjMode,'interp_tps')
[379]295    set(handles.num_RangeY_2,'Visible','off')
[204]296else
[379]297    set(handles.num_RangeY_2,'Visible','on')
[204]298end
299if strcmp(ObjectStyle,'rectangle')||strcmp(ObjectStyle,'ellipse')
[379]300    set(handles.num_RangeX_2,'Visible','on')
[204]301else
[379]302   set(handles.num_RangeX_2,'Visible','off')
[204]303end
[379]304set(handles.num_RangeZ_1,'Visible','off')
305set(handles.num_RangeZ_2,'Visible','off')
306set(handles.num_DX,'Visible','off')
307set(handles.num_DY,'Visible','off')
308set(handles.num_DZ,'Visible','off')
[204]309
310switch ObjectStyle
311    case 'points'
[675]312        set(handles.num_RangeY_2,'TooltipString','num_RangeY_2: range of projection around each point')
[379]313%         set(handles.XObject,'TooltipString','XObject: set of x coordinates of the points')
314%         set(handles.YObject,'TooltipString','YObject: set of y coordinates of the points')
315%         set(handles.ZObject,'TooltipString','ZObject: set of z coordinates of the points')
[204]316    case {'line','polyline','polygon'}
[675]317        set(handles.num_RangeY_2,'TooltipString','num_RangeY_2: range of projection around the line')
[379]318         set(handles.Coord,'TooltipString','Coord: table of x,y, z coordinates defining the line')
319%         set(handles.YObject,'TooltipString','YObject: set of y coordinates defining the line')
320%         set(handles.ZObject,'TooltipString','ZObject: set of z coordinates defining the line')
[581]321        if isequal(ProjMode,'interp_lin')|| isequal(ProjMode,'interp_tps')
[379]322            set(handles.num_DX,'Visible','on')
323            set(handles.num_DX,'TooltipString','num_DX: mesh for the interpolated field along the line')
[204]324        end       
325    case {'rectangle','ellipse'}
[675]326        set(handles.num_RangeX_2,'TooltipString',['num_RangeX_2: half length of the ' ObjectStyle])
327        set(handles.num_RangeY_2,'TooltipString',['num_RangeY_2: half width of the ' ObjectStyle])
[204]328    case {'plane'} 
[379]329        set(handles.num_Angle_3,'Visible','on')
330        set(handles.num_RangeX_1,'Visible','on')
331        set(handles.num_RangeX_2,'Visible','on')
332        set(handles.num_RangeY_1,'Visible','on')
333        set(handles.num_RangeY_2,'Visible','on')
334        set(handles.num_RangeZ_2,'TooltipString','num_ZMax: range of projection normal to the plane')
[204]335        if test3D
[379]336            set(handles.num_Angle_2,'Visible','on')
337            set(handles.num_Angle_1,'Visible','on')
338            set(handles.num_RangeZ_2,'Visible','on')
[204]339        end
[581]340        if isequal(ProjMode,'interp_lin')|| isequal(ProjMode,'interp_tps')
[379]341            set(handles.num_DX,'Visible','on')
342            set(handles.num_DY,'Visible','on')
[204]343        else
[379]344            set(handles.num_DX,'Visible','off')
345            set(handles.num_DY,'Visible','off')
[204]346        end
[581]347        if  isequal(ProjMode,'interp_lin')
[379]348            set(handles.num_DZ,'Visible','on') 
[204]349        end
350     case {'volume'} 
[379]351        set(handles.num_RangeX_1,'Visible','on')
352        set(handles.num_RangeX_2,'Visible','on')
353        set(handles.num_RangeY_1,'Visible','on')
354        set(handles.num_RangeY_2,'Visible','on')
[204]355        set(handles.XObject,'TooltipString',['XObject:  x coordinate of the axis origin for the ' ObjectStyle])
356        set(handles.YObject,'TooltipString',['YObject:  y coordinate of the axis origin for the ' ObjectStyle])
[379]357        set(handles.num_Angle_1,'Visible','on')
358        set(handles.num_Angle_2,'Visible','on')
359        set(handles.num_Angle_3,'Visible','on')
360        set(handles.num_RangeZ_1,'Visible','on')
361        set(handles.num_RangeZ_2,'Visible','on')
[581]362        if isequal(ProjMode,'interp_lin')|| isequal(ProjMode,'interp_tps')
[379]363            set(handles.num_DX,'Visible','on')
364            set(handles.num_DY,'Visible','on')
365            set(handles.num_DZ,'Visible','on')
[204]366        else
[379]367            set(handles.num_DX,'Visible','off')
368            set(handles.num_DY,'Visible','off')
369            set(handles.num_DZ,'Visible','off')
[204]370        end
371end
[864]372% set default values read in the plot of uvmat to initiate the mesh
373if isequal(ProjMode,'interp_lin')|| isequal(ProjMode,'interp_tps')
374    if isempty(str2num(get(handles.num_DX,'String')))||isempty(str2num(get(handles.num_DY,'String')));     
375        huvmat=findobj('Tag','uvmat');%find the current uvmat interface handle
376        UvData=get(huvmat,'UserData');%Data associated to the current uvmat interface
377        Field=UvData.Field;
378        if  isfield(UvData.Field,'CoordMesh')&&~isempty(UvData.Field.CoordMesh)
379            set(handles.num_DX,'String',num2str(UvData.Field.CoordMesh))
380            set(handles.num_DY,'String',num2str(UvData.Field.CoordMesh))
381            set(handles.num_RangeX_1,'String',num2str(UvData.Field.XMin))
382            set(handles.num_RangeX_2,'String',num2str(UvData.Field.XMax))
383            set(handles.num_RangeY_1,'String',num2str(UvData.Field.YMin))
384            set(handles.num_RangeY_2,'String',num2str(UvData.Field.YMax))
385        end
386        if isempty(get(handles.CoordUnit,'String'))
387            set(handles.CoordUnit,'String',Field.CoordUnit)
388        end       
389    end
390end
391
[204]392%------------------------------------------------------------------------
393
394%------------------------------------------------------------------------
[379]395function num_Angle_1_Callback(hObject, eventdata, handles)
[204]396update_slider(hObject, eventdata,handles)
397%------------------------------------------------------------------------
398%------------------------------------------------------------------------
[379]399function num_Angle_2_Callback(hObject, eventdata, handles)
[204]400update_slider(hObject, eventdata,handles)
401%------------------------------------------------------------------------
402function update_slider(hObject, eventdata,handles)
403%rotation angles
[379]404PlaneAngle(1)=str2num(get(handles.num_Angle_1,'String'));%first  angle in degrees
405PlaneAngle(2)=str2num(get(handles.num_Angle_2,'String'));%second  angle in degrees
406PlaneAngle(3)=str2num(get(handles.num_Angle_3,'String'));%second  angle in degrees
[212]407om=norm(PlaneAngle);%norm of rotation angle in radians
408OmAxis=PlaneAngle/om; %unit vector marking the rotation axis
409cos_om=cos(pi*om/180);
410sin_om=sin(pi*om/180);
411coeff=OmAxis(3)*(1-cos_om);
412%components of the unity vector norm_plane normal to the projection plane
413norm_plane(1)=OmAxis(1)*coeff+OmAxis(2)*sin_om;
414norm_plane(2)=OmAxis(2)*coeff-OmAxis(1)*sin_om;
415norm_plane(3)=OmAxis(3)*coeff+cos_om;
[204]416huvmat=findobj('Tag','uvmat');%find the current uvmat interface handle
417UvData=get(huvmat,'UserData');%Data associated to the current uvmat interface
418if isfield(UvData,'X') & isfield(UvData,'Y') & isfield(UvData,'Z')
[212]419    Z=norm_plane(1)*(UvData.X)+norm_plane(2)*(UvData.Y)+norm_plane(3)*(UvData.Z);
[204]420    set(handles.z_slider,'Min',min(Z))
421    set(handles.z_slider,'Max',max(Z))
422    ZMax_Callback(hObject, eventdata, handles)
423end
424%------------------------------------------------------------------------
[379]425function num_DX_Callback(hObject, eventdata, handles)
[204]426%------------------------------------------------------------------------
427%------------------------------------------------------------------------
[379]428function num_DY_Callback(hObject, eventdata, handles)
[204]429%------------------------------------------------------------------------
430%------------------------------------------------------------------------
[379]431function num_DZ_Callback(hObject, eventdata, handles)
[204]432%------------------------------------------------------------------------
433
434
435%------------------------------------------------------------------------
[680]436% --- Executes on button press in REFRESH: refresh the current object , refresh the object and its projected field
437%------------------------------------------------------------------------
438function REFRESH_Callback(hObject, eventdata, handles)
[204]439
[681]440set(handles.REFRESH,'BackgroundColor',[1 1 0])% indicate activation of REFRESH
[508]441drawnow
442
[630]443%% update the object in the GUI series if relevant
444if strcmp(get(handles.set_object,'Name'),'edit_object_series')
445    hseries=findobj(allchild(0),'Tag','series');
446    if ~isempty(hseries)
447        SeriesData=get(hseries,'UserData');
448    SeriesData.ProjObject=read_GUI(handles.set_object);%read the parameters defining the object in the GUI set_object
449    set(hseries,'UserData',SeriesData);
450    end
[680]451    set(handles.REFRESH,'BackgroundColor',[1 0 0])
[630]452    return
453end
454
[410]455%% read the object parameters in the GUI set_object
[388]456ObjectData=read_GUI(handles.set_object);%read the parameters defining the object in the GUI set_object
[397]457if iscell(ObjectData.Coord)%check for empty line
458    ObjectData.Coord=[0 0 0];
459    hhset_object=guidata(handles.set_object);
460    set(hhset_object.Coord,'Data',ObjectData.Coord)
461end
[388]462checknan=isnan(sum(ObjectData.Coord,2));%check for NaN lines
463if ~isempty(checknan)
464    ObjectData.Coord(checknan,:)=[];%remove the NaN lines
465end
[402]466ObjectName=ObjectData.Name;%name of the current object defined in set_object
[309]467if isempty(ObjectName)
[432]468     ObjectName=ObjectData.Type;% name the object by the object type type by default
[204]469end
[413]470
[432]471%% read the current object selection in the GUI uvmat
[413]472huvmat=findobj('tag','uvmat');%find the current uvmat GUI handle
473UvData=get(huvmat,'UserData');%Data associated to the GUI uvmat
474hhuvmat=guidata(huvmat);%handles of the objects children of the  GUI uvmat
[432]475ListObject=get(hhuvmat.ListObject,'String');% list of objects displayed in uvmat
[625]476
[650]477if isequal(get(hhuvmat.CheckEditObject,'Value'),0) %we append a new object
[625]478    ListObject=[ListObject;{''}];
479    IndexObj=length(ListObject);
480    set(hhuvmat.ListObject,'String',ListObject)
481    set(hhuvmat.ListObject,'Value',IndexObj)
482    UvData.ProjObject{IndexObj}=[]; %create a new empty object
483    UvData.ProjObject{IndexObj}.DisplayHandle.uvmat=hhuvmat.PlotAxes; % axes for plot_object
484    UvData.ProjObject{IndexObj}.DisplayHandle.view_field=[]; %no plot handle before plot_field operation
[630]485else   
[625]486    IndexObj=get(hhuvmat.ListObject,'Value');% index of the selected object for display in uvmat
487end
488
[630]489%set or modify(edit mode) the name of the currently selected object
[445]490detectname=1;
491ObjectNameNew=ObjectName;
492vers=0;% index of the name
493ListOther=ListObject;
494ListOther(IndexObj)=[];
495while ~isempty(detectname)
496    detectname=find(strcmp(ObjectNameNew,ListOther),1);%test the existence of the proposed name in the list
497    if detectname% if the object name already exists
498        indstr=regexp(ObjectNameNew,'\D');%indices of non number characters
499        if indstr(end)<length(ObjectNameNew) %object name ends by a number
500            vers=str2double(ObjectNameNew(indstr(end)+1:end))+1;
501            ObjectNameNew=[ObjectNameNew(1:indstr(end)) num2str(vers)];
502        else
503            vers=vers+1;
504            ObjectNameNew=[ObjectNameNew(1:indstr(end)) '_' num2str(vers)];
[204]505        end
506    end
[445]507end
508ObjectName=ObjectNameNew;
509set(handles.Name,'String',ObjectName)% display the default name in set_object
510ListObject{IndexObj}=ObjectName;
511set(hhuvmat.ListObject,'String',ListObject);%complement the object list
[625]512set(hhuvmat.ListObject_1,'String',ListObject);%complement the object list
[650]513set(hhuvmat.CheckViewObject,'Value',1)% indicate that the currently selected objected is viewed on set_object
[622]514check_handle=isfield(UvData.ProjObject{IndexObj},'DisplayHandle') && isfield(UvData.ProjObject{IndexObj}.DisplayHandle,'uvmat')...
515    && ~isempty(UvData.ProjObject{IndexObj}.DisplayHandle.uvmat) && ishandle(UvData.ProjObject{IndexObj}.DisplayHandle.uvmat);
[445]516if check_handle
[622]517    obj_handle=UvData.ProjObject{IndexObj}.DisplayHandle.uvmat;
[445]518end
[622]519UvData.ProjObject{IndexObj}=ObjectData;%record the current object properties in uvmat
[445]520if check_handle
[622]521    UvData.ProjObject{IndexObj}.DisplayHandle.uvmat=obj_handle; %preserve the object plot handle if valid
[445]522else
[622]523    UvData.ProjObject{IndexObj}.DisplayHandle.uvmat=hhuvmat.PlotAxes; %axes taken as object display handle by defualt
[204]524end
[309]525
[680]526%% refresh the field projected on the object
[432]527hview_field=[];%default
[445]528IndexObj_1=get(hhuvmat.ListObject_1,'Value');
[388]529if strcmp(ObjectData.ProjMode,'mask_inside')||strcmp(ObjectData.ProjMode,'mask_outside')||strcmp(ObjectData.ProjMode,'none')
530    PlotType='text';
531else
[581]532    % create tps coeff if needed for ProjMode 'interp_tps'
[586]533    if strcmp(ObjectData.ProjMode,'interp_tps')&&~isfield(UvData.Field,'Coord_tps')
534        %UvData.Field=calc_tps(UvData.Field,1);
535        [UvData.Field,errormsg]=tps_coeff_field(UvData.Field,1);
536        if ~isempty(errormsg)
537            msgbox_uvmat('ERROR', ['set_object/tps_coeff_field/' errormsg])
[680]538            set(handles.REFRESH,'enable','on')
[586]539            return
540        end
[508]541    end
[388]542    [ProjData,errormsg]= proj_field(UvData.Field,ObjectData);%project the current field of uvmat on ObjectData
543    if ~isempty(errormsg)
[586]544        msgbox_uvmat('ERROR', ['set_object/proj_field/' errormsg])
[680]545        set(handles.REFRESH,'enable','on')
[388]546        return
[586]547    end
[445]548    if isequal(IndexObj_1,IndexObj) % if  the projection is in uvmat
[586]549        PlotType=plot_field(ProjData,hhuvmat.PlotAxes,read_GUI(get(hhuvmat.PlotAxes,'parent')));%update the current uvmat plot
[491]550    else  % if the projection is in view_field
[413]551        hview_field=findobj(allchild(0),'tag','view_field');
552        if isempty(hview_field)
[429]553            hview_field=view_field(ProjData); %open the view_field GUI for plot
[413]554        else
555            hhview_field=guidata(hview_field);
[511]556            [PlotType,PlotParam]=plot_field(ProjData,hhview_field.PlotAxes,read_GUI(hview_field));%update an existing  plot in view_field
[591]557            errormsg=fill_GUI(PlotParam,hview_field);
[569]558            if ~isempty(errormsg)
559                msgbox_uvmat('ERROR',errormsg)
560                return
561            end
[586]562            %     write_plot_param(hhview_field,PlotParam); %update the display of plotting parameters for the current object
[413]563        end
[429]564        haxes=findobj(hview_field,'tag','axes3');
565        Data=get(hview_field,'UserData');
566        if strcmp(get(haxes,'Visible'),'off')%sempty(PlotParam.Coordinates)% case of no plot display (pure text table)
567            h_TableDisplay=findobj(hview_field,'tag','TableDisplay');
568            pos_table=get(h_TableDisplay,'Position');
569            pos=get(hview_field,'Position');
570            set(hview_field,'Position',[pos(1)+pos(3)-pos_table(3) pos(2)+pos(4)-pos_table(4) pos_table(3) pos_table(4)])
571            drawnow
572            set(hview_field,'UserData',Data);% restore the previously stored GUI position after GUI resizing
573        else
574            set(hview_field,'Position',Data.GUISize)
575        end
[388]576    end
[227]577end
[429]578
[680]579%% update the object refresh
[622]580hobject=UvData.ProjObject{IndexObj}.DisplayHandle.uvmat;
[445]581% if we are editing the object used for projection in uvmat
582if isequal(IndexObj_1,IndexObj)
583    %update the representation of the current object for projection field represented in view_field
[622]584    for iobj=1:numel(UvData.ProjObject)
585        UvData.ProjObject{iobj}.DisplayHandle.uvmat=...
586            plot_object(UvData.ProjObject{iobj},UvData.ProjObject{IndexObj_1},UvData.ProjObject{iobj}.DisplayHandle.uvmat,'b');
[432]587    end
[445]588else %  we are editing the object used for projection field represented in view_field
589    %update the representation of the current object in uvmat
[622]590    UvData.ProjObject{IndexObj}.DisplayHandle.uvmat=...
591             plot_object(UvData.ProjObject{IndexObj},UvData.ProjObject{IndexObj_1},UvData.ProjObject{IndexObj}.DisplayHandle.uvmat,'m');
[680]592    %indicate the object index in the user data of the object refresh (needed for further mouse editing)
[622]593    ObjectInfo=get(UvData.ProjObject{IndexObj}.DisplayHandle.uvmat,'UserData');
[445]594    ObjectInfo.IndexObj=IndexObj;
[622]595    set(UvData.ProjObject{IndexObj}.DisplayHandle.uvmat,'UserData',ObjectInfo)
[445]596    % update the representation of all objects in view_field
[622]597    for iobj=1:numel(UvData.ProjObject)
[729]598        if isfield(UvData.ProjObject{iobj},'DisplayHandle') && isfield(UvData.ProjObject{iobj}.DisplayHandle,'view_field')
[622]599            UvData.ProjObject{iobj}.DisplayHandle.view_field=...
600                plot_object(UvData.ProjObject{iobj},UvData.ProjObject{iobj},UvData.ProjObject{iobj}.DisplayHandle.view_field,'b');
[445]601        end
[432]602    end
603end
604set(huvmat,'UserData',UvData)
605
[204]606%% update the GUI uvmat
[650]607set(hhuvmat.CheckEditObject,'Value',1) % set uvmat to object edit mode to allow further object update
608set(hhuvmat.CheckViewField,'Value',1)
[681]609
[680]610set(handles.REFRESH,'BackgroundColor',[1 0 0])
[681]611set(handles.num_RangeY_2,'BackgroundColor',[1 1 1])
612
[772]613% --- Executes on button press in DisplayCoord.
614function DisplayCoord_Callback(hObject, eventdata, handles)
[738]615global Coord
616Coord=get(handles.Coord,'Data');
617evalin('base','global Coord')%make Coord global in the workspace
618display('object coordinates:')
619evalin('base','Coord') %display Coord in the workspace
620commandwindow; %brings the Matlab command window to the front
621
[204]622%----------------------------------------------------
[379]623function num_RangeY_1_Callback(hObject, eventdata, handles)
[204]624%------------------------------------------------------------------------
625
[379]626function num_RangeZ_1_Callback(hObject, eventdata, handles)
[204]627%------------------------------------------------------------------------
628
[379]629function num_RangeZ_2_Callback(hObject, eventdata, handles)
630DZ=str2num(get(handles.num_RangeZ_2,'String'));
[204]631ZMin=get(handles.z_slider,'Min');
632ZMax=get(handles.z_slider,'Max');
633if ~isequal(ZMax-ZMin,0)
634    rel_step(1)=DZ/(ZMax-ZMin);
635    rel_step(2)=0.2;
636    set(handles.z_slider,'SliderStep',rel_step)
637end
638%------------------------------------------------------------------------
[379]639function num_RangeY_2_Callback(hObject, eventdata, handles)
[204]640%------------------------------------------------------------------------
641
[379]642function num_RangeX_1_Callback(hObject, eventdata, handles)
[204]643%------------------------------------------------------------------------
644
[379]645function num_RangeX_2_Callback(hObject, eventdata, handles)
[681]646
[204]647%------------------------------------------------------------------------
648%------------------------------------------------------------------------
649function SAVE_Callback(hObject, eventdata, handles)
650% ------------------------------------------------------
[379]651Object=read_GUI(handles.set_object);
[204]652huvmat=findobj('Tag','uvmat');
653if isempty(huvmat)
654    huvmat=findobj(allchild(0),'Name','series');
655end
656hchild=get(huvmat,'Children');
657hrootpath=findobj(hchild,'Tag','RootPath');
658if isempty(hrootpath)
659    RootPath='';
660else
661    RootPath=get(hrootpath,'String');
662    if iscell(RootPath)
663        RootPath=RootPath{1};
664    end
665end
666title={'object name'};
[667]667dir_save=uigetfile_uvmat('select the folder for the new xml object file:',RootPath,'uigetdir');
668if ~isempty(dir_save)
669    ObjectName=get(handles.Name,'String');
670    if ~isempty(ObjectName)&&~strcmp(ObjectName,'')
[675]671        def=[ObjectName '.xml'];
[667]672    else
[675]673        def=[Object.Style '.xml'];
[667]674    end
[705]675    displ_txt={['save object in' dir_save]; 'with file name (.xml):'};%default display
[667]676    menu=get(handles.ProjMode,'String');
677    value=get(handles.ProjMode,'Value');
678    ProjMode=menu{value};
679    if strcmp(ProjMode,'mask_inside')||strcmp(ProjMode,'mask_outside')
[745]680        displ_txt=[displ_txt; '(note: to create a mask image, use ''Tools/make mask'' on the upper bar menu of uvmat)'];
[667]681    end
[675]682    answer=msgbox_uvmat('INPUT_TXT',displ_txt,def);
[705]683    if ischar(answer)
684        FullName=fullfile(dir_save,answer);
[667]685        t=struct2xml(Object);
686        t=set(t,1,'name','ProjObject');
[675]687        save(t,FullName)
688        msgbox_uvmat('CONFIRMATION',[FullName  ' saved'])
[667]689    end
[204]690end
[227]691
[204]692%------------------------------------------------------------------------
693% --- Executes on slider movement.
694function z_slider_Callback(hObject, eventdata, handles)
695%---------------------------------------------------------
696Z_value=get(handles.z_slider,'Value');
697%rotation angles
[227]698PlaneAngle=[0 0 0];
699norm_plane=[0 0 1];
700cos_om=1;
701sin_om=0;
702
[379]703PlaneAngle(1)=str2double(get(handles.num_Angle_1,'String'));%first  angle in degrees
704PlaneAngle(2)=str2double(get(handles.num_Angle_2,'String'));%second  angle in degrees
705PlaneAngle(3)=str2double(get(handles.num_Angle_3,'String'));%second  angle in degrees
[227]706PlaneAngle=(pi/180)*PlaneAngle;
[212]707om=norm(PlaneAngle);%norm of rotation angle in radians
[227]708if isequal(om,0)
709    norm_plane=[0 0 1];
710else
711    OmAxis=PlaneAngle/om; %unit vector marking the rotation axis
712    cos_om=cos(om);
713    sin_om=sin(om);
714    coeff=OmAxis(3)*(1-cos_om);
715    %components of the unity vector norm_plane normal to the projection plane
716    norm_plane(1)=OmAxis(1)*coeff+OmAxis(2)*sin_om;
717    norm_plane(2)=OmAxis(2)*coeff-OmAxis(1)*sin_om;
718    norm_plane(3)=OmAxis(3)*coeff+cos_om;
719end
[747]720Coord=get(handles.Coord,'Data');
721Coord(3)=Z_value;
722set(handles.Coord,'Data',Coord)
[204]723
[388]724% update graph
[680]725REFRESH_Callback(hObject, eventdata, handles)
[227]726
[204]727%------------------------------------------------------------------------
728% --- Executes on button press in HELP.
729function HELP_Callback(hObject, eventdata, handles)
[227]730%------------------------------------------------------------------------
[773]731web('http://servforge.legi.grenoble-inp.fr/projects/soft-uvmat/wiki/UvmatHelp#ProjObject')
[204]732
[850]733%------------------------------------------------------------------------
[674]734% --- Executes when selected cell(s) is changed in ListCoord.
[850]735%------------------------------------------------------------------------
[674]736function Coord_CellSelectionCallback(hObject, eventdata, handles)
737if ~isempty(eventdata.Indices)
738    iline=eventdata.Indices(1);% selected line number
[850]739    set(handles.CoordLine,'String',num2str(iline))
[674]740end
[379]741
[681]742
743function num_Angle_3_Callback(hObject, eventdata, handles)
744
[674]745%------------------------------------------------------------------------
[681]746% --- Executes on key press with selection of a uicontrol
[674]747%------------------------------------------------------------------------
[681]748function KeyPressFcn(hObject, eventdata, handles)
749set(handles.REFRESH,'BackgroundColor',[1 0 1])% se REFRESH to magenta color, indicates that refresh needs to be done
[684]750if strcmp(get(hObject,'Tag'),'num_Angle_3')
751    set(handles.num_RangeX_1,'String','')
752    set(handles.num_RangeX_2,'String','')
753        set(handles.num_RangeY_1,'String','')
754    set(handles.num_RangeY_2,'String','')
755end
[674]756
[681]757%------------------------------------------------------------------------
758% --- Executes on key press with focus on Coord and none of its controls.
759%------------------------------------------------------------------------
760function Coord_KeyPressFcn(hObject, eventdata, handles)
761
762set(handles.REFRESH,'BackgroundColor',[1 0 1])
[674]763xx=double(get(handles.set_object,'CurrentCharacter')); %get the keyboard character
[680]764if ismember(xx,[127 31])% delete, or downward
[674]765    Coord=get(handles.Coord,'Data');
[850]766    iline=str2double(get(handles.CoordLine,'String'));
[674]767            if isequal(xx, 31)
768                if isequal(iline,size(Coord,1))% arrow downward
769                Coord=[Coord;zeros(1,size(Coord,2))];
770                end
771            else
772    Coord(iline,:)=[];% suppress the current line
773            end
774    set(handles.Coord,'Data',Coord);
775end
[850]776
777%------------------------------------------------------------------------
778% --- Executes on button press in clear_line.
779%------------------------------------------------------------------------
780function clear_line_Callback(hObject, eventdata, handles)
781
782Coord=get(handles.Coord,'Data');
783iline=str2double(get(handles.CoordLine,'String'));
784if isempty(iline)
785    msgbox_uvmat('WARNING','no line suppressed, select a line in the table')
786else
787    Coord(iline,:)=[];
788    set(handles.REFRESH,'BackgroundColor',[1 0 1])
789    set(handles.Coord,'Data',Coord);
790    set(handles.CoordLine,'String','')
791    REFRESH_Callback(hObject,eventdata,handles)
792end
Note: See TracBrowser for help on using the repository browser.