]> Creatis software - CreaPhase.git/blob - octave_packages/image-1.0.15/mat2gray.m
Add a useful package (from Source forge) for octave
[CreaPhase.git] / octave_packages / image-1.0.15 / mat2gray.m
1 ## Copyright (C) 1999,2000  Kai Habel
2 ##
3 ## This program is free software; you can redistribute it and/or modify
4 ## it under the terms of the GNU General Public License as published by
5 ## the Free Software Foundation; either version 2 of the License, or
6 ## (at your option) any later version.
7 ##
8 ## This program is distributed in the hope that it will be useful,
9 ## but WITHOUT ANY WARRANTY; without even the implied warranty of
10 ## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
11 ## GNU General Public License for more details.
12 ##
13 ## You should have received a copy of the GNU General Public License
14 ## along with this program; If not, see <http://www.gnu.org/licenses/>.
15
16 ## -*- texinfo -*-
17 ## @deftypefn {Function File} @var{I}= mat2gray (@var{M},[min max])
18 ## Converts a matrix to a intensity image.
19 ## @end deftypefn
20
21 ## Author:      Kai Habel <kai.habel@gmx.de>
22 ## Date:        22/03/2000
23
24 function I = mat2gray (M, scale)
25
26   if (nargin < 1|| nargin > 2)
27     usage ("mat2gray(...) number of arguments must be 1 or 2");
28   endif
29
30   if (!ismatrix (M))
31     usage ("mat2gray(M,...) M must be a matrix");
32   endif
33
34   if (nargin == 1)
35     Mmin = min (min (M));
36     Mmax = max (max (M));
37   else 
38     if (isvector (scale))
39       Mmin = min (scale (1), scale (2));
40       Mmax = max (scale (1), scale (2));
41     endif
42   endif
43
44   I = (M < Mmin) .* 0;
45   I = I + (M >= Mmin & M < Mmax) .* (1 / (Mmax - Mmin) * (M - Mmin));
46   I = I + (M >= Mmax);
47
48 endfunction