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           

注：新手學習，歡迎交流！