【注】经试验,本程序要在ENVI5.3 SP1以上版本运行,编写过程参考了ENVI-IDL官方博客http://blog.sina.com.cn/s/blog_764b1e9d0102xxrk.html的示例代码
去年7月提供了FLAASH的接口,官博给出了利用接口处理多光谱数据Landsat 8以及高光谱数据AVIRIS的示例代码,由于个人需要,我尝试了试用此接口处理MODIS数据,MODIS数据下载网址为https://ladsweb.modaps.eosdis.nasa.gov/search/
程序使用MYD02QKM Level1B数据,地面分辨率为250m,包含modis的前两个波段:红波段(band1)、近红外波段(band2)。要注意,FLAASH工具以及接口均要求BIL或者BIL存储格式的数据,经MCTK工具定标后的MODIS数据为BSQ格式,需要借助ENVI工具或者IDL实现转换,IDL转换可直接利用enviRaster.export,URL,Format[, keywords = value]进行转换,keywords中的interleave可实现bsq,bil,bip三种输出格式。
eg:enviRaster.export,outputfile,'dat', interleave = 'BIL'
核心代码如下:
【注】若为landsat 8数据,光谱响应函数位置有所不同