图像处理-基本算法之幂次变换

幂次变换的基本表达式为：y=cx^r+b

    其中c、r均为正数。与对数变换相同，幂次变换将部分灰度区域映射到更宽的区域中。当r=1时，幂次变换转变为线性变换。

   (1)  当r<0时，变换函数曲线在正比函数上方。此时扩展低灰度级，压缩高灰度级，使图像变亮。这一点与对数变换十分相似。

   (2)  当r>0时，变换函数曲线在正比函数下方。此时扩展高灰度级，压缩低灰度级，使图像变暗。

   代码如下：

[cpp] view plaincopyprint?

1. /******************************************************************************     
2. *   作用:     幂次变换函数
3. *   参数: pDst     输出图像的像素数组
4. *   参数: pSrc     原始图像的像素数组
5. *   参数: nWidth   原始图像宽度
6. *   参数: nHeight  原始图像高度
7. *   参数: b          控制参数,表示曲线的上下偏移量
8. *   参数: c          控制参数,表示曲线的弯曲程度
9. *   参数: r          控制参数,表示函数的幂次
10. *   备注: 此函数对于彩色图同样适用
11. ******************************************************************************/ 
12. int PowerTrans(BYTE* pDst, BYTE* pSrc, int nWidth, int nHeight, double b, double c, double r) 
13. { 
14.     if (!pSrc || !pDst) 
15.     { 
16.         return EXIT_FAILURE; 
17.     } 
18.  
19.     // 映射表,用于256种灰度变换后的值 
20.     BYTE map[256]; 
21.  
22.     // 保存运算后的临时值 
23.     double dTemp; 
24.     int i, j; 
25.     for (i = 0; i < 256; i++) 
26.     { 
27.         // 计算当前像素变换后的值 
28.         dTemp = c * pow(i / 255.0, r) * 255 + b; 
29.  
30.         // 如果超界则修改其值 
31.         if (dTemp < 0) 
32.             dTemp = 0.0; 
33.         else if (dTemp > 255) 
34.             dTemp = 255; 
35.  
36.         // 四舍五入 
37.         map[i] = int(dTemp + 0.5); 
38.     } 
39.  
40.     for (i = 0; i < nWidth * nHeight; i++) 
41.     {    
42.         for (j = 0; j < 4; j++) 
43.             pDst[i*4+j] = map[ pSrc[i*4+j] ]; 
44.     } 
45.  
46.     return EXIT_SUCCESS; 
47. }