IDL_植被变化监测.txt

PRO ChangeDetection

  COMPILE_OPT idl2
  e=ENVI()

  ;第一和第二时相数据
  inputFile1 = 'D:\ENVI_outputdata\LC8_20160712_Subset'
  inputFile2 = 'D:\ENVI_outputdata\LC8_20180429_Subset'

  ;输出文件路径
  outputFile = e.GetTemporaryFilename()
  PRINT, outputFile

  Raster1 = e.OPENRASTER(inputFile1)
  Raster2 = e.OPENRASTER(inputFile2)

  View = e.GetView()
  Layer1 = View.CreateLayer(Raster1)
  Layer2 = View.CreateLayer(Raster2)
  View.Zoom, 1, /FULL_EXTENT

  ;快速大气校正
  QUACRaster1 = ENVIQUACRASTER(Raster1)
  QUACRaster2 = ENVIQUACRASTER(Raster2)

  ;计算NDVI
  Index = 'NDVI'
  Subset1 = ENVISpectralIndexRaster(Raster1, Index)
  Subset2 = ENVISpectralIndexRaster(Raster2, Index)