[11] | 1 | %'mouse_motion': permanently called by mouse motion over a figure (Callback for 'WindowButtonMotionFcn' of the figure) |
---|
| 2 | %----------------------------------------------------------------------- |
---|
| 3 | % |
---|
| 4 | % function mouse_motion(hObject,eventdata,handles) |
---|
| 5 | % activated by the command: |
---|
| 6 | % set(hObject,'WindowButtonMotionFcn',{'mouse_motion',handles}) |
---|
| 7 | % where hObject is the handle of the figure |
---|
[809] | 8 | |
---|
| 9 | %======================================================================= |
---|
[1061] | 10 | % Copyright 2008-2019, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France |
---|
[809] | 11 | % http://www.legi.grenoble-inp.fr |
---|
| 12 | % Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr |
---|
[11] | 13 | % |
---|
| 14 | % This file is part of the toolbox UVMAT. |
---|
[809] | 15 | % |
---|
[11] | 16 | % UVMAT is free software; you can redistribute it and/or modify |
---|
[809] | 17 | % it under the terms of the GNU General Public License as published |
---|
| 18 | % by the Free Software Foundation; either version 2 of the license, |
---|
| 19 | % or (at your option) any later version. |
---|
| 20 | % |
---|
[11] | 21 | % UVMAT is distributed in the hope that it will be useful, |
---|
| 22 | % but WITHOUT ANY WARRANTY; without even the implied warranty of |
---|
| 23 | % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|
[809] | 24 | % GNU General Public License (see LICENSE.txt) for more details. |
---|
| 25 | %======================================================================= |
---|
[11] | 26 | |
---|
| 27 | function mouse_motion(hObject,eventdata,handles) |
---|
[78] | 28 | |
---|
[187] | 29 | FigData=get(hObject,'UserData'); |
---|
| 30 | if ishandle(FigData)% case of a zoom plot, the handle of the parent rectangle is stored in UserData, its parent is the plotting axes of the rectangle |
---|
[625] | 31 | hCurrentFig=get(get(FigData,'parent'),'parent'); |
---|
[187] | 32 | else |
---|
[625] | 33 | hCurrentFig=hObject;%usual plot |
---|
[187] | 34 | end |
---|
[714] | 35 | if strcmp(get(hCurrentFig,'Pointer'),'watch') |
---|
| 36 | return % no action if a calculation is running |
---|
| 37 | end |
---|
[625] | 38 | hhCurrentFig=guidata(hCurrentFig);%handles of the elements in the GUI containing the current figure (uvmat or view_field) |
---|
[622] | 39 | CheckZoom=get(hhCurrentFig.CheckZoom,'Value');% check for zoom on mode |
---|
| 40 | CheckZoomFig=get(hhCurrentFig.CheckZoomFig,'Value');% check for zoom sub fig creation mode |
---|
[625] | 41 | hPlotAxes=hhCurrentFig.PlotAxes';% handles of the main plot axes |
---|
| 42 | AxeData=get(hPlotAxes,'UserData');% data attached to the axis |
---|
[627] | 43 | htext_display(1)=handles.text_display; |
---|
| 44 | if isfield(AxeData,'htext_display')&&ishandle(AxeData.htext_display) |
---|
| 45 | htext_display(2)=AxeData.htext_display; |
---|
| 46 | end |
---|
[625] | 47 | test_draw=0;%test for mouse drawing of object, =0 by default |
---|
| 48 | if isfield(AxeData,'Drawing')&& ~isempty(AxeData.Drawing) |
---|
| 49 | test_draw=~isequal(AxeData.Drawing,'off');%=1 if mouse drawing of object is active |
---|
| 50 | end |
---|
[622] | 51 | test_zoom_draw=0; |
---|
[150] | 52 | test_object=0; %test for object editing or creation |
---|
[701] | 53 | test_create_object=0; |
---|
[681] | 54 | test_edit_object=0;% edit test for mouse shape: an arrow |
---|
[150] | 55 | test_ruler=0;%test for active ruler |
---|
[663] | 56 | test_transform=0; |
---|
[150] | 57 | huvmat=findobj(allchild(0),'tag','uvmat');%find the uvmat interface handle |
---|
[71] | 58 | if ~isempty(huvmat) |
---|
[432] | 59 | hhuvmat=guidata(huvmat);%handles of the elements in uvma |
---|
[701] | 60 | test_create_object=strcmp(get(hhuvmat.MenuObject,'checked'),'on'); |
---|
[650] | 61 | test_edit_object=get(hhuvmat.CheckEditObject,'Value'); |
---|
[156] | 62 | test_ruler=isequal(get(hhuvmat.MenuRuler,'checked'),'on'); |
---|
[664] | 63 | test_transform=~isequal(get(hhuvmat.TransformName,'Value'),1); |
---|
[150] | 64 | end |
---|
[231] | 65 | test_piv=0; |
---|
[856] | 66 | if isfield(FigData,'CivHandle')% look for handle of the civ_input GUI |
---|
[231] | 67 | if ~ishandle(FigData.CivHandle) |
---|
| 68 | delete(hObject) |
---|
[682] | 69 | set(hCurrentFig,'Pointer','arrow'); |
---|
[231] | 70 | return |
---|
| 71 | end |
---|
[856] | 72 | hhciv=guidata(FigData.CivHandle);% list of handles in the GUI civ_input |
---|
[231] | 73 | test_piv=1; |
---|
[150] | 74 | end |
---|
[432] | 75 | %find the current axe 'CurrentAxes' and display the current mouse position or uicontrol tag |
---|
[11] | 76 | text_displ_1=''; |
---|
| 77 | text_displ_2=''; |
---|
| 78 | text_displ_3=''; |
---|
| 79 | text_displ_4=''; |
---|
| 80 | |
---|
[625] | 81 | % AxeData=[];%default |
---|
[67] | 82 | xy=[];%default |
---|
[741] | 83 | set(hCurrentFig,'Units','normalized') |
---|
[622] | 84 | xy_fig=get(hObject,'CurrentPoint');% current point of the current figure (gcbo) |
---|
[11] | 85 | pointershape='arrow';% default pointer is an arrow |
---|
[150] | 86 | |
---|
[432] | 87 | %% loop on all the objects in the current figure, detect whether the mouse is over a plot axes |
---|
[622] | 88 | hchildren=get(hObject,'Children');%handles of all objects in the current figure |
---|
| 89 | check_visible=strcmp(get(hchildren,'Visible'),'on');% if visible='on', =0 otherwise |
---|
| 90 | hchildren=hchildren(check_visible); %kkep only the visible children |
---|
| 91 | PosChildren=get(hchildren,'Position');% set of object positions |
---|
| 92 | if iscell(PosChildren)% only one child |
---|
| 93 | PosLength=cellfun('length',PosChildren);% set of vector lengths for object positions |
---|
| 94 | hchildren=hchildren(PosLength==4);% keep only objects with position defined by a 4 element vector |
---|
| 95 | PosChildren=cell2mat(PosChildren(PosLength==4));% convert cells to matrix of positions |
---|
| 96 | end |
---|
| 97 | if size(PosChildren,2)~=4 |
---|
[682] | 98 | set(hCurrentFig,'Pointer','arrow'); |
---|
[622] | 99 | return |
---|
| 100 | end |
---|
| 101 | xy_fig_mat=ones(size(PosChildren,1),1)*xy_fig;% mouse position set to a matrix |
---|
| 102 | check_pos=xy_fig_mat >= PosChildren(:,1:2) & xy_fig_mat <= PosChildren(:,1:2)+PosChildren(:,3:4);% compare object to mouse position |
---|
| 103 | ind_object=find(check_pos(:,1) & check_pos(:,2),1);% select the index of the (first) object under the mouse |
---|
| 104 | hchild=hchildren(ind_object);% corresponding object handle |
---|
[432] | 105 | CurrentAxes=[]; |
---|
[622] | 106 | |
---|
| 107 | %if the mouse is over an axis, look at the data |
---|
[682] | 108 | htype=get(hchild,'Type'); |
---|
| 109 | if strcmp(htype,'axes') |
---|
[622] | 110 | CurrentAxes=hchild; |
---|
| 111 | xy=get(CurrentAxes,'CurrentPoint');%xy(1,1),xy(1,2): current x,y positions in axes coordinates |
---|
| 112 | test_zoom_draw=test_draw && isequal(AxeData.Drawing,'zoom')&& isfield(AxeData,'CurrentOrigin') && isequal(get(gcf,'SelectionType'),'normal'); |
---|
| 113 | test_object=test_draw && isfield(AxeData,'CurrentObject') && ~isempty(AxeData.CurrentObject) && ishandle(AxeData.CurrentObject); |
---|
[701] | 114 | if CheckZoom |
---|
| 115 | pointershape='zoom'; |
---|
| 116 | elseif CheckZoomFig |
---|
| 117 | pointershape='zoomfig'; |
---|
| 118 | elseif ~test_edit_object && ~test_ruler |
---|
[699] | 119 | if CheckZoom |
---|
| 120 | pointershape='zoom'; |
---|
[701] | 121 | elseif test_draw|| test_create_object |
---|
| 122 | pointershape='crosshair';%set pointer with cross shape |
---|
[699] | 123 | else |
---|
[886] | 124 | pointershape='crosshair';%set pointer with large cross (default when mouse is over an axis) |
---|
[699] | 125 | end |
---|
[622] | 126 | end |
---|
[859] | 127 | FigData=get(hCurrentFig,'UserData'); % user data of the current figure |
---|
| 128 | tagaxes=get(CurrentAxes,'tag'); % tag name of the current axes, for instance PlotAxes |
---|
[622] | 129 | if isfield(FigData,tagaxes) |
---|
[859] | 130 | Field=FigData.(tagaxes); % the current field is sought as substructure of FigData with key tagaxes |
---|
[1009] | 131 | |
---|
[622] | 132 | if isfield(Field,'ListVarName') |
---|
| 133 | [CellInfo,NbDimArray]=find_field_cells(Field);%analyse the physical fields contained in Field |
---|
| 134 | text_displ_1=''; |
---|
| 135 | text_displ_2=''; |
---|
| 136 | text_displ_3=''; |
---|
| 137 | text_displ_4=''; |
---|
| 138 | text_displ_5=''; |
---|
| 139 | ivec=[]; |
---|
| 140 | xName=''; |
---|
| 141 | z=[]; |
---|
| 142 | for icell=1:numel(CellInfo)%look for all physical fields |
---|
| 143 | if NbDimArray(icell)>=2 % select 2D field |
---|
| 144 | if isfield(Field,'CoordMesh') && ~isempty(Field.CoordMesh)&& strcmp(CellInfo{icell}.CoordType,'scattered')%case of unstructured data |
---|
| 145 | X=Field.(Field.ListVarName{CellInfo{icell}.CoordIndex(end)}); |
---|
| 146 | Y=Field.(Field.ListVarName{CellInfo{icell}.CoordIndex(end-1)}); |
---|
| 147 | flag_vec=(X<(xy(1,1)+Field.CoordMesh/3) & X>(xy(1,1)-Field.CoordMesh/3)) & ...%flagx=1 for the vectors with x position selected by the mouse |
---|
| 148 | (Y<(xy(1,2)+Field.CoordMesh/3) & Y>(xy(1,2)-Field.CoordMesh/3));%f |
---|
| 149 | ivec=find(flag_vec,1);% search the (first) selected vector index ivec |
---|
| 150 | hhh=findobj(CurrentAxes,'Tag','vector_marker'); |
---|
| 151 | if ~isempty(ivec) |
---|
| 152 | % mark the vectors with a circle in the absence of other operations |
---|
| 153 | if ~test_object && ~test_edit_object && ~test_ruler && ~CheckZoomFig |
---|
| 154 | pointershape='arrow'; %mouse indicates the detection of a vector |
---|
| 155 | if isempty(hhh) |
---|
[625] | 156 | set(0,'CurrentFigure',hCurrentFig) |
---|
| 157 | set(hCurrentFig,'CurrentAxes',CurrentAxes) |
---|
[622] | 158 | rectangle('Curvature',[1 1],... |
---|
| 159 | 'Position',[X(ivec)-Field.CoordMesh/2 Y(ivec)-Field.CoordMesh/2 Field.CoordMesh Field.CoordMesh],'EdgeColor','m',... |
---|
| 160 | 'LineStyle','-','Tag','vector_marker'); |
---|
[251] | 161 | else |
---|
[622] | 162 | set(hhh,'Visible','on') |
---|
| 163 | set(hhh,'Position',[X(ivec)-Field.CoordMesh/2 Y(ivec)-Field.CoordMesh/2 Field.CoordMesh Field.CoordMesh]) |
---|
[251] | 164 | end |
---|
[622] | 165 | end |
---|
| 166 | %display the field values |
---|
| 167 | for ivar=1:numel(CellInfo{icell}.VarIndex) |
---|
| 168 | VarName=Field.ListVarName{CellInfo{icell}.VarIndex(ivar)}; |
---|
| 169 | VarVal=Field.(VarName)(ivec); |
---|
| 170 | var_text=[VarName '=' num2str(VarVal,4) ',']; |
---|
| 171 | if isequal(ivar,CellInfo{icell}.CoordIndex(end))||isequal(ivar,CellInfo{icell}.CoordIndex(end-1))||isequal(ivar,CellInfo{icell}.CoordIndex(1)) |
---|
| 172 | text_displ_1=[text_displ_1 var_text]; |
---|
| 173 | elseif (isfield(CellInfo{icell},'VarIndex_vector_x') && isequal(ivar,CellInfo{icell}.VarIndex_vector_x))||isequal(ivar,CellInfo{icell}.VarIndex_vector_y)||... |
---|
| 174 | (isfield(CellInfo{icell},'VarIndex_vector_z') && isequal(ivar,CellInfo{icell}.VarIndex_vector_z)) |
---|
| 175 | text_displ_4=[text_displ_4 var_text]; |
---|
| 176 | else |
---|
| 177 | text_displ_5=[text_displ_5 var_text]; |
---|
[105] | 178 | end |
---|
[61] | 179 | end |
---|
[622] | 180 | else |
---|
| 181 | if ~isempty(hhh) |
---|
| 182 | set(hhh,'Visible','off') |
---|
| 183 | end |
---|
[71] | 184 | end |
---|
[622] | 185 | elseif strcmp(CellInfo{icell}.CoordType,'grid') %structured coordinates |
---|
| 186 | yName=Field.ListVarName{CellInfo{icell}.CoordIndex(1)}; |
---|
| 187 | xName=Field.ListVarName{CellInfo{icell}.CoordIndex(2)}; |
---|
| 188 | y=Field.(yName); |
---|
| 189 | x=Field.(xName); |
---|
| 190 | VarName=Field.ListVarName{CellInfo{icell}.VarIndex(1)}; |
---|
| 191 | nxy=size(Field.(VarName)); |
---|
| 192 | MaxAY=max(y(1),y(end)); |
---|
| 193 | MinAY=min(y(1),y(end)); |
---|
| 194 | if (xy(1,1)>x(1))&(xy(1,1)<x(end))&(xy(1,2)<MaxAY)&(xy(1,2)>MinAY) |
---|
| 195 | indx0=1+round((nxy(2)-1)*(xy(1,1)-x(1))/(x(end)-x(1))); % index x of pixel |
---|
| 196 | indy0=1+round((nxy(1)-1)*(xy(1,2)-y(1))/(y(end)-y(1))); % index y of pixel |
---|
| 197 | if indx0>=1 & indx0<=nxy(2) & indy0>=1 & indy0<=nxy(1) |
---|
| 198 | text_displ_2=['i=' num2str(indx0) ',j=' num2str(indy0) ',']; |
---|
| 199 | for ivar=1:numel(CellInfo{icell}.VarIndex) |
---|
| 200 | VarName=Field.ListVarName{CellInfo{icell}.VarIndex(ivar)}; |
---|
| 201 | VarVal=Field.(VarName)(indy0,indx0,:); |
---|
| 202 | var_text=[VarName '=' num2str(VarVal) ',']; |
---|
| 203 | text_displ_4=[text_displ_4 var_text]; |
---|
| 204 | end |
---|
| 205 | end |
---|
[251] | 206 | end |
---|
| 207 | end |
---|
[622] | 208 | end |
---|
| 209 | end |
---|
| 210 | % display the current x,y plot coordinates in the absence of detected vector |
---|
| 211 | if isempty(ivec) |
---|
| 212 | if isempty(xName) |
---|
| 213 | xName='x'; |
---|
| 214 | yName='y'; |
---|
| 215 | end |
---|
| 216 | text_displ_1=[xName '=' num2str(xy(1,1),4) ', ' yName '=' num2str(xy(1,2),4) ',']; |
---|
| 217 | end |
---|
[970] | 218 | %display the original phys x,y,z coordinates if defined by the projection plane |
---|
[741] | 219 | if isfield(Field,'ProjObjectType') && strcmp(Field.ProjObjectType,'plane') && isfield(Field,'ProjObjectCoord') && length(Field.ProjObjectCoord)>=3 |
---|
[970] | 220 | pos=[xy(1,1) xy(1,2) 0];%coordinates on the graph |
---|
[741] | 221 | if isfield(Field,'ProjObjectAngle')&&~isequal(Field.ProjObjectAngle,[0 0 0]) |
---|
| 222 | om=norm(Field.ProjObjectAngle);%norm of rotation angle in radians |
---|
| 223 | OmAxis=Field.ProjObjectAngle/om; %unit vector marking the rotation axis |
---|
[622] | 224 | cos_om=cos(pi*om/180); |
---|
| 225 | sin_om=sin(pi*om/180); |
---|
| 226 | pos=[xy(1,1) xy(1,2) 0]; |
---|
[972] | 227 | %pos=cos_om*pos+sin_om*cross(OmAxis,pos)+(1-cos_om)*(OmAxis*pos')*OmAxis; |
---|
| 228 | coeff=OmAxis(3)*(1-cos_om); |
---|
| 229 | norm_plane(1)=OmAxis(1)*coeff+OmAxis(2)*sin_om; |
---|
| 230 | norm_plane(2)=OmAxis(2)*coeff-OmAxis(1)*sin_om; |
---|
| 231 | norm_plane(3)=OmAxis(3)*coeff+cos_om; |
---|
| 232 | %Z0=norm_plane*Field.ProjObjectCoord'/norm_plane(3); |
---|
| 233 | pos(3)=-(norm_plane(1)*pos(1)-norm_plane(2)*pos(2))/norm_plane(3); |
---|
[622] | 234 | end |
---|
[972] | 235 | pos=pos+Field.ProjObjectCoord; |
---|
[622] | 236 | text_displ_3=[text_displ_3 'x,y,z=' num2str(pos,4)]; |
---|
| 237 | end |
---|
| 238 | % case of PIV correlation display |
---|
| 239 | if test_piv |
---|
[856] | 240 | [dd,ind_pt]=min(abs(Field.X-xy(1,1))+abs(Field.Y-xy(1,2))); |
---|
[862] | 241 | if isfield(hhciv,'TestCiv2') && strcmp(get(hhciv.Civ2,'Visible'),'on')&& ... |
---|
| 242 | strcmp(get(hhciv.TestCiv2,'Visible'),'on')&& get(hhciv.TestCiv2,'Value')% if TestCiv2 is activated |
---|
[856] | 243 | CivOption='Civ2'; |
---|
| 244 | Param.CheckCiv1=0; |
---|
| 245 | par_civ=read_GUI(hhciv.Civ2);%read the Civ2 panel in civ_input |
---|
[862] | 246 | par_civ.Civ1_SubRange=Field.Civ1_SubRange;% get the subranges used for patch |
---|
| 247 | par_civ.Civ1_NbCentres=Field.Civ1_NbCentres;% get the nbre of civ1 data points in each subrange |
---|
| 248 | par_civ.Civ1_Coord_tps=Field.Civ1_Coord_tps;% get the positions of the centres (Civ1 data) |
---|
| 249 | par_civ.Civ1_U_tps=Field.Civ1_U_tps;% get the tps coefficients for spline, U component |
---|
| 250 | par_civ.Civ1_V_tps=Field.Civ1_V_tps;% get the tps coefficients for spline, V component |
---|
[856] | 251 | par_civ.Civ1_Dt=Field.Civ1_Dt; |
---|
[858] | 252 | par_civ.Civ2_Dt=Field.Civ2_Dt; |
---|
[862] | 253 | shiftx=Field.ShiftX(ind_pt);% get field info stored in the cuirrent figure |
---|
[856] | 254 | shifty=Field.ShiftY(ind_pt); |
---|
| 255 | else |
---|
| 256 | CivOption='Civ1'; |
---|
| 257 | Param.CheckCiv2=0; |
---|
| 258 | par_civ=read_GUI(hhciv.Civ1);%read the Civ1 panel in civ_input |
---|
| 259 | shiftx=par_civ.SearchBoxShift(1); |
---|
| 260 | shifty=par_civ.SearchBoxShift(2); |
---|
| 261 | end |
---|
[622] | 262 | xround=Field.X(ind_pt); |
---|
| 263 | yround=Field.Y(ind_pt); |
---|
[856] | 264 | |
---|
[859] | 265 | % mark the correlation box with a rectangle in view_field |
---|
[856] | 266 | ibx2=floor((par_civ.CorrBoxSize(1)-1)/2); |
---|
| 267 | iby2=floor((par_civ.CorrBoxSize(2)-1)/2); |
---|
| 268 | isx2=floor((par_civ.SearchBoxSize(1)-1)/2); |
---|
| 269 | isy2=floor((par_civ.SearchBoxSize(2)-1)/2); |
---|
[622] | 270 | hhh=findobj(CurrentAxes,'Tag','PIV_box_marker'); |
---|
| 271 | hhhh=findobj(CurrentAxes,'Tag','PIV_search_marker'); |
---|
| 272 | if isempty(hhh) |
---|
[625] | 273 | set(0,'CurrentFigure',hCurrentFig) |
---|
| 274 | set(hCurrentFig,'CurrentAxes',CurrentAxes) |
---|
[622] | 275 | rectangle('Curvature',[0 0],... |
---|
| 276 | 'Position',[xround-ibx2 yround-iby2 2*ibx2 2*iby2],'EdgeColor','m',... |
---|
| 277 | 'LineStyle','-','Tag','PIV_box_marker'); |
---|
| 278 | rectangle('Curvature',[0 0],... |
---|
| 279 | 'Position',[xround-isx2+shiftx yround-isy2+shifty 2*isx2 2*isy2],'EdgeColor','m',... |
---|
| 280 | 'LineStyle','- -','Tag','PIV_search_marker'); |
---|
| 281 | else |
---|
| 282 | set(hhh,'Position',[xround-ibx2 yround-iby2 2*ibx2 2*iby2]) |
---|
| 283 | set(hhhh,'Position',[xround-isx2+shiftx yround-isy2+shifty 2*isx2 2*isy2]) |
---|
| 284 | end |
---|
[856] | 285 | |
---|
| 286 | % perform the PIV calculation as the single point xround yround |
---|
| 287 | Param.CheckFix1=0; |
---|
| 288 | Param.CheckPatch1=0;%desactivate all calculations except Civ2 or Civ1 |
---|
[855] | 289 | Param.Action.RUN=1; |
---|
| 290 | Param.ActionInput.ListCompareMode='PIV'; |
---|
[856] | 291 | par_civ.ImageA=Field.A; |
---|
| 292 | par_civ.ImageB=Field.B; |
---|
| 293 | par_civ.ImageHeight=size(par_civ.ImageA,1); |
---|
| 294 | par_civ.ImageWidth=size(par_civ.ImageA,2); |
---|
| 295 | par_civ.Grid=[xround yround];% PIV calculation with a single point |
---|
[858] | 296 | Param.ActionInput.(CivOption)=par_civ; |
---|
[855] | 297 | [Data,errormsg,result_conv]= civ_series(Param); |
---|
[622] | 298 | if ~isempty(errormsg) |
---|
| 299 | text_displ_5=errormsg; |
---|
| 300 | else |
---|
| 301 | rangx(1)=-(isx2-ibx2)+shiftx; |
---|
| 302 | rangx(2)=isx2-ibx2+shiftx; |
---|
| 303 | rangy(1)=-(isy2-iby2)-shifty; |
---|
| 304 | rangy(2)=(isy2-iby2)-shifty; |
---|
| 305 | hcorr=[]; |
---|
| 306 | if isfield(AxeData,'CurrentCorrImage') |
---|
| 307 | hcorr=AxeData.CurrentCorrImage; |
---|
| 308 | if ~ishandle(hcorr) |
---|
| 309 | hcorr=[]; |
---|
[207] | 310 | end |
---|
[11] | 311 | end |
---|
[622] | 312 | if isempty(hcorr) |
---|
| 313 | corrfig=findobj(allchild(0),'tag','corrfig'); |
---|
| 314 | if ~isempty(corrfig) |
---|
| 315 | set(0,'CurrentFigure',corrfig(1)) |
---|
[862] | 316 | AxeData.CurrentCorrImage=imagesc(rangx,-rangy,result_conv,[0 1]);% plot the correlation map, with full colormap from 0 to 1 |
---|
[858] | 317 | AxeData.CurrentVector=line([0 Data.([CivOption '_U'])],[0 Data.([CivOption '_V'])],'Tag','vector'); |
---|
[862] | 318 | AxeData.TitleHandle=title(['[x,y]= ' num2str(par_civ.Grid) ' px']); |
---|
[622] | 319 | colorbar |
---|
| 320 | set(CurrentAxes,'UserData',AxeData) |
---|
[862] | 321 | ParentAxe=get(AxeData.CurrentCorrImage,'parent'); |
---|
| 322 | set(ParentAxe,'YDir','normal') |
---|
| 323 | set(ParentAxe,'DataAspectRatio',[1 1 1]); %equal axis scale |
---|
[150] | 324 | end |
---|
[622] | 325 | else |
---|
| 326 | set(AxeData.CurrentCorrImage,'CData',result_conv) |
---|
| 327 | set(AxeData.CurrentCorrImage,'XData',rangx) |
---|
| 328 | set(AxeData.CurrentCorrImage,'YData',-rangy) |
---|
[862] | 329 | ParentAxe=get(AxeData.CurrentCorrImage,'parent'); |
---|
| 330 | rangx(1)=min(rangx(1),0); |
---|
| 331 | rangx(2)=max(rangx(2),0); |
---|
| 332 | rangy(2)=max(rangy(2),0); |
---|
| 333 | rangy(1)=min(rangy(1),0); |
---|
| 334 | set(ParentAxe,'XLim',rangx) |
---|
| 335 | set(ParentAxe,'YLim',[-rangy(2) -rangy(1)]) |
---|
[858] | 336 | set(AxeData.CurrentVector,'XData',[0 Data.([CivOption '_U'])],'YData',[0 Data.([CivOption '_V'])]); |
---|
[862] | 337 | set(AxeData.TitleHandle,'String',['[x,y]= ' num2str(par_civ.Grid) ' px']) |
---|
[11] | 338 | end |
---|
| 339 | end |
---|
[251] | 340 | end |
---|
[11] | 341 | end |
---|
| 342 | end |
---|
| 343 | end |
---|
[292] | 344 | if ~isempty(text_displ_1) |
---|
[542] | 345 | text_displ=[{text_displ_1};{text_displ_2};{text_displ_3};{text_displ_4};{text_displ_5}]; |
---|
| 346 | ind_blank=find(strcmp('',text_displ)); |
---|
| 347 | if ~isempty(ind_blank) |
---|
| 348 | text_displ(ind_blank)=[]; |
---|
| 349 | end |
---|
[627] | 350 | %set(handles.text_display,'String',text_displ) |
---|
| 351 | set(htext_display,'String',text_displ) |
---|
[292] | 352 | else |
---|
[627] | 353 | %set(handles.text_display,'String',get(handles.text_display,'UserData')) |
---|
| 354 | set(htext_display,'String',get(handles.text_display,'UserData')) |
---|
[292] | 355 | end |
---|
[11] | 356 | |
---|
[71] | 357 | %%%%%%%%%%%%% |
---|
[622] | 358 | %% draw a zoom rectangle if checkZoomFig has been selected |
---|
[71] | 359 | if test_zoom_draw |
---|
[622] | 360 | xy_rect=AxeData.CurrentOrigin;% mark the previous position from mouse down |
---|
[71] | 361 | if ~isempty(xy_rect) |
---|
| 362 | rect(1)=min(xy(1,1),xy_rect(1));%origin rectangle, x coordinate |
---|
| 363 | rect(2)=min(xy(1,2),xy_rect(2));%origin rectangle, y coordinate |
---|
| 364 | rect(3)=abs(xy(1,1)-xy_rect(1));%rectangle width |
---|
| 365 | rect(4)=abs(xy(1,2)-xy_rect(2));%rectangle height |
---|
[622] | 366 | if rect(3)>0 && rect(4)>0 |
---|
| 367 | if isfield(AxeData,'CurrentRectZoom')&& ~isempty(AxeData.CurrentRectZoom) && ishandle(AxeData.CurrentRectZoom) |
---|
[71] | 368 | set(AxeData.CurrentRectZoom,'Position',rect);%update the rectangle position |
---|
| 369 | else |
---|
[622] | 370 | AxeData.CurrentRectZoom=rectangle('Position',rect,'Tag','rect_zoom','EdgeColor','b'); |
---|
[432] | 371 | set(CurrentAxes,'UserData',AxeData) |
---|
[71] | 372 | end |
---|
| 373 | end |
---|
| 374 | end |
---|
| 375 | end |
---|
| 376 | |
---|
[11] | 377 | %%%%%%%%%%%%%%%%% |
---|
[150] | 378 | %% create or modify an object |
---|
[682] | 379 | if strcmp(htype,'axes') && ~isempty(huvmat) && test_object |
---|
[432] | 380 | UvData=get(huvmat,'UserData'); |
---|
[11] | 381 | PlotData=get(AxeData.CurrentObject,'UserData'); |
---|
[432] | 382 | if ~isfield(PlotData,'IndexObj') |
---|
[682] | 383 | set(hCurrentFig,'Pointer','arrow'); |
---|
[432] | 384 | return |
---|
| 385 | end |
---|
[966] | 386 | if numel(UvData.ProjObject)<PlotData.IndexObj |
---|
| 387 | return |
---|
| 388 | end |
---|
[622] | 389 | ObjectData=UvData.ProjObject{PlotData.IndexObj}; |
---|
[432] | 390 | if isequal(hObject,huvmat)% if the mouse ifs over the GUI uvmat |
---|
[966] | 391 | if numel(UvData.ProjObject)<get(hhuvmat.ListObject_1,'Value') |
---|
| 392 | return |
---|
| 393 | end |
---|
[622] | 394 | ProjObject=UvData.ProjObject{get(hhuvmat.ListObject_1,'Value')}; |
---|
[432] | 395 | else |
---|
[966] | 396 | if numel(UvData.ProjObject)<get(hhuvmat.ListObject,'Value') |
---|
| 397 | return |
---|
| 398 | end |
---|
[622] | 399 | ProjObject=UvData.ProjObject{get(hhuvmat.ListObject,'Value')}; |
---|
[432] | 400 | end |
---|
| 401 | XYData=AxeData.CurrentOrigin; |
---|
| 402 | if isequal(AxeData.Drawing,'create') && isfield(AxeData,'CurrentOrigin') && ~isempty(AxeData.CurrentOrigin) |
---|
[954] | 403 | switch ObjectData.Type |
---|
| 404 | case {'line','polyline','polygon','points','plane_z'} |
---|
| 405 | ObjectData.Coord=[ObjectData.Coord ;xy(1,1:2)]; |
---|
| 406 | % ObjectData.Coord(end,:)=xy(1,:); |
---|
| 407 | case {'rectangle','ellipse','volume'} |
---|
| 408 | ObjectData.Coord=(AxeData.CurrentOrigin+xy(1,1:2))/2;% keep only the first point coordinate |
---|
[622] | 409 | ObjectData.RangeX=abs(ObjectData.Coord(1,1)-xy(1,1));%rectangle width |
---|
[954] | 410 | ObjectData.RangeY=abs(ObjectData.Coord(1,2)-xy(1,2));%rectangle height |
---|
| 411 | case 'plane' %case of 'plane' |
---|
| 412 | DX=(xy(1,1)-ObjectData.Coord(1,1)); |
---|
| 413 | DY=(xy(1,2)-ObjectData.Coord(1,2)); |
---|
| 414 | ObjectData.Phi=(angle(DX+i*DY))*180/pi;%rectangle widt |
---|
| 415 | if isfield(ObjectData,'RangeX') |
---|
| 416 | XMax=sqrt(DX*DX+DY*DY); |
---|
| 417 | if XMax>max(ObjectData.RangeX) |
---|
| 418 | ObjectData.RangeX=[min(ObjectData.RangeX) XMax]; |
---|
| 419 | end |
---|
[11] | 420 | end |
---|
| 421 | end |
---|
[432] | 422 | plot_object(ObjectData,ProjObject,AxeData.CurrentObject,'m'); |
---|
| 423 | pointershape='crosshair'; |
---|
[681] | 424 | elseif test_edit_object && isequal(AxeData.Drawing,'translate') |
---|
[432] | 425 | DX=xy(1,1)-XYData(1);%translation from initial position |
---|
| 426 | DY=xy(1,2)-XYData(2); |
---|
| 427 | ObjectData.Coord(:,1)=ObjectData.Coord(:,1)+DX; |
---|
| 428 | ObjectData.Coord(:,2)=ObjectData.Coord(:,2)+DY; |
---|
| 429 | plot_object(ObjectData,ProjObject,AxeData.CurrentObject,'m'); |
---|
| 430 | pointershape='fleur'; |
---|
[681] | 431 | elseif test_edit_object && isequal(AxeData.Drawing,'deform') |
---|
[432] | 432 | ind_move=AxeData.CurrentIndex; |
---|
| 433 | ObjectData.Coord(ind_move,1)=xy(1,1); |
---|
| 434 | ObjectData.Coord(ind_move,2)=xy(1,2); |
---|
| 435 | plot_object(ObjectData,ProjObject,AxeData.CurrentObject,'m'); |
---|
| 436 | pointershape='circle'; |
---|
[11] | 437 | end |
---|
[432] | 438 | end |
---|
[67] | 439 | |
---|
[150] | 440 | %% detect calibration points if the GUI geometry_calib is opened |
---|
[67] | 441 | h_geometry_calib=findobj(allchild(0),'Name','geometry_calib'); %find the geomterty_calib GUI |
---|
[682] | 442 | if strcmp(htype,'axes') && ~CheckZoom && ~isempty(h_geometry_calib) |
---|
[67] | 443 | pointershape='crosshair';%default for geometry_calib: ready to create new points |
---|
| 444 | hh_geometry_calib=guidata(h_geometry_calib); |
---|
[156] | 445 | if ~isempty(xy) && isfield(hh_geometry_calib,'ListCoord') |
---|
[67] | 446 | h_ListCoord=hh_geometry_calib.ListCoord; %findobj(h_geometry_calib,'Tag','ListCoord'); |
---|
[657] | 447 | data.Coord=get(h_ListCoord,'Data'); |
---|
| 448 | if isnumeric(data.Coord)&&~isempty(data.Coord) |
---|
[663] | 449 | if test_transform |
---|
[738] | 450 | XCoord=(data.Coord(:,1)); |
---|
| 451 | YCoord=(data.Coord(:,2)); |
---|
[663] | 452 | else |
---|
[738] | 453 | XCoord=(data.Coord(:,4)); |
---|
| 454 | YCoord=(data.Coord(:,5)); |
---|
[663] | 455 | end |
---|
[432] | 456 | xy=get(CurrentAxes,'CurrentPoint');%xy(1,1),xy(1,2): current x,y positions in axes coordinates |
---|
[67] | 457 | if ~isempty(xy) |
---|
[432] | 458 | xlim=get(CurrentAxes,'XLim'); |
---|
[738] | 459 | % ind_range_x=abs((xlim(2)-xlim(1))/50); |
---|
[432] | 460 | ylim=get(CurrentAxes,'YLim'); |
---|
[738] | 461 | ind_range=max(abs(xlim(2)-xlim(1)),abs(ylim(end)-ylim(1)))/25;%defines the size of the circle marker |
---|
| 462 | % ind_range_y=abs((ylim(2)-ylim(1))/50); |
---|
| 463 | % ind_range=sqrt(ind_range_x*ind_range_y); |
---|
| 464 | index_point=find((XCoord<xy(1,1)+ind_range/2) & (XCoord>xy(1,1)-ind_range/2) & ...%flagx=1 for the vectors with x position selected by the mouse |
---|
| 465 | (YCoord<xy(1,2)+ind_range/2) & (YCoord>xy(1,2)-ind_range/2),1);%find the first calibration point in the neighborhood of the mouse |
---|
[67] | 466 | if ~isempty(index_point) |
---|
[738] | 467 | pointershape='arrow';% default pointer is an arrow |
---|
[78] | 468 | end |
---|
| 469 | hh=findobj('Tag','calib_points');%look for handle of calibration points |
---|
[738] | 470 | if ~isempty(hh) && ~isempty(get(hh,'UserData')) %&& get(hh_geometry_calib.CheckEnableMouse,'Value') |
---|
| 471 | %set(hh,'UserData',index_point) |
---|
[173] | 472 | index_point=get(hh,'UserData'); |
---|
[78] | 473 | XCoord(index_point)=xy(1,1); |
---|
| 474 | YCoord(index_point)=xy(1,2); |
---|
| 475 | set(hh,'XData',XCoord) |
---|
| 476 | set(hh,'YData',YCoord) |
---|
[738] | 477 | end |
---|
[78] | 478 | if ~isempty(index_point) |
---|
[738] | 479 | set(hh_geometry_calib.CoordLine,'String',num2str(index_point)) |
---|
| 480 | % Data=get(h_ListCoord,'Data'); |
---|
| 481 | % Data(:,6)=zeros(size(Data,1),1); |
---|
| 482 | % Data(index_point,6)=-1;%mrk the point on the GUI geometry_calib |
---|
| 483 | % set(h_ListCoord,'Data',Data); |
---|
[67] | 484 | hhh=findobj('Tag','calib_marker');%look for handle of point marker (circle) |
---|
| 485 | if ~isempty(hhh) |
---|
[71] | 486 | set(hhh,'Position',[XCoord(index_point)-ind_range/2 YCoord(index_point)-ind_range/2 ind_range ind_range]) |
---|
[657] | 487 | else |
---|
[837] | 488 | axes(hchild) |
---|
[738] | 489 | rectangle('Curvature',[1 1],... |
---|
| 490 | 'Position',[XCoord(index_point)-ind_range/2 YCoord(index_point)-ind_range/2 ind_range ind_range],'EdgeColor','m',... |
---|
| 491 | 'LineStyle','-','Tag','calib_marker'); |
---|
[67] | 492 | end |
---|
[78] | 493 | end |
---|
[67] | 494 | end |
---|
| 495 | end |
---|
| 496 | end |
---|
| 497 | end |
---|
| 498 | |
---|
[150] | 499 | %% draw ruler |
---|
[156] | 500 | if test_ruler && isfield(AxeData,'Drawing') && isequal(AxeData.Drawing,'ruler') |
---|
| 501 | if isfield(AxeData,'RulerHandle') |
---|
| 502 | pointershape='crosshair'; %give the mouse pointer a cross shape |
---|
| 503 | RulerCoord=[AxeData.RulerCoord ;xy(1,1:2)]; %coordinates defining the ruler segment |
---|
| 504 | set(AxeData.RulerHandle,'XData',RulerCoord(:,1));% updtate the x coordinates for the ruler graphic object |
---|
| 505 | set(AxeData.RulerHandle,'YData',RulerCoord(:,2));% updtate the y coordinates for the ruler graphic object |
---|
| 506 | end |
---|
[60] | 507 | end |
---|
[150] | 508 | |
---|
| 509 | %% update the mouse pointer |
---|
[699] | 510 | if strcmp(pointershape,'zoom')||strcmp(pointershape,'zoomfig') |
---|
| 511 | CData=set_pointershape(pointershape); |
---|
| 512 | set(hCurrentFig,'Pointer','custom','PointerShapeCData',CData,'PointerShapeHotSpot',[9 9]) |
---|
| 513 | else |
---|
[625] | 514 | set(hCurrentFig,'Pointer',pointershape); |
---|
[699] | 515 | end |
---|
| 516 | |
---|
| 517 | function CData=set_pointershape(pointershape) |
---|
| 518 | CData=ones(16,16); |
---|
| 519 | [ind_x,ind_y]=meshgrid([1:16],[1:16]); |
---|
| 520 | if strcmp(pointershape,'zoom') |
---|
| 521 | radius=(ind_x-9).*(ind_x-9)+(ind_y-8.5).*(ind_y-9); |
---|
| 522 | CData(radius<25 & radius>16)=2; %make white circle |
---|
| 523 | CData(radius<16 | radius>40)=NaN; %make the centre transparent |
---|
| 524 | CData(16,16)=1; CData(15,16)=2; |
---|
| 525 | CData(15,15)=1; CData(14,15)=2; |
---|
| 526 | CData(14,14)=1; CData(13,14)=2; |
---|
| 527 | CData(13,13)=1; CData(12,13)=2; |
---|
| 528 | elseif strcmp(pointershape,'zoomfig') |
---|
| 529 | CData(:,1:3)=2; |
---|
| 530 | CData(:,14:16)=2; |
---|
| 531 | CData(1:3,:)=2; |
---|
| 532 | CData(14:16,:)=2; |
---|
| 533 | CData(CData==1)=NaN; |
---|
| 534 | end |
---|