ENVI/IDL混合编程-批量影像拼接

;README
;used for batch mosaicing with sentinel-2 images (band by band)
;usage: function_mosaic(['/20200131_1/B1','/20200131_2/B1'],'/20200131/B1_mosaic')
;requirement: ENVI/IDL 5.3 or higher
;reference: the reference book of envi/idl software
;coded by Yinxiacao
;date: April 8, 2020

FUNCTION function_mosaic,files,outpath
  COMPILE_OPT IDL2
  ; Start the application
  e = ENVI()
  FOR i=0, N_ELEMENTS(files)-1 DO BEGIN
    IF FILE_TEST(files[i])  EQ 0 THEN RETURN, 0 ;IF THE INPUT FILE NOT EXIST
  ENDFOR
  ; Select input scenes 
  ;files = REVERSE(FILE_SEARCH(inpath, inname))
  scenes = !NULL
  FOR i=0, N_ELEMENTS(files)-1 DO BEGIN
    raster = e.OpenRaster(files[i])
    ;query and update the data ignore value
    metadata = raster.METADATA
    pos = WHERE(STRCMP(metadata.tags, 'data ignore value'