IDL实现经纬度转行列号影像裁切

最近在做毕业设计，其中有一步需要利用经纬度裁切区域，作为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