ENVI下DOM数据生产专题之图像增强

      多光谱影像彩色合成方法主要分为2种：自然真彩色合成和非自然假彩色合成。自然真彩色合成是指合成后的彩色影像上地物色彩与实际地物色彩接近或者一致，一般的方法就是多光谱影像的红、绿、蓝对应R/G/B合成；非自然假彩色则反之。

      遥感影像自然真彩色合成可分为以下几种方法：1、直接用多光谱影像的红、绿、蓝通道合成，一般用于高分辨率影像；2、利用其它波段加权处理，重新生成红、绿、蓝波段，一般用于增强某种地物颜色层次，如植被、水体等；3、利用其它波段信息重新生成某一波段，一般用于缺少波段的传感器，如SPOT、Aster等。

1)      波段加权运算(包含R、G、B、NIR波段的数据)

     最常见的就是增强植被信息，可使用绿色和近红外波段加权运算。如下公式：

          Band_new=a*B_gree+（1-a）*B_nir

         a是权重值，取0~1。

下面使用ENVI下的Bandmath和layerstacking工具，利用TM影像合成真彩色图像。

(1)     打开包括R、G、B、NIR波段的TM影像。

(2)     选择Basic tools->Bandmath，输入表达式：byte(b2*0.8+b4*0.2)，选择相应的波段执行运算。

(3)     在波段列表中，选择：R-红色波段，G-合成波段，B-蓝色波段显示，可以看到合成想效果。

(4)     Basictools->layerstacking，将相应三个波段组成一个三波段图像文件。

                             

另外，为了让植被之外的地物颜色更加真实，可以只对植被进行增强，这里使用NDVI对植被进行区分。首先先计算NDVI，使用以下波段运算表达式进行加强运算：

（B3 gt 0.2）*(b2*0.8+b4*0.2)+（B3 le 0.2）*b2

B3：NDVI。

2)    波段生成

常见的就是SPOT影像蓝色波段的生成，下面是SPOT常见的方法：

(1)     原来的绿波段（0.50-0.59 μm）当作蓝波段（该波段靠近蓝波段的光谱范围），红波段（0.61-0.68μm）仍采用原来的波段，绿波段用绿波段、红波段、近红外波段的算术平均值来代替。

R：XS2

G：（XS1+XS2+ XS3）/3，bandmath表达式：byte((fix(b1)+b2+b3)/3)

B：XS1

   效果如下：

 

(2)     将原来的绿波段（0.50-0.59 μm）当作蓝波段，红波段（0.61-0.68μm）仍采用原来的波段，绿波段用绿波段、近红外波段按3：1的加权算术平均值来代替。

R：XS2

G：（XS1×3＋XS3）/4，    bandmath表达式：(b1*3+b3)/4

B：XS1

效果如下：

 

以上两种由于其算法是比较相近的，所以生成影像的色彩效果区别不大。

以上方法不仅限于SPOT影像，也可用于ASTER影像，甚至用于其他包括绿波段、红波段、近红外波段的数据。