ENVI/IDL混编 批量波段合成 简洁函数

环境: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

  • 0
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值