/usr/share/octave/packages/geometry-1.7.0/shape2d/shapetransform.m is in octave-geometry 1.7.0-1build1.
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 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 | ## Copyright (C) 2012 Juan Pablo Carbajal <carbajal@ifi.uzh.ch>
##
## This program is free software; you can redistribute it and/or modify it under
## the terms of the GNU General Public License as published by the Free Software
## Foundation; either version 3 of the License, or (at your option) any later
## version.
##
## This program is distributed in the hope that it will be useful, but WITHOUT
## ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
## FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more
## details.
##
## You should have received a copy of the GNU General Public License along with
## this program; if not, see <http://www.gnu.org/licenses/>.
## -*- texinfo -*-
## @deftypefn {Function File} {@var{nshape} = } shapetransform (@var{shape}, @var{T})
## Applies transformation to a shape defined by piecewise smooth polynomials.
##
## @var{shape} is a cell where each elements is a 2-by-(poly_degree+1) matrix
## containing a pair of polynomials.
##
## Format of @var{T} can be one of :
## @example
## @group
## [c] , [a b] , [a b c] or [a b c]
## [f] [d e] [d e f] [d e f]
## [0 0 1]
## @end group
## @end example
##
## @seealso{shape2polygon, shapeplot}
## @end deftypefn
function nshape = shapetransform (shape, Trans)
if size(Trans,1) < 2
error("geometry:shapetransform:InvalidArgument", ...
"Transformation can be 2x1, 2x2, 2x3 or 3x3. See help.");
end
if ~iscell(shape)
error("geometry:shapetransform:InvalidArgument", "Shape must be a cell of 2D polynomials.");
end
A =[];
v = [];
switch size(Trans,2)
case 1
# Just translation
v = Trans;
case 2
# Just linear transformation
A = Trans;
case 3
# Affine transform
A = Trans(1:2,1:2);
v = Trans(1:2,3);
end
nshape = cellfun (@(x)polytransform (x,A,v), shape, 'UniformOutput',false);
endfunction
function np = polytransform(p,A,v)
np = p;
if ~isempty (A)
np = A*np;
end
if ~isempty (v)
np(:,end) = np(:,end) + v;
end
endfunction
%!demo
%! shape = {[-93.172 606.368 -476.054 291.429; ...
%! -431.196 637.253 11.085 163.791]; ...
%! [-75.3626 -253.2337 457.1678 328.5714; ...
%! 438.7659 -653.6278 -7.9953 380.9336]; ...
%! [-89.5841 344.9716 -275.3876 457.1429; ...
%! -170.3613 237.8858 1.0469 158.0765];...
%! [32.900 -298.704 145.804 437.143; ...
%! -243.903 369.597 -34.265 226.648]; ...
%! [-99.081 409.127 -352.903 317.143; ...
%! 55.289 -114.223 -26.781 318.076]; ...
%! [-342.231 191.266 168.108 274.286; ...
%! 58.870 -38.083 -89.358 232.362]};
%!
%! A = shapearea (shape);
%! T = eye(2)/sqrt(A);
%! shape = shapetransform (shape,T);
%! T = shapecentroid (shape)(:);
%! shape = shapetransform (shape,-T + [2; 0]);
%!
%! close
%! shapeplot (shape,'-r','linewidth',2);
%! hold on
%! for i = 1:9
%! T = createRotation (i*pi/5)(1:2,1:2)/exp(0.3*i);
%! shapeplot (shapetransform(shape, T), 'color',rand(1,3),'linewidth',2);
%! end
%! hold off
%! axis tight
%! axis square
|