图像处理系列——图像融合之色彩变换3（圆柱体IHS）

IHS变换融合算法相关的基础知识和原理请参考博文--色彩变换1。这里直接续上讲其变种模型的实现。

1、圆柱体IHS变换（间接IHS）

//进行正变换//G、B比较
if(dRGB[1]>=dRGB[2]) C = 0;
if(dRGB[1]<dRGB[2]) C = PI;
fI = (1/sqrt(3.0))*(dRGB[2]+dRGB[1]+dRGB[0]);
v1 = (dRGB[2]+dRGB[1]-2*dRGB[0])/(sqrt(6.0));
v2 = (dRGB[2]-dRGB[1])/(sqrt(2.0));
fH = atan(v1/v2) + C;
fS = sqrt(v1*v1+v2*v2);

dRGBt[0] = fI;dRGBt[1] = fS;dRGBt[2] = fH;
//用高空间分辨率的图像替换Intensity亮度分量
dRGBt[0] = pPANImage[i*width+j];

//进行逆变换
br = (1/sqrt(3.0))*dRGBt[0]-(2/sqrt(6.0))*v1;
bg = (1/sqrt(3.0))*dRGBt[0]+(1/sqrt(6.0))*v1-(1/sqrt(2.0))*v2;
bb = (1/sqrt(3.0))*dRGBt[0]+(1/sqrt(6.0))*v1+(1/sqrt(2.0))*v2;		
br *= 1.5;
bg *= 1.5;
bb *= 1.5;

2、圆柱体IHS变换结果图

3、如果稍微改动其逆变换，如下

//进行逆变换
br = (1/sqrt(3.0))*dRGBt[0]-(2/sqrt(6.0))*v1;
bg = (1/sqrt(3.0))*dRGBt[0]+(1/sqrt(6.0))*v1-(1/sqrt(2.0))*v2;
bb = (1/sqrt(3.0))*dRGBt[0]+(1/sqrt(6.0))*v1+(1/sqrt(2.0))*v2;		
//br *= 1.5;
//bg *= 1.5;
//bb *= 1.5;

此时的结果图如下：