/usr/share/freemat/toolbox/graph/clabel.m is in freemat-data 4.0-5.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 | % CLABEL CLABEL Add Labels To Contour Plot
%
% Usage
%
% The clabel function adds labels to a contour plot
% Generate contour labels for a contour plot. The syntax
% for its use is either:
%
% handles = clabel(contourhandle,property,value,property,value,...)
%
% which labels all of the contours in the plot, or
%
% handles = clabel(contourhandle,vals,property,value,property,value,...)
%
% which only labels those contours indicated by the vector vals.
% The contourhandle must be the handle to a contour plot object.
% The remaining property/value pairs are passed to the text function
% to control the parameters of the generated text labels. See the
% text properties for the details on what can be used in those labels.
% Copyright (c) 2002-2007 Samit Basu
% Licensed under the GPL
function handles = clabel(contourhandle,varargin)
if (~strcmp(get(contourhandle,'type'),'contour'))
error('First argument to clabel must be the handle of a contour.');
end
C = get(contourhandle,'contourmatrix');
vals_to_label = get(contourhandle,'levellist');
if ((numel(varargin) > 0) && (~isstr(varargin{1})))
vals_to_label = varargin{1};
varargin = varargin(2:end);
end
% Loop through the contours
colptr = 1;
while (colptr < size(C,2))
% Get the next line segment
linestart = colptr+1;
linestop = C(2,colptr) + linestart - 1;
lineval = C(1,colptr);
if (any(close_enough(vals_to_label-lineval)))
xvals = C(1,linestart:linestop);
yvals = C(2,linestart:linestop);
% Compute the svec - the distance along the line
dxvals = xvals(2:end)-xvals(1:end-1);
dyvals = yvals(2:end)-yvals(1:end-1);
angs = atan2(dyvals,dxvals)*180/pi;
svals = sqrt(dxvals.^2 + dyvals.^2);
sdist = cumsum(svals);
sndx = 10;
hlist = [];
while (sndx < numel(sdist)-10)
h = text(xvals(sndx),yvals(sndx),sprintf('%g',lineval), ...
'rotation',angs(sndx),varargin{:});
hlist = [hlist,h];
[a,sndx] = find(sdist > sdist(sndx)+max(sdist)/3,1);
end
end
colptr = linestop + 1;
end
if (nargout > 0)
handles = hlist;
end
function u = close_enough(d)
u = abs(d) < 100*eps;
|