plotIGES.m

function handlePlot=plotIGES(ParameterData,srf,fignr,subd,holdoff_flag,fine_flag,plotNonOriginalCrv)
% PLOTIGES plots surfaces, curves and points from IGES-file
%
% Simple usage:
%
% handlePlot=plotIGES(ParameterData)
%
% or
%
% handlePlot=plotIGES(ParameterData,1)
%
% Ordinary usage:
%
% plotIGES(ParameterData,srf,fignr,subd,holdoff_flag)
%
% Input:
%
% ParameterData - Parameter data from IGES file. ParameterData
%                 is the output from IGES2MATLAB
% srf - Flag for surface plotting. (0,1,2)
%        0 (default), no surfaces are plotted,
%        1 the surface is plotted as triangular patches,
%        2 the surface is plotted as triangular mesh
% fignr - Figure number of the plot. 1 default
% subd - Nuber of subdivisions when plotting curves
%        subd is nubmer of subdivisions for each parameter when
%        plotting surfaces. 100 default
% holdoff_flag - Bolean value (1/0). If 1 then hold off the plot
%                when the plot is done. 1 default
% fine_flag - Bolean value (1/0). If 0 the surface will be rough
%             and if 1 the surface will be finer. 0 default
% plotNonOriginalCrv - Flag for curve plotting (1/0).
%        0, no non original curves are plotted,
%        1 (default),
%
% Output:
%
% handlePlot - plothandle
%
% m-file can be downloaded at
% http://www.mathworks.com/matlabcentral/fileexchange/13253-iges-toolbox
%
% written by Per Bergstrm 2012-01-09
%

if nargin<7
    plotNonOriginalCrv=1;
    if nargin<6
        fine_flag=0;
        if nargin<5
            holdoff_flag=1;
            if nargin<4
                subd=100;
                if nargin<3
                    fignr=1;
                    if nargin<2
                        srf=0;
                        if nargin<1
                            error('plotIGES must have an input');
                        end
                    end
                end
            end
        end
    end
end

if isempty(ParameterData)
    error('Empty ParameterData');
elseif not(iscell(ParameterData))
    error('Invalid ParameterData. Must be a cell array!');
end

if isempty(srf)
    srf=0;
elseif not(srf==0 | srf==1 | srf==2)
    srf=0;
end

srf0not=not(srf==0);
srf1=srf==1;

if isempty(fignr)
    fignr=1;
else
    fignr=round(fignr);
end

if isempty(subd)
    subd=100;
else
    subd=round(subd);
end

if subd<5
    subd=5;
end

if isempty(holdoff_flag)
    holdoff_flag=1;
elseif not(holdoff_flag==0 | holdoff_flag==1)
    holdoff_flag=1;
end

if isempty(fine_flag)
    fine_flag=0;
elseif not(fine_flag==0 | fine_flag==1)
    fine_flag=0;
end

if isempty(plotNonOriginalCrv)
    plotNonOriginalCrv=1;
elseif not(plotNonOriginalCrv==0 | plotNonOriginalCrv==1)
    plotNonOriginalCrv=1;
end

subd=subd+1;  % subd now number of points, not number of subintervals

siz=length(ParameterData);

if srf1
    clr=[0.53 0.24 0.24];
    lclrf=0.3;
    for i=siz:-1:1
        if ParameterData{i}.type==314
            clr=[ParameterData{i}.cc1 ParameterData{i}.cc2 ParameterData{i}.cc3]/100;
            break
        end
    end
else
    clr=[0 0 1];
    lclrf=0.5;
end

if fignr>0
    figure(fignr),hold on
end

if nargout>0
    
    handlePlot=cell(1,3);
    
    for i=1:siz
        
        [P,isSCP,isSup]=retSrfCrvPnt(2,ParameterData,1,i,subd,3);
        
        if and(isSCP,not(isSup))
            
            if plotNonOriginalCrv
                hl=plot3(P(1,:),P(2,:),P(3,:),'Color',lclrf*clr,'LineWidth',1);
            else
                if ParameterData{i}.original
                    hl=plot3(P(1,:),P(2,:),P(3,:),'Color',lclrf*clr,'LineWidth',1);
                end
            end
            
            handlePlot{1}=[handlePlot{1};hl];
            
        elseif not(isSCP)
            
            [P,isSCP]=retSrfCrvPnt(3,ParameterData,1,i);
            
            if isSCP
                
                hp=plot3(ParameterData{i}.x,ParameterData{i}.y,ParameterData{i}.z,'r.');
                
                handlePlot{2}=[handlePlot{2};hp];
                
            elseif srf0not
                
                if fine_flag
                    [P,isSCP,isSup,TRI]=retSrfCrvPnt(1,ParameterData,1,i,subd,1);
                else
                    [P,isSCP,isSup,TRI]=retSrfCrvPnt(1,ParameterData,1,i,subd);
                end
                
                if and(isSCP,not(isSup))
                    if srf1
                        hs=patch('faces',TRI,'vertices',P','FaceColor',clr,'EdgeColor','none');
                    else
                        hs=patch('faces',TRI,'vertices',P','FaceColor','none','EdgeColor','b');
                    end
                    
                    handlePlot{3}=[handlePlot{3};hs];
                end
                
            end
        end
        
    end
    
else
    
    for i=1:siz
        
        [P,isSCP,isSup]=retSrfCrvPnt(2,ParameterData,1,i,subd,3);
        
        if and(isSCP,not(isSup))
            
            if plotNonOriginalCrv
                plot3(P(1,:),P(2,:),P(3,:),'Color',lclrf*clr,'LineWidth',1);
            else
                if ParameterData{i}.original
                    plot3(P(1,:),P(2,:),P(3,:),'Color',lclrf*clr,'LineWidth',1);
                end
            end
            
        elseif not(isSCP)
            
            [P,isSCP]=retSrfCrvPnt(3,ParameterData,1,i);
            
            if isSCP
                
                plot3(ParameterData{i}.x,ParameterData{i}.y,ParameterData{i}.z,'r.');
                
            elseif srf0not
                
                if fine_flag
                    [P,isSCP,isSup,TRI]=retSrfCrvPnt(1,ParameterData,1,i,subd,1);
                else
                    [P,isSCP,isSup,TRI]=retSrfCrvPnt(1,ParameterData,1,i,subd);
                end
                
                if and(isSCP,not(isSup))
                    if srf1
                        patch('faces',TRI,'vertices',P','FaceColor',clr,'EdgeColor','none');
                    else
                        patch('faces',TRI,'vertices',P','FaceColor','none','EdgeColor','b');
                    end
                end
                
            end
        end
        
    end
    
end

axis image

sc=0.2;

xl=xlim;
dx=xl(2)-xl(1);
xl(1)=xl(1)-sc*dx;
xl(2)=xl(2)+sc*dx;
xlim(xl);

yl=ylim;
dy=yl(2)-yl(1);
yl(1)=yl(1)-sc*dy;
yl(2)=yl(2)+sc*dy;
ylim(yl);

zl=zlim;
dz=zl(2)-zl(1);
zl(1)=zl(1)-sc*dz;
zl(2)=zl(2)+sc*dz;
zlim(zl);

if holdoff_flag
    hold off;
end