source: trunk/src/mouse_up.m @ 660

Last change on this file since 660 was 660, checked in by sommeria, 7 years ago

a few bugs corrected in geometry_calib

File size: 20.1 KB
Line 
1%'mouse_up': function  activated when the mouse button is released
2%------------------------------------------------------------------------
3% function mouse_up(hObject,eventdata,handles)
4% activated by the command:
5% set(hObject,'WindowButtonUpFcn',{'mouse_up'}),
6% where hObject is the handle of the figure
7
8%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
9%  Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.
10%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
11%     This file is part of the toolbox UVMAT.
12%
13%     UVMAT is free software; you can redistribute it and/or modify
14%     it under the terms of the GNU General Public License as published by
15%     the Free Software Foundation; either version 2 of the License, or
16%     (at your option) any later version.
17%
18%     UVMAT is distributed in the hope that it will be useful,
19%     but WITHOUT ANY WARRANTY; without even the implied warranty of
20%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
21%     GNU General Public License (file UVMAT/COPYING.txt) for more details.
22%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
23
24function mouse_up(hObject,eventdata,handles)
25
26test_ruler=0;%default
27AxeData=get(gca,'UserData');
28if isfield(AxeData,'ParentAxes')% case of a zoom plot as current axis
29    hcurrentaxes=AxeData.ParentAxes;
30    AxeData=get(hcurrentaxes,'UserData');
31    hcurrentfig=get(hcurrentaxes,'parent');%handles of the GUI parent of the zoom plot
32    testsubplot=1;% mouse selection is on a zoom subplot
33else
34    hcurrentfig=hObject;
35    hcurrentaxes=gca; %store the current axes handle
36    testsubplot=0;
37end
38CurrentOrigin=[];
39if isfield(AxeData,'CurrentOrigin')
40    CurrentOrigin=AxeData.CurrentOrigin;
41end
42FigTag=get(hcurrentfig,'tag');
43hhcurrentfig=guidata(hcurrentfig);%the current figure is a GUI (uvmat or view_field)
44CheckZoom=get(hhcurrentfig.CheckZoom,'Value');
45CheckZoomFig=get(hhcurrentfig.CheckZoomFig,'Value');%exclusive to CheckZoom
46huvmat=findobj(allchild(0),'tag','uvmat');%find the uvmat interface handle
47if ~isempty(huvmat)
48    hhuvmat=guidata(huvmat);
49    UvData=get(huvmat,'UserData');
50   test_ruler=~CheckZoom && isequal(get(hhuvmat.MenuRuler,'checked'),'on');%test for ruler  action, second priority
51end
52test_drawing=0;%default, =1 to allow drawing by further mouse action
53xy=get(gca,'CurrentPoint');%xy(1,1),xy(1,2): current x,y positions in axes coordinates
54
55
56%% proceed with the creation or editing (translation/deformation) of an object
57if ~isempty(huvmat) && isfield(AxeData,'Drawing') && ~isequal(AxeData.Drawing,'off') && isfield(AxeData,'CurrentObject')...
58        && ~isempty(AxeData.CurrentObject) && ishandle(AxeData.CurrentObject)
59    %     xy=get(currentaxes,'CurrentPoint');%xy(1,1),xy(1,2): current x,y positions in axes coordinates
60    PlotData=get(AxeData.CurrentObject,'UserData');%get data attached to the current projection object
61    IndexObj=PlotData.IndexObj;
62    ObjectData=UvData.ProjObject{IndexObj};
63    check_multiple=0;
64   
65    h_set_object=findobj(allchild(0),'Tag','set_object');
66    hh_set_object=guidata(h_set_object);
67    ObjectData.Coord=get(hh_set_object.Coord,'Data');
68   
69    % ending translation
70    if isequal(AxeData.Drawing,'translate')
71        XYData=AxeData.CurrentOrigin;
72        DX=xy(1,1)-XYData(1);%translation from initial position
73        DY=xy(1,2)-XYData(2);
74        ObjectData.Coord(:,1)=ObjectData.Coord(:,1)+DX;
75        ObjectData.Coord(:,2)=ObjectData.Coord(:,2)+DY;
76       
77        %ending object deformation
78    elseif isequal(AxeData.Drawing,'deform')
79        ind_move=AxeData.CurrentIndex;
80        ObjectData.Coord(ind_move,1)=xy(1,1);
81        ObjectData.Coord(ind_move,2)=xy(1,2);
82       
83        %creating object
84    else
85        switch ObjectData.Type
86            case {'line'}
87                if size(ObjectData.Coord,1)==1 % this is the mouse up for the first point, continue until next click
88                    check_multiple=1;
89                else
90                    %ObjectData.Coord=[ObjectData.Coord ;CurrentOrigin];% append the second point of the line (the last pointed position during mouse down)
91                end
92            case {'rectangle','ellipse','volume'}
93                %                  if size(ObjectData.Coord,1)==1 % this is the mouse up for the first point, continue until next click
94                %                     check_multiple=1;
95                %                  else
96                ObjectData.Coord=(CurrentOrigin+xy(1,1:2))/2;% keep only the first point coordinate
97                ObjectData.RangeX=abs(ObjectData.Coord(1,1)-xy(1,1));%rectangle width
98                ObjectData.RangeY=abs(ObjectData.Coord(1,2)-xy(1,2));%rectangle height
99                if isequal(ObjectData.RangeX,0)||isequal(ObjectData.RangeY,0)
100                    check_multiple=1;% pass to next mous up if width of height=0
101                end
102                %                 ObjectData.Coord(1,1)=(xy(1,1)+XYData(1))/2;%origin rectangle, x coordinate
103                %                 ObjectData.Coord(1,2)=(xy(1,2)+XYData(2))/2;
104                %                 ObjectData.RangeX=abs(xy(1,1)-XYData(1))/2;%rectangle width
105                %                 ObjectData.RangeY=abs(xy(1,2)-XYData(2))/2;%rectangle height
106                %                  end
107            case 'plane' %case of 'plane'
108                DX=(xy(1,1)-ObjectData.Coord(1,1));
109                DY=(xy(1,2)-ObjectData.Coord(1,2));
110                ObjectData.Phi=(angle(DX+i*DY))*180/pi;%rectangle width
111                if isfield(ObjectData,'RangeX')
112                    XMax=sqrt(DX*DX+DY*DY);
113                    if XMax>max(ObjectData.RangeX)
114                        ObjectData.RangeX=[min(ObjectData.RangeX) XMax];
115                    end
116                end
117            otherwise
118                check_multiple=1;
119        end
120    end
121   
122    %show object coordinates in the GUI set_object
123    %     h_set_object=findobj(allchild(0),'Tag','set_object');
124    %     hh_set_object=guidata(h_set_object);
125    %     set(hh_set_object.Coord,'Data',ObjectData.Coord);
126    if strcmp(ObjectData.Type,'rectangle')||strcmp(ObjectData.Type,'ellipse')
127        set(hh_set_object.num_RangeX_2,'String',num2str(ObjectData.RangeX,4));
128        set(hh_set_object.num_RangeY_2,'String',num2str(ObjectData.RangeY,4));
129    end
130   
131    %% stop drawing and plot the projected field if the object manipulation is finished
132    if check_multiple==0  || isequal(get(hcurrentfig,'SelectionType'),'alt')
133        AxeData.CurrentOrigin=[]; %suppress the current origin
134        hobject=UvData.ProjObject{IndexObj}.DisplayHandle.(FigTag);
135        if ~isempty(hObject)
136            ProjObject=UvData.ProjObject{get(hhuvmat.ListObject_1,'Value')};
137            AxeData.CurrentObject=plot_object(ObjectData,ProjObject,hobject,'m');%draw the object and its handle becomes AxeData.CurrentObject
138        end
139        %%
140        if  ~isempty(ObjectData)
141            % plot the field projected on the object
142            ProjData= proj_field(UvData.Field,ObjectData);%project the current interface field on ObjectData
143            if ~isempty(ProjData)
144                if strcmp(FigTag,'uvmat')% uvmat plot selected, projection plot seen in view_field
145                    hview_field=findobj(allchild(0),'tag','view_field');
146                    if isempty(hview_field)
147                        hview_field=view_field(ProjData); %open the view_field GUI for plot
148                    else
149                        hhview_field=guidata(hview_field);
150                        [PlotType,PlotParam]=plot_field(ProjData,hhview_field.PlotAxes,read_GUI(hview_field));%update an existing  plot in view_field
151                        errormsg=fill_GUI(PlotParam,hview_field);
152                    end
153                    ViewFieldData=get(hview_field,'UserData');
154                    haxes=findobj(hview_field,'tag','axes3');
155                    if strcmp(get(haxes,'Visible'),'off')%sempty(PlotParam.Coordinates)% case of no plot display (pure text table)
156                        h_TableDisplay=findobj(hview_field,'tag','TableDisplay');
157                        pos_table=get(h_TableDisplay,'Position');
158                        pos=get(hview_field,'Position');
159                        set(hview_field,'Position',[pos(1)+pos(3)-pos_table(3) pos(2)+pos(4)-pos_table(4) pos_table(3) pos_table(4)])
160                        drawnow
161                        set(hview_field,'UserData',ViewFieldData);% restore the previously stored GUI position after GUI resizing
162                    else
163                        set(hview_field,'Position',ViewFieldData.GUISize)
164                    end
165                else
166                    UvData.PlotAxes=ProjData;
167                    [PlotType,PlotParam]=plot_field(ProjData,hhuvmat.PlotAxes,read_GUI(huvmat));%update an existing field plot
168                    errormsg=fill_GUI(PlotParam,huvmat);
169                end
170            end
171            set(hhuvmat.CheckViewField,'Value',1);%
172            set(hhuvmat.CheckEditObject,'Value',1);%           
173            set(hhuvmat.CheckEditObject,'Enable','on');%
174            set(get(h_set_object,'children'),'Enable','on')
175        end
176        UvData.ProjObject{IndexObj}=ObjectData;
177        if isfield(UvData.ProjObject{IndexObj},'CreateMode')
178            UvData.ProjObject{IndexObj}=rmfield(UvData.ProjObject{IndexObj},'CreateMode');%remove createMode to mark the object as finished
179        end
180    else
181        test_drawing=1;%allow continuation of drawing object
182        AxeData.CurrentOrigin=[xy(1,1) xy(1,2)]; %the current point becomes the next current origin
183    end
184    %     UvData.ProjObject{IndexObj}=ObjectData;
185    hother=findobj('Tag','deformpoint');%find all the deformpoints
186    set(hother,'Color','b');%reset all the deformpoints in 'blue'
187end
188
189%% creation or update of a  zoom sub-plot
190if CheckZoomFig && isequal(get(hcurrentfig,'SelectionType'),'normal')&&...%if left button has been pressed
191     ~isempty(CurrentOrigin) && ~isequal(CurrentOrigin(1),xy(1,1)) && ~isequal(CurrentOrigin(2),xy(1,2))%if mouse moved in x and y since presed down
192    hparentfig=hcurrentfig;
193    %open or update a new zoom figure if a rectangle has been drawn
194    if ishandle(hcurrentaxes);
195        if isfield(AxeData,'CurrentRectZoom') && ~isempty(AxeData.CurrentRectZoom) && ishandle(AxeData.CurrentRectZoom)
196            %PosRect=get(AxeData.CurrentRectZoom,'Position');
197            if isfield(AxeData,'CurrentVec') && ~isempty(AxeData.CurrentVec) && ishandle(AxeData.CurrentVec)
198                delete(AxeData.CurrentVec)
199            end
200            if ~testsubplot% if we are not already on a zoom plot
201                hfig2=findobj(allchild(0),'Tag','zoom_fig');
202                if isempty(hfig2)% create zoom sub plot if absent
203                    hfig2=figure('name',['zoom_' FigTag],'tag','zoom_fig');%create new figure (unit='pixels' by default)
204                    set(0,'Unit','pixels')
205                    FigPos=get(hfig2,'Position');%get the standard width and height of the fig
206                    ScreenSize=get(0,'ScreenSize');% get the size of the screen, to put the fig on the upper right
207                    Left=ScreenSize(3)- FigPos(3)-40; %right edge close to the right, with margin=40
208                    Bottom=ScreenSize(4)-FigPos(4)-40; %put fig at top right
209                    FigPos(1:2)=[Left Bottom];
210                    set(hfig2,'Position',FigPos);% put the zoom fig close to the upper right of the screen
211                    map=colormap(hcurrentaxes);
212                    colormap(map);%transmit the current colormap to the zoom fig
213                    set(hfig2,'KeyPressFcn',{@keyboard_callback,handles})%set keyboard action function
214                    set(hfig2,'WindowButtonMotionFcn',{@mouse_motion,handles})%set mouse action function
215                    set(hfig2,'WindowButtonDownFcn',{@mouse_down})%set mouse click action function
216                    set(hfig2,'WindowButtonUpFcn',{@mouse_up,handles})
217                else
218                    zoom_axes=findobj(hfig2,'Type','axes');%delete existing axes
219                    axes(zoom_axes);%make the zoom axes apparent
220                    delete(zoom_axes)
221                end
222                set(hfig2,'DeleteFcn',{@close_fig,AxeData.CurrentRectZoom})
223                set(hfig2,'UserData',AxeData.CurrentRectZoom)% record the parent object (zoom rectangle) in the new fig   
224                AxeData.ZoomAxes=copyobj(hcurrentaxes,hfig2); %copy the current graph axes to the zoom figure
225                hrect_zoom=findobj(AxeData.ZoomAxes,'Tag','rect_zoom');%find and delete the copy of the rect_zoom rectangle
226                delete(hrect_zoom)
227                hcol=findobj(hparentfig,'Tag','Colorbar'); %look for colorbar axes
228                if ~isempty(hcol)
229                    hcol_new=colorbar;
230                    YTick=get(hcol,'YTick');
231                    YTicklabel=get(hcol,'Yticklabel');
232                    colbarlim=get(hcol,'YLim');
233                    newcolbarlim=get(hcol_new,'YLim');
234                    scale_bar=(newcolbarlim(2)-newcolbarlim(1))/(colbarlim(2)-colbarlim(1));
235                    YTick_rescaled=newcolbarlim(1)+scale_bar*(YTick-colbarlim(1));
236                    set(hcol_new,'YTick',YTick_rescaled);
237                    set(hcol_new,'Yticklabel',YTicklabel);
238                end
239                ChildAxeData=get(AxeData.ZoomAxes,'UserData');
240            end 
241            ChildAxeData.CurrentOrigin=[];% forget the current origin
242            ChildAxeData.CurrentRectZoom=[]; % no rect zoom in the new window
243            ChildAxeData.Drawing='off';
244            ChildAxeData.ParentAxes=hcurrentaxes;
245            ChildAxeData.ParentRect=AxeData.CurrentRectZoom;%set the rectangle drawing as a 'parent' associated to the new axe
246            %PosRect=CurrentOrigin;
247           % xy=get(hcurrentaxes,'CurrentPoint');%xy(1,1),xy(1,2): current x,y positions in axes coordinates
248            if xy(1,1)>CurrentOrigin(1)
249            set(AxeData.ZoomAxes,'Xlim',[CurrentOrigin(1) xy(1,1)])
250            else
251                set(AxeData.ZoomAxes,'Xlim',[xy(1,1) CurrentOrigin(1)])
252            end
253            if xy(1,2)>CurrentOrigin(2)
254            set(AxeData.ZoomAxes,'Ylim',[CurrentOrigin(2) xy(1,2)])
255            else
256                set(AxeData.ZoomAxes,'Ylim',[xy(1,2) CurrentOrigin(2)])
257            end
258            set(AxeData.ZoomAxes,'UserData',ChildAxeData);%update the AxeData of the new axes
259        end
260    end
261end
262
263%% zoom in or out by a factor 2 if no new figure is created
264if CheckZoom
265    if testsubplot
266        haxes=gca;% zoom on a zoom sub-plot
267    else
268        haxes=hcurrentaxes;% zoom on the main plot
269    end
270   % xy=get(haxes,'CurrentPoint');%xy(1,1),xy(1,2): current x,y positions in axes coordinates
271    xlim=get(haxes,'XLim');
272    ylim=get(haxes,'YLim');
273    % if left mouse button has been pressed, zoom in by a factor of 2
274    if  isequal(get(gcf,'SelectionType'),'normal');%if left button has been pressed, zoom in by a factor of 2
275        PlotBoxAspectRatio=get(haxes,'PlotBoxAspectRatio');
276        yoverx=PlotBoxAspectRatio(2)/PlotBoxAspectRatio(1);
277        if yoverx <2
278            xlim(1)=0.5*xy(1,1)+0.5*xlim(1);
279            xlim(2)=0.5*xy(1,1)+0.5*xlim(2);%double the field whith the middle at the selected points
280            set(haxes,'XLim',xlim)
281        end
282        if yoverx >0.5
283            ylim(2)=0.5*xy(1,2)+0.5*ylim(2);
284            ylim(1)=0.5*xy(1,2)+0.5*ylim(1);
285            set(haxes,'YLim',ylim)
286        end
287       
288        % if right mouse button has been pressed, zoom out by a factor of 2
289    else
290        xlim(1)=2*xlim(1)-xy(1,1);% reverse of the zoom on action
291        xlim(2)=2*xlim(2)-xy(1,1);
292        ylim(1)=2*ylim(1)-xy(1,2);
293        ylim(2)=2*ylim(2)-xy(1,2);
294        % adjust the zoom out to the available field
295        if ~testsubplot && isfield(AxeData,'RangeX')&& isfield(AxeData,'RangeY')
296            xlim(1)=max(AxeData.RangeX(1),xlim(1));
297            xlim(2)=min(AxeData.RangeX(2),xlim(2));
298            ylim(1)=max(AxeData.RangeY(1),ylim(1));
299            ylim(2)=min(AxeData.RangeY(2),ylim(2));
300            if ylim(1)>=ylim(2)|| xlim(1)>=xlim(2)
301                xlim=AxeData.RangeX;
302                ylim=AxeData.RangeY;
303            end
304            % desactivate the zoom if the full field is visible within the axes
305            if isequal(xlim,AxeData.RangeX) && isequal(ylim,AxeData.RangeY)
306                set(hhuvmat.CheckZoom,'Value',0)
307               % set(hhuvmat.CheckZoom,'BackgroundColor',[0.7 0.7 0.7])
308                set(hhuvmat.CheckFixLimits,'Value',0)
309              %  set(hhuvmat.CheckFixLimits,'BackgroundColor',[0.7 0.7 0.7])
310            end
311        end
312        set(haxes,'XLim',xlim)
313        set(haxes,'YLim',ylim)
314        %test whther zoom out is operating (to inactivate AxedAta
315        if ~isfield(AxeData,'CurrentXLim')|| ~isequal(xlim,AxeData.CurrentXLim)
316            AxeData.CurrentXLim=xlim;%
317        end
318    end
319    %if isfield(AxeData,'LimEditBox')&& AxeData.LimEditBox% update display of the GUI containing the axis (uvmat or view_field)
320    if testsubplot
321        set(AxeData.CurrentRectZoom,'Position',[xlim(1) ylim(1) xlim(2)-xlim(1) ylim(2)-ylim(1)])
322    else
323        set(hhcurrentfig.num_MinX,'String',num2str(xlim(1)))
324        set(hhcurrentfig.num_MaxX,'String',num2str(xlim(2)))
325        set(hhcurrentfig.num_MinY,'String',num2str(ylim(1)))
326        set(hhcurrentfig.num_MaxY,'String',num2str(ylim(2)))
327    end
328end
329
330%% editing calibration point
331if ~CheckZoom && isfield(AxeData,'Drawing') && isequal(AxeData.Drawing,'calibration')
332    h_geometry_calib=findobj(allchild(0),'tag','geometry_calib'); %find the geomterty_calib GUI
333    if ~isempty(h_geometry_calib)
334        hh_geometry_calib=guidata(h_geometry_calib);
335        edit_test=get(hh_geometry_calib.CheckEnableMouse,'Value');
336        hh=findobj(hcurrentaxes,'tag','calib_points');%look for handle of calibration points
337        if ~isempty(hh) && edit_test
338            index_point=get(hh,'UserData');
339            set(hh,'UserData',[])%remove edit mode
340            h_ListCoord=hh_geometry_calib.ListCoord; %handles of the coordinate list
341            Coord=get(h_ListCoord,'Data');
342           % data=read_geometry_calib(Coord);
343            %         val=get(h_ListCoord,'Value');
344          %  xy=get(hcurrentaxes,'CurrentPoint');%xy(1,1),xy(1,2): current x,y positions in axes coordinates
345            Coord(index_point,4)=xy(1,1);
346            Coord(index_point,5)=xy(1,2);
347            set(h_ListCoord,'Data',Coord)
348%             for ipoint=1:size(Coord,1)
349%                 for jcoord=1:5
350%                     Coord_cell{ipoint,jcoord}=num2str(data.Coord(ipoint,jcoord),4);%display coordiantes with 4 digits
351%                 end
352%             end
353%             Tabchar=cell2tab(Coord_cell,' | ');
354%             Tabchar=[Tabchar ;{'......'}];
355%             set(h_ListCoord,'String',Tabchar)
356            set(hh,'XData',Coord(:,4))
357            set(hh,'YData',Coord(:,5))
358        end
359    end
360end
361
362%% finalising ruler
363if test_ruler && ~isempty(xy)
364    set(hhuvmat.MenuRuler,'checked','off')%desable the ruler option in uvmat
365    xy=get(hcurrentaxes,'CurrentPoint');% get the current mouse coordinates
366    RulerCoord=[AxeData.RulerCoord ;xy(1,1:2)];% append the recorded ruler origin to the current mouse coordinates
367    RulerCoord=diff(RulerCoord,1);% coordiante difference between segment end and beginning
368    RulerCoord=RulerCoord(1)+i*RulerCoord(2);
369    distance=abs(RulerCoord);
370    azimuth=(180/pi)*angle(RulerCoord);
371    msgbox_uvmat('RULER','',['length: ' num2str(distance,3) ',  angle(degrees): ' num2str(azimuth,3)])
372    delete(AxeData.RulerHandle)%delete the ruler graphic object
373    AxeData=rmfield(AxeData,'RulerHandle');%remove the ruler handle in AxeData
374    AxeData.Drawing='off';%exit the ruler drawing mode
375end
376
377
378%% update
379if test_drawing==0
380        AxeData.Drawing='off';%stop current drawing action
381end
382set(hcurrentaxes,'UserData',AxeData);
383if ~isempty(huvmat)
384    set(huvmat,'UserData',UvData);
385end
386
387%------------------------------------------------------------------------   
388% --- 'close_fig': function  activated when a zoom figure is closed
389%------------------------------------------------------------------------
390function close_fig(ggg,eventdata,hparent)
391
392hfig=get(get(hparent,'parent'),'parent');
393hbutton=findobj(hfig,'Tag','CheckZoomFig');
394if ~isempty(hbutton)
395    set(hbutton,'Value',0)% desactivate the zoom fig option
396end
397delete(hparent)  % delete the rectangle showing the zoom graph in the parent fig
398
Note: See TracBrowser for help on using the repository browser.