+ itk::Vector< double, 3 > rgb;
+ rgb[ 0 ] = double( r );
+ rgb[ 1 ] = double( b );
+ rgb[ 2 ] = double( g );
+ rgb = iM * rgb;
+
+ P out;
+ out[ 0 ] = TValue( rgb[ 0 ] );
+ out[ 1 ] = TValue( rgb[ 1 ] );
+ out[ 2 ] = TValue( rgb[ 2 ] );
+ return( out );