gmviz.m

算断裂的

function gmviz(obj, color, vizdim)
%% GMVIZ: plot a brep or simplicial complex in matlab or VRML
%% gmviz(obj {, color {, vizdim }});
%% Displays an object.  The object is either a brep or simplicial
%% complex.  The object is either displayed in the
%% matlab graphics window or is sent to a file in VRML1.0 format, 
%% and may be automatically loaded by a web browser.  The
%% selection of the graphics engine is controlled by gmvizgui.
%% Argument color is either a string with
%% one or more color letters, e.g., 'rw' for
%% red white, or a matrix with a 4-element color spec per row.
%% vizdim is the dimension or dimensions to plot.  The number
%% of colors specified must be the same as the number of dimensions
%% specified.  See the documentation for more information.
%% A second calling format is 
%% gmviz({simpcomp,brep} {, vizcolor {, vizdim}})
%% in which case the simplicial complex colors are derived
%% from the brep color specifications.


if nargin < 1
  error('gmviz requires at least one argument');
end

[l1,l2] = size(obj);
if l1 * l2 == 2
  pair = 1;
  vizobj = obj{1};
  obj2 = obj{2};
else
  pair = 0;
  vizobj = obj;
end

if ~isa(vizobj,'zba')
  error('Object passed to gmviz not a zba');
end

if length(vizobj) < 1
  error('Object passed to gmviz not an simpcomp or brep');
end


global GM_BREP_TYPE_CODE
global GM_SIMPCOMP_TYPE_CODE
typecode = double(vizobj{0});
if strcmpi(typecode, GM_BREP_TYPE_CODE)
  objtype = 1;
elseif  strcmpi(typecode, GM_SIMPCOMP_TYPE_CODE)
  objtype = 2;
else
  error('Object passed to gmviz not an simpcomp or brep');
end

if pair & objtype == 1
  error('In calling format gmviz({a,b}), a must be a simpcomplex');
end
errmsg2 = 'In calling format gmviz({a,b}), b must be a brep';
if pair == 1
  if ~isa(obj2, 'zba')
    error(errmsg2)
  end
  if length(obj2) < 1
    error(errmsg2)
  end
  if ~strcmpi(double(obj2{0}), GM_BREP_TYPE_CODE)
    error(errmsg2)
  end
end

di = double(vizobj{2});
gdim = double(vizobj{1});
if gdim < 0
  error('Object to plot has no faces')
end
if di ~= 2 & di ~= 3
  error('Object to gmviz must have dimension 2 or 3')
end
bbox = 1e307 * [ones(di,1),-ones(di,1)];

[engine, filename, filestat, notify, format, ...
 thickness, default_color, beziersub] = gmvizgui_get_state(di);

if nargin < 2
  color = default_color;
end
if nargin < 3
  vizdim = min(gdim, di - 1);
end

for p = 1 : length(vizdim)
  if vizdim(p) < 0 | vizdim(p) > gdim
    error('Viz dimension out of range')
  end
end


if (isstr(color))
  colorn = [];
  for j = 1 : length(color)
    if color(j) == 'y'
      color1 = [1,1,0,1];
    elseif color(j) == 'm'
      color1 = [1,0,1,1];
    elseif color(j) == 'c'
      color1 = [0,1,1,1];
    elseif color(j) == 'r'
      color1 = [1,0,0,1];
    elseif color(j) == 'g'
      color1 = [0,1,0,1];
    elseif color(j) == 'b'
      color1 = [0,0,1,1];
    elseif color(j) == 'w'
      color1 = [1,1,1,1];
    elseif color(j) == 'k'
      color1 = [0,0,0,1];
    elseif color(j) == 'i'
      color1 = [0,0,0,0];
    else
      error('Unknown color name')
    end
    colorn = [colorn; color1];
  end
else
  colorn = color;
end

[colrow,colcol] = size(colorn);

if colcol ~= 4
   error('color spec must have 4 entries')
end

if colrow ~= length(vizdim)
  error('Number of colors must equal number of dimensions')
end


for k = 1 : length(vizdim)
  if pair
    [pts,simps,cols] = gm_vizp(vizobj, obj2, colorn(k,:), ...
                               vizdim(k), beziersub);
  else
    [pts,simps,cols] = gm_vizp(vizobj,  colorn(k,:), ...
                               vizdim(k), beziersub);
  end
  bbox = gmviz_update_bbox(bbox, pts);
end

if bbox(1,1) > bbox(1,2)
  error('No objects to plot')
end


if strcmpi(engine,'local')
  handle1 = gmviz_init_matlabplot(bbox);
else
  handle1 = gmviz_init_vrmlplot(bbox, filename, filestat);
end


for k = 1 : length(vizdim)
  if pair
    [pts,simps,cols] = gm_vizp(vizobj, obj2, colorn(k,:), ...
                               vizdim(k), beziersub);
  else
    [pts,simps,cols] = gm_vizp(vizobj,  colorn(k,:), ...
                               vizdim(k), beziersub);
  end
  if strcmpi(engine,'local')
    gmviz_matlabplot(handle1, bbox, pts, simps, cols, vizdim(k), thickness);
  else
    gmviz_vrmlplot(handle1, bbox, pts, simps, cols, vizdim(k), thickness);
  end
end

if strcmpi(engine,'local')
  gmviz_terminate_matlabplot(handle1);
else
  gmviz_terminate_vrmlplot(handle1, filename, notify);
end







% ------------------------------------------------------------------
% Copyright (c) 1999 by Cornell University.  All rights reserved.
% See the accompanying file 'Copyright' for authorship information,
% the terms of the license governing this software, and disclaimers
% concerning this software.
% ------------------------------------------------------------------
% This file is part of the QMG software.  
% Version 2.0 of QMG, release date September 3, 1999
% ------------------------------------------------------------------