&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;data.X(:,i+(j-1)*10)=&nbsp;x;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=eval>feval</span>(options.fun,data);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>else</span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=graph>figure</span>;&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;imshow(bmp,[]);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br><span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;figure(7);&nbsp;</span><br><span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;imshow(bmp,[]);</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=label>case</span>&nbsp;<span class=quotes>'alt'</span>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=comment>%&nbsp;right</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Cla</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;cells&nbsp;=&nbsp;getappdata(&nbsp;<span class=graph>gcf</span>,&nbsp;<span class=quotes>'cells'</span>&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x&nbsp;=&nbsp;<span class=graph>get</span>(&nbsp;<span class=graph>gca</span>,&nbsp;<span class=quotes>'currentpoint'</span>&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;[r&nbsp;c]&nbsp;=&nbsp;index(x([1&nbsp;3]));<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>if</span>&nbsp;~isempty(cells{r,c})&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=graph>set</span>(cells{r,c},&nbsp;<span class=quotes>'erasemode'</span>,<span class=quotes>'normal'</span>);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;delete(cells{r,c});<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;cells{r,c}&nbsp;=&nbsp;[];<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;setappdata(&nbsp;<span class=graph>gcf</span>,&nbsp;<span class=quotes>'cells'</span>,&nbsp;cells&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br><br><span class=defun_kw>function</span>&nbsp;<span class=defun_name>Cla</span>()<br>&nbsp;&nbsp;<span class=graph>cla</span>;<br>&nbsp;&nbsp;<span class=graph>plot</span>(&nbsp;[&nbsp;0&nbsp;0&nbsp;1&nbsp;1&nbsp;0&nbsp;],&nbsp;[&nbsp;0&nbsp;.5&nbsp;.5&nbsp;0&nbsp;0&nbsp;]&nbsp;);<br>&nbsp;&nbsp;hold&nbsp;on;<br>&nbsp;&nbsp;<span class=keyword>for</span>&nbsp;i&nbsp;=&nbsp;1:9,&nbsp;<span class=graph>plot</span>(&nbsp;[i/10&nbsp;i/10],[0&nbsp;.5]&nbsp;);&nbsp;<span class=keyword>end</span><br>&nbsp;&nbsp;<span class=keyword>for</span>&nbsp;i&nbsp;=&nbsp;1:4,&nbsp;<span class=graph>plot</span>(&nbsp;[0&nbsp;1],[i/10&nbsp;i/10]&nbsp;);&nbsp;&nbsp;<span class=keyword>end</span><br><br>&nbsp;&nbsp;<span class=graph>axis</span>&nbsp;equal;<br>&nbsp;&nbsp;<br>&nbsp;&nbsp;<span class=graph>set</span>(&nbsp;<span class=graph>gca</span>,&nbsp;<span class=quotes>'drawmode'</span>,&nbsp;<span class=quotes>'fast'</span>&nbsp;);<br>&nbsp;&nbsp;<span class=graph>set</span>(&nbsp;<span class=graph>gca</span>,&nbsp;<span class=quotes>'interruptible'</span>,&nbsp;<span class=quotes>'off'</span>&nbsp;);<br>&nbsp;&nbsp;<span class=graph>set</span>(&nbsp;<span class=graph>gca</span>,&nbsp;<span class=quotes>'xlimmode'</span>,&nbsp;<span class=quotes>'manual'</span>,&nbsp;<span class=quotes>'ylimmode'</span>,&nbsp;<span class=quotes>'manual'</span>,&nbsp;<span class=quotes>'zlimmode'</span>,&nbsp;<span class=quotes>'manual'</span>&nbsp;);<br><span class=comment>%&nbsp;axis&nbsp;off</span><br><br>&nbsp;&nbsp;title(<span class=quotes>'Control:&nbsp;left\_button&nbsp;-&nbsp;draw,&nbsp;middle\_button&nbsp;-&nbsp;classify,&nbsp;right\_button&nbsp;-&nbsp;erase.'</span>);<br><br><span class=defun_kw>function</span>&nbsp;<span class=defun_name>Plot</span>(<span class=defun_in>varargin</span>)<br>&nbsp;&nbsp;&nbsp;x&nbsp;=<span class=graph>get</span>(&nbsp;<span class=graph>gca</span>,&nbsp;<span class=quotes>'currentpoint'</span>&nbsp;);<br>&nbsp;&nbsp;&nbsp;<span class=keyword>if</span>(&nbsp;x(1)&nbsp;&gt;&nbsp;0&nbsp;&&nbsp;x(1)&nbsp;&lt;&nbsp;1&nbsp;&&nbsp;x(3)&nbsp;&gt;&nbsp;0&nbsp;&&nbsp;x(3)&nbsp;&lt;&nbsp;1&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;l&nbsp;=&nbsp;getappdata(<span class=graph>gcf</span>,&nbsp;<span class=quotes>'last'</span>);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>if</span>(&nbsp;isempty(&nbsp;l&nbsp;)&nbsp;),<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;l&nbsp;=&nbsp;<span class=graph>plot</span>(&nbsp;x(1),&nbsp;x(3),&nbsp;<span class=quotes>'.-'</span>&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=graph>set</span>(&nbsp;l,&nbsp;<span class=quotes>'erasemode'</span>,&nbsp;<span class=quotes>'none'</span>,&nbsp;<span class=quotes>'tag'</span>,&nbsp;<span class=quotes>'brush_stoke'</span>,&nbsp;<span class=quotes>'color'</span>,&nbsp;[0.5&nbsp;0&nbsp;0]&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;setappdata(<span class=graph>gcf</span>,&nbsp;<span class=quotes>'last'</span>,&nbsp;l&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>else</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;X&nbsp;=&nbsp;<span class=graph>get</span>(&nbsp;l,&nbsp;<span class=quotes>'xdata'</span>&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Y&nbsp;=&nbsp;<span class=graph>get</span>(&nbsp;l,&nbsp;<span class=quotes>'ydata'</span>&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=graph>set</span>(&nbsp;l,&nbsp;<span class=quotes>'xdata'</span>,&nbsp;[X&nbsp;x(1)],&nbsp;<span class=quotes>'ydata'</span>,&nbsp;[Y&nbsp;x(3)]&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br><span class=keyword>end</span><br><br><br><span class=defun_kw>function</span>&nbsp;<span class=defun_out>bmp&nbsp;</span>=&nbsp;<span class=defun_name>plot2bmp</span>(<span class=defun_in>&nbsp;handles&nbsp;</span>)<br>&nbsp;&nbsp;&nbsp;options=getappdata(&nbsp;<span class=graph>gcf</span>,&nbsp;<span class=quotes>'options'</span>);<br><br>&nbsp;&nbsp;&nbsp;Width&nbsp;=&nbsp;options.width*10;<br>&nbsp;&nbsp;&nbsp;Height&nbsp;=&nbsp;options.height*5;<br>&nbsp;&nbsp;&nbsp;bmp&nbsp;=&nbsp;zeros(&nbsp;Height,&nbsp;Width&nbsp;);<br>&nbsp;&nbsp;&nbsp;<br>&nbsp;&nbsp;&nbsp;<span class=keyword>for</span>&nbsp;i&nbsp;=&nbsp;1:length(handles&nbsp;),<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br><span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;X&nbsp;=&nbsp;get(&nbsp;handles(i),&nbsp;'xdata');</span><br><span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Y&nbsp;=&nbsp;get(&nbsp;handles(i),&nbsp;'ydata');
</span><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;points&nbsp;=&nbsp;<span class=graph>get</span>(&nbsp;handles(i),&nbsp;<span class=quotes>'Userdata'</span>);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;X&nbsp;=&nbsp;points.xdata;
<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Y&nbsp;=&nbsp;points.ydata;<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x1&nbsp;=&nbsp;min(fix(X(1)*Width)+1,Width);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y1&nbsp;=&nbsp;min(fix(2*Y(1)*Height)+1,Height);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>for</span>&nbsp;j=1:length(&nbsp;X&nbsp;)<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x2&nbsp;=&nbsp;min(fix(X(j)*Width)+1,Width);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y2&nbsp;=&nbsp;min(fix(2*Y(j)*Height)+1,Height);<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;n&nbsp;=&nbsp;max(&nbsp;ceil(&nbsp;max(&nbsp;abs(x2-x1),&nbsp;abs(y2-y1)&nbsp;)&nbsp;*&nbsp;2&nbsp;),&nbsp;1&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;a&nbsp;=&nbsp;[0:n]/n;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x&nbsp;=&nbsp;round(&nbsp;x1&nbsp;*&nbsp;a&nbsp;+&nbsp;x2&nbsp;*&nbsp;(1-a)&nbsp;);<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;y&nbsp;=&nbsp;Height&nbsp;-&nbsp;round(&nbsp;y1&nbsp;*&nbsp;a&nbsp;+&nbsp;y2&nbsp;*&nbsp;(1-a)&nbsp;)&nbsp;+&nbsp;1;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;bmp(&nbsp;y&nbsp;+&nbsp;(x&nbsp;-&nbsp;1)&nbsp;*&nbsp;Height&nbsp;)&nbsp;=&nbsp;1;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;x1=x2;&nbsp;y1=y2;<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br>&nbsp;&nbsp;&nbsp;<span class=keyword>end</span><br>&nbsp;&nbsp;<br><span class=defun_kw>function</span>&nbsp;<span class=defun_name>normalize</span>(<span class=defun_in>&nbsp;corner,&nbsp;sz,&nbsp;h&nbsp;</span>)
<br>&nbsp;&nbsp;
<br>x&nbsp;=&nbsp;<span class=graph>get</span>(&nbsp;h,&nbsp;<span class=quotes>'xdata'</span>&nbsp;);<br>y&nbsp;=&nbsp;<span class=graph>get</span>(&nbsp;h,&nbsp;<span class=quotes>'ydata'</span>&nbsp;);<br><span class=keyword>if</span>(&nbsp;iscell(x)&nbsp;),&nbsp;x&nbsp;=&nbsp;[x{:}];&nbsp;<span class=keyword>end</span><br><span class=keyword>if</span>(&nbsp;iscell(y)&nbsp;),&nbsp;y&nbsp;=&nbsp;[y{:}];&nbsp;<span class=keyword>end</span><br><br>mx&nbsp;=&nbsp;min(&nbsp;x&nbsp;);<br>Mx&nbsp;=&nbsp;max(&nbsp;x&nbsp;);<br>sx&nbsp;=&nbsp;Mx&nbsp;-&nbsp;mx;<br>my&nbsp;=&nbsp;min(&nbsp;y&nbsp;);<br>My&nbsp;=&nbsp;max(&nbsp;y&nbsp;);<br>sy&nbsp;=&nbsp;My&nbsp;-&nbsp;my;<br><br>centerx&nbsp;=&nbsp;(mx&nbsp;+&nbsp;Mx)&nbsp;/&nbsp;2;<br>centery&nbsp;=&nbsp;(my&nbsp;+&nbsp;My)&nbsp;/&nbsp;2;<br>center&nbsp;=&nbsp;corner&nbsp;+&nbsp;sz/2;<br><br><span class=keyword>if</span>(&nbsp;sy/sx&nbsp;&gt;&nbsp;&nbsp;sz(1)/sz(2)&nbsp;)<br>&nbsp;&nbsp;scale&nbsp;=&nbsp;sz(1)&nbsp;/&nbsp;sy;<br><span class=keyword>else</span><br>&nbsp;&nbsp;scale&nbsp;=&nbsp;sz(2)&nbsp;/&nbsp;sx;<br><span class=keyword>end</span><br>
<br><span class=keyword>for</span>&nbsp;hnd&nbsp;=&nbsp;h<br><span class=comment>%&nbsp;&nbsp;set(&nbsp;hnd,&nbsp;'erasemode',&nbsp;'normal'&nbsp;);</span><br><span class=comment>%&nbsp;&nbsp;set(&nbsp;hnd,&nbsp;'xdata',&nbsp;...</span><br><span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;(&nbsp;get(&nbsp;hnd,&nbsp;'xdata'&nbsp;)&nbsp;-&nbsp;centerx&nbsp;)&nbsp;*&nbsp;scale&nbsp;+&nbsp;center(2),&nbsp;...</span><br><span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;'ydata',&nbsp;...</span><br><span class=comment>%&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;(&nbsp;get(&nbsp;hnd,&nbsp;'ydata'&nbsp;)&nbsp;-&nbsp;centery&nbsp;)&nbsp;*&nbsp;scale&nbsp;+&nbsp;center(1)&nbsp;);
</span><br>
<br>&nbsp;&nbsp;points.xdata&nbsp;=&nbsp;(&nbsp;<span class=graph>get</span>(&nbsp;hnd,&nbsp;<span class=quotes>'xdata'</span>&nbsp;)&nbsp;-&nbsp;centerx&nbsp;)&nbsp;*&nbsp;scale&nbsp;+&nbsp;center(2);
<br>&nbsp;&nbsp;points.ydata&nbsp;=&nbsp;(&nbsp;<span class=graph>get</span>(&nbsp;hnd,&nbsp;<span class=quotes>'ydata'</span>&nbsp;)&nbsp;-&nbsp;centery&nbsp;)&nbsp;*&nbsp;scale&nbsp;+&nbsp;center(1);
<br>&nbsp;&nbsp;<span class=graph>set</span>(&nbsp;hnd,&nbsp;<span class=quotes>'userdata'</span>,&nbsp;points&nbsp;);<br><span class=keyword>end</span><br>
<br>
<br><span class=defun_kw>function</span>&nbsp;<span class=defun_out>[r,&nbsp;c]&nbsp;</span>=&nbsp;<span class=defun_name>index</span>(<span class=defun_in>&nbsp;x&nbsp;</span>)<br>r&nbsp;=&nbsp;min(&nbsp;floor(&nbsp;x(2)&nbsp;*&nbsp;10&nbsp;)&nbsp;+&nbsp;1,&nbsp;5&nbsp;);<br>c&nbsp;=&nbsp;min(&nbsp;floor(&nbsp;x(1)&nbsp;*&nbsp;10&nbsp;)&nbsp;+&nbsp;1,&nbsp;10&nbsp;);<br></code>