source: trunk/src/set_object.m @ 1072

Last change on this file since 1072 was 1072, checked in by sommeria, 3 years ago

LIF updated and bug corrections

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