环境:ENVI5.3+IDL8.5
代码:
此段代码主要用于sentinel-2波段批量合成
以合成 band2,3,4,8,8a,11为例:
1. 定义函数
FUNCTION Layerstack, inpath, inname,outpath,outname
;used for layer stacking
outpath_name=FILEPATH(outname,ROOT_DIR=outpath)
IF FILE_TEST(outpath_name) EQ 1 THEN RETURN, 0
e = ENVI()
; test whether file exist
FOR i=0, size(inname,/N_ELEMENTS)-1 DO BEGIN
localpath=FILEPATH(inname[i],ROOT_DIR=inpath)
IF FILE_TEST(localpath) EQ 0 THEN RETURN, 0 ;IF THE INPUT FILE NOT EXIST
ENDFOR
; Select a Landsat TM scene from 1985
localpath=FILEPATH(inname[0],ROOT_DIR=inpath)
INPUT_RASTER=e.OpenRaster(localpath)
FOR i=1, SIZE(inname,/N_ELEMENTS)-1 DO BEGIN
localpath=FILEPATH(inname[i],ROOT_DIR=inpath)
INPUT_RASTER=[INPUT_RASTER,e.OpenRaster(localpath)]
ENDFOR
; Ge