堆叠多个Landsat图像

时间:2018-05-19 10:07:00

标签: python bash gdal landsat rasterio

我想在特定format中堆叠大约700张Landsat图像作为YATSM连续变化检测和分类(CCDC)算法的输入。我尝试使用此bash脚本4_stack_TIFF.sh。在此脚本中,使用此文件landsat_stack.py完成堆叠。我有两个错误:

  • 首先与缺少模块module load batch_landsat有关。我可以找到这个共享的存储库

  • 其次与landsat_stack.py有关,我得到以下内容: Stacking to extent: Traceback (most recent call last): File "./landsat_stack.py", line 936, in <module> sys.exit(main()) File "./landsat_stack.py", line 920, in main utm, resume, fformat, creation_opts)) File "./landsat_stack.py", line 757, in landsat_stack print('\tUpper Left: {ulx},{uly}'.format(ulx=extent[0], uly=extent[1])) TypeError: 'NoneType' object is not subscriptable Done!

问题:

  1. 是否有人使用过YATSM并在地理数据预处理中遇到同样的问题?
  2. 如果没有;你能建议我使用Rasterio或Rio吗?
  3. 格式:

    8 band images: 
       Band 1 SR (SR * 10000)
       Band 2 SR (SR * 10000)
       Band 3 SR (SR * 10000)
       Band 4 SR (SR * 10000)
       Band 5 SR (SR * 10000)
       Band 7 SR (SR * 10000)
       Band 6 Thermal Brightness (C * 100)
       Fmask
           0 - clear land
           1 - clear water
           2 - cloud
           3 - snow
           4 - shadow
           255 - NoData
    

0 个答案:

没有答案