source: trunk/src/plot_object.m @ 965

Last change on this file since 965 was 955, checked in by sommeria, 8 years ago

various

File size: 22.2 KB
Line 
1%'plot_object': draws a projection object (points, line, plane...)
2%-------------------------------------------------------------------
3% function hh=plot_object(ObjectData,ProjObject,hplot,col)
4%
5%OUTPUT
6%             hh: handles of the graphic object (core part)
7%
8%INPUT:
9%
10% ObjectData: structure representing the object properties:
11%        .Type : style of projection object
12%        .Coord: set of coordinates defining the object position;
13%        .ProjMode=type of projection ;
14%       .ProjAngle=angle of projection;
15%       .DX,.DY,.DZ=increments;
16%       .YMax,YMin: min and max Y
17% ProjObject: projection object corresponding to the current plot (e. g. plane)
18% hplot: handle of the object plot to modify or if it is an axis, the axis
19%            where the object must be plotted, or if it is a figure the plotting figure
20% col: color of the plot, e;g; 'm', 'b' ..;
21
22%=======================================================================
23% Copyright 2008-2016, LEGI UMR 5519 / CNRS UGA G-INP, Grenoble, France
24%   http://www.legi.grenoble-inp.fr
25%   Joel.Sommeria - Joel.Sommeria (A) legi.cnrs.fr
26%
27%     This file is part of the toolbox UVMAT.
28%
29%     UVMAT is free software; you can redistribute it and/or modify
30%     it under the terms of the GNU General Public License as published
31%     by the Free Software Foundation; either version 2 of the license,
32%     or (at your option) any later version.
33%
34%     UVMAT is distributed in the hope that it will be useful,
35%     but WITHOUT ANY WARRANTY; without even the implied warranty of
36%     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
37%     GNU General Public License (see LICENSE.txt) for more details.
38%=======================================================================
39
40function [hh]=plot_object(ObjectData,ProjObject,hplot,col)
41
42%% default output
43hh=[];%default output
44% object representation is canceled if the field is not projected on a plane or is the same as the represented object
45if ~isfield(ObjectData,'Type')|| isequal(ProjObject,ObjectData)||~isfield(ProjObject,'Type')|| ~strcmp(ProjObject.Type,'plane')
46    if ~isempty(hplot) && ishandle(hplot) && ~strcmp(get(hplot,'Type'),'axes')
47        ObjectPlotData=get(hplot,'UserData');
48        if isfield(ObjectPlotData,'SubObject') & ishandle(ObjectPlotData.SubObject)
49            delete(ObjectPlotData.SubObject);
50        end
51        if isfield(ObjectPlotData,'DeformPoint') & ishandle(ObjectPlotData.DeformPoint)
52            delete(ObjectPlotData.DeformPoint);
53        end
54        delete(hplot)
55    end
56    return
57end
58XMin=0;%default range for the graph
59XMax=0;
60YMin=0;
61YMax=0;
62ZMin=0;
63ZMax=0;
64XMinRange=[];%default range set by input
65XMaxRange=[];
66YMinRange=[];
67YMaxRange=[];
68ZMinRange=[];
69ZMaxRange=[];
70
71%% determine the plotting axes (with handle 'haxes')
72test_newobj=1;
73if ishandle(hplot)
74    if isequal(get(hplot,'Tag'),'proj_object')% hplot is the handle of an object representation 
75        test_newobj=0;
76        haxes=get(hplot,'parent');
77        currentfig=get(haxes,'parent');
78    elseif isequal(get(hplot,'Type'),'axes')% hplot is the handle of an axis
79        haxes=hplot;
80        currentfig=get(hplot,'parent');
81    elseif isequal(get(hplot,'Type'),'figure')% hplot is the handle of a figure
82        set(0,'CurrentFigure',hplot);%set the input figure as the current one
83        haxes=findobj(hplot,'Type','axes');%look for axes in the figure
84        haxes=haxes(1);
85        currentfig=hplot;
86    else
87        currentfig=figure; %create new figure
88        hplot=axes;%create new axes
89        haxes=hplot;
90    end
91else
92    currentfig=figure; %create new figure
93    hplot=axes;%create new axes
94    haxes=hplot;
95end
96set(0,'CurrentFigure',currentfig)%set the currentfigure as the current one
97set(currentfig,'CurrentAxes',haxes);%set the current axes in the current figure
98
99%% default input parameters
100if ~isfield(ObjectData,'ProjMode')||isempty(ObjectData.ProjMode)
101     ObjectData.ProjMode='projection';%default
102end
103if ~isfield(ObjectData,'Coord')||isempty(ObjectData.Coord)
104     ObjectData.Coord=[0 0 0];%default
105end
106if isfield(ObjectData,'RangeX') && ~isempty(ObjectData.RangeX)
107    XMax=max(ObjectData.RangeX);
108    XMin=min(ObjectData.RangeX);
109        XMaxRange=max(ObjectData.RangeX);
110        if numel(ObjectData.RangeX)==2
111    XMinRange=min(ObjectData.RangeX);
112        end
113end
114if isfield(ObjectData,'RangeY')&&~isempty(ObjectData.RangeY)
115    YMax=max(ObjectData.RangeY);
116    YMin=min(ObjectData.RangeY);
117        YMaxRange=max(ObjectData.RangeY);
118        if numel(ObjectData.RangeY)==2
119    YMinRange=min(ObjectData.RangeY);
120        end
121end
122if isfield(ObjectData,'RangeZ')&&~isempty(ObjectData.RangeZ)
123    ZMax=max(ObjectData.RangeZ);
124    ZMin=min(ObjectData.RangeZ);
125    ZMaxRange=max(ObjectData.RangeZ);
126    ZMinRange=min(ObjectData.RangeZ);
127end
128if strcmp(ObjectData.Type,'points') && strcmp(ObjectData.ProjMode,'projection')
129    YMax=max(XMax,YMax);
130    YMax=max(YMax,ZMax);
131end
132sizcoord=size(ObjectData.Coord);
133
134%% determine the coordinates xline, yline,xsup,xinf, yinf,ysup determining the new object plot
135%test_line= isequal(ObjectData.Type,'points')||isequal(ObjectData.Type,'line')||...
136%   isequal(ObjectData.Type,'polyline')||isequal(ObjectData.Type,'polygon')|| isequal(ObjectData.Type,'plane')|| isequal(ObjectData.Type,'volume');
137%test_patch=isequal(ObjectData.ProjMode,'inside')||isequal(ObjectData.ProjMode,'outside')||isequal(ObjectData.Type,'volume')...
138%    ||isequal(ObjectData.ProjMode,'mask_inside')||isequal(ObjectData.ProjMode,'mask_outside');
139test_line=ismember(ObjectData.Type,{'points','line','polyline','polygon','plane','plane_z','volume'});
140test_patch=ismember(ObjectData.ProjMode,{'inside','outside','mask_inside','mask_outside'});
141if test_line
142    xline=ObjectData.Coord(:,1);
143    yline=ObjectData.Coord(:,2);
144    nbpoints=numel(xline);
145    switch ObjectData.Type
146        case 'line_x'
147            xline=[xline; ObjectData.RangeX(2)];%creating the line
148            yline=[yline; ObjectData.RangeY(2)];%creating the line
149        case 'polygon'
150            xline=[xline; ObjectData.Coord(1,1)];%closing the line
151            yline=[yline; ObjectData.Coord(1,2)];
152        case {'plane','plane_z','volume'}
153            if ~isfield(ObjectData,'Angle')
154                ObjectData.Angle=[0 0 0];
155            end
156            if strcmp(ObjectData.Type,'plane_z')
157                Angle_1=ObjectData.Angle(1);
158                Angle_2=ObjectData.Angle(2);
159                norm_angle=sqrt(Angle_1*Angle_1+Angle_2*Angle_2);
160                cosphi=Angle_1/norm_angle;%angle in radians
161                sinphi=Angle_2/norm_angle;%angle in radians
162            else
163                cosphi=cos(ObjectData.Angle(3)*pi/180);%angle in radians
164                sinphi=sin(ObjectData.Angle(3)*pi/180);%angle in radians
165            end
166            x0=xline(1); y0=yline(1);
167            xlim=get(haxes,'XLim');
168            ylim=get(haxes,'YLim');
169            graph_scale=max(abs(xlim(2)-xlim(1)),abs(ylim(2)-ylim(1)))/2;% estimate the length of axes plots
170            XMax=graph_scale;
171            YMax=graph_scale;
172            XMin=-graph_scale;
173            YMin=-graph_scale;
174            if  ~isempty(XMaxRange)
175                XMax=XMaxRange;
176            end
177            if  ~isempty(XMinRange)
178                XMin=XMinRange;
179            end
180            if  ~isempty(YMaxRange)
181                YMax=YMaxRange;
182            end
183            if  ~isempty(YMinRange)
184                YMin=YMinRange;
185            end
186            % axes lines
187            xline=NaN(1,13);
188            xline(1)=x0+min(0,XMin)*cosphi; % min end of the x axes
189            yline(1)=y0+min(0,XMin)*sinphi;
190            xline(2)=x0+XMax*cosphi;% max end of the x axes
191            yline(2)=y0+XMax*sinphi;
192            xline(8)=x0-min(0,YMin)*sinphi;% min end of the y axes
193            yline(8)=y0+min(0,YMin)*cosphi;
194            xline(9)=x0-YMax*sinphi;% max end of the y axes
195            yline(9)=y0+YMax*cosphi;
196           
197            %arrows on x axis
198            arrow_scale=graph_scale/20;
199            phi=acos(cosphi);
200            xline(3)=xline(2)-arrow_scale*cos(phi-pi/8);
201            yline(3)=yline(2)-arrow_scale*sin(phi-pi/8);
202            xline(5)=xline(2);
203            yline(5)=yline(2);
204            xline(6)=xline(2)-arrow_scale*cos(phi+pi/8);
205            yline(6)=yline(2)-arrow_scale*sin(phi+pi/8);
206           
207            %arrows on y axis
208            xline(10)=xline(9)-arrow_scale*cos(phi+pi/2-pi/8);
209            yline(10)=yline(9)-arrow_scale*sin(phi+pi/2-pi/8);
210            xline(12)=xline(9);
211            yline(12)=yline(9);
212            xline(13)=xline(9)-arrow_scale*cos(phi+pi/2+pi/8);
213            yline(13)=yline(9)-arrow_scale*sin(phi+pi/2+pi/8);
214            %xline=[Xbeg_x Xend_x NaN Ybeg_x Yend_x];
215            %yline=[Xbeg_y Xend_y NaN Ybeg_y Yend_y];
216            %  dashed lines indicating bounds
217            xsup=NaN(1,5);
218            ysup=NaN(1,5);
219            if ~isempty(XMaxRange)
220                xsup(1)=xline(2)-YMin*sin(phi);
221                ysup(1)=yline(2)+YMin*cos(phi);
222                xsup(2)=xline(2)-YMax*sin(phi);
223                ysup(2)=yline(2)+YMax*cos(phi);
224            end
225            if ~isempty(YMaxRange)
226                xsup(2)=xline(2)-YMax*sin(phi);
227                ysup(2)=yline(2)+YMax*cos(phi);
228                xsup(3)=xline(9)+XMin*cos(phi);
229                ysup(3)=yline(9)+XMin*sin(phi);
230            end
231            if ~isempty(XMinRange)
232                xsup(3)=xline(9)+XMin*cos(phi);
233                ysup(3)=yline(9)+XMin*sin(phi);
234                xsup(4)=x0+XMin*cos(phi)-YMin*sin(phi);
235                ysup(4)=y0+XMin*sin(phi)+YMin*cos(phi);
236            end
237            if ~isempty(YMinRange)
238                xsup(4)=x0+XMin*cos(phi)-YMin*sin(phi);
239                ysup(4)=y0+XMin*sin(phi)+YMin*cos(phi);
240                xsup(5)=xline(8)-YMin*sin(phi);
241                ysup(5)=yline(8)+YMin*cos(phi);
242            end
243    end
244end
245SubLineStyle='none';%default
246if isfield(ObjectData,'ProjMode')
247    if isequal(ObjectData.ProjMode,'projection')
248        SubLineStyle='--'; %range of projection marked by dash
249        if isfield (ObjectData,'DX')
250            ObjectData=rmfield(ObjectData,'DX');
251        end
252        if isfield (ObjectData,'DY')
253            ObjectData=rmfield(ObjectData,'DY');
254        end
255    elseif isequal(ObjectData.ProjMode,'filter')
256        SubLineStyle=':';%range of projection not visible
257    end
258end
259if ismember(ObjectData.Type,{'line','polyline','polygon'})
260    if length(xline)<2
261        theta=0;
262    else
263        theta=angle(diff(xline)+1i*diff(yline));
264        theta(length(xline))=theta(length(xline)-1);
265    end
266    xsup(1)=xline(1)+YMax*sin(theta(1));
267    xinf(1)=xline(1)-YMax*sin(theta(1));
268    ysup(1)=yline(1)-YMax*cos(theta(1));
269    yinf(1)=yline(1)+YMax*cos(theta(1));
270    for ip=2:length(xline)
271        xsup(ip)=xline(ip)+YMax*sin((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
272        xinf(ip)=xline(ip)-YMax*sin((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
273        ysup(ip)=yline(ip)-YMax*cos((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
274        yinf(ip)=yline(ip)+YMax*cos((theta(ip)+theta(ip-1))/2)/cos((theta(ip-1)-theta(ip))/2);
275    end
276end
277
278%% shading image
279if test_patch
280    npMx=512;
281    npMy=512; 
282    flag=zeros(npMy,npMx);
283    if isequal(ObjectData.Type,'ellipse')
284        XimaMin=ObjectData.Coord(1,1)-XMax;
285        XimaMax=ObjectData.Coord(1,1)+XMax;
286        YimaMin=ObjectData.Coord(1,2)-YMax;
287        YimaMax=ObjectData.Coord(1,2)+YMax;
288        xlim=[1.2*XimaMin-0.2*XimaMax 1.2*XimaMax-0.2*XimaMin];%create an image around the ellipse
289        ylim=[1.2*YimaMin-0.2*YimaMax 1.2*YimaMax-0.2*YimaMin];
290        scale_x=2*1.4*XMax/npMx;
291        scale_y=2*1.4*YMax/npMy;
292        xi=(0.5:npMx-0.5)*scale_x+xlim(1);
293        yi=(0.5:npMy-0.5)*scale_y+ylim(1);
294        [Xi,Yi]=meshgrid(xi,yi);
295        X2Max=XMax*XMax;
296        Y2Max=YMax*YMax;
297        distX=(Xi-ObjectData.Coord(1,1));
298        distY=(Yi-ObjectData.Coord(1,2));
299        flag=(distX.*distX/X2Max+distY.*distY/Y2Max)<1;
300    elseif isequal(ObjectData.Type,'rectangle')||isequal(ObjectData.Type,'volume')
301        XimaMin=ObjectData.Coord(1,1)-XMax;
302        XimaMax=ObjectData.Coord(1,1)+XMax;
303        YimaMin=ObjectData.Coord(1,2)-YMax;
304        YimaMax=ObjectData.Coord(1,2)+YMax;
305        xlim=[1.2*XimaMin-0.2*XimaMax 1.2*XimaMax-0.2*XimaMin];%create an image around the ellipse
306        ylim=[1.2*YimaMin-0.2*YimaMax 1.2*YimaMax-0.2*YimaMin];
307        scale_x=2*1.4*XMax/npMx;
308        scale_y=2*1.4*YMax/npMy;
309        xi=(0.5:npMx-0.5)*scale_x+xlim(1);
310        yi=(0.5:npMy-0.5)*scale_y+ylim(1);
311        [Xi,Yi]=meshgrid(xi,yi);
312        distX=abs(Xi-ObjectData.Coord(1,1));
313        distY=abs(Yi-ObjectData.Coord(1,2));
314        flag=distX<XMax & distY< YMax;
315    elseif isequal(ObjectData.Type,'polygon')
316        XimaMin=min(ObjectData.Coord(:,1));
317        XimaMax=max(ObjectData.Coord(:,1));
318        YimaMin=min(ObjectData.Coord(:,2));
319        YimaMax=max(ObjectData.Coord(:,2));
320        xlim=[1.2*XimaMin-0.2*XimaMax 1.2*XimaMax-0.2*XimaMin];
321        ylim=[1.2*YimaMin-0.2*YimaMax 1.2*YimaMax-0.2*YimaMin];
322        [Xlim,Ylim]=meshgrid(linspace(xlim(1),xlim(2),npMx),linspace(ylim(1),ylim(2),npMy));
323        %flag=roipoly(xlim,ylim,flag,ObjectData.Coord(:,1),ObjectData.Coord(:,2));%=1 inside the polygon, 0 outsid
324        flag=inpolygon(Xlim,Ylim,ObjectData.Coord(:,1),ObjectData.Coord(:,2));%=1 inside the polygon, 0 outsid
325    end
326    if isequal(ObjectData.ProjMode,'outside')||isequal(ObjectData.ProjMode,'mask_outside')
327        flag=~flag;
328    end
329    imflag=zeros(npMx,npMy,3);
330    imflag(:,:,3)=flag; % blue color
331    if isequal(col,'m')
332         imflag(:,:,1)=flag; % magenta color
333    end
334    dx=(xlim(2)-xlim(1))/npMx;
335    dy=(ylim(2)-ylim(1))/npMy;
336end
337
338PlotData=[];%default
339
340%% MODIFY AN EXISTING OBJECT PLOT
341if test_newobj==0;
342    hh=hplot;
343    PlotData=get(hplot,'UserData');
344    if test_line
345        set(hplot,'XData',xline)
346        set(hplot,'YData',yline)
347        %modify subobjects
348        if isfield(PlotData,'SubObject')
349            if isequal(ObjectData.Type,'points')
350                if ~isequal(YMax,0)
351                    for ipt=1:min(length(PlotData.SubObject),size(ObjectData.Coord,1))
352                        set(PlotData.SubObject(ipt),'Position',[ObjectData.Coord(ipt,1)-YMax ObjectData.Coord(ipt,2)-YMax 2*YMax 2*YMax])
353                    end
354                    %complement missing points
355                    if length(PlotData.SubObject)>nbpoints% fpoints in excess on the graph
356                        for ii=nbpoints+1: length(PlotData.SubObject);
357                            if ishandle(PlotData.SubObject(ii))
358                                delete(PlotData.SubObject(ii))
359                            end
360                        end
361                    end
362                    %                NbDeformPoint=nbpoints;
363                   
364                    if nbpoints>length(PlotData.SubObject)
365                        for ipt=length(PlotData.SubObject)+1:nbpoints
366                            PlotData.SubObject(ipt)=rectangle('Curvature',[1 1],...
367                                'Position',[ObjectData.Coord(ipt,1)-YMax ObjectData.Coord(ipt,2)-YMax 2*YMax 2*YMax],'EdgeColor',col,...
368                                'LineStyle',SubLineStyle,'Tag','proj_object');
369                        end
370                    end
371                end
372            elseif length(PlotData.SubObject)==2
373                set(PlotData.SubObject(1),'XData',xinf);
374                set(PlotData.SubObject(1),'YData',yinf);
375                set(PlotData.SubObject(2),'XData',xsup);
376                set(PlotData.SubObject(2),'YData',ysup);
377            elseif length(PlotData.SubObject)==1
378                set(PlotData.SubObject(1),'XData',xsup);
379                set(PlotData.SubObject(1),'YData',ysup);
380            end
381        end
382        if isfield(PlotData,'DeformPoint')
383            NbDeformPoint=length(PlotData.DeformPoint);
384            if NbDeformPoint>nbpoints% fpoints in excess on the graph
385                for ii=nbpoints+1:NbDeformPoint;
386                    if ishandle(PlotData.DeformPoint(ii))
387                        delete(PlotData.DeformPoint(ii))
388                    end
389                end
390                NbDeformPoint=nbpoints;
391            end
392            for ipt=1:NbDeformPoint
393                if ishandle(PlotData.DeformPoint(ipt))
394                    if nbpoints>=ipt
395                        set(PlotData.DeformPoint(ipt),'XData',xline(ipt),'YData',yline(ipt));
396                    end
397                end
398            end
399            if nbpoints>length(PlotData.DeformPoint)
400                for ipt=length(PlotData.DeformPoint)+1:nbpoints
401                    PlotData.DeformPoint(ipt)=line(xline(ipt),yline(ipt),'Color',col,'LineStyle','-','Tag','DeformPoint',...
402                        'SelectionHighlight','off','UserData',hplot);
403                end
404                set(hplot,'UserData',PlotData)
405            end
406        end
407    elseif (isequal(ObjectData.Type,'rectangle')||isequal(ObjectData.Type,'ellipse'))&&XMax>0 && YMax>0
408        set(hplot,'Position',[ObjectData.Coord(1,1)-XMax ObjectData.Coord(1,2)-YMax 2*XMax 2*YMax])
409    end
410    if test_patch
411        if isfield(PlotData,'SubObject')
412        for iobj=1:length(PlotData.SubObject)
413            if ~ishandle(PlotData.SubObject(iobj))
414                hold on
415                hhh=image([xlim(1)+dx/2 xlim(2)-dx/2],[ylim(1)+dy/2 ylim(2)-dy/2],imflag,'Tag','proj_object','HitTest','off');
416                set(hhh,'AlphaData',(flag)*0.2)% set partial transparency to the filling color
417                PlotData.SubObject(iobj)=hhh;
418            else
419                objtype=get(PlotData.SubObject(iobj),'Type');
420                if isequal(objtype,'image')
421                    set(PlotData.SubObject(iobj),'CData',imflag,'AlphaData',(flag)*0.2)
422                    set(PlotData.SubObject(iobj),'XData',[xlim(1)+dx/2 xlim(2)-dx/2])
423                    set(PlotData.SubObject(iobj),'YData',[ylim(1)+dy/2 ylim(2)-dy/2])
424                end
425            end
426        end
427        end
428    else% no patch image requested, erase existing ones
429        if isfield(PlotData,'SubObject')
430            for iobj=1:length(PlotData.SubObject)
431                if ishandle(PlotData.SubObject(iobj)) && strcmp(get(PlotData.SubObject(iobj),'Type'),'image')
432                    delete(PlotData.SubObject(iobj))
433                end
434            end
435        end
436    end
437end
438
439%% create the object
440if test_newobj
441%     axes(haxes)
442    hother=findobj('Tag','proj_object');%find all the proj objects
443    for iobj=1:length(hother)
444        if strcmp(get(hother(iobj),'Type'),'rectangle')|| strcmp(get(hother(iobj),'Type'),'patch')
445            set(hother(iobj),'EdgeColor','b')
446            if isequal(get(hother(iobj),'FaceColor'),'m')
447                set(hother(iobj),'FaceColor','b')
448            end
449        elseif isequal(get(hother(iobj),'Type'),'image')
450               Acolor=get(hother(iobj),'CData');
451               Acolor(:,:,1)=zeros(size(Acolor,1),size(Acolor,2));
452               set(hother(iobj),'CData',Acolor);
453        else
454             set(hother(iobj),'Color','b')
455        end
456        set(hother(iobj),'Selected','off')
457    end
458    hother=findobj('Tag','DeformPoint');
459    set(hother,'Color','b');
460    set(hother,'Selected','off') 
461    switch ObjectData.Type
462        case 'points'
463            hh=line(ObjectData.Coord(:,1),ObjectData.Coord(:,2),'Color',col,'LineStyle','none','Marker','+');
464            for ipt=1:length(xline)
465                PlotData.DeformPoint(ipt)=line(ObjectData.Coord(ipt,1),ObjectData.Coord(ipt,2),'Color',...
466                    col,'LineStyle','none','SelectionHighlight','off','UserData',hh,'Tag','DeformPoint');
467                %create circle around each point
468                if ~isequal(YMax,0)
469                    PlotData.SubObject(ipt)=rectangle('Curvature',[1 1],...
470                        'Position',[ObjectData.Coord(ipt,1)-YMax ObjectData.Coord(ipt,2)-YMax 2*YMax 2*YMax],'EdgeColor',col,...
471                        'LineStyle',SubLineStyle,'Tag','proj_object');
472                end
473            end
474        case {'line','polyline','polygon'}
475            hh=line(xline,yline,'Color',col);
476                PlotData.SubObject(1)=line(xinf,yinf,'Color',col,'LineStyle',SubLineStyle,'Tag','proj_object');%draw sub-lines
477                PlotData.SubObject(2)=line(xsup,ysup,'Color',col,'LineStyle',SubLineStyle,'Tag','proj_object');
478                for ipt=1:sizcoord(1)
479                    PlotData.DeformPoint(ipt)=line(ObjectData.Coord(ipt,1),ObjectData.Coord(ipt,2),'Color',...
480                        col,'LineStyle','none','Marker','.','Tag','DeformPoint','SelectionHighlight','off','UserData',hh);
481                end
482        case {'plane','volume','plane_z'}
483            hh=line(xline,yline,'Color',col);
484            PlotData.SubObject(1)=line(xsup,ysup,'Color',col,'LineStyle',SubLineStyle,'Tag','proj_object');
485        case 'rectangle'
486            hh=rectangle('Position',[ObjectData.Coord(1,1)-XMax ObjectData.Coord(1,2)-YMax 2*XMax 2*YMax],'LineWidth',2,'EdgeColor',col);
487        case 'ellipse'
488            hh=rectangle('Curvature',[1 1],'Position',[ObjectData.Coord(1,1)-XMax ObjectData.Coord(1,2)-YMax 2*XMax 2*YMax],'EdgeColor',col,'LineWidth',2);
489        otherwise
490            msgbox_uvmat('ERROR','unknown ObjectData.Type in plot_object.m')
491            return
492    end
493    set(hh,'Tag','proj_object')
494     if test_patch
495         hold on
496        hhh=image([xlim(1)+dx/2 xlim(2)-dx/2],[ylim(1)+dy/2 ylim(2)-dy/2],imflag,'Tag','proj_object','HitTest','off');
497       set(hhh,'AlphaData',(flag)*0.2)% set partial transparency to the filling color
498         PlotData.SubObject=hhh;   
499     end
500    if isfield(PlotData,'SubObject')
501        set(PlotData.SubObject,'UserData',hh)%record the parent handles in the SubObjects
502    end
503    if isfield(PlotData,'DeformPoint')
504        for ipt=1:sizcoord(1)
505            set(PlotData.DeformPoint(ipt),'UserData',hh);%record the parent handles in the SubObjects
506        end
507        set(PlotData.DeformPoint,'UserData',hh)%record the parent handles in the SubObjects
508    end
509end
510set(hh,'UserData',PlotData)
Note: See TracBrowser for help on using the repository browser.