1 /*=========================================================================
2 *
3 * Copyright Insight Software Consortium
4 *
5 * Licensed under the Apache License, Version 2.0 (the "License");
6 * you may not use this file except in compliance with the License.
7 * You may obtain a copy of the License at
8 *
9 * http://www.apache.org/licenses/LICENSE-2.0.txt
10 *
11 * Unless required by applicable law or agreed to in writing, software
12 * distributed under the License is distributed on an "AS IS" BASIS,
13 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
14 * See the License for the specific language governing permissions and
15 * limitations under the License.
16 *
17 *=========================================================================*/
18 #ifndef itkHotColormapFunction_hxx
19 #define itkHotColormapFunction_hxx
20
21 #include "itkHotColormapFunction.h"
22
23 namespace itk
24 {
25 namespace Function
26 {
27 template< typename TScalar, typename TRGBPixel >
28 typename HotColormapFunction< TScalar, TRGBPixel >::RGBPixelType
29 HotColormapFunction< TScalar, TRGBPixel >
operator ()(const TScalar & v) const30 ::operator()(const TScalar & v) const
31 {
32 // Map the input scalar between [0, 1].
33 RealType value = this->RescaleInputValue(v);
34
35 // Apply the color mapping.
36 RealType red = 63.0 / 26.0 * value - 1.0 / 13.0;
37
38 red = std::max(0.0, red);
39 red = std::min(1.0, red);
40
41 RealType green = 63.0 / 26.0 * value - 11.0 / 13.0;
42 green = std::max(0.0, green);
43 green = std::min(1.0, green);
44
45 RealType blue = 4.5 * value - 3.5;
46 blue = std::max(0.0, blue);
47 blue = std::min(1.0, blue);
48
49 // Set the rgb components after rescaling the values.
50 RGBPixelType pixel;
51 NumericTraits<TRGBPixel>::SetLength(pixel, 3);
52
53 pixel[0] = this->RescaleRGBComponentValue(red);
54 pixel[1] = this->RescaleRGBComponentValue(green);
55 pixel[2] = this->RescaleRGBComponentValue(blue);
56
57 return pixel;
58 }
59 } // end namespace Function
60 } // end namespace itk
61
62 #endif
63