1 | <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" |
---|
2 | "http://www.w3.org/TR/REC-html40/loose.dtd"> |
---|
3 | <html> |
---|
4 | <head> |
---|
5 | <title>Description of geometry_calib</title> |
---|
6 | <meta name="keywords" content="geometry_calib"> |
---|
7 | <meta name="description" content="'geometry_calib': performs geometric calibration from a set of reference points"> |
---|
8 | <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1"> |
---|
9 | <meta name="generator" content="m2html © 2003 Guillaume Flandin"> |
---|
10 | <meta name="robots" content="index, follow"> |
---|
11 | <link type="text/css" rel="stylesheet" href="../m2html.css"> |
---|
12 | </head> |
---|
13 | <body> |
---|
14 | <a name="_top"></a> |
---|
15 | <div><a href="../index.html">Home</a> > <a href="index.html">.</a> > geometry_calib.m</div> |
---|
16 | |
---|
17 | <!--<table width="100%"><tr><td align="left"><a href="../index.html"><img alt="<" border="0" src="../left.png"> Master index</a></td> |
---|
18 | <td align="right"><a href="index.html">Index for . <img alt=">" border="0" src="../right.png"></a></td></tr></table>--> |
---|
19 | |
---|
20 | <h1>geometry_calib |
---|
21 | </h1> |
---|
22 | |
---|
23 | <h2><a name="_name"></a>PURPOSE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
24 | <div class="box"><strong>'geometry_calib': performs geometric calibration from a set of reference points</strong></div> |
---|
25 | |
---|
26 | <h2><a name="_synopsis"></a>SYNOPSIS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
27 | <div class="box"><strong>function varargout = geometry_calib(varargin) </strong></div> |
---|
28 | |
---|
29 | <h2><a name="_description"></a>DESCRIPTION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
30 | <div class="fragment"><pre class="comment">'geometry_calib': performs geometric calibration from a set of reference points |
---|
31 | |
---|
32 | function varargout = geometry_calib(varargin) |
---|
33 | |
---|
34 | A%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
---|
35 | Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org. |
---|
36 | AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA |
---|
37 | This file is part of the toolbox UVMAT. |
---|
38 | |
---|
39 | UVMAT is free software; you can redistribute it and/or modify |
---|
40 | it under the terms of the GNU General Public License as published by |
---|
41 | the Free Software Foundation; either version 2 of the License, or |
---|
42 | (at your option) any later version. |
---|
43 | |
---|
44 | UVMAT is distributed in the hope that it will be useful, |
---|
45 | but WITHOUT ANY WARRANTY; without even the implied warranty of |
---|
46 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the |
---|
47 | GNU General Public License (file UVMAT/COPYING.txt) for more details. |
---|
48 | AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</pre></div> |
---|
49 | |
---|
50 | <!-- crossreference --> |
---|
51 | <h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
52 | This function calls: |
---|
53 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
54 | <li><a href="ListDir.html" class="code" title="function [ListDevices,ListRecords,ListXml,List]=ListDir(CurrentPath,ListExperiments,ListDevices_in,ListRecords_in)">ListDir</a> 'ListDir': scan the structure of the directory tree (for dataview.m)</li><li><a href="cell2tab.html" class="code" title="function Tabchar=cell2tab(Tabcell,separator)">cell2tab</a> 'cell2tab': transform a Matlab cell in a character array suitable for display in a table</li><li><a href="dataview.html" class="code" title="function varargout = dataview(varargin)">dataview</a> TEST 'dataview': function for scanning directories in a campaign (TEST)</li><li><a href="delete_object.html" class="code" title="function delete_object(hObject)">delete_object</a> 'delete_object': delete a projection object, defined by its index in the Uvmat list or by its graphic handle</li><li><a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a> 'msgbox_uvmat': associated with GUI msgbox_uvmat.fig to display message boxes, for error, warning or input calls</li><li><a href="name2display.html" class="code" title="function [RootPath,RootFile,field_count,str2,str_a,str_b,ext,nom_type,subdir]=name2display(fileinput)">name2display</a> 'name2display': extracts the root name and field numbers from an input filename</li><li><a href="plot_field.html" class="code" title="function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)">plot_field</a> 'plot_field': plot any field with the structure defined in the uvmat package</li><li><a href="px.html" class="code" title="function [DataOut,DataOut_1]=px(Data,CalibData,Data_1,CalibData_1)%DataIn,Calib)">px</a> 'px': transform fields from physical coordinates (phys) to image (px) coordinates</li><li><a href="px_XYZ.html" class="code" title="function [X,Y]=px_XYZ(Calib,Xphys,Yphys,Zphys)">px_XYZ</a> 'px_XYZ': transform phys coordinates to image coordinates (px)</li><li><a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a> 'read_geometry_calib': read data on the GUI geometry_calib</li><li><a href="struct2xml.html" class="code" title="function t=struct2xml(Object,t,root_uid)">struct2xml</a> 'struct2xml': transform a matlab structure to a xml tree.</li><li><a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a> 'uvmat': function associated with the GUI 'uvmat.fig' for images and data field visualization</li><li><a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a> 'warndlg_uvmat': display warning message (error, warning, confirmation) in a given figure</li></ul> |
---|
55 | This function is called by: |
---|
56 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
57 | <li><a href="mouse_down.html" class="code" title="function xy=mouse_down(hObject,eventdata)">mouse_down</a> 'mouse_down': function activated when the mouse button is pressed on a figure (callback for 'WindowButtonDownFcn'</li><li><a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a> 'uvmat': function associated with the GUI 'uvmat.fig' for images and data field visualization</li></ul> |
---|
58 | <!-- crossreference --> |
---|
59 | |
---|
60 | <h2><a name="_subfunctions"></a>SUBFUNCTIONS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
61 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
62 | <li><a href="#_sub1" class="code">function geometry_calib_OpeningFcn(hObject, eventdata, handles, handles_uvmat,pos,inputfile)</a></li><li><a href="#_sub2" class="code">function varargout = geometry_calib_OutputFcn(hObject, eventdata, handles)</a></li><li><a href="#_sub3" class="code">function Phi_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub4" class="code">function OPEN_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub5" class="code">function loadfile(handles,fileinput)</a></li><li><a href="#_sub6" class="code">function closefcn(gcbo,eventdata,handles_uvmat)</a></li><li><a href="#_sub7" class="code">function PLOT_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub8" class="code">function MenuCoord_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub9" class="code">function delete_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub10" class="code">function APPLY_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub11" class="code">function REPLICATE_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub12" class="code">function update_imadoc(GeometryCalib,outputfile)</a></li><li><a href="#_sub13" class="code">function GeometryCalib=calib_rescale(Coord)</a></li><li><a href="#_sub14" class="code">function GeometryCalib=calib_linear(Coord)</a></li><li><a href="#_sub15" class="code">function GeometryCalib=calib_tsai(Coord)</a></li><li><a href="#_sub16" class="code">function translation_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub17" class="code">function T_x_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub18" class="code">function T_y_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub19" class="code">function T_z_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub20" class="code">function rotation_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub21" class="code">function XImage_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub22" class="code">function YImage_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub23" class="code">function XObject_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub24" class="code">function YObject_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub25" class="code">function ZObject_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub26" class="code">function update_list(hObject, eventdata, handles)</a></li><li><a href="#_sub27" class="code">function ListCoord_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub28" class="code">function translation_plus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub29" class="code">function translation_minus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub30" class="code">function translation(handles,T)</a></li><li><a href="#_sub31" class="code">function rotation_plus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub32" class="code">function rotation_minus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub33" class="code">function rotation(handles,Phi)</a></li><li><a href="#_sub34" class="code">function O_x_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub35" class="code">function O_y_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub36" class="code">function O_z_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub37" class="code">function edit_append_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub38" class="code">function NEW_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub39" class="code">function HELP_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub40" class="code">function key_press_fcn(hObject,eventdata,handles)</a></li><li><a href="#_sub41" class="code">function append_point_Callback(hObject, eventdata, handles)</a></li></ul> |
---|
63 | <h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
64 | <div class="fragment"><pre>0001 <span class="comment">%'geometry_calib': performs geometric calibration from a set of reference points</span> |
---|
65 | 0002 <span class="comment">%</span> |
---|
66 | 0003 <span class="comment">% function varargout = geometry_calib(varargin)</span> |
---|
67 | 0004 <span class="comment">%</span> |
---|
68 | 0005 <span class="comment">%A%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
69 | 0006 <span class="comment">% Copyright Joel Sommeria, 2008, LEGI / CNRS-UJF-INPG, sommeria@coriolis-legi.org.</span> |
---|
70 | 0007 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
71 | 0008 <span class="comment">% This file is part of the toolbox UVMAT.</span> |
---|
72 | 0009 <span class="comment">%</span> |
---|
73 | 0010 <span class="comment">% UVMAT is free software; you can redistribute it and/or modify</span> |
---|
74 | 0011 <span class="comment">% it under the terms of the GNU General Public License as published by</span> |
---|
75 | 0012 <span class="comment">% the Free Software Foundation; either version 2 of the License, or</span> |
---|
76 | 0013 <span class="comment">% (at your option) any later version.</span> |
---|
77 | 0014 <span class="comment">%</span> |
---|
78 | 0015 <span class="comment">% UVMAT is distributed in the hope that it will be useful,</span> |
---|
79 | 0016 <span class="comment">% but WITHOUT ANY WARRANTY; without even the implied warranty of</span> |
---|
80 | 0017 <span class="comment">% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the</span> |
---|
81 | 0018 <span class="comment">% GNU General Public License (file UVMAT/COPYING.txt) for more details.</span> |
---|
82 | 0019 <span class="comment">%AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA</span> |
---|
83 | 0020 |
---|
84 | 0021 <a name="_sub0" href="#_subfunctions" class="code">function varargout = geometry_calib(varargin)</a> |
---|
85 | 0022 <span class="comment">% GEOMETRY_CALIB M-file for geometry_calib.fig</span> |
---|
86 | 0023 <span class="comment">% GEOMETRY_CALIB, by itself, creates a MenuCoord GEOMETRY_CALIB or raises the existing</span> |
---|
87 | 0024 <span class="comment">% singleton*.</span> |
---|
88 | 0025 <span class="comment">%</span> |
---|
89 | 0026 <span class="comment">% H = GEOMETRY_CALIB returns the handle to a MenuCoord GEOMETRY_CALIB or the handle to</span> |
---|
90 | 0027 <span class="comment">% the existing singleton*.</span> |
---|
91 | 0028 <span class="comment">%</span> |
---|
92 | 0029 <span class="comment">% GEOMETRY_CALIB('CALLBACK',hObject,eventData,handles,...) calls the local</span> |
---|
93 | 0030 <span class="comment">% function named CALLBACK in GEOMETRY_CALIB.M with the given input arguments.</span> |
---|
94 | 0031 <span class="comment">%</span> |
---|
95 | 0032 <span class="comment">% GEOMETRY_CALIB('Property','Value',...) creates a MenuCoord GEOMETRY_CALIB or raises the</span> |
---|
96 | 0033 <span class="comment">% existing singleton*. Starting from the left, property value pairs are</span> |
---|
97 | 0034 <span class="comment">% applied to the GUI before geometry_calib_OpeningFunction gets called. An</span> |
---|
98 | 0035 <span class="comment">% unrecognized property name or invalid value makes property application</span> |
---|
99 | 0036 <span class="comment">% stop. All inputs are passed to geometry_calib_OpeningFcn via varargin.</span> |
---|
100 | 0037 <span class="comment">%</span> |
---|
101 | 0038 <span class="comment">% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one</span> |
---|
102 | 0039 <span class="comment">% instance to run (singleton)".</span> |
---|
103 | 0040 <span class="comment">%</span> |
---|
104 | 0041 <span class="comment">% See also: GUIDE, GUIDATA, GUIHANDLES</span> |
---|
105 | 0042 |
---|
106 | 0043 <span class="comment">% Edit the above text to modify the response to help geometry_calib</span> |
---|
107 | 0044 |
---|
108 | 0045 <span class="comment">% Last Modified by GUIDE v2.5 07-Nov-2009 21:53:58</span> |
---|
109 | 0046 |
---|
110 | 0047 <span class="comment">% Begin initialization code - DO NOT edit</span> |
---|
111 | 0048 gui_Singleton = 1; |
---|
112 | 0049 gui_State = struct(<span class="string">'gui_Name'</span>, mfilename, <span class="keyword">...</span> |
---|
113 | 0050 <span class="string">'gui_Singleton'</span>, gui_Singleton, <span class="keyword">...</span> |
---|
114 | 0051 <span class="string">'gui_OpeningFcn'</span>, @<a href="#_sub1" class="code" title="subfunction geometry_calib_OpeningFcn(hObject, eventdata, handles, handles_uvmat,pos,inputfile)">geometry_calib_OpeningFcn</a>, <span class="keyword">...</span> |
---|
115 | 0052 <span class="string">'gui_OutputFcn'</span>, @<a href="#_sub2" class="code" title="subfunction varargout = geometry_calib_OutputFcn(hObject, eventdata, handles)">geometry_calib_OutputFcn</a>, <span class="keyword">...</span> |
---|
116 | 0053 <span class="string">'gui_LayoutFcn'</span>, [] , <span class="keyword">...</span> |
---|
117 | 0054 <span class="string">'gui_Callback'</span>, []); |
---|
118 | 0055 <span class="keyword">if</span> nargin & isstr(varargin{1}) |
---|
119 | 0056 gui_State.gui_Callback = str2func(varargin{1}); |
---|
120 | 0057 <span class="keyword">end</span> |
---|
121 | 0058 |
---|
122 | 0059 <span class="keyword">if</span> nargout |
---|
123 | 0060 [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); |
---|
124 | 0061 <span class="keyword">else</span> |
---|
125 | 0062 gui_mainfcn(gui_State, varargin{:}); |
---|
126 | 0063 <span class="keyword">end</span> |
---|
127 | 0064 <span class="comment">% End initialization code - DO NOT edit</span> |
---|
128 | 0065 |
---|
129 | 0066 |
---|
130 | 0067 <span class="comment">% --- Executes just before geometry_calib is made visible.</span> |
---|
131 | 0068 <span class="comment">%INPUT:</span> |
---|
132 | 0069 <span class="comment">%handles: handles of the geometry_calib interface elements</span> |
---|
133 | 0070 <span class="comment">% PlotHandles: set of handles of the elements contolling the plotting</span> |
---|
134 | 0071 <span class="comment">% parameters on the uvmat interface (obtained by 'get_plot_handle.m')</span> |
---|
135 | 0072 <a name="_sub1" href="#_subfunctions" class="code">function geometry_calib_OpeningFcn(hObject, eventdata, handles, handles_uvmat,pos,inputfile)</a> |
---|
136 | 0073 set(handles.Phi,<span class="string">'TooltipString'</span>,<span class="string">'Phi: rotation angle of the physical point coordiantes (in degrees)'</span>)<span class="comment">% TO PUT IN GUIDE</span> |
---|
137 | 0074 <span class="comment">% Choose default command line output for geometry_calib</span> |
---|
138 | 0075 handles.output = hObject; |
---|
139 | 0076 |
---|
140 | 0077 <span class="comment">% Update handles structure</span> |
---|
141 | 0078 guidata(hObject, handles); |
---|
142 | 0079 movegui(hObject,<span class="string">'east'</span>);<span class="comment">% position the GUI ton the right of the screen</span> |
---|
143 | 0080 <span class="keyword">if</span> exist(<span class="string">'handles_uvmat'</span>,<span class="string">'var'</span>) <span class="comment">%& isfield(data,'ParentButton')</span> |
---|
144 | 0081 set(hObject,<span class="string">'DeleteFcn'</span>,{@<a href="#_sub6" class="code" title="subfunction closefcn(gcbo,eventdata,handles_uvmat)">closefcn</a>,handles_uvmat})<span class="comment">%</span> |
---|
145 | 0082 <span class="keyword">end</span> |
---|
146 | 0083 <span class="comment">%set the position of the interface</span> |
---|
147 | 0084 <span class="keyword">if</span> exist(<span class="string">'pos'</span>,<span class="string">'var'</span>)& length(pos)>2 |
---|
148 | 0085 pos_gui=get(hObject,<span class="string">'Position'</span>); |
---|
149 | 0086 pos_gui(1)=pos(1); |
---|
150 | 0087 pos_gui(2)=pos(2); |
---|
151 | 0088 set(hObject,<span class="string">'Position'</span>,pos_gui); |
---|
152 | 0089 <span class="keyword">end</span> |
---|
153 | 0090 inputxml=<span class="string">''</span>; |
---|
154 | 0091 <span class="keyword">if</span> exist(<span class="string">'inputfile'</span>,<span class="string">'var'</span>)& ~isempty(inputfile) |
---|
155 | 0092 [Path,Name,ext]=fileparts(inputfile); |
---|
156 | 0093 form=imformats(ext([2:end])); |
---|
157 | 0094 <span class="keyword">if</span> ~isempty(form)<span class="comment">% if the input file is an image</span> |
---|
158 | 0095 set(hObject,<span class="string">'UserData'</span>,inputfile) |
---|
159 | 0096 [Pathsub,RootFile,field_count,str2,str_a,str_b,ext,nom_type,subdir]=<a href="name2display.html" class="code" title="function [RootPath,RootFile,field_count,str2,str_a,str_b,ext,nom_type,subdir]=name2display(fileinput)">name2display</a>(inputfile); |
---|
160 | 0097 inputxml=[fullfile(Pathsub,RootFile) <span class="string">'.xml'</span>]; |
---|
161 | 0098 <span class="keyword">end</span> |
---|
162 | 0099 <span class="keyword">end</span> |
---|
163 | 0100 set(handles.ListCoord,<span class="string">'String'</span>,{<span class="string">''</span>}) |
---|
164 | 0101 <span class="keyword">if</span> exist(inputxml,<span class="string">'file'</span>) |
---|
165 | 0102 <a href="#_sub5" class="code" title="subfunction loadfile(handles,fileinput)">loadfile</a>(handles,inputxml)<span class="comment">% load the point coordiantes existing in the xml file</span> |
---|
166 | 0103 <span class="keyword">end</span> |
---|
167 | 0104 |
---|
168 | 0105 set(handles.ListCoord,<span class="string">'KeyPressFcn'</span>,{@<a href="#_sub40" class="code" title="subfunction key_press_fcn(hObject,eventdata,handles)">key_press_fcn</a>,handles})<span class="comment">%set keyboard action function</span> |
---|
169 | 0106 <span class="comment">%set(hObject,'KeyPressFcn',{'keyboard_callback',handles})%set keyboard action function on uvmat interface when geometry_calib is on top</span> |
---|
170 | 0107 <span class="comment">%htable=uitable(10,5)</span> |
---|
171 | 0108 <span class="comment">%set(htable,'ColumnNames',{'x','y','z','X(pixels)','Y(pixels)'})</span> |
---|
172 | 0109 |
---|
173 | 0110 <span class="comment">% --- Outputs from this function are returned to the command line.</span> |
---|
174 | 0111 <a name="_sub2" href="#_subfunctions" class="code">function varargout = geometry_calib_OutputFcn(hObject, eventdata, handles)</a> |
---|
175 | 0112 <span class="comment">% varargout cell array for returning output args (see VARARGOUT);</span> |
---|
176 | 0113 <span class="comment">% hObject handle to figure</span> |
---|
177 | 0114 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
178 | 0115 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
179 | 0116 |
---|
180 | 0117 <span class="comment">% Get default command line output from handles structure</span> |
---|
181 | 0118 varargout{1} = handles.output; |
---|
182 | 0119 varargout{2}=handles; |
---|
183 | 0120 |
---|
184 | 0121 <span class="comment">%------------</span> |
---|
185 | 0122 <a name="_sub3" href="#_subfunctions" class="code">function Phi_Callback(hObject, eventdata, handles)</a> |
---|
186 | 0123 |
---|
187 | 0124 |
---|
188 | 0125 <span class="comment">%-----------------------------------------------------</span> |
---|
189 | 0126 <span class="comment">% --- Executes on button press in OPEN.</span> |
---|
190 | 0127 <a name="_sub4" href="#_subfunctions" class="code">function OPEN_Callback(hObject, eventdata, handles)</a> |
---|
191 | 0128 <span class="comment">%get the object file</span> |
---|
192 | 0129 huvmat=findobj(allchild(0),<span class="string">'Name'</span>,<span class="string">'uvmat'</span>); |
---|
193 | 0130 UvData=get(huvmat,<span class="string">'UserData'</span>); |
---|
194 | 0131 hchild=get(huvmat,<span class="string">'Children'</span>); |
---|
195 | 0132 hrootpath=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'RootPath'</span>); |
---|
196 | 0133 oldfile=get(hrootpath,<span class="string">'String'</span>); |
---|
197 | 0134 <span class="keyword">if</span> isempty(oldfile) |
---|
198 | 0135 oldfile=<span class="string">''</span>; |
---|
199 | 0136 <span class="keyword">end</span> |
---|
200 | 0137 <span class="comment">%[FileName,PathName] = uigetfile('*.civ','Select a .civ file',oldfile)</span> |
---|
201 | 0138 [FileName, PathName, filterindex] = uigetfile( <span class="keyword">...</span> |
---|
202 | 0139 {<span class="string">'*.xml;*.mat'</span>, <span class="string">' (*.xml,*.mat)'</span>; |
---|
203 | 0140 <span class="string">'*.xml'</span>, <span class="string">'.xml files '</span>; <span class="keyword">...</span> |
---|
204 | 0141 <span class="string">'*.mat'</span>, <span class="string">'.mat matlab files '</span>}, <span class="keyword">...</span> |
---|
205 | 0142 <span class="string">'Pick a file'</span>,oldfile); |
---|
206 | 0143 fileinput=[PathName FileName];<span class="comment">%complete file name</span> |
---|
207 | 0144 testblank=findstr(fileinput,<span class="string">' '</span>);<span class="comment">%look for blanks</span> |
---|
208 | 0145 <span class="keyword">if</span> ~isempty(testblank) |
---|
209 | 0146 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'forbidden input file name or path: no blank character allowed'</span>,<span class="string">'ERROR'</span>) |
---|
210 | 0147 <span class="keyword">return</span> |
---|
211 | 0148 <span class="keyword">end</span> |
---|
212 | 0149 sizf=size(fileinput); |
---|
213 | 0150 <span class="keyword">if</span> (~ischar(fileinput)|~isequal(sizf(1),1)),<span class="keyword">return</span>;<span class="keyword">end</span> |
---|
214 | 0151 <a href="#_sub5" class="code" title="subfunction loadfile(handles,fileinput)">loadfile</a>(handles,fileinput) |
---|
215 | 0152 |
---|
216 | 0153 <span class="comment">%--------------------------------------------------</span> |
---|
217 | 0154 <span class="comment">%read input xml file and update the edit boxes</span> |
---|
218 | 0155 <a name="_sub5" href="#_subfunctions" class="code">function loadfile(handles,fileinput)</a> |
---|
219 | 0156 |
---|
220 | 0157 <span class="comment">%read the input xml file</span> |
---|
221 | 0158 t=xmltree(fileinput); |
---|
222 | 0159 s=convert(t);<span class="comment">%convert to matlab structure</span> |
---|
223 | 0160 <span class="comment">%read data currently displayed on the interface</span> |
---|
224 | 0161 PointCoord=[]; |
---|
225 | 0162 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
226 | 0163 data=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
227 | 0164 <span class="comment">%data=read_geometry_calib(handles);</span> |
---|
228 | 0165 Coord=[]; <span class="comment">%default</span> |
---|
229 | 0166 <span class="keyword">if</span> isfield(data,<span class="string">'Coord'</span>) |
---|
230 | 0167 Coord=data.Coord; |
---|
231 | 0168 <span class="keyword">end</span> |
---|
232 | 0169 TabChar_0=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
233 | 0170 nbcoord_0=size(TabChar_0,1); |
---|
234 | 0171 <span class="keyword">if</span> isequal(get(handles.edit_append,<span class="string">'Value'</span>),2) <span class="comment">%edit mode A REVOIR</span> |
---|
235 | 0172 val=get(handles.ListCoord,<span class="string">'Value'</span>)-1; |
---|
236 | 0173 <span class="keyword">else</span> |
---|
237 | 0174 val=length(TabChar_0); |
---|
238 | 0175 <span class="keyword">end</span> |
---|
239 | 0176 nbcoord=0; |
---|
240 | 0177 |
---|
241 | 0178 <span class="comment">%case of calibration (ImaDoc) input file</span> |
---|
242 | 0179 <span class="keyword">if</span> isfield(s,<span class="string">'Heading'</span>) |
---|
243 | 0180 CalibData=s.Heading; |
---|
244 | 0181 <span class="keyword">end</span> |
---|
245 | 0182 CalibData.XmlInput=fileinput; |
---|
246 | 0183 hcalib=get(handles.OPEN,<span class="string">'parent'</span>); |
---|
247 | 0184 set(hcalib,<span class="string">'UserData'</span>,CalibData);<span class="comment">%store the heading in the interface 'UserData'</span> |
---|
248 | 0185 <span class="keyword">if</span> isfield(s,<span class="string">'GeometryCalib'</span>) |
---|
249 | 0186 Calib=s.GeometryCalib; |
---|
250 | 0187 <span class="keyword">if</span> isfield(Calib,<span class="string">'CalibrationType'</span>) |
---|
251 | 0188 CalibrationType=Calib.CalibrationType; |
---|
252 | 0189 <span class="keyword">switch</span> CalibrationType |
---|
253 | 0190 <span class="keyword">case</span> <span class="string">'linear'</span> |
---|
254 | 0191 set(handles.calib_type,<span class="string">'Value'</span>,2) |
---|
255 | 0192 <span class="keyword">case</span> <span class="string">'tsai'</span> |
---|
256 | 0193 set(handles.calib_type,<span class="string">'Value'</span>,3) |
---|
257 | 0194 <span class="keyword">end</span> |
---|
258 | 0195 <span class="keyword">end</span> |
---|
259 | 0196 <span class="keyword">if</span> isfield(Calib,<span class="string">'SourceCalib'</span>) |
---|
260 | 0197 <span class="keyword">if</span> isfield(Calib.SourceCalib,<span class="string">'PointCoord'</span>) |
---|
261 | 0198 PointCoord=Calib.SourceCalib.PointCoord; |
---|
262 | 0199 <span class="keyword">end</span> |
---|
263 | 0200 <span class="comment">% if isfield(Calib.SourceCalib,'ImageCalib')</span> |
---|
264 | 0201 <span class="comment">% hcalib=get(handles.OPEN,'parent');</span> |
---|
265 | 0202 <span class="comment">% set(hcalib,'UserData',Calib.SourceCalib.ImageCalib);%store the source image name in the interface 'UserData'</span> |
---|
266 | 0203 <span class="comment">% end</span> |
---|
267 | 0204 <span class="keyword">end</span> |
---|
268 | 0205 nbcoord=length(PointCoord); |
---|
269 | 0206 <span class="keyword">if</span> ~isfield(Calib,<span class="string">'ErrorRms'</span>)&~isfield(Calib,<span class="string">'ErrorMax'</span>) <span class="comment">%old convention of Gauthier (cord in mm)</span> |
---|
270 | 0207 <span class="keyword">for</span> i=1:length(PointCoord) |
---|
271 | 0208 line=str2num(PointCoord{i}); |
---|
272 | 0209 Coord(i+val,4:5)=line(4:5);<span class="comment">%px x</span> |
---|
273 | 0210 Coord(i+val,1:3)=line(1:3)/10;<span class="comment">%phys x</span> |
---|
274 | 0211 <span class="keyword">end</span> |
---|
275 | 0212 <span class="keyword">else</span> |
---|
276 | 0213 <span class="keyword">for</span> i=1:length(PointCoord) |
---|
277 | 0214 line=str2num(PointCoord{i}); |
---|
278 | 0215 Coord(i,4:5)=line(4:5);<span class="comment">%px x</span> |
---|
279 | 0216 Coord(i,1:3)=line(1:3);<span class="comment">%phys x</span> |
---|
280 | 0217 <span class="keyword">end</span> |
---|
281 | 0218 <span class="keyword">end</span> |
---|
282 | 0219 <span class="keyword">end</span> |
---|
283 | 0220 <span class="comment">%case of xml files of points</span> |
---|
284 | 0221 <span class="keyword">if</span> isfield(s,<span class="string">'Coord'</span>) |
---|
285 | 0222 PointCoord=s.Coord; |
---|
286 | 0223 nbcoord=length(PointCoord); |
---|
287 | 0224 <span class="comment">%case of image coordinates</span> |
---|
288 | 0225 <span class="keyword">if</span> isfield(s,<span class="string">'CoordType'</span>)& isequal(s.CoordType,<span class="string">'px'</span>) |
---|
289 | 0226 <span class="keyword">for</span> i=1:nbcoord |
---|
290 | 0227 line=str2num(PointCoord{i}); |
---|
291 | 0228 Coord(i+val,4:5)=line(1:2); |
---|
292 | 0229 <span class="keyword">end</span> |
---|
293 | 0230 <span class="comment">%case of physical coordinates</span> |
---|
294 | 0231 <span class="keyword">else</span> |
---|
295 | 0232 <span class="keyword">for</span> i=1:nbcoord |
---|
296 | 0233 line=str2num(PointCoord{i}); |
---|
297 | 0234 Coord(i+val,1:3)=line(1:3); |
---|
298 | 0235 nbcolumn=size(Coord,2); |
---|
299 | 0236 <span class="keyword">if</span> nbcolumn<5 |
---|
300 | 0237 Coord(i+val,nbcolumn+1:5)=zeros(1,5-nbcolumn); |
---|
301 | 0238 <span class="keyword">end</span> |
---|
302 | 0239 <span class="keyword">end</span> |
---|
303 | 0240 <span class="keyword">end</span> |
---|
304 | 0241 <span class="keyword">end</span> |
---|
305 | 0242 CoordCell={}; |
---|
306 | 0243 <span class="keyword">for</span> iline=1:size(Coord,1) |
---|
307 | 0244 <span class="keyword">for</span> j=1:5 |
---|
308 | 0245 CoordCell{iline,j}=num2str(Coord(iline,j),4); |
---|
309 | 0246 <span class="keyword">end</span> |
---|
310 | 0247 <span class="keyword">end</span> |
---|
311 | 0248 Tabchar=<a href="cell2tab.html" class="code" title="function Tabchar=cell2tab(Tabcell,separator)">cell2tab</a>(CoordCell,<span class="string">' | '</span>);<span class="comment">%transform cells into table ready for display</span> |
---|
312 | 0249 set(handles.ListCoord,<span class="string">'Value'</span>,1) |
---|
313 | 0250 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
314 | 0251 |
---|
315 | 0252 |
---|
316 | 0253 <span class="comment">%----------------------------------------------------</span> |
---|
317 | 0254 <span class="comment">% executed when closing: set the parent interface button to value 0</span> |
---|
318 | 0255 <a name="_sub6" href="#_subfunctions" class="code">function closefcn(gcbo,eventdata,handles_uvmat)</a> |
---|
319 | 0256 huvmat=findobj(allchild(0),<span class="string">'Name'</span>,<span class="string">'uvmat'</span>); |
---|
320 | 0257 <span class="keyword">if</span> exist(<span class="string">'handles_uvmat'</span>,<span class="string">'var'</span>) |
---|
321 | 0258 set(handles_uvmat.cal,<span class="string">'Value'</span>,0) |
---|
322 | 0259 <a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a>(<span class="string">'cal_Callback'</span>,huvmat,[],handles_uvmat); |
---|
323 | 0260 <span class="comment">% set(parent_button,'Value',0)%put unactivated buttons to green</span> |
---|
324 | 0261 <span class="comment">% set(parent_button,'BackgroundColor',[0 1 0]);</span> |
---|
325 | 0262 <span class="keyword">end</span> |
---|
326 | 0263 |
---|
327 | 0264 <span class="comment">%-----------------------------------------------------------------------</span> |
---|
328 | 0265 <span class="comment">% --- Executes on button press in PLOT (plot): PLOT the defined object and its projected field</span> |
---|
329 | 0266 <a name="_sub7" href="#_subfunctions" class="code">function PLOT_Callback(hObject, eventdata, handles)</a> |
---|
330 | 0267 |
---|
331 | 0268 huvmat=findobj(allchild(0),<span class="string">'Name'</span>,<span class="string">'uvmat'</span>);<span class="comment">%find the current uvmat interface handle</span> |
---|
332 | 0269 UvData=get(huvmat,<span class="string">'UserData'</span>);<span class="comment">%Data associated to the current uvmat interface</span> |
---|
333 | 0270 hhuvmat=guidata(huvmat); <span class="comment">%handles of GUI elements in uvmat</span> |
---|
334 | 0271 hplot=findobj(huvmat,<span class="string">'Tag'</span>,<span class="string">'axes3'</span>);<span class="comment">%main plotting axis of uvmat</span> |
---|
335 | 0272 h_menu_coord=findobj(huvmat,<span class="string">'Tag'</span>,<span class="string">'menu_coord'</span>); |
---|
336 | 0273 menu=get(h_menu_coord,<span class="string">'String'</span>); |
---|
337 | 0274 choice=get(h_menu_coord,<span class="string">'Value'</span>); |
---|
338 | 0275 <span class="keyword">if</span> iscell(menu) |
---|
339 | 0276 option=menu{choice}; |
---|
340 | 0277 <span class="keyword">else</span> |
---|
341 | 0278 option=<span class="string">'px'</span>; <span class="comment">%default</span> |
---|
342 | 0279 <span class="keyword">end</span> |
---|
343 | 0280 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
344 | 0281 ObjectData=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
345 | 0282 <span class="comment">%ObjectData=read_geometry_calib(handles);%read the interface input parameters defining the object</span> |
---|
346 | 0283 <span class="keyword">if</span> isequal(option,<span class="string">'phys'</span>) |
---|
347 | 0284 ObjectData.Coord=ObjectData.Coord(:,[1:3]); |
---|
348 | 0285 <span class="keyword">elseif</span> isequal(option,<span class="string">'px'</span>)||isequal(option,<span class="string">''</span>) |
---|
349 | 0286 ObjectData.Coord=ObjectData.Coord(:,[4:5]); |
---|
350 | 0287 <span class="keyword">else</span> |
---|
351 | 0288 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'the choice in menu_coord of uvmat must be px or phys '</span>,<span class="string">'ERROR'</span>) |
---|
352 | 0289 <span class="keyword">end</span> |
---|
353 | 0290 axes(hhuvmat.axes3) |
---|
354 | 0291 hh=findobj(<span class="string">'Tag'</span>,<span class="string">'calib_points'</span>); |
---|
355 | 0292 <span class="keyword">if</span> isempty(hh) |
---|
356 | 0293 hh=line(ObjectData.Coord(:,1),ObjectData.Coord(:,2),<span class="string">'Color'</span>,<span class="string">'m'</span>,<span class="string">'Tag'</span>,<span class="string">'calib_points'</span>,<span class="string">'LineStyle'</span>,<span class="string">'.'</span>,<span class="string">'Marker'</span>,<span class="string">'+'</span>); |
---|
357 | 0294 <span class="keyword">else</span> |
---|
358 | 0295 set(hh,<span class="string">'XData'</span>,ObjectData.Coord(:,1)) |
---|
359 | 0296 set(hh,<span class="string">'YData'</span>,ObjectData.Coord(:,2)) |
---|
360 | 0297 <span class="keyword">end</span> |
---|
361 | 0298 |
---|
362 | 0299 <span class="comment">% --- Executes on button press in MenuCoord.</span> |
---|
363 | 0300 <a name="_sub8" href="#_subfunctions" class="code">function MenuCoord_Callback(hObject, eventdata, handles)</a> |
---|
364 | 0301 |
---|
365 | 0302 |
---|
366 | 0303 <span class="comment">% --- Executes on button press in delete.</span> |
---|
367 | 0304 <a name="_sub9" href="#_subfunctions" class="code">function delete_Callback(hObject, eventdata, handles)</a> |
---|
368 | 0305 SetData=get(gcbf,<span class="string">'UserData'</span>);<span class="comment">%get the interface data</span> |
---|
369 | 0306 IndexObj=SetData.IndexObj; |
---|
370 | 0307 <a href="delete_object.html" class="code" title="function delete_object(hObject)">delete_object</a>(IndexObj); |
---|
371 | 0308 |
---|
372 | 0309 |
---|
373 | 0310 <span class="comment">%------------------------------------------------------------------</span> |
---|
374 | 0311 <span class="comment">% --- Executes on button press in calibrate_lin.</span> |
---|
375 | 0312 <a name="_sub10" href="#_subfunctions" class="code">function APPLY_Callback(hObject, eventdata, handles)</a> |
---|
376 | 0313 <span class="comment">%------------------------------------------------------------------</span> |
---|
377 | 0314 calib_cell=get(handles.calib_type,<span class="string">'String'</span>); |
---|
378 | 0315 val=get(handles.calib_type,<span class="string">'Value'</span>); |
---|
379 | 0316 calib_type=calib_cell{val}; |
---|
380 | 0317 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
381 | 0318 Object=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
382 | 0319 |
---|
383 | 0320 <span class="keyword">if</span> isequal(calib_type,<span class="string">'rescale'</span>) |
---|
384 | 0321 GeometryCalib=<a href="#_sub13" class="code" title="subfunction GeometryCalib=calib_rescale(Coord)">calib_rescale</a>(Object.Coord); |
---|
385 | 0322 <span class="keyword">elseif</span> isequal(calib_type,<span class="string">'linear'</span>) |
---|
386 | 0323 GeometryCalib=<a href="#_sub14" class="code" title="subfunction GeometryCalib=calib_linear(Coord)">calib_linear</a>(Object.Coord); |
---|
387 | 0324 <span class="keyword">elseif</span> isequal(calib_type,<span class="string">'tsai'</span>) |
---|
388 | 0325 GeometryCalib=<a href="#_sub15" class="code" title="subfunction GeometryCalib=calib_tsai(Coord)">calib_tsai</a>(Object.Coord); |
---|
389 | 0326 <span class="keyword">end</span> |
---|
390 | 0327 unitlist=get(handles.CoordUnit,<span class="string">'String'</span>); |
---|
391 | 0328 unit=unitlist{get(handles.CoordUnit,<span class="string">'value'</span>)}; |
---|
392 | 0329 GeometryCalib.CoordUnit=unit; |
---|
393 | 0330 |
---|
394 | 0331 huvmat=findobj(allchild(0),<span class="string">'Name'</span>,<span class="string">'uvmat'</span>); |
---|
395 | 0332 hhuvmat=guidata(huvmat);<span class="comment">%handles of elements in the GUI uvmat</span> |
---|
396 | 0333 RootPath=<span class="string">''</span>; |
---|
397 | 0334 RootFile=<span class="string">''</span>; |
---|
398 | 0335 <span class="keyword">if</span> ~isempty(hhuvmat.RootPath)& ~isempty(hhuvmat.RootFile) |
---|
399 | 0336 testhandle=1; |
---|
400 | 0337 RootPath=get(hhuvmat.RootPath,<span class="string">'String'</span>); |
---|
401 | 0338 RootFile=get(hhuvmat.RootFile,<span class="string">'String'</span>); |
---|
402 | 0339 filebase=fullfile(RootPath,RootFile); |
---|
403 | 0340 outputfile=[filebase <span class="string">'.xml'</span>]; |
---|
404 | 0341 <span class="keyword">else</span> |
---|
405 | 0342 question={<span class="string">'save the calibration data and point coordinates in'</span>}; |
---|
406 | 0343 def={fullfile(RootPath,[<span class="string">'ObjectCalib.xml'</span>])}; |
---|
407 | 0344 options.Resize=<span class="string">'on'</span>; |
---|
408 | 0345 answer=inputdlg(question,<span class="string">'save average in a new file'</span>,1,def,options); |
---|
409 | 0346 outputfile=answer{1}; |
---|
410 | 0347 <span class="keyword">end</span> |
---|
411 | 0348 testappend=0; |
---|
412 | 0349 <span class="keyword">if</span> exist(outputfile,<span class="string">'file'</span>);<span class="comment">%=1 if the output file already exists, 0 else</span> |
---|
413 | 0350 t=xmltree(outputfile); <span class="comment">%read the file</span> |
---|
414 | 0351 uid=find(t,<span class="string">'ImaDoc'</span>); |
---|
415 | 0352 <span class="keyword">if</span> ~isequal(uid,1)<span class="comment">%if the xml file is not ImaDoc, delete it (after backup)</span> |
---|
416 | 0353 backupfile=outputfile; |
---|
417 | 0354 testexist=2; |
---|
418 | 0355 <span class="keyword">while</span> testexist==2 |
---|
419 | 0356 backupfile=[backupfile <span class="string">'~'</span>]; |
---|
420 | 0357 testexist=exist(backupfile,<span class="string">'file'</span>); |
---|
421 | 0358 <span class="keyword">end</span> |
---|
422 | 0359 [success,message]=copyfile(outputfile,backupfile);<span class="comment">%make backup</span> |
---|
423 | 0360 <span class="keyword">if</span> isequal(success,1) |
---|
424 | 0361 delete(outputfile) |
---|
425 | 0362 <span class="keyword">else</span> |
---|
426 | 0363 <span class="keyword">return</span> |
---|
427 | 0364 <span class="keyword">end</span> |
---|
428 | 0365 <span class="keyword">else</span> |
---|
429 | 0366 uid_calib=find(t,<span class="string">'ImaDoc/GeometryCalib'</span>); |
---|
430 | 0367 <span class="keyword">if</span> ~isempty(uid) <span class="comment">%if GeometryCalib already exists, delete its content</span> |
---|
431 | 0368 backupfile=outputfile; |
---|
432 | 0369 testexist=2; |
---|
433 | 0370 <span class="keyword">while</span> testexist==2 |
---|
434 | 0371 backupfile=[backupfile <span class="string">'~'</span>]; |
---|
435 | 0372 testexist=exist(backupfile,<span class="string">'file'</span>); |
---|
436 | 0373 <span class="keyword">end</span> |
---|
437 | 0374 [success,message]=copyfile(outputfile,backupfile);<span class="comment">%make backup</span> |
---|
438 | 0375 <span class="keyword">if</span> isequal(success,1) |
---|
439 | 0376 delete(outputfile) |
---|
440 | 0377 <span class="keyword">else</span> |
---|
441 | 0378 <span class="keyword">return</span> |
---|
442 | 0379 <span class="keyword">end</span> |
---|
443 | 0380 uid_child=children(t,uid_calib); |
---|
444 | 0381 t=delete(t,uid_child); |
---|
445 | 0382 testappend=1; |
---|
446 | 0383 <span class="keyword">end</span> |
---|
447 | 0384 <span class="keyword">end</span> |
---|
448 | 0385 <span class="keyword">end</span> |
---|
449 | 0386 <span class="keyword">if</span> ~testappend |
---|
450 | 0387 t=xmltree; |
---|
451 | 0388 t=set(t,1,<span class="string">'name'</span>,<span class="string">'ImaDoc'</span>); |
---|
452 | 0389 [t,uid_calib]=add(t,1,<span class="string">'element'</span>,<span class="string">'GeometryCalib'</span>); |
---|
453 | 0390 <span class="comment">% t=struct2xml(GeometryCalib,t,uid_calib);</span> |
---|
454 | 0391 <span class="keyword">end</span> |
---|
455 | 0392 <span class="comment">% hgrid=get(handles.REPLICATE,'parent');%read the calibration image source on the interface userdata</span> |
---|
456 | 0393 <span class="comment">% imagename=get(hgrid,'UserData');</span> |
---|
457 | 0394 <span class="comment">% if exist(imagename,'file')</span> |
---|
458 | 0395 <span class="comment">% GeometryCalib.SourceCalib.ImageCalib=imagename;</span> |
---|
459 | 0396 <span class="comment">% end</span> |
---|
460 | 0397 GeometryCalib.SourceCalib.PointCoord=Object.Coord; |
---|
461 | 0398 t=<a href="struct2xml.html" class="code" title="function t=struct2xml(Object,t,root_uid)">struct2xml</a>(GeometryCalib,t,uid_calib); |
---|
462 | 0399 save(t,outputfile); |
---|
463 | 0400 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>({[outputfile <span class="string">' updated with calibration data'</span>];<span class="keyword">...</span> |
---|
464 | 0401 [<span class="string">'Error rms (along x,y)='</span> num2str(GeometryCalib.ErrorRms) <span class="string">' pixels'</span>];<span class="keyword">...</span> |
---|
465 | 0402 [<span class="string">'Error max (along x,y)='</span> num2str(GeometryCalib.ErrorMax) <span class="string">' pixels'</span>]},<span class="string">'CONFIRMATION'</span>) |
---|
466 | 0403 |
---|
467 | 0404 <span class="comment">%display image with new calibration in the currently opened uvmat interface</span> |
---|
468 | 0405 hhh=findobj(hhuvmat.axes3,<span class="string">'Tag'</span>,<span class="string">'calib_marker'</span>);<span class="comment">% delete calib points and markers</span> |
---|
469 | 0406 <span class="keyword">if</span> ~isempty(hhh) |
---|
470 | 0407 delete(hhh); |
---|
471 | 0408 <span class="keyword">end</span> |
---|
472 | 0409 hhh=findobj(hhuvmat.axes3,<span class="string">'Tag'</span>,<span class="string">'calib_points'</span>); |
---|
473 | 0410 <span class="keyword">if</span> ~isempty(hhh) |
---|
474 | 0411 delete(hhh); |
---|
475 | 0412 <span class="keyword">end</span> |
---|
476 | 0413 set(hhuvmat.FixedLimits,<span class="string">'Value'</span>,0)<span class="comment">% put FixedLimits option to 'off'</span> |
---|
477 | 0414 set(hhuvmat.FixedLimits,<span class="string">'BackgroundColor'</span>,[0.7 0.7 0.7]) |
---|
478 | 0415 <a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a>(<span class="string">'RootPath_Callback'</span>,hObject,eventdata,hhuvmat); <span class="comment">%file input with xml reading in uvmat</span> |
---|
479 | 0416 |
---|
480 | 0417 <span class="comment">%------------------------------------------------------------------</span> |
---|
481 | 0418 <span class="comment">% --- Executes on button press in calibrate_lin.</span> |
---|
482 | 0419 <a name="_sub11" href="#_subfunctions" class="code">function REPLICATE_Callback(hObject, eventdata, handles)</a> |
---|
483 | 0420 <span class="comment">%------------------------------------------------------------------</span> |
---|
484 | 0421 calib_cell=get(handles.calib_type,<span class="string">'String'</span>); |
---|
485 | 0422 val=get(handles.calib_type,<span class="string">'Value'</span>); |
---|
486 | 0423 calib_type=calib_cell{val}; |
---|
487 | 0424 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
488 | 0425 Object=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
489 | 0426 |
---|
490 | 0427 <span class="keyword">if</span> isequal(calib_type,<span class="string">'rescale'</span>) |
---|
491 | 0428 GeometryCalib=<a href="#_sub13" class="code" title="subfunction GeometryCalib=calib_rescale(Coord)">calib_rescale</a>(Object.Coord); |
---|
492 | 0429 <span class="keyword">elseif</span> isequal(calib_type,<span class="string">'linear'</span>) |
---|
493 | 0430 GeometryCalib=<a href="#_sub14" class="code" title="subfunction GeometryCalib=calib_linear(Coord)">calib_linear</a>(Object.Coord); |
---|
494 | 0431 <span class="keyword">elseif</span> isequal(calib_type,<span class="string">'tsai'</span>) |
---|
495 | 0432 GeometryCalib=<a href="#_sub15" class="code" title="subfunction GeometryCalib=calib_tsai(Coord)">calib_tsai</a>(Object.Coord); |
---|
496 | 0433 <span class="keyword">end</span> |
---|
497 | 0434 <span class="comment">% %record image source</span> |
---|
498 | 0435 <span class="comment">% hgrid=get(handles.REPLICATE,'parent');%read the calibration image source on the interface userdata</span> |
---|
499 | 0436 <span class="comment">% imagename=get(hgrid,'UserData');</span> |
---|
500 | 0437 <span class="comment">% if exist(imagename,'file')</span> |
---|
501 | 0438 <span class="comment">% GeometryCalib.SourceCalib.ImageCalib=imagename;</span> |
---|
502 | 0439 <span class="comment">% end</span> |
---|
503 | 0440 GeometryCalib.SourceCalib.PointCoord=Object.Coord; |
---|
504 | 0441 |
---|
505 | 0442 <span class="comment">%open and read the dataview GUI</span> |
---|
506 | 0443 h_dataview=findobj(allchild(0),<span class="string">'name'</span>,<span class="string">'dataview'</span>); |
---|
507 | 0444 Device=[];<span class="comment">%default</span> |
---|
508 | 0445 <span class="keyword">if</span> isempty(h_dataview) |
---|
509 | 0446 h_dataview=<a href="dataview.html" class="code" title="function varargout = dataview(varargin)">dataview</a>; |
---|
510 | 0447 hhdataview=guidata(h_dataview); |
---|
511 | 0448 drawnow |
---|
512 | 0449 hGUI=get(handles.REPLICATE,<span class="string">'parent'</span>);<span class="comment">%read the calibration image source on the interface userdata</span> |
---|
513 | 0450 Heading=get(hGUI,<span class="string">'UserData'</span>); |
---|
514 | 0451 <span class="keyword">if</span> isfield(Heading,<span class="string">'XmlInput'</span>) |
---|
515 | 0452 XmlInput=fileparts(Heading.XmlInput); |
---|
516 | 0453 [XmlInput,filename,ext]=fileparts(XmlInput); |
---|
517 | 0454 <span class="keyword">end</span> |
---|
518 | 0455 <span class="keyword">if</span> isfield(Heading,<span class="string">'Record'</span>) && isequal([filename ext],Heading.Record) |
---|
519 | 0456 [XmlInput,filename,ext]=fileparts(XmlInput); |
---|
520 | 0457 <span class="keyword">end</span> |
---|
521 | 0458 <span class="keyword">if</span> isfield(Heading,<span class="string">'Device'</span>) && isequal([filename ext],Heading.Device) |
---|
522 | 0459 [XmlInput,filename,ext]=fileparts(XmlInput); |
---|
523 | 0460 Device=Heading.Device; |
---|
524 | 0461 <span class="keyword">end</span> |
---|
525 | 0462 <span class="keyword">if</span> isfield(Heading,<span class="string">'Experiment'</span>) && isequal([filename ext],Heading.Experiment) |
---|
526 | 0463 [PP,filename,ext]=fileparts(XmlInput); |
---|
527 | 0464 <span class="keyword">end</span> |
---|
528 | 0465 testinput=0; |
---|
529 | 0466 <span class="keyword">if</span> isfield(Heading,<span class="string">'SubCampaign'</span>) && isequal([filename ext],Heading.SubCampaign) |
---|
530 | 0467 set(hhdataview.RootDirectory,<span class="string">'String'</span>,XmlInput) |
---|
531 | 0468 set(hhdataview.SubCampaignTest,<span class="string">'Value'</span>,1) |
---|
532 | 0469 testinput=1; |
---|
533 | 0470 <span class="keyword">elseif</span> isfield(Heading,<span class="string">'Campaign'</span>) && isequal([filename ext],Heading.Campaign) |
---|
534 | 0471 set(hhdataview.RootDirectory,<span class="string">'String'</span>,XmlInput) |
---|
535 | 0472 set(hhdataview.SubCampaignTest,<span class="string">'Value'</span>,0) |
---|
536 | 0473 testinput=1; |
---|
537 | 0474 <span class="keyword">end</span> |
---|
538 | 0475 <span class="keyword">if</span> testinput |
---|
539 | 0476 <a href="dataview.html" class="code" title="function varargout = dataview(varargin)">dataview</a>(<span class="string">'RootDirectory_Callback'</span>,hObject,eventdata,hhdataview) |
---|
540 | 0477 ListDevices=get(hhdataview.ListDevices,<span class="string">'String'</span>); |
---|
541 | 0478 <span class="keyword">for</span> ilist=1:length(ListDevices) |
---|
542 | 0479 <span class="keyword">if</span> isequal(ListDevices{ilist},Device) |
---|
543 | 0480 set(hhdataview.ListDevices,<span class="string">'Value'</span>,ilist) |
---|
544 | 0481 <a href="dataview.html" class="code" title="function varargout = dataview(varargin)">dataview</a>(<span class="string">'ListDevices_Callback'</span>,hObject,eventdata,hhdataview) |
---|
545 | 0482 <span class="keyword">break</span> |
---|
546 | 0483 <span class="keyword">end</span> |
---|
547 | 0484 <span class="keyword">end</span> |
---|
548 | 0485 <span class="keyword">end</span> |
---|
549 | 0486 <span class="keyword">return</span> |
---|
550 | 0487 <span class="keyword">end</span> |
---|
551 | 0488 |
---|
552 | 0489 hhdataview=guidata(h_dataview); |
---|
553 | 0490 CurrentPath=get(hhdataview.RootDirectory,<span class="string">'String'</span>); |
---|
554 | 0491 ListExperiments=get(hhdataview.ListExperiments,<span class="string">'String'</span>); |
---|
555 | 0492 Value=get(hhdataview.ListExperiments,<span class="string">'Value'</span>); |
---|
556 | 0493 <span class="keyword">if</span> ~isequal(Value,1) |
---|
557 | 0494 ListExperiments=ListExperiments(Value); |
---|
558 | 0495 <span class="keyword">end</span> |
---|
559 | 0496 ListDevices=get(hhdataview.ListDevices,<span class="string">'String'</span>); |
---|
560 | 0497 Value=get(hhdataview.ListDevices,<span class="string">'Value'</span>); |
---|
561 | 0498 <span class="keyword">if</span> isequal(Value,1) |
---|
562 | 0499 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'manually select in the GUI dataview the device being calibrated'</span>,<span class="string">'ERROR'</span>) |
---|
563 | 0500 <span class="keyword">return</span> |
---|
564 | 0501 <span class="keyword">else</span> |
---|
565 | 0502 ListDevices=ListDevices(Value); |
---|
566 | 0503 <span class="keyword">end</span> |
---|
567 | 0504 ListRecords=get(hhdataview.ListRecords,<span class="string">'String'</span>); |
---|
568 | 0505 Value=get(hhdataview.ListRecords,<span class="string">'Value'</span>); |
---|
569 | 0506 <span class="keyword">if</span> ~isequal(Value,1) |
---|
570 | 0507 ListRecords=ListRecords(Value); |
---|
571 | 0508 <span class="keyword">end</span> |
---|
572 | 0509 [ListDevices,ListRecords,ListXml,List]=<a href="ListDir.html" class="code" title="function [ListDevices,ListRecords,ListXml,List]=ListDir(CurrentPath,ListExperiments,ListDevices_in,ListRecords_in)">ListDir</a>(CurrentPath,ListExperiments,ListDevices,ListRecords); |
---|
573 | 0510 ListXml=get(hhdataview.ListXml,<span class="string">'String'</span>); |
---|
574 | 0511 Value=get(hhdataview.ListXml,<span class="string">'Value'</span>); |
---|
575 | 0512 <span class="keyword">if</span> isequal(Value,1) |
---|
576 | 0513 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'you need to select in the GUI dataview the xml files to edit'</span>,<span class="string">'ERROR'</span>) |
---|
577 | 0514 <span class="keyword">return</span> |
---|
578 | 0515 <span class="keyword">else</span> |
---|
579 | 0516 ListXml=ListXml(Value); |
---|
580 | 0517 <span class="keyword">end</span> |
---|
581 | 0518 |
---|
582 | 0519 <span class="comment">%update all the selected xml files</span> |
---|
583 | 0520 answer=<a href="msgbox_uvmat.html" class="code" title="function varargout = msgbox_uvmat(varargin)">msgbox_uvmat</a>(<span class="string">'INPUT_Y-N'</span>,[num2str(length(Value)) <span class="string">' xml files for device '</span> ListDevices{1} <span class="string">' will be refreshed with '</span> calib_type <span class="string">' calibration data'</span>]); |
---|
584 | 0521 <span class="keyword">if</span> ~isequal(answer{1},<span class="string">'OK'</span>) |
---|
585 | 0522 <span class="keyword">return</span> |
---|
586 | 0523 <span class="keyword">end</span> |
---|
587 | 0524 <span class="comment">% 'TEST'</span> |
---|
588 | 0525 <span class="comment">% List</span> |
---|
589 | 0526 <span class="comment">% return</span> |
---|
590 | 0527 <span class="keyword">for</span> iexp=1:length(List.Experiment) |
---|
591 | 0528 ExpName=List.Experiment{iexp}.name; |
---|
592 | 0529 <span class="keyword">if</span> isfield(List.Experiment{iexp},<span class="string">'Device'</span>) |
---|
593 | 0530 <span class="keyword">for</span> idevice=1:length(List.Experiment{iexp}.Device) |
---|
594 | 0531 DeviceName=List.Experiment{iexp}.Device{idevice}.name; |
---|
595 | 0532 <span class="keyword">if</span> isfield(List.Experiment{iexp}.Device{idevice},<span class="string">'xmlfile'</span>) |
---|
596 | 0533 <span class="keyword">for</span> ixml=1:length(List.Experiment{iexp}.Device{idevice}.xmlfile) |
---|
597 | 0534 FileName=List.Experiment{iexp}.Device{idevice}.xmlfile{ixml}; |
---|
598 | 0535 <span class="keyword">for</span> ilistxml=1:length(ListXml) |
---|
599 | 0536 <span class="keyword">if</span> isequal(FileName,ListXml{ilistxml}) |
---|
600 | 0537 set(hhdataview.ListXml,<span class="string">'Value'</span>,Value(ilistxml)) |
---|
601 | 0538 drawnow |
---|
602 | 0539 xmlfullname=fullfile(CurrentPath,ExpName,DeviceName,FileName); |
---|
603 | 0540 <a href="#_sub12" class="code" title="subfunction update_imadoc(GeometryCalib,outputfile)">update_imadoc</a>(GeometryCalib,xmlfullname) |
---|
604 | 0541 <span class="keyword">break</span> |
---|
605 | 0542 <span class="keyword">end</span> |
---|
606 | 0543 <span class="keyword">end</span> |
---|
607 | 0544 <span class="keyword">end</span> |
---|
608 | 0545 <span class="keyword">elseif</span> isfield(List.Experiment{iexp}.Device{idevice},<span class="string">'Record'</span>) |
---|
609 | 0546 <span class="keyword">for</span> irecord=1:length(List.Experiment{iexp}.Device{idevice}.Record) |
---|
610 | 0547 RecordName=List.Experiment{iexp}.Device{idevice}.Record{irecord}.name; |
---|
611 | 0548 <span class="keyword">if</span> isfield(List.Experiment{iexp}.Device{idevice}.Record{irecord},<span class="string">'xmlfile'</span>) |
---|
612 | 0549 <span class="keyword">for</span> ixml=1:length(List.Experiment{iexp}.Device{idevice}.Record{irecord}.xmlfile) |
---|
613 | 0550 FileName=List.Experiment{iexp}.Device{idevice}.Record{irecord}.xmlfile{ixml}; |
---|
614 | 0551 <span class="keyword">for</span> ilistxml=1:length(ListXml) |
---|
615 | 0552 <span class="keyword">if</span> isequal(FileName,ListXml{ilistxml}) |
---|
616 | 0553 set(hhdataview.ListXml,<span class="string">'Value'</span>,Value(ilistxml)) |
---|
617 | 0554 drawnow |
---|
618 | 0555 xmlfullname=fullfile(CurrentPath,ExpName,DeviceName,RecordName,FileName); |
---|
619 | 0556 <a href="#_sub12" class="code" title="subfunction update_imadoc(GeometryCalib,outputfile)">update_imadoc</a>(GeometryCalib,xmlfullname) |
---|
620 | 0557 <span class="keyword">break</span> |
---|
621 | 0558 <span class="keyword">end</span> |
---|
622 | 0559 <span class="keyword">end</span> |
---|
623 | 0560 <span class="keyword">end</span> |
---|
624 | 0561 <span class="keyword">end</span> |
---|
625 | 0562 <span class="keyword">end</span> |
---|
626 | 0563 <span class="keyword">end</span> |
---|
627 | 0564 <span class="keyword">end</span> |
---|
628 | 0565 <span class="keyword">end</span> |
---|
629 | 0566 <span class="keyword">end</span> |
---|
630 | 0567 set(hhdataview.ListXml,<span class="string">'Value'</span>,Value) |
---|
631 | 0568 |
---|
632 | 0569 |
---|
633 | 0570 <span class="comment">%-------------------------------------------------------------</span> |
---|
634 | 0571 <a name="_sub12" href="#_subfunctions" class="code">function update_imadoc(GeometryCalib,outputfile)</a> |
---|
635 | 0572 testappend=0; |
---|
636 | 0573 <span class="keyword">if</span> exist(outputfile,<span class="string">'file'</span>);<span class="comment">%=1 if the output file already exists, 0 else</span> |
---|
637 | 0574 t=xmltree(outputfile); <span class="comment">%read the file</span> |
---|
638 | 0575 uid=find(t,<span class="string">'ImaDoc'</span>); |
---|
639 | 0576 <span class="keyword">if</span> isequal(uid,1)<span class="comment">%if the xml file is ImaDoc</span> |
---|
640 | 0577 uid_calib=find(t,<span class="string">'ImaDoc/GeometryCalib'</span>); |
---|
641 | 0578 <span class="keyword">if</span> ~isempty(uid) <span class="comment">%if GeometryCalib already exists, delete its content</span> |
---|
642 | 0579 backupfile=outputfile; |
---|
643 | 0580 testexist=2; |
---|
644 | 0581 <span class="keyword">while</span> testexist==2 |
---|
645 | 0582 backupfile=[backupfile <span class="string">'~'</span>]; |
---|
646 | 0583 testexist=exist(backupfile,<span class="string">'file'</span>); |
---|
647 | 0584 <span class="keyword">end</span> |
---|
648 | 0585 [success,message]=copyfile(outputfile,backupfile);<span class="comment">%make backup</span> |
---|
649 | 0586 <span class="keyword">if</span> isequal(success,1) |
---|
650 | 0587 delete(outputfile) |
---|
651 | 0588 <span class="keyword">else</span> |
---|
652 | 0589 <span class="keyword">return</span> |
---|
653 | 0590 <span class="keyword">end</span> |
---|
654 | 0591 uid_child=children(t,uid_calib); |
---|
655 | 0592 t=delete(t,uid_child); |
---|
656 | 0593 testappend=1; |
---|
657 | 0594 <span class="keyword">end</span> |
---|
658 | 0595 <span class="keyword">end</span> |
---|
659 | 0596 <span class="keyword">end</span> |
---|
660 | 0597 <span class="keyword">if</span> ~testappend |
---|
661 | 0598 t=xmltree; |
---|
662 | 0599 t=set(t,1,<span class="string">'name'</span>,<span class="string">'ImaDoc'</span>); |
---|
663 | 0600 [t,uid_calib]=add(t,1,<span class="string">'element'</span>,<span class="string">'GeometryCalib'</span>); |
---|
664 | 0601 <span class="comment">% t=struct2xml(GeometryCalib,t,uid_calib);</span> |
---|
665 | 0602 <span class="keyword">end</span> |
---|
666 | 0603 |
---|
667 | 0604 t=<a href="struct2xml.html" class="code" title="function t=struct2xml(Object,t,root_uid)">struct2xml</a>(GeometryCalib,t,uid_calib); |
---|
668 | 0605 save(t,outputfile); |
---|
669 | 0606 |
---|
670 | 0607 |
---|
671 | 0608 <span class="comment">%-----------------------------------------------------------------</span> |
---|
672 | 0609 <span class="comment">% determine the parameters for a calibration by an affine function (rescaling and offset, no rotation)</span> |
---|
673 | 0610 <a name="_sub13" href="#_subfunctions" class="code">function GeometryCalib=calib_rescale(Coord)</a> |
---|
674 | 0611 <span class="comment">%------------------------------------------------------------------</span> |
---|
675 | 0612 |
---|
676 | 0613 X=Coord(:,1); |
---|
677 | 0614 Y=Coord(:,2); |
---|
678 | 0615 x_ima=Coord(:,4); |
---|
679 | 0616 y_ima=Coord(:,5); |
---|
680 | 0617 [<a href="px.html" class="code" title="function [DataOut,DataOut_1]=px(Data,CalibData,Data_1,CalibData_1)%DataIn,Calib)">px</a>,sx]=polyfit(X,x_ima,1); |
---|
681 | 0618 [py,sy]=polyfit(Y,y_ima,1); |
---|
682 | 0619 T_x=<a href="px.html" class="code" title="function [DataOut,DataOut_1]=px(Data,CalibData,Data_1,CalibData_1)%DataIn,Calib)">px</a>(2); |
---|
683 | 0620 T_y=py(2); |
---|
684 | 0621 GeometryCalib.CalibrationType=<span class="string">'rescale'</span>; |
---|
685 | 0622 GeometryCalib.focal=1; |
---|
686 | 0623 GeometryCalib.CoordUnit=[];<span class="comment">% default value, to be updated by the calling function</span> |
---|
687 | 0624 GeometryCalib.Tx_Ty_Tz=[T_x T_y 1]; |
---|
688 | 0625 GeometryCalib.R=[<a href="px.html" class="code" title="function [DataOut,DataOut_1]=px(Data,CalibData,Data_1,CalibData_1)%DataIn,Calib)">px</a>(1),0,0;0,py(1),0;0,0,1]; |
---|
689 | 0626 |
---|
690 | 0627 <span class="comment">%check error</span> |
---|
691 | 0628 Calib.dpx=1; |
---|
692 | 0629 Calib.dpy=1; |
---|
693 | 0630 Calib.sx=1; |
---|
694 | 0631 Calib.Cx=0; |
---|
695 | 0632 Calib.Cy=0; |
---|
696 | 0633 Calib.Tz=1; |
---|
697 | 0634 Calib.kappa1=0; |
---|
698 | 0635 Calib.f=GeometryCalib.focal; |
---|
699 | 0636 Calib.Tx=T_x; |
---|
700 | 0637 Calib.Ty=T_y; |
---|
701 | 0638 Calib.R=GeometryCalib.R; |
---|
702 | 0639 [Xpoints,Ypoints]=<a href="px_XYZ.html" class="code" title="function [X,Y]=px_XYZ(Calib,Xphys,Yphys,Zphys)">px_XYZ</a>(Calib,X,Y,0); |
---|
703 | 0640 GeometryCalib.ErrorRms(1)=sqrt(mean((Xpoints-x_ima).*(Xpoints-x_ima))); |
---|
704 | 0641 GeometryCalib.ErrorMax(1)=max(abs(Xpoints-x_ima)); |
---|
705 | 0642 GeometryCalib.ErrorRms(2)=sqrt(mean((Ypoints-y_ima).*(Ypoints-y_ima))); |
---|
706 | 0643 GeometryCalib.ErrorMax(2)=max(abs(Ypoints-y_ima)); |
---|
707 | 0644 |
---|
708 | 0645 |
---|
709 | 0646 <span class="comment">%------------------------------------------------------------------</span> |
---|
710 | 0647 <span class="comment">% determine the parameters for a calibration by a linear transform matrix (rescale and rotation)</span> |
---|
711 | 0648 <a name="_sub14" href="#_subfunctions" class="code">function GeometryCalib=calib_linear(Coord)</a> |
---|
712 | 0649 <span class="comment">%------------------------------------------------------------------</span> |
---|
713 | 0650 X=Coord(:,1); |
---|
714 | 0651 Y=Coord(:,2); |
---|
715 | 0652 x_ima=Coord(:,4); |
---|
716 | 0653 y_ima=Coord(:,5); |
---|
717 | 0654 XY_mat=[ones(size(X)) X Y]; |
---|
718 | 0655 a_X1=XY_mat\x_ima; <span class="comment">%transformation matrix for X</span> |
---|
719 | 0656 x1=XY_mat*a_X1;<span class="comment">%reconstruction</span> |
---|
720 | 0657 err_X1=max(abs(x1-x_ima));<span class="comment">%error</span> |
---|
721 | 0658 a_Y1=XY_mat\y_ima;<span class="comment">%transformation matrix for X</span> |
---|
722 | 0659 y1=XY_mat*a_Y1; |
---|
723 | 0660 err_Y1=max(abs(y1-y_ima));<span class="comment">%error</span> |
---|
724 | 0661 T_x=a_X1(1); |
---|
725 | 0662 T_y=a_Y1(1); |
---|
726 | 0663 GeometryCalib.CalibrationType=<span class="string">'linear'</span>; |
---|
727 | 0664 GeometryCalib.focal=1; |
---|
728 | 0665 GeometryCalib.CoordUnit=[];<span class="comment">% default value, to be updated by the calling function</span> |
---|
729 | 0666 GeometryCalib.Tx_Ty_Tz=[T_x T_y 1]; |
---|
730 | 0667 GeometryCalib.R=[a_X1(2),a_X1(3),0;a_Y1(2),a_Y1(3),0;0,0,1]; |
---|
731 | 0668 |
---|
732 | 0669 <span class="comment">%check error</span> |
---|
733 | 0670 GeometryCalib.ErrorRms(1)=sqrt(mean((x1-x_ima).*(x1-x_ima))); |
---|
734 | 0671 GeometryCalib.ErrorMax(1)=max(abs(x1-x_ima)); |
---|
735 | 0672 GeometryCalib.ErrorRms(2)=sqrt(mean((y1-y_ima).*(y1-y_ima))); |
---|
736 | 0673 GeometryCalib.ErrorMax(2)=max(abs(y1-y_ima)); |
---|
737 | 0674 |
---|
738 | 0675 |
---|
739 | 0676 |
---|
740 | 0677 |
---|
741 | 0678 <span class="comment">%------------------------------------------------------------------</span> |
---|
742 | 0679 <a name="_sub15" href="#_subfunctions" class="code">function GeometryCalib=calib_tsai(Coord)</a> |
---|
743 | 0680 <span class="comment">%------------------------------------------------------------------</span> |
---|
744 | 0681 <span class="comment">%TSAI</span> |
---|
745 | 0682 <span class="comment">% 'calibration_lin' provides a linear transform on coordinates,</span> |
---|
746 | 0683 path_uvmat=which(<span class="string">'uvmat'</span>);<span class="comment">% check the path detected for source file uvmat</span> |
---|
747 | 0684 path_UVMAT=fileparts(path_uvmat); <span class="comment">%path to UVMAT</span> |
---|
748 | 0685 <span class="keyword">if</span> isunix |
---|
749 | 0686 <span class="comment">%fid = fopen(fullfile(path_UVMAT,'PARAM_LINUX.txt'),'r');%open the file with civ binary names</span> |
---|
750 | 0687 xmlfile=fullfile(path_UVMAT,<span class="string">'PARAM_LINUX.xml'</span>); |
---|
751 | 0688 <span class="keyword">if</span> exist(xmlfile,<span class="string">'file'</span>) |
---|
752 | 0689 t=xmltree(xmlfile); |
---|
753 | 0690 sparam=convert(t); |
---|
754 | 0691 <span class="keyword">end</span> |
---|
755 | 0692 <span class="keyword">else</span> |
---|
756 | 0693 <span class="comment">%fid = fopen(fullfile(path_UVMAT,'PARAM_WIN.txt'),'r');%open the file with civ binary names</span> |
---|
757 | 0694 xmlfile=fullfile(path_UVMAT,<span class="string">'PARAM_WIN.xml'</span>); |
---|
758 | 0695 <span class="keyword">if</span> exist(xmlfile,<span class="string">'file'</span>) |
---|
759 | 0696 t=xmltree(xmlfile); |
---|
760 | 0697 sparam=convert(t); |
---|
761 | 0698 <span class="keyword">end</span> |
---|
762 | 0699 <span class="keyword">end</span> |
---|
763 | 0700 <span class="keyword">if</span> ~isfield(sparam,<span class="string">'GeometryCalib_exe'</span>) |
---|
764 | 0701 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'calibration program <GeometryCalib_exe> undefined in parameter file '</span> xmlfile],<span class="string">'ERROR'</span>) |
---|
765 | 0702 <span class="keyword">return</span> |
---|
766 | 0703 <span class="keyword">end</span> |
---|
767 | 0704 Tsai_exe=sparam.GeometryCalib_exe; |
---|
768 | 0705 <span class="keyword">if</span> ~exist(Tsai_exe,<span class="string">'file'</span>) |
---|
769 | 0706 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([<span class="string">'calibration program '</span> Tsai_exe <span class="string">' does not exist'</span>],<span class="string">'ERROR'</span>) |
---|
770 | 0707 <span class="keyword">return</span> |
---|
771 | 0708 <span class="keyword">end</span> |
---|
772 | 0709 |
---|
773 | 0710 textcoord=num2str(Coord,4); |
---|
774 | 0711 dlmwrite(<span class="string">'t.txt'</span>,textcoord,<span class="string">''</span>); |
---|
775 | 0712 <span class="comment">% ['!' Tsai_exe ' -f1 0 -f2 t.txt']</span> |
---|
776 | 0713 eval([<span class="string">'!'</span> Tsai_exe <span class="string">' -f t.txt > tsaicalib.log'</span>]); |
---|
777 | 0714 <span class="keyword">if</span> ~exist(<span class="string">'calib.dat'</span>,<span class="string">'file'</span>) |
---|
778 | 0715 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'no output from calibration program Tsai_exe: possibly too few points'</span>,<span class="string">'ERROR'</span>) |
---|
779 | 0716 <span class="keyword">end</span> |
---|
780 | 0717 calibdat=dlmread(<span class="string">'calib.dat'</span>); |
---|
781 | 0718 GeometryCalib.CalibrationType=<span class="string">'tsai'</span>; |
---|
782 | 0719 GeometryCalib.focal=calibdat(10); |
---|
783 | 0720 GeometryCalib.dpx_dpy=[calibdat(5) calibdat(6)]; |
---|
784 | 0721 GeometryCalib.Cx_Cy=[calibdat(7) calibdat(8)]; |
---|
785 | 0722 GeometryCalib.sx=calibdat(9); |
---|
786 | 0723 GeometryCalib.kappa1=calibdat(11); |
---|
787 | 0724 GeometryCalib.CoordUnit=[];<span class="comment">% default value, to be updated by the calling function</span> |
---|
788 | 0725 GeometryCalib.Tx_Ty_Tz=[calibdat(12) calibdat(13) calibdat(14)]; |
---|
789 | 0726 Rx_Ry_Rz=calibdat([15:17]); |
---|
790 | 0727 sa = sin(Rx_Ry_Rz(1)) ; |
---|
791 | 0728 ca=cos(Rx_Ry_Rz(1)); |
---|
792 | 0729 sb=sin(Rx_Ry_Rz(2)); |
---|
793 | 0730 cb =cos(Rx_Ry_Rz(2)); |
---|
794 | 0731 sg =sin(Rx_Ry_Rz(3)); |
---|
795 | 0732 cg =cos(Rx_Ry_Rz(3)); |
---|
796 | 0733 r1 = cb * cg; |
---|
797 | 0734 r2 = cg * sa * sb - ca * sg; |
---|
798 | 0735 r3 = sa * sg + ca * cg * sb; |
---|
799 | 0736 r4 = cb * sg; |
---|
800 | 0737 r5 = sa * sb * sg + ca * cg; |
---|
801 | 0738 r6 = ca * sb * sg - cg * sa; |
---|
802 | 0739 r7 = -sb; |
---|
803 | 0740 r8 = cb * sa; |
---|
804 | 0741 r9 = ca * cb; |
---|
805 | 0742 <span class="comment">%EN DEDUIRE MATRICE R ??</span> |
---|
806 | 0743 GeometryCalib.R=[r1,r2,r3;r4,r5,r6;r7,r8,r9]; |
---|
807 | 0744 <span class="comment">%erreur a caracteriser?</span> |
---|
808 | 0745 <span class="comment">%check error</span> |
---|
809 | 0746 Calib.dpx=GeometryCalib.dpx_dpy(1); |
---|
810 | 0747 Calib.dpy=GeometryCalib.dpx_dpy(2); |
---|
811 | 0748 Calib.sx=GeometryCalib.sx; |
---|
812 | 0749 Calib.Cx=GeometryCalib.Cx_Cy(1); |
---|
813 | 0750 Calib.Cy=GeometryCalib.Cx_Cy(2); |
---|
814 | 0751 Calib.kappa1=GeometryCalib.kappa1; |
---|
815 | 0752 Calib.f=GeometryCalib.focal; |
---|
816 | 0753 Calib.Tx=GeometryCalib.Tx_Ty_Tz(1); |
---|
817 | 0754 Calib.Ty=GeometryCalib.Tx_Ty_Tz(2); |
---|
818 | 0755 Calib.Tz=GeometryCalib.Tx_Ty_Tz(3); |
---|
819 | 0756 Calib.R=GeometryCalib.R; |
---|
820 | 0757 X=Coord(:,1); |
---|
821 | 0758 Y=Coord(:,2); |
---|
822 | 0759 Z=Coord(:,3); |
---|
823 | 0760 x_ima=Coord(:,4); |
---|
824 | 0761 y_ima=Coord(:,5); |
---|
825 | 0762 [Xpoints,Ypoints]=<a href="px_XYZ.html" class="code" title="function [X,Y]=px_XYZ(Calib,Xphys,Yphys,Zphys)">px_XYZ</a>(Calib,X,Y,Z); |
---|
826 | 0763 |
---|
827 | 0764 GeometryCalib.ErrorRms(1)=sqrt(mean((Xpoints-x_ima).*(Xpoints-x_ima))); |
---|
828 | 0765 GeometryCalib.ErrorMax(1)=max(abs(Xpoints-x_ima)); |
---|
829 | 0766 GeometryCalib.ErrorRms(2)=sqrt(mean((Ypoints-y_ima).*(Ypoints-y_ima))); |
---|
830 | 0767 GeometryCalib.ErrorMax(2)=max(abs(Ypoints-y_ima)); |
---|
831 | 0768 <span class="comment">% Nfx</span> |
---|
832 | 0769 <span class="comment">% dx</span> |
---|
833 | 0770 <span class="comment">% dy</span> |
---|
834 | 0771 <span class="comment">% 5 dpx</span> |
---|
835 | 0772 <span class="comment">% 6 dpy</span> |
---|
836 | 0773 <span class="comment">% cx</span> |
---|
837 | 0774 <span class="comment">% cy</span> |
---|
838 | 0775 <span class="comment">% sx</span> |
---|
839 | 0776 <span class="comment">% f</span> |
---|
840 | 0777 <span class="comment">% kappa1</span> |
---|
841 | 0778 <span class="comment">% tx</span> |
---|
842 | 0779 <span class="comment">% ty</span> |
---|
843 | 0780 <span class="comment">% tz</span> |
---|
844 | 0781 <span class="comment">% rx</span> |
---|
845 | 0782 <span class="comment">% ry</span> |
---|
846 | 0783 <span class="comment">% rz</span> |
---|
847 | 0784 <span class="comment">% p1</span> |
---|
848 | 0785 <span class="comment">% p2</span> |
---|
849 | 0786 |
---|
850 | 0787 <span class="comment">%calibcoeff=str2num(calibdat)</span> |
---|
851 | 0788 |
---|
852 | 0789 <span class="comment">% --- Executes on button press in translation.</span> |
---|
853 | 0790 <a name="_sub16" href="#_subfunctions" class="code">function translation_Callback(hObject, eventdata, handles)</a> |
---|
854 | 0791 |
---|
855 | 0792 |
---|
856 | 0793 <a name="_sub17" href="#_subfunctions" class="code">function T_x_Callback(hObject, eventdata, handles)</a> |
---|
857 | 0794 <span class="comment">% hObject handle to T_x (see GCBO)</span> |
---|
858 | 0795 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
859 | 0796 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
860 | 0797 |
---|
861 | 0798 <span class="comment">% Hints: get(hObject,'String') returns contents of T_x as text</span> |
---|
862 | 0799 <span class="comment">% str2double(get(hObject,'String')) returns contents of T_x as a double</span> |
---|
863 | 0800 |
---|
864 | 0801 |
---|
865 | 0802 |
---|
866 | 0803 |
---|
867 | 0804 |
---|
868 | 0805 <a name="_sub18" href="#_subfunctions" class="code">function T_y_Callback(hObject, eventdata, handles)</a> |
---|
869 | 0806 <span class="comment">% hObject handle to T_y (see GCBO)</span> |
---|
870 | 0807 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
871 | 0808 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
872 | 0809 |
---|
873 | 0810 <span class="comment">% Hints: get(hObject,'String') returns contents of T_y as text</span> |
---|
874 | 0811 <span class="comment">% str2double(get(hObject,'String')) returns contents of T_y as a double</span> |
---|
875 | 0812 |
---|
876 | 0813 |
---|
877 | 0814 <a name="_sub19" href="#_subfunctions" class="code">function T_z_Callback(hObject, eventdata, handles)</a> |
---|
878 | 0815 <span class="comment">% hObject handle to T_z (see GCBO)</span> |
---|
879 | 0816 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
880 | 0817 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
881 | 0818 |
---|
882 | 0819 <span class="comment">% Hints: get(hObject,'String') returns contents of T_z as text</span> |
---|
883 | 0820 <span class="comment">% str2double(get(hObject,'String')) returns contents of T_z as a double</span> |
---|
884 | 0821 |
---|
885 | 0822 |
---|
886 | 0823 <span class="comment">% --- Executes on button press in rotation.</span> |
---|
887 | 0824 <a name="_sub20" href="#_subfunctions" class="code">function rotation_Callback(hObject, eventdata, handles)</a> |
---|
888 | 0825 angle_rot=(pi/180)*str2num(get(handles.Phi,<span class="string">'String'</span>)); |
---|
889 | 0826 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
890 | 0827 data=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
891 | 0828 data.Coord(:,1)=cos(angle_rot)*data.Coord(:,1)+sin(angle_rot)*data.Coord(:,2); |
---|
892 | 0829 data.Coord(:,1)=-sin(angle_rot)*data.Coord(:,1)+cos(angle_rot)*data.Coord(:,2); |
---|
893 | 0830 set(handles.XObject,<span class="string">'String'</span>,num2str(data.Coord(:,1),4)); |
---|
894 | 0831 set(handles.YObject,<span class="string">'String'</span>,num2str(data.Coord(:,2),4)); |
---|
895 | 0832 |
---|
896 | 0833 |
---|
897 | 0834 <a name="_sub21" href="#_subfunctions" class="code">function XImage_Callback(hObject, eventdata, handles)</a> |
---|
898 | 0835 <a href="#_sub26" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
899 | 0836 |
---|
900 | 0837 <a name="_sub22" href="#_subfunctions" class="code">function YImage_Callback(hObject, eventdata, handles)</a> |
---|
901 | 0838 <a href="#_sub26" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
902 | 0839 |
---|
903 | 0840 <a name="_sub23" href="#_subfunctions" class="code">function XObject_Callback(hObject, eventdata, handles)</a> |
---|
904 | 0841 <a href="#_sub26" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
905 | 0842 |
---|
906 | 0843 <a name="_sub24" href="#_subfunctions" class="code">function YObject_Callback(hObject, eventdata, handles)</a> |
---|
907 | 0844 <a href="#_sub26" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
908 | 0845 |
---|
909 | 0846 <a name="_sub25" href="#_subfunctions" class="code">function ZObject_Callback(hObject, eventdata, handles)</a> |
---|
910 | 0847 <a href="#_sub26" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
911 | 0848 |
---|
912 | 0849 <a name="_sub26" href="#_subfunctions" class="code">function update_list(hObject, eventdata, handles)</a> |
---|
913 | 0850 str4=get(handles.XImage,<span class="string">'String'</span>); |
---|
914 | 0851 str5=get(handles.YImage,<span class="string">'String'</span>); |
---|
915 | 0852 str1=get(handles.XObject,<span class="string">'String'</span>); |
---|
916 | 0853 tt=double(str1); |
---|
917 | 0854 str2=get(handles.YObject,<span class="string">'String'</span>); |
---|
918 | 0855 str3=get(handles.ZObject,<span class="string">'String'</span>); |
---|
919 | 0856 <span class="keyword">if</span> ~isempty(str1) & ~isequal(double(str1),32) & (isempty(str3)|isequal(double(str3),32)) |
---|
920 | 0857 str3=<span class="string">'0'</span>;<span class="comment">%put z to 0 by default</span> |
---|
921 | 0858 <span class="keyword">end</span> |
---|
922 | 0859 strline=[str1 <span class="string">' | '</span> str2 <span class="string">' | '</span> str3 <span class="string">' | '</span> str4 <span class="string">' | '</span> str5]; |
---|
923 | 0860 Coord=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
924 | 0861 val=get(handles.ListCoord,<span class="string">'Value'</span>); |
---|
925 | 0862 Coord{val}=strline; |
---|
926 | 0863 set(handles.ListCoord,<span class="string">'String'</span>,Coord) |
---|
927 | 0864 |
---|
928 | 0865 <span class="comment">%--------------------------------------------------------------------</span> |
---|
929 | 0866 <span class="comment">% --- Executes on selection change in ListCoord.</span> |
---|
930 | 0867 <span class="comment">%--------------------------------------------------------------------</span> |
---|
931 | 0868 <a name="_sub27" href="#_subfunctions" class="code">function ListCoord_Callback(hObject, eventdata, handles)</a> |
---|
932 | 0869 <span class="comment">% hObject handle to ListCoord (see GCBO)</span> |
---|
933 | 0870 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
934 | 0871 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
935 | 0872 |
---|
936 | 0873 <span class="comment">% Hints: contents = get(hObject,'String') returns ListCoord contents as cell array</span> |
---|
937 | 0874 <span class="comment">% contents{get(hObject,'Value')} returns selected item from ListCoord</span> |
---|
938 | 0875 <span class="comment">%set(handles.edit_append,'Value',2); %set to edit mode</span> |
---|
939 | 0876 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
940 | 0877 val=get(handles.ListCoord,<span class="string">'Value'</span>); |
---|
941 | 0878 <span class="keyword">if</span> length(Coord_cell)>0 |
---|
942 | 0879 coord_str=Coord_cell{val}; |
---|
943 | 0880 k=findstr(<span class="string">'|'</span>,coord_str); |
---|
944 | 0881 <span class="keyword">if</span> isempty(k) |
---|
945 | 0882 <span class="keyword">return</span> |
---|
946 | 0883 <span class="keyword">end</span> |
---|
947 | 0884 set(handles.XObject,<span class="string">'String'</span>,coord_str(1:k(1)-5)) |
---|
948 | 0885 set(handles.YObject,<span class="string">'String'</span>,coord_str(k(1)+5:k(2)-5)) |
---|
949 | 0886 set(handles.ZObject,<span class="string">'String'</span>,coord_str(k(2)+5:k(3)-5)) |
---|
950 | 0887 set(handles.XImage,<span class="string">'String'</span>,coord_str(k(3)+5:k(4)-5)) |
---|
951 | 0888 set(handles.YImage,<span class="string">'String'</span>,coord_str(k(4)+5:end)) |
---|
952 | 0889 huvmat=findobj(allchild(0),<span class="string">'Name'</span>,<span class="string">'uvmat'</span>);<span class="comment">%find the current uvmat interface handle</span> |
---|
953 | 0890 hplot=findobj(huvmat,<span class="string">'Tag'</span>,<span class="string">'axes3'</span>);<span class="comment">%main plotting axis of uvmat</span> |
---|
954 | 0891 h_menu_coord=findobj(huvmat,<span class="string">'Tag'</span>,<span class="string">'menu_coord'</span>); |
---|
955 | 0892 menu=get(h_menu_coord,<span class="string">'String'</span>); |
---|
956 | 0893 choice=get(h_menu_coord,<span class="string">'Value'</span>); |
---|
957 | 0894 <span class="keyword">if</span> iscell(menu) |
---|
958 | 0895 option=menu{choice}; |
---|
959 | 0896 <span class="keyword">else</span> |
---|
960 | 0897 option=<span class="string">'px'</span>; <span class="comment">%default</span> |
---|
961 | 0898 <span class="keyword">end</span> |
---|
962 | 0899 <span class="keyword">if</span> isequal(option,<span class="string">'phys'</span>) |
---|
963 | 0900 XCoord=str2num(coord_str(1:k(1)-5)); |
---|
964 | 0901 YCoord=str2num(coord_str(k(1)+5:k(2)-5)); |
---|
965 | 0902 <span class="keyword">elseif</span> isequal(option,<span class="string">'px'</span>)|| isequal(option,<span class="string">''</span>) |
---|
966 | 0903 XCoord=str2num(coord_str(k(3)+5:k(4)-5)); |
---|
967 | 0904 YCoord=str2num(coord_str(k(4)+5:end)); |
---|
968 | 0905 <span class="keyword">else</span> |
---|
969 | 0906 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'the choice in menu_coord of uvmat must be px or phys '</span>,<span class="string">'ERROR'</span>) |
---|
970 | 0907 <span class="keyword">end</span> |
---|
971 | 0908 huvmat=findobj(allchild(0),<span class="string">'Name'</span>,<span class="string">'uvmat'</span>);<span class="comment">%find the current uvmat interface handle</span> |
---|
972 | 0909 hplot=findobj(huvmat,<span class="string">'Tag'</span>,<span class="string">'axes3'</span>);<span class="comment">%main plotting axis of uvmat</span> |
---|
973 | 0910 hhh=findobj(hplot,<span class="string">'Tag'</span>,<span class="string">'calib_marker'</span>); |
---|
974 | 0911 <span class="keyword">if</span> isempty(hhh) |
---|
975 | 0912 axes(hplot) |
---|
976 | 0913 line(XCoord,YCoord,<span class="string">'Color'</span>,<span class="string">'m'</span>,<span class="string">'Tag'</span>,<span class="string">'calib_marker'</span>,<span class="string">'LineStyle'</span>,<span class="string">'.'</span>,<span class="string">'Marker'</span>,<span class="string">'o'</span>,<span class="string">'MarkerSize'</span>,20); |
---|
977 | 0914 <span class="keyword">else</span> |
---|
978 | 0915 set(hhh,<span class="string">'XData'</span>,XCoord) |
---|
979 | 0916 set(hhh,<span class="string">'YData'</span>,YCoord) |
---|
980 | 0917 <span class="keyword">end</span> |
---|
981 | 0918 <span class="keyword">end</span> |
---|
982 | 0919 <span class="comment">%------------------------------------------------------</span> |
---|
983 | 0920 <span class="comment">% --- Executes on button press in translation_plus.</span> |
---|
984 | 0921 <a name="_sub28" href="#_subfunctions" class="code">function translation_plus_Callback(hObject, eventdata, handles)</a> |
---|
985 | 0922 |
---|
986 | 0923 T=[0 0 0]; |
---|
987 | 0924 T_x=get(handles.T_x,<span class="string">'String'</span>); |
---|
988 | 0925 T_y=get(handles.T_y,<span class="string">'String'</span>); |
---|
989 | 0926 T_z=get(handles.T_z,<span class="string">'String'</span>); |
---|
990 | 0927 <span class="keyword">if</span> ~isempty(T_x) |
---|
991 | 0928 T(1)=str2num(T_x); |
---|
992 | 0929 <span class="keyword">end</span> |
---|
993 | 0930 <span class="keyword">if</span> ~isempty(T_y) |
---|
994 | 0931 T(2)=str2num(T_y); |
---|
995 | 0932 <span class="keyword">end</span> |
---|
996 | 0933 <span class="keyword">if</span> ~isempty(T_z) |
---|
997 | 0934 T(3)=str2num(T_z); |
---|
998 | 0935 <span class="keyword">end</span> |
---|
999 | 0936 <a href="#_sub30" class="code" title="subfunction translation(handles,T)">translation</a>(handles,T) |
---|
1000 | 0937 |
---|
1001 | 0938 |
---|
1002 | 0939 |
---|
1003 | 0940 <span class="comment">% --- Executes on button press in translation_minus.</span> |
---|
1004 | 0941 <a name="_sub29" href="#_subfunctions" class="code">function translation_minus_Callback(hObject, eventdata, handles)</a> |
---|
1005 | 0942 |
---|
1006 | 0943 T=[0 0 0]; |
---|
1007 | 0944 T_x=get(handles.T_x,<span class="string">'String'</span>); |
---|
1008 | 0945 T_y=get(handles.T_y,<span class="string">'String'</span>); |
---|
1009 | 0946 T_z=get(handles.T_z,<span class="string">'String'</span>); |
---|
1010 | 0947 <span class="keyword">if</span> ~isempty(T_x) |
---|
1011 | 0948 T(1)=-str2num(T_x); |
---|
1012 | 0949 <span class="keyword">end</span> |
---|
1013 | 0950 <span class="keyword">if</span> ~isempty(T_y) |
---|
1014 | 0951 T(2)=-str2num(T_y); |
---|
1015 | 0952 <span class="keyword">end</span> |
---|
1016 | 0953 <span class="keyword">if</span> ~isempty(T_z) |
---|
1017 | 0954 T(3)=-str2num(T_z); |
---|
1018 | 0955 <span class="keyword">end</span> |
---|
1019 | 0956 <a href="#_sub30" class="code" title="subfunction translation(handles,T)">translation</a>(handles,T) |
---|
1020 | 0957 |
---|
1021 | 0958 |
---|
1022 | 0959 <span class="comment">%%%--------------------------------------</span> |
---|
1023 | 0960 <a name="_sub30" href="#_subfunctions" class="code">function translation(handles,T)</a> |
---|
1024 | 0961 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
1025 | 0962 data=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
1026 | 0963 <span class="comment">% data=read_geometry_calib(handles);</span> |
---|
1027 | 0964 data.Coord(:,1)=T(1)+data.Coord(:,1); |
---|
1028 | 0965 data.Coord(:,2)=T(2)+data.Coord(:,2); |
---|
1029 | 0966 data.Coord(:,3)=T(3)+data.Coord(:,3); |
---|
1030 | 0967 data.Coord(:,[4 5])=data.Coord(:,[4 5]); |
---|
1031 | 0968 <span class="keyword">for</span> i=1:size(data.Coord,1) |
---|
1032 | 0969 <span class="keyword">for</span> j=1:5 |
---|
1033 | 0970 Coord{i,j}=num2str(data.Coord(i,j),4);<span class="comment">%phys x,y,z</span> |
---|
1034 | 0971 <span class="keyword">end</span> |
---|
1035 | 0972 <span class="keyword">end</span> |
---|
1036 | 0973 Tabchar=<a href="cell2tab.html" class="code" title="function Tabchar=cell2tab(Tabcell,separator)">cell2tab</a>(Coord,<span class="string">' | '</span>); |
---|
1037 | 0974 set(handles.ListCoord,<span class="string">'Value'</span>,1) |
---|
1038 | 0975 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
1039 | 0976 |
---|
1040 | 0977 <span class="comment">%----------------------------------------------------</span> |
---|
1041 | 0978 <span class="comment">% --- Executes on button press in rotation_plus.</span> |
---|
1042 | 0979 <a name="_sub31" href="#_subfunctions" class="code">function rotation_plus_Callback(hObject, eventdata, handles)</a> |
---|
1043 | 0980 Phi=0; |
---|
1044 | 0981 Phi=get(handles.Phi,<span class="string">'String'</span>); |
---|
1045 | 0982 <span class="keyword">if</span> ~isempty(Phi) |
---|
1046 | 0983 Phi=str2num(Phi); |
---|
1047 | 0984 <span class="keyword">end</span> |
---|
1048 | 0985 <a href="#_sub33" class="code" title="subfunction rotation(handles,Phi)">rotation</a>(handles,Phi) |
---|
1049 | 0986 |
---|
1050 | 0987 <span class="comment">%-------------------------------------------------</span> |
---|
1051 | 0988 <span class="comment">% --- Executes on button press in rotation_minus.</span> |
---|
1052 | 0989 <a name="_sub32" href="#_subfunctions" class="code">function rotation_minus_Callback(hObject, eventdata, handles)</a> |
---|
1053 | 0990 Phi=0; |
---|
1054 | 0991 Phi=get(handles.Phi,<span class="string">'String'</span>); |
---|
1055 | 0992 <span class="keyword">if</span> ~isempty(Phi) |
---|
1056 | 0993 Phi=-str2num(Phi); |
---|
1057 | 0994 <span class="keyword">end</span> |
---|
1058 | 0995 <a href="#_sub33" class="code" title="subfunction rotation(handles,Phi)">rotation</a>(handles,Phi) |
---|
1059 | 0996 |
---|
1060 | 0997 <span class="comment">%-----------------------------------------------------</span> |
---|
1061 | 0998 <span class="comment">%rotation</span> |
---|
1062 | 0999 <a name="_sub33" href="#_subfunctions" class="code">function rotation(handles,Phi)</a> |
---|
1063 | 1000 O_x=str2num(get(handles.O_x,<span class="string">'String'</span>)); |
---|
1064 | 1001 O_y=str2num(get(handles.O_y,<span class="string">'String'</span>)); |
---|
1065 | 1002 <span class="keyword">if</span> isempty(O_x) |
---|
1066 | 1003 O_x=0;<span class="comment">%default</span> |
---|
1067 | 1004 <span class="keyword">end</span> |
---|
1068 | 1005 <span class="keyword">if</span> isempty(O_y) |
---|
1069 | 1006 O_y=0;<span class="comment">%default</span> |
---|
1070 | 1007 <span class="keyword">end</span> |
---|
1071 | 1008 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
1072 | 1009 data=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
1073 | 1010 <span class="comment">%data=read_geometry_calib(handles);</span> |
---|
1074 | 1011 r1=cos(pi*Phi/180); |
---|
1075 | 1012 r2=-sin(pi*Phi/180); |
---|
1076 | 1013 r3=sin(pi*Phi/180); |
---|
1077 | 1014 r4=cos(pi*Phi/180); |
---|
1078 | 1015 x=data.Coord(:,1); |
---|
1079 | 1016 y=data.Coord(:,2); |
---|
1080 | 1017 data.Coord(:,1)=r1*x+r2*y; |
---|
1081 | 1018 data.Coord(:,2)=r3*x+r4*y; |
---|
1082 | 1019 <span class="comment">% data.Coord(:,[4 5])=data.Coord(:,[4 5]);</span> |
---|
1083 | 1020 <span class="keyword">for</span> i=1:size(data.Coord,1) |
---|
1084 | 1021 <span class="keyword">for</span> j=1:5 |
---|
1085 | 1022 Coord{i,j}=num2str(data.Coord(i,j),4);<span class="comment">%phys x,y,z</span> |
---|
1086 | 1023 <span class="keyword">end</span> |
---|
1087 | 1024 <span class="keyword">end</span> |
---|
1088 | 1025 Tabchar=<a href="cell2tab.html" class="code" title="function Tabchar=cell2tab(Tabcell,separator)">cell2tab</a>(Coord,<span class="string">' | '</span>); |
---|
1089 | 1026 set(handles.ListCoord,<span class="string">'Value'</span>,1) |
---|
1090 | 1027 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
1091 | 1028 |
---|
1092 | 1029 |
---|
1093 | 1030 <a name="_sub34" href="#_subfunctions" class="code">function O_x_Callback(hObject, eventdata, handles)</a> |
---|
1094 | 1031 |
---|
1095 | 1032 |
---|
1096 | 1033 <a name="_sub35" href="#_subfunctions" class="code">function O_y_Callback(hObject, eventdata, handles)</a> |
---|
1097 | 1034 |
---|
1098 | 1035 |
---|
1099 | 1036 <a name="_sub36" href="#_subfunctions" class="code">function O_z_Callback(hObject, eventdata, handles)</a> |
---|
1100 | 1037 |
---|
1101 | 1038 |
---|
1102 | 1039 <span class="comment">% --- Executes on selection change in edit_append.</span> |
---|
1103 | 1040 <a name="_sub37" href="#_subfunctions" class="code">function edit_append_Callback(hObject, eventdata, handles)</a> |
---|
1104 | 1041 <span class="comment">% val=get(handles.PLOT_append,'Value');</span> |
---|
1105 | 1042 <span class="comment">% if isequal(val,2); %append mode</span> |
---|
1106 | 1043 <span class="comment">% %appeler mouse</span> |
---|
1107 | 1044 <span class="comment">% end</span> |
---|
1108 | 1045 choice=get(handles.edit_append,<span class="string">'Value'</span>); |
---|
1109 | 1046 <span class="keyword">if</span> choice==1 |
---|
1110 | 1047 Coord=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
1111 | 1048 val=length(Coord); |
---|
1112 | 1049 <span class="keyword">if</span> val>=1 & isequal(Coord{val},<span class="string">''</span>) |
---|
1113 | 1050 val=val-1; <span class="comment">%do not take into account blank</span> |
---|
1114 | 1051 <span class="keyword">end</span> |
---|
1115 | 1052 Coord{val+1}=<span class="string">''</span>; |
---|
1116 | 1053 set(handles.ListCoord,<span class="string">'String'</span>,Coord) |
---|
1117 | 1054 set(handles.ListCoord,<span class="string">'Value'</span>,val+1) |
---|
1118 | 1055 <span class="keyword">end</span> |
---|
1119 | 1056 |
---|
1120 | 1057 |
---|
1121 | 1058 <span class="comment">%A REVOIR</span> |
---|
1122 | 1059 <span class="comment">% if choice==2</span> |
---|
1123 | 1060 <span class="comment">% %display image with px coordinates</span> |
---|
1124 | 1061 <span class="comment">% hrootpath=findobj(huvmat,'Tag','RootPath');</span> |
---|
1125 | 1062 <span class="comment">% hrootfile=findobj(huvmat,'Tag','RootFile');</span> |
---|
1126 | 1063 <span class="comment">% RootPath='';</span> |
---|
1127 | 1064 <span class="comment">% RootFile='';</span> |
---|
1128 | 1065 <span class="comment">% % if ~isempty(hrootpath)& ~isempty(hrootfile)</span> |
---|
1129 | 1066 <span class="comment">% testhandle=1;</span> |
---|
1130 | 1067 <span class="comment">% RootPath=get(hrootpath,'String');</span> |
---|
1131 | 1068 <span class="comment">% RootFile=get(hrootfile,'String');</span> |
---|
1132 | 1069 <span class="comment">% % filebase=fullfile(RootPath,RootFile);</span> |
---|
1133 | 1070 <span class="comment">% % outputfile=[filebase '.xml'];</span> |
---|
1134 | 1071 <span class="comment">% Indices=get(findobj(huvmat,'Tag','FileIndex'),'String');</span> |
---|
1135 | 1072 <span class="comment">% Ext=get(findobj(huvmat,'Tag','FileExt'),'String');</span> |
---|
1136 | 1073 <span class="comment">% imagename=[fullfile(RootPath,RootFile) Indices Ext];</span> |
---|
1137 | 1074 <span class="comment">% % input.menu_coord=1;</span> |
---|
1138 | 1075 <span class="comment">% h_menu_coord=findobj(huvmat,'Tag','menu_coord');</span> |
---|
1139 | 1076 <span class="comment">% set(h_menu_coord,'Value',3)</span> |
---|
1140 | 1077 <span class="comment">% huvmat=uvmat(imagename);%open uvmat, set phys coord (Value 1)</span> |
---|
1141 | 1078 <span class="comment">%</span> |
---|
1142 | 1079 <span class="comment">% % end</span> |
---|
1143 | 1080 <span class="comment">% end</span> |
---|
1144 | 1081 |
---|
1145 | 1082 <a name="_sub38" href="#_subfunctions" class="code">function NEW_Callback(hObject, eventdata, handles)</a> |
---|
1146 | 1083 <span class="comment">%A METTRE SOUS UN BOUTON</span> |
---|
1147 | 1084 huvmat=findobj(allchild(0),<span class="string">'Name'</span>,<span class="string">'uvmat'</span>); |
---|
1148 | 1085 hchild=get(huvmat,<span class="string">'children'</span>); |
---|
1149 | 1086 hcoord=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'menu_coord'</span>); |
---|
1150 | 1087 coordtype=get(hcoord,<span class="string">'Value'</span>); |
---|
1151 | 1088 haxes=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'axes3'</span>); |
---|
1152 | 1089 AxeData=get(haxes,<span class="string">'UserData'</span>); |
---|
1153 | 1090 <span class="keyword">if</span> ~isequal(hcoord,2) |
---|
1154 | 1091 set(hcoord,<span class="string">'Value'</span>,2) |
---|
1155 | 1092 huvmat=<a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a>(AxeData); |
---|
1156 | 1093 <span class="string">'relancer uvmat'</span> |
---|
1157 | 1094 <span class="keyword">end</span> |
---|
1158 | 1095 <span class="keyword">if</span> ~isfield(AxeData,<span class="string">'ZoomAxes'</span>) |
---|
1159 | 1096 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'first draw a window around a grid marker'</span>,<span class="string">'ERRROR'</span>) |
---|
1160 | 1097 <span class="keyword">return</span> |
---|
1161 | 1098 <span class="keyword">end</span> |
---|
1162 | 1099 XLim=get(AxeData.ZoomAxes,<span class="string">'XLim'</span>); |
---|
1163 | 1100 YLim=get(AxeData.ZoomAxes,<span class="string">'YLim'</span>); |
---|
1164 | 1101 np=size(AxeData.A); |
---|
1165 | 1102 ind_sub_x=round(XLim) |
---|
1166 | 1103 ind_sub_y=np(1)-round(YLim) |
---|
1167 | 1104 Mfiltre=AxeData.A([ind_sub_y(2):ind_sub_y(1)] ,ind_sub_x,:); |
---|
1168 | 1105 Mfiltre_norm=double(Mfiltre); |
---|
1169 | 1106 Mfiltre_norm=Mfiltre_norm/sum(sum(Mfiltre_norm)); |
---|
1170 | 1107 Mfiltre_norm=100*(Mfiltre_norm-mean(mean(Mfiltre_norm))); |
---|
1171 | 1108 Atype=class(AxeData.A) |
---|
1172 | 1109 Data.NbDim=2; |
---|
1173 | 1110 Data.A=filter2(Mfiltre_norm,double(AxeData.A)); |
---|
1174 | 1111 Data.A=feval(Atype,Data.A); |
---|
1175 | 1112 Data.AName=<span class="string">'image'</span>; |
---|
1176 | 1113 Data.AX=AxeData.AX; |
---|
1177 | 1114 Data.AY=AxeData.AY; |
---|
1178 | 1115 Data.CoordType=<span class="string">'px'</span>; |
---|
1179 | 1116 <a href="plot_field.html" class="code" title="function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)">plot_field</a>(Data) |
---|
1180 | 1117 |
---|
1181 | 1118 |
---|
1182 | 1119 <span class="comment">% --- Executes on button press in HELP.</span> |
---|
1183 | 1120 <a name="_sub39" href="#_subfunctions" class="code">function HELP_Callback(hObject, eventdata, handles)</a> |
---|
1184 | 1121 <span class="comment">% hObject handle to HELP (see GCBO)</span> |
---|
1185 | 1122 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
1186 | 1123 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
1187 | 1124 path_to_uvmat=which (<span class="string">'uvmat'</span>);<span class="comment">% check the path of uvmat</span> |
---|
1188 | 1125 pathelp=fileparts(path_to_uvmat); |
---|
1189 | 1126 helpfile=fullfile(pathelp,<span class="string">'UVMAT_DOC'</span>,<span class="string">'uvmat_doc.html'</span>); |
---|
1190 | 1127 <span class="keyword">if</span> isempty(dir(helpfile)), <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>(<span class="string">'The help file uvmat_doc.html needs to be put in the directory UVMAT/UVMAT_DOC'</span>,<span class="string">'ERROR'</span>) |
---|
1191 | 1128 <span class="keyword">else</span> |
---|
1192 | 1129 web([helpfile <span class="string">'#geometry_calib'</span>]) |
---|
1193 | 1130 <span class="keyword">end</span> |
---|
1194 | 1131 |
---|
1195 | 1132 <span class="comment">%'key_press_fcn:' function activated when a key is pressed on the keyboard</span> |
---|
1196 | 1133 <span class="comment">%-----------------------------------</span> |
---|
1197 | 1134 <a name="_sub40" href="#_subfunctions" class="code">function key_press_fcn(hObject,eventdata,handles)</a> |
---|
1198 | 1135 hh=get(hObject,<span class="string">'parent'</span>); |
---|
1199 | 1136 xx=double(get(hh,<span class="string">'CurrentCharacter'</span>)) <span class="comment">%get the keyboard character</span> |
---|
1200 | 1137 |
---|
1201 | 1138 <span class="keyword">if</span> isequal(xx,8)<span class="comment">%move arrow right</span> |
---|
1202 | 1139 Coord_cell=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
1203 | 1140 data=<a href="read_geometry_calib.html" class="code" title="function data=read_geometry_calib(Coord_cell)">read_geometry_calib</a>(Coord_cell); |
---|
1204 | 1141 Coord=[]; <span class="comment">%default</span> |
---|
1205 | 1142 <span class="keyword">if</span> isfield(data,<span class="string">'Coord'</span>) |
---|
1206 | 1143 Coord=data.Coord; |
---|
1207 | 1144 <span class="keyword">end</span> |
---|
1208 | 1145 val=get(handles.ListCoord,<span class="string">'Value'</span>); |
---|
1209 | 1146 Coord(val,:)=[];<span class="comment">%suppress the selected item in the list</span> |
---|
1210 | 1147 CoordCell={}; |
---|
1211 | 1148 <span class="keyword">for</span> iline=1:size(Coord,1) |
---|
1212 | 1149 <span class="keyword">for</span> j=1:5 |
---|
1213 | 1150 CoordCell{iline,j}=num2str(Coord(iline,j),4); |
---|
1214 | 1151 <span class="keyword">end</span> |
---|
1215 | 1152 <span class="keyword">end</span> |
---|
1216 | 1153 Tabchar=<a href="cell2tab.html" class="code" title="function Tabchar=cell2tab(Tabcell,separator)">cell2tab</a>(CoordCell,<span class="string">' | '</span>);<span class="comment">%transform cells into table ready for display</span> |
---|
1217 | 1154 val=min(size(Coord,1),val); |
---|
1218 | 1155 set(handles.ListCoord,<span class="string">'Value'</span>,max(val,1)) |
---|
1219 | 1156 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
1220 | 1157 <a href="#_sub27" class="code" title="subfunction ListCoord_Callback(hObject, eventdata, handles)">ListCoord_Callback</a>(hObject, eventdata, handles) |
---|
1221 | 1158 <span class="keyword">end</span> |
---|
1222 | 1159 |
---|
1223 | 1160 |
---|
1224 | 1161 <span class="comment">% --- Executes on button press in append_point.</span> |
---|
1225 | 1162 <a name="_sub41" href="#_subfunctions" class="code">function append_point_Callback(hObject, eventdata, handles)</a> |
---|
1226 | 1163 |
---|
1227 | 1164 Coord=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
1228 | 1165 val=length(Coord); |
---|
1229 | 1166 <span class="keyword">if</span> val>=1 & isequal(Coord{val},<span class="string">''</span>) |
---|
1230 | 1167 val=val-1; <span class="comment">%do not take into account blank</span> |
---|
1231 | 1168 <span class="keyword">end</span> |
---|
1232 | 1169 Coord{val+1}=<span class="string">''</span>; |
---|
1233 | 1170 set(handles.ListCoord,<span class="string">'String'</span>,Coord) |
---|
1234 | 1171 set(handles.ListCoord,<span class="string">'Value'</span>,val+1) |
---|
1235 | 1172 |
---|
1236 | 1173 |
---|
1237 | 1174</pre></div> |
---|
1238 | <hr><address>Generated on Fri 13-Nov-2009 11:17:03 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> © 2003</address> |
---|
1239 | </body> |
---|
1240 | </html> |
---|