annotate scripts/image/cmpermute.m @ 15515:4beb3a4bd440

rgbplot.m, cmpermute.m, cmunique.m, ind2rgb.m, imwrite.m: use core iscolormap function
author Carnë Draug <carandraug+dev@gmail.com>
date Fri, 12 Oct 2012 23:11:08 +0200
parents 0ba7be7fed1c
children d1285ebe60ca
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
14896
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
1 ## Copyright (C) 2004 Josep Mones i Teixidor
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
2 ## Copyright (C) 2012 Rik Wehbring
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
3 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
4 ## This file is part of Octave.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
5 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
6 ## Octave is free software; you can redistribute it and/or modify it
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
7 ## under the terms of the GNU General Public License as published by
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
8 ## the Free Software Foundation; either version 3 of the License, or (at
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
9 ## your option) any later version.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
10 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
11 ## Octave is distributed in the hope that it will be useful, but
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
12 ## WITHOUT ANY WARRANTY; without even the implied warranty of
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
13 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
14 ## General Public License for more details.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
15 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
16 ## You should have received a copy of the GNU General Public License
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
17 ## along with Octave; see the file COPYING. If not, see
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
18 ## <http://www.gnu.org/licenses/>.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
19
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
20 ## -*- texinfo -*-
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
21 ## @deftypefn {Function File} {[@var{Y}, @var{newmap}] =} cmpermute (@var{X}, @var{map})
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
22 ## @deftypefnx {Function File} {[@var{Y}, @var{newmap}] =} cmpermute (@var{X}, @var{map}, @var{index})
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
23 ## Reorder colors in a colormap.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
24 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
25 ## When called with only two arguments, @code{cmpermute} randomly rearranges
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
26 ## the colormap @var{map} and returns a new colormap @var{newmap}. It also
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
27 ## returns the indexed image @var{Y} which is the equivalent of the original
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
28 ## input image @var{X} when displayed using @var{newmap}. The input image
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
29 ## @var{X} must be an indexed image of class uint8 or double.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
30 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
31 ## When called with an optional third argument the order of colors in the
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
32 ## new colormap is defined by @var{index}.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
33 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
34 ## @strong{Caution:} @code{index} should not have repeated elements or the
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
35 ## function will fail.
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
36 ##
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
37 ## @end deftypefn
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
38
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
39 ## Author: Josep Mones i Teixidor <jmones@puntbarra.com>
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
40
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
41 function [Y, newmap] = cmpermute (X, map, index)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
42
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
43 if (nargin < 2 || nargin > 3)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
44 print_usage ();
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
45 endif
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
46
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
47 ## FIXME: Matlab only accepts 2 types. Expand to uint16 & single??
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
48 if (! (isa (X, "uint8") || isa (X, "double")))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
49 error ("cmpermute: X must be of class uint8 or double");
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
50 endif
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
51
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
52 if (! isreal (X) || issparse (X)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
53 || (isfloat (X) && (any (X(:) < 1 || any (X(:) != fix (X(:)))))))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
54 error ("cmpermute: X must be an indexed image");
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
55 endif
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
56
15515
4beb3a4bd440 rgbplot.m, cmpermute.m, cmunique.m, ind2rgb.m, imwrite.m: use core iscolormap function
Carnë Draug <carandraug+dev@gmail.com>
parents: 14896
diff changeset
57 if (! iscolormap (map))
14896
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
58 error ("cmpermute: MAP must be a valid colormap");
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
59 endif
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
60
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
61 if (nargin < 3)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
62 index = randperm (rows (map));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
63 elseif (! isvector (index) || length (index) != rows (map))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
64 error ("cmpermute: invalid parameter INDEX");
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
65 endif
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
66
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
67 ## new colormap
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
68 newmap = map(index,:);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
69
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
70 ## build reverse index
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
71 rindex = zeros (size (index));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
72 rindex(index) = 1:length (index);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
73
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
74 ## adapt indices
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
75 if (isa (X, "uint8"))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
76 rindex = uint8 (rindex-1);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
77 ## 0-based indices
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
78 Y = rindex(double (X) + 1);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
79 else
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
80 Y = rindex(X);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
81 endif
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
82
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
83 endfunction
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
84
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
85
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
86 %!demo
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
87 %! [Y, newmap] = cmpermute ([1:4], hot (4), 4:-1:1)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
88 %! ## colormap will be arranged in reverse order (so will image)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
89
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
90 %!shared X, map
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
91 %! X = uint8 (magic (16));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
92 %! [X, map] = cmunique (X);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
93
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
94 %!test # random permutation, 0-based index
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
95 %! [Y, newmap] = cmpermute (X, map);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
96 %! ## test we didn't lose colors
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
97 %! assert (sort (map), sortrows (newmap));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
98 %! ## test if images are equal
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
99 %! assert (map(double (X)+1), newmap(double (Y)+1));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
100
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
101 %!test # reverse map, 0-based index
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
102 %! [Y, newmap] = cmpermute (X, map, rows (map):-1:1);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
103 %! ## we expect a reversed colormap
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
104 %! assert (flipud (newmap), map);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
105 %! ## we expect reversed indices in image
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
106 %! assert (X, max (Y(:)) - Y);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
107
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
108 %!shared X,map
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
109 %! X = uint16 (magic (20));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
110 %! [X, map] = cmunique (X);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
111
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
112 %!test # random permutation, 1-based index
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
113 %! [Y, newmap] = cmpermute (X, map);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
114 %! ## test we didn't lose colors
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
115 %! assert (sort (map), sortrows (newmap));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
116 %! ## test if images are equal
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
117 %! assert (map(X), newmap(Y));
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
118
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
119 %!test # reverse map, 1-based index
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
120 %! [Y, newmap] = cmpermute (X, map, rows (map):-1:1);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
121 %! ## we expect a reversed colormap
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
122 %! assert (newmap (rows (newmap):-1:1,:), map);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
123 %! ## we expect reversed indices in image
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
124 %! assert (X, max (Y(:)) + 1 - Y);
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
125
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
126 ## Test input validation
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
127 %!error cmpermute ()
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
128 %!error cmpermute (1,2,3,4)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
129 %!error <X must be of class uint8> cmpermute (uint16 (magic (16)), jet (256))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
130 %!error <X must be an indexed image> cmpermute (1+i, jet (256))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
131 %!error <X must be an indexed image> cmpermute (sparse (1), jet (256))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
132 %!error <X must be an indexed image> cmpermute (0, jet (256))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
133 %!error <X must be an indexed image> cmpermute (1.5, jet (256))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
134 %!error <MAP must be a valid colormap> cmpermute (1, "a")
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
135 %!error <MAP must be a valid colormap> cmpermute (1, i)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
136 %!error <MAP must be a valid colormap> cmpermute (1, ones (3,3,3))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
137 %!error <MAP must be a valid colormap> cmpermute (1, ones (3,2))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
138 %!error <MAP must be a valid colormap> cmpermute (1, [-1 1 1])
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
139 %!error <MAP must be a valid colormap> cmpermute (1, [2 1 1])
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
140 %!error <invalid parameter INDEX> cmpermute (1, [0 1 0;1 0 1], ones (3))
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
141 %!error <invalid parameter INDEX> cmpermute (1, [0 1 0;1 0 1], 1:3)
0ba7be7fed1c Add new functions cmpermute(), cmunque().
Rik <rik@octave.org>
parents:
diff changeset
142