source: trunk/src/set_grid.m @ 604

Last change on this file since 604 was 591, checked in by sommeria, 11 years ago

various updates, in particular modification of series to do calculations in the cluster

File size: 14.2 KB
Line 
1%'set_grid':produce grid for PIV with one or two images (stereo case)
2%------------------------------------------------------------------------
3% function varargout = set_grid(varargin)
4% associated with the GUI set_grid.fig
5%
6%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
7%  Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.
8%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
9%     This file is part of the toolbox UVMAT.
10%
11%     UVMAT is free software; you can redistribute it and/or modify
12%     it under the terms of the GNU General Public License as published by
13%     the Free Software Foundation; either version 2 of the License, or
14%     (at your option) any later version.
15%
16%     UVMAT is distributed in the hope that it will be useful,
17%     but WITHOUT ANY WARRANTY; without even the implied warranty of
18%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
19%     GNU General Public License (file UVMAT/COPYING.txt) for more details.
20%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA
21
22function varargout = set_grid(varargin)
23
24% Last Modified by GUIDE v2.5 01-Mar-2013 22:41:43
25
26% Begin initialization code - DO NOT PLOT
27gui_Singleton = 1;
28gui_State = struct('gui_Name',       mfilename, ...
29                   'gui_Singleton',  gui_Singleton, ...
30                   'gui_OpeningFcn', @set_grid_OpeningFcn, ...
31                   'gui_OutputFcn',  @set_grid_OutputFcn, ...
32                   'gui_LayoutFcn',  [] , ...
33                   'gui_Callback',   []);
34               
35if nargin && ischar(varargin{1})
36    gui_State.gui_Callback = str2func(varargin{1});
37end
38if nargout
39    [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
40else
41    gui_mainfcn(gui_State, varargin{:});
42end
43
44%-------------------------------------------------------------------
45% --- Executes just before set_grid is made visible.
46%INPUT:
47% handles: handles of the set_grid interface elements
48%'IndexObj': index of the object (on the UvData list) that set_grid will modify
49%        if =[] or absent: index still undefined (create mode in uvmat)
50%        if=0; no associated object (used for series), the button 'PLOT' is  then unvisible
51%'data': read from an existing object selected in the interface
52%      .TITLE : class of object ('POINTS','LINE',....)
53%      .num_DX,num_DY,DZ; meshes for regular grids
54%      .Coord: object position coordinates
55%      .ParentButton: handle of the uicontrol object calling the interface
56% PlotHandles: set of handles of the elements contolling the plotting of the projected field:
57%  if =[] or absent, no plot (mask mode in uvmat)
58% parameters on the uvmat interface (obtained by 'get_plot_handle.m')
59function set_grid_OpeningFcn(hObject, eventdata, handles,InputFile,InputField)
60
61% Choose default command line output for set_grid
62handles.output = hObject;
63
64% Update handles structure
65guidata(hObject, handles);
66
67%default
68set(hObject,'DeleteFcn',@closefcn)
69set(hObject,'WindowButtonDownFcn',{'mouse_down'})%set mouse click action function
70set(handles.CoordType,'ListboxTop',1)
71set(handles.CoordType,'Value',1);
72set(handles.CoordType,'String',{'phys';'px'});
73if exist('InputFile','var')
74   set(handles.ImageA,'String',InputFile)
75end
76
77%% use InputField input from uvmat
78check_pixel=0;
79if exist('InputField','var')
80    if strcmp(InputField.CoordUnit,'pixel')
81        set(handles.CoordType,'Value',2)
82        set(handles.TxtWarning,'Visible','on')
83        Mesh=20;%default mesh in pixel
84        check_pixel=1;
85    else
86        set(handles.CoordType,'Value',1)
87        InputField.CoordMesh=20*InputField.CoordMesh; % about 20 pixels
88        % adjust the mesh to a value 1, 2 , 5 *10^n
89        ord=10^(floor(log10(InputField.CoordMesh)));%order of magnitude
90        if InputField.CoordMesh/ord>=5
91            Mesh=5*ord;
92        elseif InputField.CoordMesh/ord>=2
93            Mesh=2*ord;
94        else
95            Mesh=ord;
96        end
97    end
98    Input.DX=Mesh;
99    Input.DY=Mesh;
100    Input.XMin=(Mesh/2)*ceil(InputField.XMin/(Mesh/2))-0.5*check_pixel;
101    Input.XMax=Input.XMin+Mesh*floor((InputField.XMax-Input.XMin)/Mesh)-0.5*check_pixel;
102    Input.YMin=(Mesh/2)*ceil(InputField.YMin/(Mesh/2))-0.5*check_pixel;
103    Input.YMax=Input.YMin+Mesh*floor((InputField.YMax-Input.YMin)/Mesh)-0.5*check_pixel;
104    errormsg=fill_GUI(Input,handles.set_grid);
105end
106
107% --- Outputs from this function are returned to the command line.
108function varargout = set_grid_OutputFcn(hObject, eventdata, handles)
109% Get default command line output from handles structure
110varargout{1} = handles.output;
111varargout{2}=handles;
112
113%----------------------------------------------------
114% executed when closing: set the parent interface button to value 0
115function closefcn(gcbo,eventdata)
116huvmat=findobj(allchild(0),'Name','uvmat');%find the current uvmat interface handle
117parent_button=findobj(huvmat,'Tag','grid');
118if ~isempty(parent_button)
119    set(parent_button,'Value',0)%put unactivated buttons to green
120    tag=get(parent_button,'Tag');
121    if isequal(tag,'edit')
122        set(parent_button,'BackgroundColor',[0.7 0.7 0.7]);
123    else
124        set(parent_button,'BackgroundColor',[0 1 0]);
125    end
126end
127
128%-----------------------------------------------------------------------
129% --- Executes on button press in plot: PLOT the defined object and its projected field
130function plot_Callback(hObject, eventdata, handles)
131[grid_pix_A,grid_pix_B,grid_phys]=get_grid(read_GUI(handles.set_grid));
132huvmat=findobj(allchild(0),'tag','uvmat');
133hhuvmat=guidata(huvmat);
134axes(hhuvmat.PlotAxes);
135hold on
136UvData=get(huvmat,'UserData');
137if isfield(UvData.Field, 'CoordUnit')&& strcmp(UvData.Field.CoordUnit,'pixel')
138    plot(grid_pix_A(:,1),grid_pix_A(:,2),'.')
139else
140    plot(grid_phys(:,1),grid_phys(:,2),'.')
141end
142
143%% display grid in second image defiend
144if ~isempty(grid_pix_B)
145    hviewfield=view_field(get(handles.imageB,'String'));
146    hhviewfield=guidata(hviewfield);
147    axes(hhviewfield.PlotAxes);
148    hold on
149    if isfield(UvData.Field, 'CoordUnit')&& strcmp(UvData.Field.CoordUnit,'pixel')
150        plot(grid_pix_B(:,1),grid_pix_B(:,2),'.')
151    else
152        plot(grid_phys(:,1),grid_phys(:,2),'.')
153    end
154end
155
156%------------------------------------------------------------------------
157% --- Executes on button press in CoordType.
158function CoordType_Callback(hObject, eventdata, handles)
159%------------------------------------------------------------------------
160set(handles.num_XMin,'String','')
161set(handles.num_XMax,'String','')
162set(handles.num_DX,'String','')
163set(handles.num_YMin,'String','')
164set(handles.num_YMax,'String','')
165set(handles.num_DY,'String','')
166set(handles.num_Z,'String','')
167if isequal(get(handles.CoordType,'Value'),2)
168    set(handles.TxtWarning,'visible','on')
169else
170    set(handles.TxtWarning,'visible','on')
171end
172
173%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
174% --- Executes on button press in delete.
175function delete_Callback(hObject, eventdata, handles)
176
177%SetData=get(gcbf,'UserData');%get the interface data
178%IndexObj=SetData.IndexObj;
179huvmat=findobj('Name','uvmat');%find the current uvmat interface handle
180UvData=get(huvmat,'UserData');%Data associated to the current uvmat interface
181if isfield(UvData,'CurrentObjectIndex')
182    IndexObj=UvData.CurrentObjectIndex;
183else
184    IndexObj=[];
185end
186delete_object(IndexObj);
187
188
189% ------------------------------------------------------
190function Save_Callback(hObject, eventdata, handles)
191% ------------------------------------------------------
192[grid_pix_A,grid_pix_B]=get_grid(read_GUI(handles.set_grid));
193
194
195 %ECRIRE FICHIERS
196nbpointsA=size(grid_pix_A);
197XA=grid_pix_A(:,1);%index=position+0.5 rounded at the nearest integer value
198YA=grid_pix_A(:,2);
199unitcolumn=32*ones(size(XA));
200Xchar=num2str(XA);
201blanc=char(unitcolumn);
202Ychar=num2str(YA);
203tete=['1 ' num2str(nbpointsA(1))];
204txt=[Xchar blanc Ychar];
205textgrid={tete;txt};
206textout=char(textgrid);
207imageA=get(handles.ImageA,'String');
208RootPath=fileparts_uvmat(imageA);
209Answer = msgbox_uvmat('INPUT_TXT','grid file name (*.grid)',fullfile(RootPath,'gridA.grid'));
210dlmwrite(Answer,textout,'');
211msgbox_uvmat('CONFIRMATION',[Answer ' written as ASCII text file']);
212if ~isempty(grid_pix_B)
213    nbpointsB=size(grid_pix_B);
214    XB=round(grid_pix_B(:,1)+0.5);%index=position+0.5 rounded at the nearest integer value
215    YB=round(grid_pix_B(:,2)+0.5);
216    unitcolumn=32*ones(size(XB));
217    Xchar=num2str(XB);
218    blanc=char(unitcolumn);
219    Ychar=num2str(YB);
220    tete=['1 ' num2str(nbpointsB(1))];
221    txt=[Xchar blanc Ychar];
222    textgrid={tete;txt};
223    textout=char(textgrid);
224    Answer = msgbox_uvmat('INPUT_TXT','grid file name (*.grid)',fullfile(RootPath,'gridB.grid'));
225    dlmwrite(Answer,textout,'');
226    msgbox_uvmat('CONFIRMATION',[Answer ' written as ASCII text file']);
227end
228
229
230%------------------------------------------------------------------------
231function [grid_pix_A,grid_pix_B,grid_phys]=get_grid(GUI)
232%------------------------------------------------------------------------
233grid_pix_B=[];%default
234array_x=GUI.XMin:GUI.DX:GUI.XMax;% array of x values
235array_y=GUI.YMin:GUI.DY:GUI.YMax;% array of y values
236[grid_x,grid_y]=meshgrid(array_x,array_y);% matrices of x and y values
237grid_x=reshape(grid_x,[],1); %matrix of x  values reshaped in line
238grid_y=reshape(grid_y,[],1);%matrix of y values reshaped in line
239% grid_z=zeros(nx_patch*ny_patch,1);% plane coordinates (TODO: 3D grids)
240
241%% check the input image A
242if ~exist(GUI.ImageA,'file')
243    msgbox_uvmat('ERROR',['input image file' imageA 'does not exist'])
244    return
245end
246[FileType,tild,VideoObject]=get_file_type(GUI.ImageA);
247switch FileType
248    case {'image','multimage','video','mmreader'}% case of input image or movie OK
249    otherwise
250        msgbox_uvmat('ERROR',['error: ' GUI.ImageA ' is not an image type recognized by Matlab '])
251        return
252end
253[RootPath,SubDir,RootFile,tild,tild,tild,tild,FileExt]=fileparts_uvmat(GUI.ImageA);
254
255%% transform to pixels if the grid is defined in phys coordinates
256grid_x_imaA=grid_x;%default grid in image A coordinates
257grid_y_imaA=grid_y;
258% MenuCoord=get(handles.CoordType,'String');% type of coordinates for grid definition, phys or pixel
259if strcmp(GUI.CoordType,'phys')
260    fileAxml=fullfile(RootPath,[SubDir '.xml']);% new convention for xml name
261    if ~exist(fileAxml,'file')
262        fileAxml=[fullfile(RootPath,RootFile) '.xml'];% old convention for xml name
263    end
264    tsaiA=[];%default
265    if exist(fileAxml,'file')
266        [XmlDataA,errormsg]=imadoc2struct(fileAxml);
267        if ~isempty(errormsg)
268            msgbox_uvmat('ERROR',['error in ' fileAxml ': ' errormsg])
269            return
270        end
271        if isfield(XmlDataA,'GeometryCalib')
272            tsaiA=XmlDataA.GeometryCalib;
273        end
274    end
275    if isempty(tsaiA)
276        msgbox_uvmat('WARNING','no geometric calibration available for image A, phys =pixel')
277    else
278        [grid_x_imaA,grid_y_imaA]=px_XYZ(tsaiA,grid_x,grid_y,GUI.Z);
279    end
280end
281
282%% detect the grid points which are inside image A
283A=read_image(GUI.ImageA,FileType,VideoObject,1);
284npxA=size(A,2);
285npyA=size(A,1);
286flag=grid_x_imaA>=1 & grid_x_imaA<=npxA & grid_y_imaA>=1 & grid_y_imaA<=npyA;% ='true' inside the image
287
288%% detect the grid points which are inside image B if relevant (use for stereo PIV)
289if isfield(GUI,'ImageB')
290    if ~exist(imageB,'file')
291        msgbox_uvmat('ERROR',['input image file' GUI.ImageB 'does not exist'])
292        return
293    end
294    [RootPathB,SubDirB,RootFileB,tild,tild,tild,tild,FileExt]=fileparts_uvmat(GUI.ImageB);
295    fileBxml=fullfile(RootPathB,[SubDirB '.xml']);% new convention for xml name
296    if ~exist(fileBxml,'file')
297        fileBxml=[fullfile(RootPathB,RootFileB) '.xml'];% old convention for xml name
298    end
299    tsaiB=[];%default
300    if exist(fileBxml,'file')
301        [XmlDataB,errormsg]=imadoc2struct(fileBxml);
302        if ~isempty(errormsg)
303            msgbox_uvmat('ERROR',['error in ' fileAxml ': ' errormsg])
304            return
305        end
306        if isfield(XmlDataB,'GeometryCalib')
307            tsaiB=XmlDataB.GeometryCalib;
308        end
309    end
310    if isempty(tsaiB)
311        msgbox_uvmat('WARNING','no geometric calibration available for image B, phys =pixel')
312        grid_x_imaB=grid_x;
313        grid_y_imaB=grid_y;
314    else
315        [grid_x_imaB,grid_y_imaB]=px_XYZ(tsaiB,grid_x,grid_y,GUI.Z);
316    end
317    B=imread(GUI.ImageB);
318    npxB=size(B,2);
319    npyB=size(B,1);
320    flagB=grid_x_imaB>=1 & grid_x_imaB<=npxB & grid_y_imaB>=1 & grid_y_imaB<=npyB;
321    flag=flagA & flagB;
322    grid_pix_B(:,1)=round(grid_x_imaB(flag));
323    grid_pix_B(:,2)=round(grid_y_imaB(flag));
324end
325
326grid_x_imaA=grid_x_imaA(flag);
327grid_y_imaA=grid_y_imaA(flag);
328grid_pix_A=[grid_x_imaA grid_y_imaA];
329grid_x=grid_x(flag);
330grid_y=grid_y(flag);
331grid_phys=[grid_x grid_y];
332
333
334function GetImageB_Callback(hObject, eventdata, handles)
335if isequal(get(handles.GetImageB,'Value'),1)
336    set(handles.ImageB,'Visible','on')
337    [FileName, PathName, filterindex] = uigetfile( ...
338            {'*.*', 'All Files (*.*)'}, ...
339            'Pick the second image file',fileparts(fileparts(get(handles.ImageA,'String'))));
340        ImageB=fullfile(PathName,FileName);
341        [FileType,tild,VideoObject]=get_file_type(ImageB);
342    switch FileType
343        case {'image','multimage','video','mmreader'}% case of input image or movie OK
344            set(handles.ImageB,'String',ImageB)
345        otherwise
346            msgbox_uvmat('ERROR',['error: ' imageB ' is not an image type recognized by Matlab '])
347            return
348    end
349else
350    set(handles.ImageB,'Visible','off')
351end
352
353
354%------------------------------------------------------------------------
355% --- Executes on button press in HELP.
356function HELP_Callback(hObject, eventdata, handles)
357%------------------------------------------------------------------------
358path_to_uvmat=which ('uvmat');% check the path of uvmat
359pathelp=fileparts(path_to_uvmat);
360helpfile=fullfile(pathelp,'uvmat_doc','uvmat_doc.html');
361if isempty(dir(helpfile)), errordlg('Please put the help file uvmat_doc.html in the directory UVMAT/UVMAT_DOC')
362else
363web([helpfile '#set_grid'])   
364end
365
366
367
368function ImageA_Callback(hObject, eventdata, handles)
369% hObject    handle to ImageA (see GCBO)
370% eventdata  reserved - to be defined in a future version of MATLAB
371% handles    structure with handles and user data (see GUIDATA)
372
373% Hints: get(hObject,'String') returns contents of ImageA as text
374%        str2double(get(hObject,'String')) returns contents of ImageA as a double
Note: See TracBrowser for help on using the repository browser.