ENVI/IDL波段运算

完成对波段的一下运算：

1、将数组中小于0的赋值为-999，代码如下：

(B1 LT 0)*(0-999)+(B1 GE 0)*B1

 

2、求三个波段的平均值，如果某波段小于0，则不参与运算，代码如下： 

(B1>0+B2>0+B3>0) / (((B1 GE 0) + (B2 GE 0)+(B3 GE 0)) >1)

 

3、不同波段两幅图像，b1中的云部分（象元大于200） 用b2代替，代码如下：

(B1 GT 200) * B2+(B1 LE 200) * B1

完成上述三次处理的代码示例： 

PRO WORK1

ENVI_OPEN_FILE,R_FID=FID
ENVI_FILE_QUERY,FID,DIMS=DIMS,NB=NB,NS=NS,NL=NL

;ENVI_GET_DATA(FID,DIMS,POS)
B1=ENVI_GET_DATA(FID=FID,DIMS=DIMS,POS=0)
B2=ENVI_GET_DATA(FID=FID,DIMS=DIMS,POS=1)
B3=ENVI_GET_DATA(FID=FID,DIMS=DIMS,POS=2)
;将数组中小于0的赋值为-999
A=(B1 LT 0)*(0-999)+(B1 GE 0)*B1

;求三个波段的平均值，如果某波段小于0，则不参与运算
B=(B1>0+B2>0+B3>0) / (((B1 GE 0) + (B2 GE 0)+(B3 GE 0)) >1)

; 两幅图像，B1中的云部分（象元大于200） 用B2代替
C=(B1 GT 200) * B2+(B1 LE 200) * B1

;输出到内存中
ENVI_ENTER_DATA, A, R_FID = FID
ENVI_ENTER_DATA, B, R_FID = FID
ENVI_ENTER_DATA, C, R_FID = FID

END

 

注：新手学习，欢迎交流！