最近在做毕业设计,其中有一步需要利用经纬度裁切区域,作为IDL小白,看过好几个博客,实验了很多次,终于成功了,特此写下人生中第一篇csdn博客,纪念一下。特别感谢一下我的一位同学,他的博客写的很出色,IDL学的也特别好,毕设期间经常麻烦他。看到的小伙伴可以去看看他的博客,干货很多。https://www.ixxin.cn/
我用的IDL版本是IDL85,ENVI5.3 SP1
不过这个代码应该是各种版本基本都可以,对于FLAASH接口,就得用5.3以上才可以了
第一步是经纬度转行列号
PRO band_subset
COMPILE_OPT IDL2
e = envi(/headless)
file = 'F:\2013_04\OutDir\MYD02QKM.A20130403_radiance_georef.dat'
dir = 'F:\2013_04\out'
name = file_basename(file)
oFile = dir+'\'+name+'_subset.dat' ;设置输出文件的名称
MULRaster = envi.OpenRaster(File)
Lon=[119,121]
Lat=[36.5,35]
MapXs = fltarr(2)
MapYs = fltarr(2)
FileXs = fltarr(2)
FileYs = fltarr(2)
;经纬度转行列号
s