[37] | 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 probe_calib</title> |
---|
| 6 | <meta name="keywords" content="probe_calib"> |
---|
| 7 | <meta name="description" content="'probe_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> > probe_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>probe_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>'probe_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 = probe_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">'probe_calib': performs geometric calibration from a set of reference points |
---|
| 31 | PROBE_CALIB M-file for probe_calib.fig |
---|
| 32 | PROBE_CALIB, by itself, creates a MenuCoord PROBE_CALIB or raises the existing |
---|
| 33 | singleton*. |
---|
| 34 | |
---|
| 35 | H = PROBE_CALIB returns the handle to a MenuCoord PROBE_CALIB or the handle to |
---|
| 36 | the existing singleton*. |
---|
| 37 | |
---|
| 38 | PROBE_CALIB('CALLBACK',hObject,eventData,handles,...) calls the local |
---|
| 39 | function named CALLBACK in PROBE_CALIB.M with the given input arguments. |
---|
| 40 | |
---|
| 41 | PROBE_CALIB('Property','Value',...) creates a MenuCoord PROBE_CALIB or raises the |
---|
| 42 | existing singleton*. Starting from the left, property value pairs are |
---|
| 43 | applied to the GUI before probe_calib_OpeningFunction gets called. An |
---|
| 44 | unrecognized property name or invalid value makes property application |
---|
| 45 | stop. All inputs are passed to probe_calib_OpeningFcn via varargin. |
---|
| 46 | |
---|
| 47 | *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one |
---|
| 48 | instance to run (singleton)". |
---|
| 49 | |
---|
| 50 | See also: GUIDE, GUIDATA, GUIHANDLES</pre></div> |
---|
| 51 | |
---|
| 52 | <!-- crossreference --> |
---|
| 53 | <h2><a name="_cross"></a>CROSS-REFERENCE INFORMATION <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 54 | This function calls: |
---|
| 55 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 56 | <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="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="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="plot_object.html" class="code" title="function [hh]=plot_object(ObjectDataIn,ProjObject,hplot,col)">plot_object</a> 'plot_object': draws a projection object (points, line, plane...)</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="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> |
---|
| 57 | This function is called by: |
---|
| 58 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 59 | </ul> |
---|
| 60 | <!-- crossreference --> |
---|
| 61 | |
---|
| 62 | <h2><a name="_subfunctions"></a>SUBFUNCTIONS <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 63 | <ul style="list-style-image:url(../matlabicon.gif)"> |
---|
| 64 | <li><a href="#_sub1" class="code">function probe_calib_OpeningFcn(hObject, eventdata, handles, data,pos,inputfile)</a></li><li><a href="#_sub2" class="code">function varargout = probe_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 import_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)</a></li><li><a href="#_sub7" class="code">function edit_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 calib_offset_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub11" class="code">function calib_lin_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub12" class="code">function translation_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub13" class="code">function T_x_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub14" class="code">function T_y_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub15" class="code">function T_z_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub16" class="code">function rotation_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub17" class="code">function XImage_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub18" class="code">function YImage_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub19" class="code">function XObject_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub20" class="code">function YObject_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub21" class="code">function ZObject_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub22" class="code">function update_list(hObject, eventdata, handles)</a></li><li><a href="#_sub23" class="code">function data=read_probe_calib(handles)</a></li><li><a href="#_sub24" class="code">function ListCoord_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub25" class="code">function translation_plus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub26" class="code">function translation_minus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub27" class="code">function translation(handles,T)</a></li><li><a href="#_sub28" class="code">function rotation_plus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub29" class="code">function rotation_minus_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub30" class="code">function rotation(handles,Phi)</a></li><li><a href="#_sub31" class="code">function O_x_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub32" class="code">function O_y_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub33" class="code">function O_z_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub34" class="code">function edit_append_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub35" class="code">function NEW_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub36" class="code">function HELP_Callback(hObject, eventdata, handles)</a></li><li><a href="#_sub37" class="code">function key_press_fcn(hObject,eventdata,handles)</a></li></ul> |
---|
| 65 | <h2><a name="_source"></a>SOURCE CODE <a href="#_top"><img alt="^" border="0" src="../up.png"></a></h2> |
---|
| 66 | <div class="fragment"><pre>0001 <span class="comment">%'probe_calib': performs geometric calibration from a set of reference points</span> |
---|
| 67 | 0002 <a name="_sub0" href="#_subfunctions" class="code">function varargout = probe_calib(varargin)</a> |
---|
| 68 | 0003 <span class="comment">% PROBE_CALIB M-file for probe_calib.fig</span> |
---|
| 69 | 0004 <span class="comment">% PROBE_CALIB, by itself, creates a MenuCoord PROBE_CALIB or raises the existing</span> |
---|
| 70 | 0005 <span class="comment">% singleton*.</span> |
---|
| 71 | 0006 <span class="comment">%</span> |
---|
| 72 | 0007 <span class="comment">% H = PROBE_CALIB returns the handle to a MenuCoord PROBE_CALIB or the handle to</span> |
---|
| 73 | 0008 <span class="comment">% the existing singleton*.</span> |
---|
| 74 | 0009 <span class="comment">%</span> |
---|
| 75 | 0010 <span class="comment">% PROBE_CALIB('CALLBACK',hObject,eventData,handles,...) calls the local</span> |
---|
| 76 | 0011 <span class="comment">% function named CALLBACK in PROBE_CALIB.M with the given input arguments.</span> |
---|
| 77 | 0012 <span class="comment">%</span> |
---|
| 78 | 0013 <span class="comment">% PROBE_CALIB('Property','Value',...) creates a MenuCoord PROBE_CALIB or raises the</span> |
---|
| 79 | 0014 <span class="comment">% existing singleton*. Starting from the left, property value pairs are</span> |
---|
| 80 | 0015 <span class="comment">% applied to the GUI before probe_calib_OpeningFunction gets called. An</span> |
---|
| 81 | 0016 <span class="comment">% unrecognized property name or invalid value makes property application</span> |
---|
| 82 | 0017 <span class="comment">% stop. All inputs are passed to probe_calib_OpeningFcn via varargin.</span> |
---|
| 83 | 0018 <span class="comment">%</span> |
---|
| 84 | 0019 <span class="comment">% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one</span> |
---|
| 85 | 0020 <span class="comment">% instance to run (singleton)".</span> |
---|
| 86 | 0021 <span class="comment">%</span> |
---|
| 87 | 0022 <span class="comment">% See also: GUIDE, GUIDATA, GUIHANDLES</span> |
---|
| 88 | 0023 |
---|
| 89 | 0024 <span class="comment">% Edit the above text to modify the response to help probe_calib</span> |
---|
| 90 | 0025 |
---|
| 91 | 0026 <span class="comment">% Last Modified by GUIDE v2.5 04-Feb-2008 15:46:42</span> |
---|
| 92 | 0027 |
---|
| 93 | 0028 <span class="comment">% Begin initialization code - DO NOT EDIT</span> |
---|
| 94 | 0029 gui_Singleton = 1; |
---|
| 95 | 0030 gui_State = struct(<span class="string">'gui_Name'</span>, mfilename, <span class="keyword">...</span> |
---|
| 96 | 0031 <span class="string">'gui_Singleton'</span>, gui_Singleton, <span class="keyword">...</span> |
---|
| 97 | 0032 <span class="string">'gui_OpeningFcn'</span>, @<a href="#_sub1" class="code" title="subfunction probe_calib_OpeningFcn(hObject, eventdata, handles, data,pos,inputfile)">probe_calib_OpeningFcn</a>, <span class="keyword">...</span> |
---|
| 98 | 0033 <span class="string">'gui_OutputFcn'</span>, @<a href="#_sub2" class="code" title="subfunction varargout = probe_calib_OutputFcn(hObject, eventdata, handles)">probe_calib_OutputFcn</a>, <span class="keyword">...</span> |
---|
| 99 | 0034 <span class="string">'gui_LayoutFcn'</span>, [] , <span class="keyword">...</span> |
---|
| 100 | 0035 <span class="string">'gui_Callback'</span>, []); |
---|
| 101 | 0036 <span class="keyword">if</span> nargin & isstr(varargin{1}) |
---|
| 102 | 0037 gui_State.gui_Callback = str2func(varargin{1}); |
---|
| 103 | 0038 <span class="keyword">end</span> |
---|
| 104 | 0039 |
---|
| 105 | 0040 <span class="keyword">if</span> nargout |
---|
| 106 | 0041 [varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:}); |
---|
| 107 | 0042 <span class="keyword">else</span> |
---|
| 108 | 0043 gui_mainfcn(gui_State, varargin{:}); |
---|
| 109 | 0044 <span class="keyword">end</span> |
---|
| 110 | 0045 <span class="comment">% End initialization code - DO NOT EDIT</span> |
---|
| 111 | 0046 |
---|
| 112 | 0047 |
---|
| 113 | 0048 <span class="comment">% --- Executes just before probe_calib is made visible.</span> |
---|
| 114 | 0049 <span class="comment">%INPUT:</span> |
---|
| 115 | 0050 <span class="comment">%handles: handles of the probe_calib interface elements</span> |
---|
| 116 | 0051 <span class="comment">% PlotHandles: set of handles of the elements contolling the plotting</span> |
---|
| 117 | 0052 <span class="comment">% parameters on the uvmat interface (obtained by 'get_plot_handle.m')</span> |
---|
| 118 | 0053 <a name="_sub1" href="#_subfunctions" class="code">function probe_calib_OpeningFcn(hObject, eventdata, handles, data,pos,inputfile)</a> |
---|
| 119 | 0054 |
---|
| 120 | 0055 <span class="comment">% Choose default command line output for probe_calib</span> |
---|
| 121 | 0056 handles.output = hObject; |
---|
| 122 | 0057 |
---|
| 123 | 0058 <span class="comment">% Update handles structure</span> |
---|
| 124 | 0059 guidata(hObject, handles); |
---|
| 125 | 0060 |
---|
| 126 | 0061 <span class="comment">%default</span> |
---|
| 127 | 0062 <span class="comment">% set(hObject,'Unit','Normalized')% set the unit normalized to the screen size</span> |
---|
| 128 | 0063 <span class="comment">% set(hObject,'Position',[0.7 0.1 0.25 0.5])%set the position of the probe_calib interface</span> |
---|
| 129 | 0064 set(hObject,<span class="string">'DeleteFcn'</span>,@<a href="#_sub6" class="code" title="subfunction closefcn(gcbo,eventdata)">closefcn</a>) |
---|
| 130 | 0065 |
---|
| 131 | 0066 <span class="comment">%set the position of the interface</span> |
---|
| 132 | 0067 <span class="keyword">if</span> exist(<span class="string">'pos'</span>,<span class="string">'var'</span>)& length(pos)>2 |
---|
| 133 | 0068 pos_gui=get(hObject,<span class="string">'Position'</span>); |
---|
| 134 | 0069 pos_gui(1)=pos(1); |
---|
| 135 | 0070 pos_gui(2)=pos(2); |
---|
| 136 | 0071 set(hObject,<span class="string">'Position'</span>,pos_gui); |
---|
| 137 | 0072 <span class="keyword">end</span> |
---|
| 138 | 0073 <span class="comment">% set(handles.XImage,'String','')</span> |
---|
| 139 | 0074 <span class="comment">% set(handles.YImage,'String','')</span> |
---|
| 140 | 0075 <span class="comment">% set(handles.XObject,'String','')</span> |
---|
| 141 | 0076 <span class="comment">% set(handles.YObject,'String','')</span> |
---|
| 142 | 0077 <span class="comment">% set(handles.ZObject,'String','')</span> |
---|
| 143 | 0078 inputxml=<span class="string">''</span>; |
---|
| 144 | 0079 <span class="keyword">if</span> exist(<span class="string">'inputfile'</span>,<span class="string">'var'</span>)& ~isempty(inputfile) |
---|
| 145 | 0080 [Path,Name,ext]=fileparts(inputfile); |
---|
| 146 | 0081 <span class="keyword">if</span> isequal(ext,<span class="string">'.png'</span>) |
---|
| 147 | 0082 set(hObject,<span class="string">'UserData'</span>,inputfile) |
---|
| 148 | 0083 [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); |
---|
| 149 | 0084 inputxml=[fullfile(Pathsub,RootFile) <span class="string">'.xml'</span>]; |
---|
| 150 | 0085 <span class="keyword">end</span> |
---|
| 151 | 0086 <span class="keyword">end</span> |
---|
| 152 | 0087 <span class="keyword">if</span> exist(inputxml,<span class="string">'file'</span>) |
---|
| 153 | 0088 <a href="#_sub5" class="code" title="subfunction loadfile(handles,fileinput)">loadfile</a>(handles,inputxml) |
---|
| 154 | 0089 <span class="keyword">end</span> |
---|
| 155 | 0090 set(handles.ListCoord,<span class="string">'KeyPressFcn'</span>,{@<a href="#_sub37" class="code" title="subfunction key_press_fcn(hObject,eventdata,handles)">key_press_fcn</a>,handles})<span class="comment">%set keyboard action function</span> |
---|
| 156 | 0091 |
---|
| 157 | 0092 |
---|
| 158 | 0093 <span class="comment">% --- Outputs from this function are returned to the command line.</span> |
---|
| 159 | 0094 <a name="_sub2" href="#_subfunctions" class="code">function varargout = probe_calib_OutputFcn(hObject, eventdata, handles)</a> |
---|
| 160 | 0095 <span class="comment">% varargout cell array for returning output args (see VARARGOUT);</span> |
---|
| 161 | 0096 <span class="comment">% hObject handle to figure</span> |
---|
| 162 | 0097 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 163 | 0098 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 164 | 0099 |
---|
| 165 | 0100 <span class="comment">% Get default command line output from handles structure</span> |
---|
| 166 | 0101 varargout{1} = handles.output; |
---|
| 167 | 0102 varargout{2}=handles; |
---|
| 168 | 0103 |
---|
| 169 | 0104 <span class="comment">%------------</span> |
---|
| 170 | 0105 <a name="_sub3" href="#_subfunctions" class="code">function Phi_Callback(hObject, eventdata, handles)</a> |
---|
| 171 | 0106 |
---|
| 172 | 0107 |
---|
| 173 | 0108 <span class="comment">%-----------------------------------------------------</span> |
---|
| 174 | 0109 <span class="comment">% --- Executes on button press in import.</span> |
---|
| 175 | 0110 <a name="_sub4" href="#_subfunctions" class="code">function import_Callback(hObject, eventdata, handles)</a> |
---|
| 176 | 0111 <span class="comment">%get the object file</span> |
---|
| 177 | 0112 huvmat=findobj(<span class="string">'Tag'</span>,<span class="string">'uvmat'</span>); |
---|
| 178 | 0113 UvData=get(huvmat,<span class="string">'UserData'</span>); |
---|
| 179 | 0114 hchild=get(huvmat,<span class="string">'Children'</span>); |
---|
| 180 | 0115 hrootpath=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'RootPath'</span>); |
---|
| 181 | 0116 oldfile=get(hrootpath,<span class="string">'String'</span>); |
---|
| 182 | 0117 <span class="keyword">if</span> isempty(oldfile) |
---|
| 183 | 0118 oldfile=<span class="string">''</span>; |
---|
| 184 | 0119 <span class="keyword">end</span> |
---|
| 185 | 0120 <span class="comment">%[FileName,PathName] = uigetfile('*.civ','Select a .civ file',oldfile)</span> |
---|
| 186 | 0121 [FileName, PathName, filterindex] = uigetfile( <span class="keyword">...</span> |
---|
| 187 | 0122 {<span class="string">'*.xml;*.mat'</span>, <span class="string">' (*.xml,*.mat)'</span>; |
---|
| 188 | 0123 <span class="string">'*.xml'</span>, <span class="string">'.xml files '</span>; <span class="keyword">...</span> |
---|
| 189 | 0124 <span class="string">'*.mat'</span>, <span class="string">'.mat matlab files '</span>}, <span class="keyword">...</span> |
---|
| 190 | 0125 <span class="string">'Pick a file'</span>,oldfile); |
---|
| 191 | 0126 fileinput=[PathName FileName];<span class="comment">%complete file name</span> |
---|
| 192 | 0127 testblank=findstr(fileinput,<span class="string">' '</span>);<span class="comment">%look for blanks</span> |
---|
| 193 | 0128 <span class="keyword">if</span> ~isempty(testblank) |
---|
| 194 | 0129 <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>) |
---|
| 195 | 0130 <span class="keyword">return</span> |
---|
| 196 | 0131 <span class="keyword">end</span> |
---|
| 197 | 0132 sizf=size(fileinput); |
---|
| 198 | 0133 <span class="keyword">if</span> (~ischar(fileinput)|~isequal(sizf(1),1)),<span class="keyword">return</span>;<span class="keyword">end</span> |
---|
| 199 | 0134 <a href="#_sub5" class="code" title="subfunction loadfile(handles,fileinput)">loadfile</a>(handles,fileinput) |
---|
| 200 | 0135 |
---|
| 201 | 0136 <span class="comment">%--------------------------------------------------</span> |
---|
| 202 | 0137 <span class="comment">%read input xml file and update the edit boxes</span> |
---|
| 203 | 0138 <a name="_sub5" href="#_subfunctions" class="code">function loadfile(handles,fileinput)</a> |
---|
| 204 | 0139 |
---|
| 205 | 0140 <span class="comment">%read the input xml file</span> |
---|
| 206 | 0141 t=xmltree(fileinput); |
---|
| 207 | 0142 s=convert(t);<span class="comment">%convert to matlab structure</span> |
---|
| 208 | 0143 |
---|
| 209 | 0144 <span class="comment">%read data currently displayed on the interface</span> |
---|
| 210 | 0145 PointCoord=[]; |
---|
| 211 | 0146 data=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles); |
---|
| 212 | 0147 Coord=[]; <span class="comment">%default</span> |
---|
| 213 | 0148 <span class="keyword">if</span> isfield(data,<span class="string">'Coord'</span>) |
---|
| 214 | 0149 Coord=data.Coord; |
---|
| 215 | 0150 <span class="keyword">end</span> |
---|
| 216 | 0151 TabChar_0=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
| 217 | 0152 nbcoord_0=size(TabChar_0,1); |
---|
| 218 | 0153 <span class="keyword">if</span> isequal(get(handles.edit_append,<span class="string">'Value'</span>),1) <span class="comment">%edit mode</span> |
---|
| 219 | 0154 val=get(handles.ListCoord,<span class="string">'Value'</span>)-1; |
---|
| 220 | 0155 <span class="keyword">else</span> |
---|
| 221 | 0156 val=length(TabChar_0); |
---|
| 222 | 0157 <span class="keyword">end</span> |
---|
| 223 | 0158 nbcoord=0; |
---|
| 224 | 0159 |
---|
| 225 | 0160 <span class="comment">%case of calibration (ImaDoc) input file</span> |
---|
| 226 | 0161 <span class="keyword">if</span> isfield(s,<span class="string">'GeometryCalib'</span>) |
---|
| 227 | 0162 Calib=s.GeometryCalib; |
---|
| 228 | 0163 <span class="keyword">if</span> isfield(Calib,<span class="string">'SourceCalib'</span>) |
---|
| 229 | 0164 <span class="keyword">if</span> isfield(Calib.SourceCalib,<span class="string">'PointCoord'</span>) |
---|
| 230 | 0165 PointCoord=Calib.SourceCalib.PointCoord; |
---|
| 231 | 0166 <span class="keyword">end</span> |
---|
| 232 | 0167 <span class="keyword">if</span> isfield(Calib.SourceCalib,<span class="string">'ImageCalib'</span>) |
---|
| 233 | 0168 hcalib=get(handles.import,<span class="string">'parent'</span>); |
---|
| 234 | 0169 set(hcalib,<span class="string">'UserData'</span>,Calib.SourceCalib.ImageCalib);<span class="comment">%store the source image name in the interface 'UserData'</span> |
---|
| 235 | 0170 <span class="keyword">end</span> |
---|
| 236 | 0171 <span class="keyword">end</span> |
---|
| 237 | 0172 nbcoord=length(PointCoord); |
---|
| 238 | 0173 <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> |
---|
| 239 | 0174 <span class="keyword">for</span> i=1:length(PointCoord) |
---|
| 240 | 0175 line=str2num(PointCoord{i}); |
---|
| 241 | 0176 Coord(i+val,4:5)=line(4:5);<span class="comment">%px x</span> |
---|
| 242 | 0177 Coord(i+val,1:3)=line(1:3)/10;<span class="comment">%phys x</span> |
---|
| 243 | 0178 <span class="keyword">end</span> |
---|
| 244 | 0179 <span class="keyword">else</span> |
---|
| 245 | 0180 <span class="keyword">for</span> i=1:length(PointCoord) |
---|
| 246 | 0181 line=str2num(PointCoord{i}); |
---|
| 247 | 0182 Coord(i,4:5)=line(4:5);<span class="comment">%px x</span> |
---|
| 248 | 0183 Coord(i,1:3)=line(1:3);<span class="comment">%phys x</span> |
---|
| 249 | 0184 <span class="keyword">end</span> |
---|
| 250 | 0185 <span class="keyword">end</span> |
---|
| 251 | 0186 <span class="keyword">end</span> |
---|
| 252 | 0187 |
---|
| 253 | 0188 <span class="comment">%case of xml files of points</span> |
---|
| 254 | 0189 <span class="keyword">if</span> isfield(s,<span class="string">'Coord'</span>) |
---|
| 255 | 0190 PointCoord=s.Coord; |
---|
| 256 | 0191 nbcoord=length(PointCoord); |
---|
| 257 | 0192 <span class="comment">%case of image coordinates</span> |
---|
| 258 | 0193 <span class="keyword">if</span> isfield(s,<span class="string">'CoordType'</span>)& isequal(s.CoordType,<span class="string">'px'</span>) |
---|
| 259 | 0194 <span class="keyword">for</span> i=1:nbcoord |
---|
| 260 | 0195 line=str2num(PointCoord{i}); |
---|
| 261 | 0196 Coord(i+val,4:5)=line(1:2); |
---|
| 262 | 0197 <span class="keyword">end</span> |
---|
| 263 | 0198 <span class="comment">%case of physical coordinates</span> |
---|
| 264 | 0199 <span class="keyword">else</span> |
---|
| 265 | 0200 <span class="keyword">for</span> i=1:nbcoord |
---|
| 266 | 0201 line=str2num(PointCoord{i}) |
---|
| 267 | 0202 Coord(i+val,1:3)=line(1:3); |
---|
| 268 | 0203 nbcolumn=size(Coord,2); |
---|
| 269 | 0204 <span class="keyword">if</span> nbcolumn<5 |
---|
| 270 | 0205 Coord(i+val,nbcolumn+1:5)=zeros(1,5-nbcolumn); |
---|
| 271 | 0206 <span class="keyword">end</span> |
---|
| 272 | 0207 <span class="keyword">end</span> |
---|
| 273 | 0208 <span class="keyword">end</span> |
---|
| 274 | 0209 <span class="keyword">end</span> |
---|
| 275 | 0210 CoordCell={}; |
---|
| 276 | 0211 <span class="keyword">for</span> iline=1:size(Coord,1) |
---|
| 277 | 0212 <span class="keyword">for</span> j=1:5 |
---|
| 278 | 0213 CoordCell{iline,j}=num2str(Coord(iline,j)); |
---|
| 279 | 0214 <span class="keyword">end</span> |
---|
| 280 | 0215 <span class="keyword">end</span> |
---|
| 281 | 0216 |
---|
| 282 | 0217 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> |
---|
| 283 | 0218 set(handles.ListCoord,<span class="string">'Value'</span>,1) |
---|
| 284 | 0219 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
| 285 | 0220 |
---|
| 286 | 0221 |
---|
| 287 | 0222 <span class="comment">%----------------------------------------------------</span> |
---|
| 288 | 0223 <span class="comment">% executed when closing: set the parent interface button to value 0</span> |
---|
| 289 | 0224 <a name="_sub6" href="#_subfunctions" class="code">function closefcn(gcbo,eventdata)</a> |
---|
| 290 | 0225 SetData=get(gcbf,<span class="string">'UserData'</span>); |
---|
| 291 | 0226 <span class="keyword">if</span> isfield(SetData,<span class="string">'ParentButton'</span>) & ishandle(SetData.ParentButton) |
---|
| 292 | 0227 set(SetData.ParentButton, <span class="string">'Value'</span>,0) |
---|
| 293 | 0228 <span class="keyword">end</span> |
---|
| 294 | 0229 |
---|
| 295 | 0230 <span class="comment">%-----------------------------------------------------------------------</span> |
---|
| 296 | 0231 <span class="comment">% --- Executes on button press in edit: PLOT the defined object and its projected field</span> |
---|
| 297 | 0232 <a name="_sub7" href="#_subfunctions" class="code">function edit_Callback(hObject, eventdata, handles)</a> |
---|
| 298 | 0233 <span class="comment">%hsetobject=get(hObject,'parent');</span> |
---|
| 299 | 0234 <span class="comment">%SetData=get(hsetobject,'UserData');%get the hidden interface data</span> |
---|
| 300 | 0235 <span class="comment">%IndexObj=SetData.IndexObj;%index of the current projection object in the list of projection objects (UvData.ProjObject)</span> |
---|
| 301 | 0236 huvmat=findobj(allchild(0),<span class="string">'name'</span>,<span class="string">'uvmat'</span>);<span class="comment">%find the current uvmat interface handle</span> |
---|
| 302 | 0237 UvData=get(huvmat,<span class="string">'UserData'</span>);<span class="comment">%Data associated to the current uvmat interface</span> |
---|
| 303 | 0238 hplot=findobj(huvmat,<span class="string">'Tag'</span>,<span class="string">'axes3'</span>); |
---|
| 304 | 0239 h_menu_coord=findobj(huvmat,<span class="string">'Tag'</span>,<span class="string">'menu_coord'</span>); |
---|
| 305 | 0240 menu=get(h_menu_coord,<span class="string">'String'</span>); |
---|
| 306 | 0241 choice=get(h_menu_coord,<span class="string">'Value'</span>); |
---|
| 307 | 0242 <span class="keyword">if</span> iscell(menu) |
---|
| 308 | 0243 option=menu{choice}; |
---|
| 309 | 0244 <span class="keyword">else</span> |
---|
| 310 | 0245 option=<span class="string">'px'</span>; <span class="comment">%default</span> |
---|
| 311 | 0246 <span class="keyword">end</span> |
---|
| 312 | 0247 <span class="comment">%get axis</span> |
---|
| 313 | 0248 <span class="comment">%get CoordType</span> |
---|
| 314 | 0249 ObjectData=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles);<span class="comment">%read the interface input parameters defining the object</span> |
---|
| 315 | 0250 <span class="keyword">if</span> isequal(option,<span class="string">'phys'</span>) |
---|
| 316 | 0251 ObjectData.Coord=ObjectData.Coord(:,[1:3]); |
---|
| 317 | 0252 <span class="keyword">elseif</span> isequal(option,<span class="string">'px'</span>) |
---|
| 318 | 0253 ObjectData.Coord=ObjectData.Coord(:,[4:5]); |
---|
| 319 | 0254 <span class="keyword">else</span> |
---|
| 320 | 0255 errordlg(<span class="string">'the choice in coord_coord must be px or phys '</span>) |
---|
| 321 | 0256 <span class="keyword">end</span> |
---|
| 322 | 0257 <span class="comment">% [UvData,IndexObj]=update_obj(UvData,IndexObj,ObjectData,SetData.PlotHandles);</span> |
---|
| 323 | 0258 <span class="comment">% SetData.IndexObj=IndexObj;</span> |
---|
| 324 | 0259 <span class="comment">% set(gcbf,'UserData',SetData)%update object index in the probe_calib interface</span> |
---|
| 325 | 0260 <span class="comment">% set(huvmat,'UserData',UvData)%update the data in the uvmat interface</span> |
---|
| 326 | 0261 ObjectData.ProjMode=<span class="string">'none'</span>; |
---|
| 327 | 0262 <a href="plot_object.html" class="code" title="function [hh]=plot_object(ObjectDataIn,ProjObject,hplot,col)">plot_object</a>(ObjectData,[],hplot,<span class="string">'b'</span>); |
---|
| 328 | 0263 |
---|
| 329 | 0264 |
---|
| 330 | 0265 <span class="comment">% --- Executes on button press in MenuCoord.</span> |
---|
| 331 | 0266 <a name="_sub8" href="#_subfunctions" class="code">function MenuCoord_Callback(hObject, eventdata, handles)</a> |
---|
| 332 | 0267 |
---|
| 333 | 0268 |
---|
| 334 | 0269 <span class="comment">% --- Executes on button press in delete.</span> |
---|
| 335 | 0270 <a name="_sub9" href="#_subfunctions" class="code">function delete_Callback(hObject, eventdata, handles)</a> |
---|
| 336 | 0271 SetData=get(gcbf,<span class="string">'UserData'</span>);<span class="comment">%get the interface data</span> |
---|
| 337 | 0272 IndexObj=SetData.IndexObj; |
---|
| 338 | 0273 <a href="delete_object.html" class="code" title="function delete_object(hObject)">delete_object</a>(IndexObj); |
---|
| 339 | 0274 |
---|
| 340 | 0275 <span class="comment">% --- Executes on button press in calibrate_lin.</span> |
---|
| 341 | 0276 <a name="_sub10" href="#_subfunctions" class="code">function calib_offset_Callback(hObject, eventdata, handles)</a> |
---|
| 342 | 0277 Object=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles); |
---|
| 343 | 0278 |
---|
| 344 | 0279 <span class="comment">%make linear calibration</span> |
---|
| 345 | 0280 <span class="comment">% 'calibration_lin' provides a linear transform on coordinates,</span> |
---|
| 346 | 0281 X=Object.Coord(:,1); |
---|
| 347 | 0282 Y=Object.Coord(:,2); |
---|
| 348 | 0283 x_ima=Object.Coord(:,4); |
---|
| 349 | 0284 y_ima=Object.Coord(:,5); |
---|
| 350 | 0285 [<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); |
---|
| 351 | 0286 [py,sy]=polyfit(Y,y_ima,1); |
---|
| 352 | 0287 <span class="comment">%err_X1=max(abs(x1-x_ima));%error</span> |
---|
| 353 | 0288 <span class="comment">%err_Y1=max(abs(y1-y_ima));%error</span> |
---|
| 354 | 0289 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); |
---|
| 355 | 0290 T_y=py(2); |
---|
| 356 | 0291 GeometryCalib.focal=1; |
---|
| 357 | 0292 GeometryCalib.Tx_Ty_Tz=[T_x T_y 1]; |
---|
| 358 | 0293 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]; |
---|
| 359 | 0294 <span class="comment">%check error</span> |
---|
| 360 | 0295 Calib.dpx=1; |
---|
| 361 | 0296 Calib.dpy=1; |
---|
| 362 | 0297 Calib.sx=1; |
---|
| 363 | 0298 Calib.Cx=0; |
---|
| 364 | 0299 Calib.Cy=0; |
---|
| 365 | 0300 Calib.Tz=1; |
---|
| 366 | 0301 Calib.kappa1=0; |
---|
| 367 | 0302 Calib.f=GeometryCalib.focal; |
---|
| 368 | 0303 Calib.Tx=T_x; |
---|
| 369 | 0304 Calib.Ty=T_y; |
---|
| 370 | 0305 Calib.R=GeometryCalib.R; |
---|
| 371 | 0306 [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); |
---|
| 372 | 0307 GeometryCalib.ErrorRms(1)=sqrt(mean((Xpoints-x_ima).*(Xpoints-x_ima))); |
---|
| 373 | 0308 GeometryCalib.ErrorMax(1)=max(abs(Xpoints-x_ima)); |
---|
| 374 | 0309 GeometryCalib.ErrorRms(2)=sqrt(mean((Ypoints-y_ima).*(Ypoints-y_ima))); |
---|
| 375 | 0310 GeometryCalib.ErrorMax(2)=max(abs(Ypoints-y_ima)) |
---|
| 376 | 0311 <span class="comment">%calibrate_lin calibration results and point coordinates</span> |
---|
| 377 | 0312 huvmat=findobj(<span class="string">'Tag'</span>,<span class="string">'uvmat'</span>); |
---|
| 378 | 0313 hchild=get(huvmat,<span class="string">'Children'</span>); |
---|
| 379 | 0314 hrootpath=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'RootPath'</span>); |
---|
| 380 | 0315 hrootfile=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'RootFile'</span>); |
---|
| 381 | 0316 RootPath=<span class="string">''</span>; |
---|
| 382 | 0317 RootFile=<span class="string">''</span>; |
---|
| 383 | 0318 <span class="keyword">if</span> ~isempty(hrootpath)& ~isempty(hrootfile) |
---|
| 384 | 0319 testhandle=1; |
---|
| 385 | 0320 RootPath=get(hrootpath,<span class="string">'String'</span>); |
---|
| 386 | 0321 RootFile=get(hrootfile,<span class="string">'String'</span>); |
---|
| 387 | 0322 filebase=fullfile(RootPath,RootFile); |
---|
| 388 | 0323 outputfile=[filebase <span class="string">'.xml'</span>] |
---|
| 389 | 0324 <span class="keyword">else</span> |
---|
| 390 | 0325 question={<span class="string">'save the calibration data and point coordinates in'</span>}; |
---|
| 391 | 0326 def={fullfile(RootPath,[<span class="string">'ObjectCalib.xml'</span>])}; |
---|
| 392 | 0327 options.Resize=<span class="string">'on'</span>; |
---|
| 393 | 0328 answer=inputdlg(question,<span class="string">'save average in a new file'</span>,1,def,options); |
---|
| 394 | 0329 outputfile=answer{1}; |
---|
| 395 | 0330 <span class="keyword">end</span> |
---|
| 396 | 0331 testappend=0; |
---|
| 397 | 0332 <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> |
---|
| 398 | 0333 t=xmltree(outputfile); <span class="comment">%read the file</span> |
---|
| 399 | 0334 uid=find(t,<span class="string">'ImaDoc'</span>); |
---|
| 400 | 0335 <span class="keyword">if</span> ~isequal(uid,1)<span class="comment">%if the xml file is not ImaDoc, delete it (after backup)</span> |
---|
| 401 | 0336 backupfile=outputfile; |
---|
| 402 | 0337 testexist=2; |
---|
| 403 | 0338 <span class="keyword">while</span> testexist==2 |
---|
| 404 | 0339 backupfile=[backupfile <span class="string">'~'</span>]; |
---|
| 405 | 0340 testexist=exist(backupfile,<span class="string">'file'</span>); |
---|
| 406 | 0341 <span class="keyword">end</span> |
---|
| 407 | 0342 [success,message]=copyfile(outputfile,backupfile);<span class="comment">%make backup</span> |
---|
| 408 | 0343 <span class="keyword">if</span> isequal(success,1) |
---|
| 409 | 0344 delete(outputfile) |
---|
| 410 | 0345 <span class="keyword">else</span> |
---|
| 411 | 0346 <span class="keyword">return</span> |
---|
| 412 | 0347 <span class="keyword">end</span> |
---|
| 413 | 0348 <span class="keyword">else</span> |
---|
| 414 | 0349 uid_calib=find(t,<span class="string">'ImaDoc/GeometryCalib'</span>); |
---|
| 415 | 0350 <span class="keyword">if</span> ~isempty(uid) <span class="comment">%if GeometryCalib already exists, delete its content</span> |
---|
| 416 | 0351 backupfile=outputfile; |
---|
| 417 | 0352 testexist=2; |
---|
| 418 | 0353 <span class="keyword">while</span> testexist==2 |
---|
| 419 | 0354 backupfile=[backupfile <span class="string">'~'</span>]; |
---|
| 420 | 0355 testexist=exist(backupfile,<span class="string">'file'</span>); |
---|
| 421 | 0356 <span class="keyword">end</span> |
---|
| 422 | 0357 [success,message]=copyfile(outputfile,backupfile)<span class="comment">%make backup</span> |
---|
| 423 | 0358 <span class="keyword">if</span> isequal(success,1) |
---|
| 424 | 0359 delete(outputfile) |
---|
| 425 | 0360 <span class="keyword">else</span> |
---|
| 426 | 0361 <span class="keyword">return</span> |
---|
| 427 | 0362 <span class="keyword">end</span> |
---|
| 428 | 0363 uid_child=children(t,uid_calib); |
---|
| 429 | 0364 t=delete(t,uid_child); |
---|
| 430 | 0365 testappend=1; |
---|
| 431 | 0366 <span class="keyword">end</span> |
---|
| 432 | 0367 <span class="keyword">end</span> |
---|
| 433 | 0368 <span class="keyword">end</span> |
---|
| 434 | 0369 <span class="keyword">if</span> ~testappend |
---|
| 435 | 0370 t=xmltree; |
---|
| 436 | 0371 t=set(t,1,<span class="string">'name'</span>,<span class="string">'ImaDoc'</span>); |
---|
| 437 | 0372 [t,uid_calib]=add(t,1,<span class="string">'element'</span>,<span class="string">'GeometryCalib'</span>); |
---|
| 438 | 0373 <span class="comment">% t=struct2xml(GeometryCalib,t,uid_calib);</span> |
---|
| 439 | 0374 <span class="keyword">end</span> |
---|
| 440 | 0375 Object.Coord(:,[1:3])=Object.Coord(:,[1:3])*10; <span class="comment">%transform in</span> |
---|
| 441 | 0376 GeometryCalib.SourceCalib.PointCoord=Object.Coord; |
---|
| 442 | 0377 t=<a href="struct2xml.html" class="code" title="function t=struct2xml(Object,t,root_uid)">struct2xml</a>(GeometryCalib,t,uid_calib); |
---|
| 443 | 0378 save(t,outputfile) |
---|
| 444 | 0379 |
---|
| 445 | 0380 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([outputfile <span class="string">'updated with linear calibration data'</span>],<span class="string">'CONFIRMATION'</span>) |
---|
| 446 | 0381 |
---|
| 447 | 0382 <span class="comment">%display image with new calibration in the currently opened uvmat interface</span> |
---|
| 448 | 0383 Indices=get(findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'FileIndex'</span>),<span class="string">'String'</span>); |
---|
| 449 | 0384 Ext=get(findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'FileExt'</span>),<span class="string">'String'</span>); |
---|
| 450 | 0385 imagename=[fullfile(RootPath,RootFile) Indices Ext]; |
---|
| 451 | 0386 <span class="comment">% input.menu_coord=1;</span> |
---|
| 452 | 0387 huvmat=<a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a>(imagename,1);<span class="comment">%open uvmat, set phys coord (Value 1)</span> |
---|
| 453 | 0388 |
---|
| 454 | 0389 |
---|
| 455 | 0390 |
---|
| 456 | 0391 <span class="comment">% --- Executes on button press in calibrate_lin.</span> |
---|
| 457 | 0392 <a name="_sub11" href="#_subfunctions" class="code">function calib_lin_Callback(hObject, eventdata, handles)</a> |
---|
| 458 | 0393 Object=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles); |
---|
| 459 | 0394 |
---|
| 460 | 0395 <span class="comment">%make linear calibration</span> |
---|
| 461 | 0396 <span class="comment">% 'calibration_lin' provides a linear transform on coordinates,</span> |
---|
| 462 | 0397 X=Object.Coord(:,1); |
---|
| 463 | 0398 Y=Object.Coord(:,2); |
---|
| 464 | 0399 x_ima=Object.Coord(:,4); |
---|
| 465 | 0400 y_ima=Object.Coord(:,5); |
---|
| 466 | 0401 XY_mat=[ones(size(X)) X Y]; |
---|
| 467 | 0402 a_X1=XY_mat\x_ima; <span class="comment">%transformation matrix for X</span> |
---|
| 468 | 0403 x1=XY_mat*a_X1;<span class="comment">%reconstruction</span> |
---|
| 469 | 0404 err_X1=max(abs(x1-x_ima));<span class="comment">%error</span> |
---|
| 470 | 0405 a_Y1=XY_mat\y_ima;<span class="comment">%transformation matrix for X</span> |
---|
| 471 | 0406 y1=XY_mat*a_Y1; |
---|
| 472 | 0407 err_Y1=max(abs(y1-y_ima));<span class="comment">%error</span> |
---|
| 473 | 0408 T_x=a_X1(1); |
---|
| 474 | 0409 T_y=a_Y1(1); |
---|
| 475 | 0410 GeometryCalib.focal=1; |
---|
| 476 | 0411 GeometryCalib.Tx_Ty_Tz=[T_x T_y 1]; |
---|
| 477 | 0412 GeometryCalib.R=[a_X1(2),a_X1(3),0;a_Y1(2),a_Y1(3),0;0,0,1]; |
---|
| 478 | 0413 |
---|
| 479 | 0414 <span class="comment">%check error</span> |
---|
| 480 | 0415 GeometryCalib.ErrorRms(1)=sqrt(mean((x1-x_ima).*(x1-x_ima))); |
---|
| 481 | 0416 GeometryCalib.ErrorMax(1)=max(abs(x1-x_ima)); |
---|
| 482 | 0417 GeometryCalib.ErrorRms(2)=sqrt(mean((y1-y_ima).*(y1-y_ima))); |
---|
| 483 | 0418 GeometryCalib.ErrorMax(2)=max(abs(y1-y_ima)) |
---|
| 484 | 0419 |
---|
| 485 | 0420 <span class="comment">%calibrate_lin calibration results and point coordinates</span> |
---|
| 486 | 0421 huvmat=findobj(<span class="string">'Tag'</span>,<span class="string">'uvmat'</span>); |
---|
| 487 | 0422 hchild=get(huvmat,<span class="string">'Children'</span>); |
---|
| 488 | 0423 hrootpath=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'RootPath'</span>); |
---|
| 489 | 0424 hrootfile=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'RootFile'</span>); |
---|
| 490 | 0425 RootPath=<span class="string">''</span>; |
---|
| 491 | 0426 RootFile=<span class="string">''</span>; |
---|
| 492 | 0427 <span class="keyword">if</span> ~isempty(hrootpath)& ~isempty(hrootfile) |
---|
| 493 | 0428 testhandle=1; |
---|
| 494 | 0429 RootPath=get(hrootpath,<span class="string">'String'</span>); |
---|
| 495 | 0430 RootFile=get(hrootfile,<span class="string">'String'</span>); |
---|
| 496 | 0431 filebase=fullfile(RootPath,RootFile); |
---|
| 497 | 0432 outputfile=[filebase <span class="string">'.xml'</span>]; |
---|
| 498 | 0433 <span class="keyword">else</span> |
---|
| 499 | 0434 question={<span class="string">'save the calibration data and point coordinates in'</span>}; |
---|
| 500 | 0435 def={fullfile(RootPath,[<span class="string">'ObjectCalib.xml'</span>])}; |
---|
| 501 | 0436 options.Resize=<span class="string">'on'</span>; |
---|
| 502 | 0437 answer=inputdlg(question,<span class="string">'save average in a new file'</span>,1,def,options); |
---|
| 503 | 0438 outputfile=answer{1}; |
---|
| 504 | 0439 <span class="keyword">end</span> |
---|
| 505 | 0440 testappend=0; |
---|
| 506 | 0441 <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> |
---|
| 507 | 0442 t=xmltree(outputfile); <span class="comment">%read the file</span> |
---|
| 508 | 0443 uid=find(t,<span class="string">'ImaDoc'</span>); |
---|
| 509 | 0444 <span class="keyword">if</span> ~isequal(uid,1)<span class="comment">%if the xml file is not ImaDoc, delete it (after backup)</span> |
---|
| 510 | 0445 backupfile=outputfile; |
---|
| 511 | 0446 testexist=2; |
---|
| 512 | 0447 <span class="keyword">while</span> testexist==2 |
---|
| 513 | 0448 backupfile=[backupfile <span class="string">'~'</span>]; |
---|
| 514 | 0449 testexist=exist(backupfile,<span class="string">'file'</span>); |
---|
| 515 | 0450 <span class="keyword">end</span> |
---|
| 516 | 0451 [success,message]=copyfile(outputfile,backupfile)<span class="comment">%make backup</span> |
---|
| 517 | 0452 <span class="keyword">if</span> isequal(success,1); |
---|
| 518 | 0453 delete(outputfile) |
---|
| 519 | 0454 <span class="keyword">else</span> |
---|
| 520 | 0455 <span class="keyword">return</span> |
---|
| 521 | 0456 <span class="keyword">end</span> |
---|
| 522 | 0457 <span class="keyword">else</span> |
---|
| 523 | 0458 uid_calib=find(t,<span class="string">'ImaDoc/GeometryCalib'</span>); |
---|
| 524 | 0459 <span class="keyword">if</span> ~isempty(uid) <span class="comment">%if GeometryCalib already exists, delete its content</span> |
---|
| 525 | 0460 backupfile=outputfile; |
---|
| 526 | 0461 testexist=2; |
---|
| 527 | 0462 <span class="keyword">while</span> testexist==2 |
---|
| 528 | 0463 backupfile=[backupfile <span class="string">'~'</span>]; |
---|
| 529 | 0464 testexist=exist(backupfile,<span class="string">'file'</span>); |
---|
| 530 | 0465 <span class="keyword">end</span> |
---|
| 531 | 0466 [success,message]=copyfile(outputfile,backupfile)<span class="comment">%make backup</span> |
---|
| 532 | 0467 <span class="keyword">if</span> isequal(success,1) |
---|
| 533 | 0468 delete(outputfile) |
---|
| 534 | 0469 <span class="keyword">else</span> |
---|
| 535 | 0470 <span class="keyword">return</span> |
---|
| 536 | 0471 <span class="keyword">end</span> |
---|
| 537 | 0472 uid_child=children(t,uid_calib); |
---|
| 538 | 0473 t=delete(t,uid_child); |
---|
| 539 | 0474 testappend=1; |
---|
| 540 | 0475 <span class="keyword">end</span> |
---|
| 541 | 0476 <span class="keyword">end</span> |
---|
| 542 | 0477 <span class="keyword">end</span> |
---|
| 543 | 0478 <span class="keyword">if</span> ~testappend |
---|
| 544 | 0479 t=xmltree; |
---|
| 545 | 0480 t=set(t,1,<span class="string">'name'</span>,<span class="string">'ImaDoc'</span>); |
---|
| 546 | 0481 [t,uid_calib]=add(t,1,<span class="string">'element'</span>,<span class="string">'GeometryCalib'</span>); |
---|
| 547 | 0482 <span class="comment">% t=struct2xml(GeometryCalib,t,uid_calib);</span> |
---|
| 548 | 0483 <span class="keyword">end</span> |
---|
| 549 | 0484 <span class="comment">% Object.Coord(:,[1:3])=Object.Coord(:,[1:3]); %transform in</span> |
---|
| 550 | 0485 GeometryCalib.SourceCalib.PointCoord=Object.Coord; |
---|
| 551 | 0486 t=<a href="struct2xml.html" class="code" title="function t=struct2xml(Object,t,root_uid)">struct2xml</a>(GeometryCalib,t,uid_calib); |
---|
| 552 | 0487 save(t,outputfile) |
---|
| 553 | 0488 |
---|
| 554 | 0489 <a href="warndlg_uvmat.html" class="code" title="function hwarn=warndlg_uvmat(warntext,title)">warndlg_uvmat</a>([outputfile <span class="string">'updated with linear calibration data'</span>],<span class="string">'CONFIRMATION'</span>) |
---|
| 555 | 0490 |
---|
| 556 | 0491 <span class="comment">%display image with new calibration in the currently opened uvmat interface</span> |
---|
| 557 | 0492 Indices=get(findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'FileIndex'</span>),<span class="string">'String'</span>); |
---|
| 558 | 0493 Ext=get(findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'FileExt'</span>),<span class="string">'String'</span>); |
---|
| 559 | 0494 imagename=[fullfile(RootPath,RootFile) Indices Ext]; |
---|
| 560 | 0495 <span class="comment">% input.menu_coord=1;</span> |
---|
| 561 | 0496 <span class="keyword">if</span> exist(imagename,<span class="string">'file'</span>) |
---|
| 562 | 0497 huvmat=<a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a>(imagename,1);<span class="comment">%open uvmat, set phys coord (Value 1)</span> |
---|
| 563 | 0498 <span class="keyword">else</span> |
---|
| 564 | 0499 huvmat=<a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a>; |
---|
| 565 | 0500 <span class="keyword">end</span> |
---|
| 566 | 0501 |
---|
| 567 | 0502 |
---|
| 568 | 0503 <span class="comment">% --- Executes on button press in translation.</span> |
---|
| 569 | 0504 <a name="_sub12" href="#_subfunctions" class="code">function translation_Callback(hObject, eventdata, handles)</a> |
---|
| 570 | 0505 |
---|
| 571 | 0506 |
---|
| 572 | 0507 <a name="_sub13" href="#_subfunctions" class="code">function T_x_Callback(hObject, eventdata, handles)</a> |
---|
| 573 | 0508 <span class="comment">% hObject handle to T_x (see GCBO)</span> |
---|
| 574 | 0509 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 575 | 0510 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 576 | 0511 |
---|
| 577 | 0512 <span class="comment">% Hints: get(hObject,'String') returns contents of T_x as text</span> |
---|
| 578 | 0513 <span class="comment">% str2double(get(hObject,'String')) returns contents of T_x as a double</span> |
---|
| 579 | 0514 |
---|
| 580 | 0515 |
---|
| 581 | 0516 |
---|
| 582 | 0517 |
---|
| 583 | 0518 |
---|
| 584 | 0519 <a name="_sub14" href="#_subfunctions" class="code">function T_y_Callback(hObject, eventdata, handles)</a> |
---|
| 585 | 0520 <span class="comment">% hObject handle to T_y (see GCBO)</span> |
---|
| 586 | 0521 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 587 | 0522 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 588 | 0523 |
---|
| 589 | 0524 <span class="comment">% Hints: get(hObject,'String') returns contents of T_y as text</span> |
---|
| 590 | 0525 <span class="comment">% str2double(get(hObject,'String')) returns contents of T_y as a double</span> |
---|
| 591 | 0526 |
---|
| 592 | 0527 |
---|
| 593 | 0528 <a name="_sub15" href="#_subfunctions" class="code">function T_z_Callback(hObject, eventdata, handles)</a> |
---|
| 594 | 0529 <span class="comment">% hObject handle to T_z (see GCBO)</span> |
---|
| 595 | 0530 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 596 | 0531 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 597 | 0532 |
---|
| 598 | 0533 <span class="comment">% Hints: get(hObject,'String') returns contents of T_z as text</span> |
---|
| 599 | 0534 <span class="comment">% str2double(get(hObject,'String')) returns contents of T_z as a double</span> |
---|
| 600 | 0535 |
---|
| 601 | 0536 |
---|
| 602 | 0537 <span class="comment">% --- Executes on button press in rotation.</span> |
---|
| 603 | 0538 <a name="_sub16" href="#_subfunctions" class="code">function rotation_Callback(hObject, eventdata, handles)</a> |
---|
| 604 | 0539 angle_rot=(pi/180)*str2num(get(handles.Phi,<span class="string">'String'</span>)) |
---|
| 605 | 0540 data=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles) |
---|
| 606 | 0541 data.Coord(:,1)=cos(angle_rot)*data.Coord(:,1)+sin(angle_rot)*data.Coord(:,2); |
---|
| 607 | 0542 data.Coord(:,1)=-sin(angle_rot)*data.Coord(:,1)+cos(angle_rot)*data.Coord(:,2); |
---|
| 608 | 0543 set(handles.XObject,<span class="string">'String'</span>,num2str(data.Coord(:,1))); |
---|
| 609 | 0544 set(handles.YObject,<span class="string">'String'</span>,num2str(data.Coord(:,2))); |
---|
| 610 | 0545 |
---|
| 611 | 0546 |
---|
| 612 | 0547 <a name="_sub17" href="#_subfunctions" class="code">function XImage_Callback(hObject, eventdata, handles)</a> |
---|
| 613 | 0548 <a href="#_sub22" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
| 614 | 0549 |
---|
| 615 | 0550 <a name="_sub18" href="#_subfunctions" class="code">function YImage_Callback(hObject, eventdata, handles)</a> |
---|
| 616 | 0551 <a href="#_sub22" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
| 617 | 0552 |
---|
| 618 | 0553 <a name="_sub19" href="#_subfunctions" class="code">function XObject_Callback(hObject, eventdata, handles)</a> |
---|
| 619 | 0554 <a href="#_sub22" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
| 620 | 0555 |
---|
| 621 | 0556 <a name="_sub20" href="#_subfunctions" class="code">function YObject_Callback(hObject, eventdata, handles)</a> |
---|
| 622 | 0557 <a href="#_sub22" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
| 623 | 0558 |
---|
| 624 | 0559 <a name="_sub21" href="#_subfunctions" class="code">function ZObject_Callback(hObject, eventdata, handles)</a> |
---|
| 625 | 0560 <a href="#_sub22" class="code" title="subfunction update_list(hObject, eventdata, handles)">update_list</a>(hObject, eventdata,handles) |
---|
| 626 | 0561 |
---|
| 627 | 0562 <a name="_sub22" href="#_subfunctions" class="code">function update_list(hObject, eventdata, handles)</a> |
---|
| 628 | 0563 str4=get(handles.XImage,<span class="string">'String'</span>); |
---|
| 629 | 0564 str5=get(handles.YImage,<span class="string">'String'</span>); |
---|
| 630 | 0565 str1=get(handles.XObject,<span class="string">'String'</span>); |
---|
| 631 | 0566 tt=double(str1); |
---|
| 632 | 0567 str2=get(handles.YObject,<span class="string">'String'</span>); |
---|
| 633 | 0568 str3=get(handles.ZObject,<span class="string">'String'</span>); |
---|
| 634 | 0569 <span class="keyword">if</span> ~isempty(str1) & ~isequal(double(str1),32) & (isempty(str3)|isequal(double(str3),32)) |
---|
| 635 | 0570 str3=<span class="string">'0'</span>;<span class="comment">%put z to 0 by default</span> |
---|
| 636 | 0571 <span class="keyword">end</span> |
---|
| 637 | 0572 strline=[str1 <span class="string">' | '</span> str2 <span class="string">' | '</span> str3 <span class="string">' | '</span> str4 <span class="string">' | '</span> str5]; |
---|
| 638 | 0573 Coord=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
| 639 | 0574 testappend=get(handles.edit_append,<span class="string">'Value'</span>); |
---|
| 640 | 0575 <span class="keyword">if</span> isequal(testappend,1); <span class="comment">%edit mode</span> |
---|
| 641 | 0576 val=get(handles.ListCoord,<span class="string">'Value'</span>); |
---|
| 642 | 0577 Coord{val}=strline; |
---|
| 643 | 0578 <span class="keyword">else</span> |
---|
| 644 | 0579 val=length(get(handles.ListCoord,<span class="string">'String'</span>)); |
---|
| 645 | 0580 Coord{val+1}=strline; |
---|
| 646 | 0581 set(handles.ListCoord,<span class="string">'Value'</span>,val+1) |
---|
| 647 | 0582 <span class="comment">% if val+1<=length(Coord)</span> |
---|
| 648 | 0583 <span class="comment">% set(handles.ListCoord,'Value',val+1)</span> |
---|
| 649 | 0584 <span class="comment">% ListCoord_Callback(hObject, eventdata, handles)</span> |
---|
| 650 | 0585 <span class="keyword">end</span> |
---|
| 651 | 0586 set(handles.ListCoord,<span class="string">'String'</span>,Coord) |
---|
| 652 | 0587 <span class="comment">%set(handles.ListCoord,'Value',val+1)</span> |
---|
| 653 | 0588 |
---|
| 654 | 0589 <span class="comment">%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</span> |
---|
| 655 | 0590 <a name="_sub23" href="#_subfunctions" class="code">function data=read_probe_calib(handles)</a> |
---|
| 656 | 0591 data_XIma=[]; |
---|
| 657 | 0592 data_YIma=[]; |
---|
| 658 | 0593 data_XObject=[]; |
---|
| 659 | 0594 data_YObject=[]; |
---|
| 660 | 0595 data_ZObject=[]; |
---|
| 661 | 0596 Coord=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
| 662 | 0597 <span class="comment">% XImage=get(handles.XImage,'String');</span> |
---|
| 663 | 0598 <span class="comment">% YImage=get(handles.YImage,'String');</span> |
---|
| 664 | 0599 <span class="comment">% XObject=get(handles.XObject,'String');</span> |
---|
| 665 | 0600 <span class="comment">% YObject=get(handles.YObject,'String');</span> |
---|
| 666 | 0601 <span class="comment">% ZObject=get(handles.ZObject,'String');</span> |
---|
| 667 | 0602 <span class="comment">% if ischar(Xcolumn)</span> |
---|
| 668 | 0603 <span class="comment">% Xcolumn={Xcolumn};</span> |
---|
| 669 | 0604 <span class="comment">% end</span> |
---|
| 670 | 0605 nb_defining_points=length(Coord); |
---|
| 671 | 0606 iline=0; |
---|
| 672 | 0607 <span class="keyword">for</span> i=1:nb_defining_points |
---|
| 673 | 0608 coord_str=Coord{i};<span class="comment">%character string of line number i</span> |
---|
| 674 | 0609 k=findstr(<span class="string">'|'</span>,coord_str);<span class="comment">%find separators '|'</span> |
---|
| 675 | 0610 data1=str2num(coord_str(1:k(1)-5)); |
---|
| 676 | 0611 data2=str2num(coord_str(k(1)+5:k(2)-5)); |
---|
| 677 | 0612 data3=str2num(coord_str(k(2)+5:k(3)-5)); |
---|
| 678 | 0613 data4=str2num(coord_str(k(3)+5:k(4)-5)); |
---|
| 679 | 0614 data5=str2num(coord_str(k(4)+5:end)); |
---|
| 680 | 0615 <span class="keyword">if</span> ~isempty(data1)|~isempty(data2)|~isempty(data3)|~isempty(data4)|~isempty(data5) |
---|
| 681 | 0616 iline=iline+1; |
---|
| 682 | 0617 <span class="keyword">if</span> ~isempty(data1) |
---|
| 683 | 0618 data.Coord(iline,1)=data1; |
---|
| 684 | 0619 <span class="keyword">end</span> |
---|
| 685 | 0620 <span class="keyword">if</span> ~isempty(data2) |
---|
| 686 | 0621 data.Coord(iline,2)=data2; |
---|
| 687 | 0622 <span class="keyword">end</span> |
---|
| 688 | 0623 <span class="keyword">if</span> ~isempty(data3) |
---|
| 689 | 0624 data.Coord(iline,3)=data3; |
---|
| 690 | 0625 <span class="keyword">end</span> |
---|
| 691 | 0626 <span class="keyword">if</span> ~isempty(data4) |
---|
| 692 | 0627 data.Coord(iline,4)=data4; |
---|
| 693 | 0628 <span class="keyword">end</span> |
---|
| 694 | 0629 <span class="keyword">if</span> isempty(data5) |
---|
| 695 | 0630 data.Coord(iline,5)=0; |
---|
| 696 | 0631 <span class="keyword">else</span> |
---|
| 697 | 0632 data.Coord(iline,5)=data5; |
---|
| 698 | 0633 <span class="keyword">end</span> |
---|
| 699 | 0634 <span class="keyword">end</span> |
---|
| 700 | 0635 <span class="keyword">end</span> |
---|
| 701 | 0636 data.Style=<span class="string">'points'</span>; |
---|
| 702 | 0637 |
---|
| 703 | 0638 |
---|
| 704 | 0639 <span class="comment">% --- Executes on selection change in ListCoord.</span> |
---|
| 705 | 0640 <a name="_sub24" href="#_subfunctions" class="code">function ListCoord_Callback(hObject, eventdata, handles)</a> |
---|
| 706 | 0641 <span class="comment">% hObject handle to ListCoord (see GCBO)</span> |
---|
| 707 | 0642 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 708 | 0643 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 709 | 0644 |
---|
| 710 | 0645 <span class="comment">% Hints: contents = get(hObject,'String') returns ListCoord contents as cell array</span> |
---|
| 711 | 0646 <span class="comment">% contents{get(hObject,'Value')} returns selected item from ListCoord</span> |
---|
| 712 | 0647 set(handles.edit_append,<span class="string">'Value'</span>,1); <span class="comment">%set to edit mode</span> |
---|
| 713 | 0648 Coord=get(handles.ListCoord,<span class="string">'String'</span>); |
---|
| 714 | 0649 val=get(handles.ListCoord,<span class="string">'Value'</span>); |
---|
| 715 | 0650 <span class="keyword">if</span> length(Coord)>0 |
---|
| 716 | 0651 coord_str=Coord{val}; |
---|
| 717 | 0652 k=findstr(<span class="string">'|'</span>,coord_str); |
---|
| 718 | 0653 set(handles.XObject,<span class="string">'String'</span>,coord_str(1:k(1)-5)) |
---|
| 719 | 0654 set(handles.YObject,<span class="string">'String'</span>,coord_str(k(1)+5:k(2)-5)) |
---|
| 720 | 0655 set(handles.ZObject,<span class="string">'String'</span>,coord_str(k(2)+5:k(3)-5)) |
---|
| 721 | 0656 set(handles.XImage,<span class="string">'String'</span>,coord_str(k(3)+5:k(4)-5)) |
---|
| 722 | 0657 set(handles.YImage,<span class="string">'String'</span>,coord_str(k(4)+5:end)) |
---|
| 723 | 0658 <span class="keyword">end</span> |
---|
| 724 | 0659 |
---|
| 725 | 0660 <span class="comment">%------------------------------------------------------</span> |
---|
| 726 | 0661 <span class="comment">% --- Executes on button press in translation_plus.</span> |
---|
| 727 | 0662 <a name="_sub25" href="#_subfunctions" class="code">function translation_plus_Callback(hObject, eventdata, handles)</a> |
---|
| 728 | 0663 |
---|
| 729 | 0664 T=[0 0 0]; |
---|
| 730 | 0665 T_x=get(handles.T_x,<span class="string">'String'</span>) |
---|
| 731 | 0666 T_y=get(handles.T_y,<span class="string">'String'</span>) |
---|
| 732 | 0667 T_z=get(handles.T_z,<span class="string">'String'</span>) |
---|
| 733 | 0668 <span class="keyword">if</span> ~isempty(T_x) |
---|
| 734 | 0669 T(1)=str2num(T_x); |
---|
| 735 | 0670 <span class="keyword">end</span> |
---|
| 736 | 0671 <span class="keyword">if</span> ~isempty(T_y) |
---|
| 737 | 0672 T(2)=str2num(T_y); |
---|
| 738 | 0673 <span class="keyword">end</span> |
---|
| 739 | 0674 <span class="keyword">if</span> ~isempty(T_z) |
---|
| 740 | 0675 T(3)=str2num(T_z); |
---|
| 741 | 0676 <span class="keyword">end</span> |
---|
| 742 | 0677 <a href="#_sub27" class="code" title="subfunction translation(handles,T)">translation</a>(handles,T) |
---|
| 743 | 0678 |
---|
| 744 | 0679 |
---|
| 745 | 0680 |
---|
| 746 | 0681 <span class="comment">% --- Executes on button press in translation_minus.</span> |
---|
| 747 | 0682 <a name="_sub26" href="#_subfunctions" class="code">function translation_minus_Callback(hObject, eventdata, handles)</a> |
---|
| 748 | 0683 |
---|
| 749 | 0684 T=[0 0 0]; |
---|
| 750 | 0685 T_x=get(handles.T_x,<span class="string">'String'</span>) |
---|
| 751 | 0686 T_y=get(handles.T_y,<span class="string">'String'</span>) |
---|
| 752 | 0687 T_z=get(handles.T_z,<span class="string">'String'</span>) |
---|
| 753 | 0688 <span class="keyword">if</span> ~isempty(T_x) |
---|
| 754 | 0689 T(1)=-str2num(T_x); |
---|
| 755 | 0690 <span class="keyword">end</span> |
---|
| 756 | 0691 <span class="keyword">if</span> ~isempty(T_y) |
---|
| 757 | 0692 T(2)=-str2num(T_y); |
---|
| 758 | 0693 <span class="keyword">end</span> |
---|
| 759 | 0694 <span class="keyword">if</span> ~isempty(T_z) |
---|
| 760 | 0695 T(3)=-str2num(T_z); |
---|
| 761 | 0696 <span class="keyword">end</span> |
---|
| 762 | 0697 <a href="#_sub27" class="code" title="subfunction translation(handles,T)">translation</a>(handles,T) |
---|
| 763 | 0698 |
---|
| 764 | 0699 |
---|
| 765 | 0700 <span class="comment">%%%--------------------------------------</span> |
---|
| 766 | 0701 <a name="_sub27" href="#_subfunctions" class="code">function translation(handles,T)</a> |
---|
| 767 | 0702 data=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles); |
---|
| 768 | 0703 data.Coord(:,1)=T(1)+data.Coord(:,1); |
---|
| 769 | 0704 data.Coord(:,2)=T(2)+data.Coord(:,2); |
---|
| 770 | 0705 data.Coord(:,3)=T(3)+data.Coord(:,3); |
---|
| 771 | 0706 data.Coord(:,[4 5])=data.Coord(:,[4 5]); |
---|
| 772 | 0707 <span class="keyword">for</span> i=1:size(data.Coord,1) |
---|
| 773 | 0708 <span class="keyword">for</span> j=1:5 |
---|
| 774 | 0709 Coord{i,j}=num2str(data.Coord(i,j));<span class="comment">%phys x,y,z</span> |
---|
| 775 | 0710 <span class="keyword">end</span> |
---|
| 776 | 0711 <span class="keyword">end</span> |
---|
| 777 | 0712 Tabchar=<a href="cell2tab.html" class="code" title="function Tabchar=cell2tab(Tabcell,separator)">cell2tab</a>(Coord,<span class="string">' | '</span>); |
---|
| 778 | 0713 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
| 779 | 0714 |
---|
| 780 | 0715 <span class="comment">%----------------------------------------------------</span> |
---|
| 781 | 0716 <span class="comment">% --- Executes on button press in rotation_plus.</span> |
---|
| 782 | 0717 <a name="_sub28" href="#_subfunctions" class="code">function rotation_plus_Callback(hObject, eventdata, handles)</a> |
---|
| 783 | 0718 Phi=0; |
---|
| 784 | 0719 Phi=get(handles.Phi,<span class="string">'String'</span>) |
---|
| 785 | 0720 <span class="keyword">if</span> ~isempty(Phi) |
---|
| 786 | 0721 Phi=str2num(Phi); |
---|
| 787 | 0722 <span class="keyword">end</span> |
---|
| 788 | 0723 <a href="#_sub30" class="code" title="subfunction rotation(handles,Phi)">rotation</a>(handles,Phi) |
---|
| 789 | 0724 |
---|
| 790 | 0725 <span class="comment">%-------------------------------------------------</span> |
---|
| 791 | 0726 <span class="comment">% --- Executes on button press in rotation_minus.</span> |
---|
| 792 | 0727 <a name="_sub29" href="#_subfunctions" class="code">function rotation_minus_Callback(hObject, eventdata, handles)</a> |
---|
| 793 | 0728 Phi=0; |
---|
| 794 | 0729 Phi=get(handles.Phi,<span class="string">'String'</span>) |
---|
| 795 | 0730 <span class="keyword">if</span> ~isempty(Phi) |
---|
| 796 | 0731 Phi=-str2num(Phi); |
---|
| 797 | 0732 <span class="keyword">end</span> |
---|
| 798 | 0733 <a href="#_sub30" class="code" title="subfunction rotation(handles,Phi)">rotation</a>(handles,Phi) |
---|
| 799 | 0734 |
---|
| 800 | 0735 <span class="comment">%-----------------------------------------------------</span> |
---|
| 801 | 0736 <span class="comment">%rotation</span> |
---|
| 802 | 0737 <a name="_sub30" href="#_subfunctions" class="code">function rotation(handles,Phi)</a> |
---|
| 803 | 0738 O_x=str2num(get(handles.O_x,<span class="string">'String'</span>)); |
---|
| 804 | 0739 O_y=str2num(get(handles.O_y,<span class="string">'String'</span>)); |
---|
| 805 | 0740 <span class="keyword">if</span> isempty(O_x) |
---|
| 806 | 0741 O_x=0;<span class="comment">%default</span> |
---|
| 807 | 0742 <span class="keyword">end</span> |
---|
| 808 | 0743 <span class="keyword">if</span> isempty(O_y) |
---|
| 809 | 0744 O_y=0;<span class="comment">%default</span> |
---|
| 810 | 0745 <span class="keyword">end</span> |
---|
| 811 | 0746 data=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles); |
---|
| 812 | 0747 r1=cos(pi*Phi/180); |
---|
| 813 | 0748 r2=-sin(pi*Phi/180); |
---|
| 814 | 0749 r3=sin(pi*Phi/180); |
---|
| 815 | 0750 r4=cos(pi*Phi/180); |
---|
| 816 | 0751 data.Coord(:,1)=r1*data.Coord(:,1)+r2*data.Coord(:,2); |
---|
| 817 | 0752 data.Coord(:,2)=r3*data.Coord(:,1)+r4*data.Coord(:,2); |
---|
| 818 | 0753 <span class="comment">% data.Coord(:,[4 5])=data.Coord(:,[4 5]);</span> |
---|
| 819 | 0754 <span class="keyword">for</span> i=1:size(data.Coord,1) |
---|
| 820 | 0755 <span class="keyword">for</span> j=1:5 |
---|
| 821 | 0756 Coord{i,j}=num2str(data.Coord(i,j));<span class="comment">%phys x,y,z</span> |
---|
| 822 | 0757 <span class="keyword">end</span> |
---|
| 823 | 0758 <span class="keyword">end</span> |
---|
| 824 | 0759 Tabchar=<a href="cell2tab.html" class="code" title="function Tabchar=cell2tab(Tabcell,separator)">cell2tab</a>(Coord,<span class="string">' | '</span>); |
---|
| 825 | 0760 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
| 826 | 0761 |
---|
| 827 | 0762 <a name="_sub31" href="#_subfunctions" class="code">function O_x_Callback(hObject, eventdata, handles)</a> |
---|
| 828 | 0763 <span class="comment">% hObject handle to O_x (see GCBO)</span> |
---|
| 829 | 0764 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 830 | 0765 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 831 | 0766 |
---|
| 832 | 0767 <span class="comment">% Hints: get(hObject,'String') returns contents of O_x as text</span> |
---|
| 833 | 0768 <span class="comment">% str2double(get(hObject,'String')) returns contents of O_x as a double</span> |
---|
| 834 | 0769 |
---|
| 835 | 0770 |
---|
| 836 | 0771 |
---|
| 837 | 0772 <a name="_sub32" href="#_subfunctions" class="code">function O_y_Callback(hObject, eventdata, handles)</a> |
---|
| 838 | 0773 <span class="comment">% hObject handle to O_y (see GCBO)</span> |
---|
| 839 | 0774 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 840 | 0775 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 841 | 0776 |
---|
| 842 | 0777 <span class="comment">% Hints: get(hObject,'String') returns contents of O_y as text</span> |
---|
| 843 | 0778 <span class="comment">% str2double(get(hObject,'String')) returns contents of O_y as a double</span> |
---|
| 844 | 0779 |
---|
| 845 | 0780 |
---|
| 846 | 0781 <a name="_sub33" href="#_subfunctions" class="code">function O_z_Callback(hObject, eventdata, handles)</a> |
---|
| 847 | 0782 <span class="comment">% hObject handle to O_z (see GCBO)</span> |
---|
| 848 | 0783 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 849 | 0784 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 850 | 0785 |
---|
| 851 | 0786 <span class="comment">% Hints: get(hObject,'String') returns contents of O_z as text</span> |
---|
| 852 | 0787 <span class="comment">% str2double(get(hObject,'String')) returns contents of O_z as a double</span> |
---|
| 853 | 0788 |
---|
| 854 | 0789 |
---|
| 855 | 0790 |
---|
| 856 | 0791 |
---|
| 857 | 0792 |
---|
| 858 | 0793 |
---|
| 859 | 0794 |
---|
| 860 | 0795 <span class="comment">% --- Executes on selection change in edit_append.</span> |
---|
| 861 | 0796 <a name="_sub34" href="#_subfunctions" class="code">function edit_append_Callback(hObject, eventdata, handles)</a> |
---|
| 862 | 0797 val=get(handles.edit_append,<span class="string">'Value'</span>); |
---|
| 863 | 0798 <span class="keyword">if</span> isequal(val,2); <span class="comment">%append mode</span> |
---|
| 864 | 0799 <span class="comment">%appeler mouse</span> |
---|
| 865 | 0800 <span class="keyword">end</span> |
---|
| 866 | 0801 |
---|
| 867 | 0802 |
---|
| 868 | 0803 <a name="_sub35" href="#_subfunctions" class="code">function NEW_Callback(hObject, eventdata, handles)</a> |
---|
| 869 | 0804 <span class="comment">%A METTRE SOUS UN BOUTON</span> |
---|
| 870 | 0805 huvmat=findobj(<span class="string">'name'</span>,<span class="string">'uvmat'</span>); |
---|
| 871 | 0806 hchild=get(huvmat,<span class="string">'children'</span>); |
---|
| 872 | 0807 hcoord=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'menu_coord'</span>) |
---|
| 873 | 0808 coordtype=get(hcoord,<span class="string">'Value'</span>) |
---|
| 874 | 0809 haxes=findobj(hchild,<span class="string">'Tag'</span>,<span class="string">'axes3'</span>); |
---|
| 875 | 0810 AxeData=get(haxes,<span class="string">'UserData'</span>); |
---|
| 876 | 0811 <span class="keyword">if</span> ~isequal(hcoord,2) |
---|
| 877 | 0812 set(hcoord,<span class="string">'Value'</span>,2) |
---|
| 878 | 0813 huvmat=<a href="uvmat.html" class="code" title="function varargout = uvmat(varargin)">uvmat</a>(AxeData) |
---|
| 879 | 0814 <span class="string">'relancer uvmat'</span> |
---|
| 880 | 0815 <span class="keyword">end</span> |
---|
| 881 | 0816 <span class="keyword">if</span> ~isfield(AxeData,<span class="string">'ZoomAxes'</span>) |
---|
| 882 | 0817 <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>) |
---|
| 883 | 0818 <span class="keyword">return</span> |
---|
| 884 | 0819 <span class="keyword">end</span> |
---|
| 885 | 0820 XLim=get(AxeData.ZoomAxes,<span class="string">'XLim'</span>); |
---|
| 886 | 0821 YLim=get(AxeData.ZoomAxes,<span class="string">'YLim'</span>); |
---|
| 887 | 0822 np=size(AxeData.A); |
---|
| 888 | 0823 ind_sub_x=round(XLim) |
---|
| 889 | 0824 ind_sub_y=np(1)-round(YLim) |
---|
| 890 | 0825 Mfiltre=AxeData.A([ind_sub_y(2):ind_sub_y(1)] ,ind_sub_x,:); |
---|
| 891 | 0826 Mfiltre_norm=double(Mfiltre); |
---|
| 892 | 0827 Mfiltre_norm=Mfiltre_norm/sum(sum(Mfiltre_norm)); |
---|
| 893 | 0828 Mfiltre_norm=100*(Mfiltre_norm-mean(mean(Mfiltre_norm))); |
---|
| 894 | 0829 Atype=class(AxeData.A) |
---|
| 895 | 0830 Data.NbDim=2; |
---|
| 896 | 0831 Data.A=filter2(Mfiltre_norm,double(AxeData.A)); |
---|
| 897 | 0832 Data.A=feval(Atype,Data.A); |
---|
| 898 | 0833 Data.AName=<span class="string">'image'</span>; |
---|
| 899 | 0834 Data.AX=AxeData.AX; |
---|
| 900 | 0835 Data.AY=AxeData.AY; |
---|
| 901 | 0836 Data.CoordType=<span class="string">'px'</span>; |
---|
| 902 | 0837 <a href="plot_field.html" class="code" title="function [PlotType,PlotParamOut,haxes]= plot_field(Data,haxes,PlotParam,KeepLim,PosColorbar)">plot_field</a>(Data) |
---|
| 903 | 0838 |
---|
| 904 | 0839 |
---|
| 905 | 0840 <span class="comment">% --- Executes on button press in HELP.</span> |
---|
| 906 | 0841 <a name="_sub36" href="#_subfunctions" class="code">function HELP_Callback(hObject, eventdata, handles)</a> |
---|
| 907 | 0842 <span class="comment">% hObject handle to HELP (see GCBO)</span> |
---|
| 908 | 0843 <span class="comment">% eventdata reserved - to be defined in a future version of MATLAB</span> |
---|
| 909 | 0844 <span class="comment">% handles structure with handles and user data (see GUIDATA)</span> |
---|
| 910 | 0845 path_to_uvmat=which (<span class="string">'uvmat'</span>);<span class="comment">% check the path of uvmat</span> |
---|
| 911 | 0846 pathelp=fileparts(path_to_uvmat); |
---|
| 912 | 0847 helpfile=fullfile(pathelp,<span class="string">'UVMAT_DOC'</span>,<span class="string">'uvmat_doc.html'</span>); |
---|
| 913 | 0848 <span class="keyword">if</span> isempty(dir(helpfile)), errordlg(<span class="string">'Please put the help file uvmat_doc.html in the directory UVMAT/UVMAT_DOC'</span>) |
---|
| 914 | 0849 <span class="keyword">else</span> |
---|
| 915 | 0850 web([helpfile <span class="string">'#probe_calib'</span>]) |
---|
| 916 | 0851 <span class="keyword">end</span> |
---|
| 917 | 0852 |
---|
| 918 | 0853 <span class="comment">%'move_key:' function activated when a key is pressed on the keyboard</span> |
---|
| 919 | 0854 <span class="comment">%-----------------------------------</span> |
---|
| 920 | 0855 <a name="_sub37" href="#_subfunctions" class="code">function key_press_fcn(hObject,eventdata,handles)</a> |
---|
| 921 | 0856 hh=get(hObject,<span class="string">'parent'</span>) |
---|
| 922 | 0857 xx=double(get(hh,<span class="string">'CurrentCharacter'</span>)) <span class="comment">%get the keyboard character</span> |
---|
| 923 | 0858 |
---|
| 924 | 0859 <span class="keyword">if</span> isequal(xx,8)<span class="comment">%move arrow right</span> |
---|
| 925 | 0860 data=<a href="#_sub23" class="code" title="subfunction data=read_probe_calib(handles)">read_probe_calib</a>(handles); |
---|
| 926 | 0861 Coord=[]; <span class="comment">%default</span> |
---|
| 927 | 0862 <span class="keyword">if</span> isfield(data,<span class="string">'Coord'</span>) |
---|
| 928 | 0863 Coord=data.Coord |
---|
| 929 | 0864 <span class="keyword">end</span> |
---|
| 930 | 0865 val=get(handles.ListCoord,<span class="string">'Value'</span>); |
---|
| 931 | 0866 Coord(val,:)=[]; |
---|
| 932 | 0867 CoordCell={}; |
---|
| 933 | 0868 <span class="keyword">for</span> iline=1:size(Coord,1) |
---|
| 934 | 0869 <span class="keyword">for</span> j=1:5 |
---|
| 935 | 0870 CoordCell{iline,j}=num2str(Coord(iline,j)); |
---|
| 936 | 0871 <span class="keyword">end</span> |
---|
| 937 | 0872 <span class="keyword">end</span> |
---|
| 938 | 0873 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> |
---|
| 939 | 0874 val=min(size(Coord,1),val); |
---|
| 940 | 0875 set(handles.ListCoord,<span class="string">'Value'</span>,max(val,1)) |
---|
| 941 | 0876 set(handles.ListCoord,<span class="string">'String'</span>,Tabchar) |
---|
| 942 | 0877 <span class="keyword">end</span></pre></div> |
---|
| 943 | <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> |
---|
| 944 | </body> |
---|
| 945 | </html> |
---|