Mercurial > octave-dspies
view scripts/image/ind2rgb.m @ 15688:4db08f52a6ed
ind2rgb: correction for indexes when image is uint class. Added tests block.
author | Carnë Draug <carandraug+dev@gmail.com> |
---|---|
date | Mon, 12 Nov 2012 04:21:52 +0000 |
parents | 4beb3a4bd440 |
children | 14b7679891dd |
line wrap: on
line source
## Copyright (C) 1994-2012 John W. Eaton ## ## This file is part of Octave. ## ## Octave 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. ## ## Octave 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 Octave; see the file COPYING. If not, see ## <http://www.gnu.org/licenses/>. ## -*- texinfo -*- ## @deftypefn {Function File} {@var{rgb} =} ind2rgb (@var{x}) ## @deftypefnx {Function File} {@var{rgb} =} ind2rgb (@var{x}, @var{map}) ## @deftypefnx {Function File} {[@var{R}, @var{G}, @var{B}] =} ind2rgb (@dots{}) ## Convert an indexed image to red, green, and blue color components. ## If the colormap doesn't contain enough colors, pad it with the ## last color in the map. ## If @var{map} is omitted, the current colormap is used for the conversion. ## @seealso{rgb2ind, ind2gray, hsv2rgb, ntsc2rgb} ## @end deftypefn ## Author: Tony Richardson <arichard@stark.cc.oh.us> ## Created: July 1994 ## Adapted-By: jwe function [R, G, B] = ind2rgb (x, map) ## Do we have the right number of inputs? if (nargin < 1 || nargin > 2) print_usage (); elseif (nargin == 1) map = colormap (); endif ## Check if X is an indexed image. if (ndims (x) != 2 || (isfloat (x) && ! isindex (x)) || ! ismember (class (x), {"double", "single", "uint8", "uint16"})) error ("ind2rgb: X must be an indexed image"); endif ## Check the color map. if (! iscolormap (map)) error ("ind2rgb: MAP must be a valid colormap"); endif ## Do we have enough colors in the color map? ## there's an offset of 1 when the indexed image is an integer class so we fix ## it now and convert it to float only if really necessary and even then only ## to single precision since its enough for both uint8 and uint16 maxidx = max (x(:)); if (isinteger (x)) if (maxidx == intmax (class (x))) x = single (x); endif x += 1; maxidx += 1; endif rm = rows (map); if (rm < maxidx) ## Pad with the last color in the map. pad = repmat (map(end,:), maxidx-rm, 1); map(end+1:maxidx, :) = pad; endif ## Compute result [hi, wi] = size (x); R = reshape (map (x(:), 1), hi, wi); G = reshape (map (x(:), 2), hi, wi); B = reshape (map (x(:), 3), hi, wi); ## Use 3D array if only one output is requested. if (nargout <= 1) R(:,:,3) = B; R(:,:,2) = G; endif endfunction %!shared img, map, ergb, rgb, r, g, b %! img = [2 4 5; 3 2 5; 1 2 4]; %! map = [0.0 0.0 0.0 %! 0.2 0.4 0.6 %! 0.4 0.4 0.5 %! 0.3 0.7 1.0 %! 0.1 0.5 0.8]; %! ergb(:,:,1) = [0.2 0.3 0.1; 0.4 0.2 0.1; 0.0 0.2 0.3]; %! ergb(:,:,2) = [0.4 0.7 0.5; 0.4 0.4 0.5; 0.0 0.4 0.7]; %! ergb(:,:,3) = [0.6 1.0 0.8; 0.5 0.6 0.8; 0.0 0.6 1.0]; %! ## test basic usage with 1 and 3 outputs %! [rgb] = ind2rgb (img, map); %! [r, g, b] = ind2rgb (img, map); %!assert (ergb, rgb); %!assert (ergb, reshape ([r(:) g(:) b(:)], [size(img) 3])); %! ## test correction for integers %! img = uint8 (img -1); %! [rgb] = ind2rgb (img, map); %!assert (ergb, rgb); %! ## check it fails when image is a float with an index value of 0 %!fail("[rgb] = ind2rgb (double(img), map)")