source: trunk/src/set_grid.m @ 1126

Last change on this file since 1126 was 1126, checked in by g7moreau, 5 months ago

Update Copyright

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