tri2list.m

物流分析工具包。Facility location: Continuous minisum facility location, alternate location-allocation (ALA)

function IJ = tri2list(T)
%TRI2LIST Convert triangle indices to arc list representation.
%    IJ = tri2list(T)
%     T = n x 3 matrix, where each row defines indices for one of n 
%         triangles
%    IJ = 2-column arc list
%
% Converts set of triangles defined by X and Y vector indices into IJ. 
% Useful for ploting results from DELAUNAY triangulation using PPLOT.
%
% See also CONVHULL and DSEARCH

% Copyright (c) 1994-2006 by Michael G. Kay
% Matlog Version 9 13-Jan-2006 (http://www.ie.ncsu.edu/kay/matlog)

% Input Error Checking ****************************************************
if size(T,2) ~= 3
   error('"T" must be an n x 3 matrix of triangle indices.')
end
% End (Input Error Checking) **********************************************

m = max(T(:));
A = sparse(m,m);
for i = 1:length(T(:,1))
   A(T(i,1),T(i,2)) = 1;
   A(T(i,2),T(i,3)) = 1;
   A(T(i,3),T(i,1)) = 1;
end
A = A | triu(A)' | tril(A)';
IJ = adj2list(A);
IJ = IJ(:,[1 2]);