GLASS数据批量处理

遥感数据处理系列

一些项目及科研中遇到的小需求，一方面记录自己的学习历程，另一方面帮助大家学习。本系列文章的开发环境为：ArcGIS 10.2.2 + Python 2.7、ENVI 5.3 + IDL 8.5

ArcPy批量计算栅格数据平均值
GLDAS数据下载及处理（NC转TIF）
ArcGIS批量裁剪栅格数据
ArcPy批量栅格重采样
ArcPy批量裁剪栅格数据

IDL多进程批处理遥感数据
ArcPy批量拼接栅格数据

---

---

前言

大佬们搞出了各种数据集，并且开放下载、全球共享、打破科研壁垒，促进研究进程。不过，真要用他们的数据集却遇到各种问题。在此记录下~

GLASS数据集的FVC数据在下载后如何批量处理呢？参考MODIS数据的处理流程进行数据处理，基于MRT工具 + IDL语言。

---

一、HDF文件整理

1. 原理简介

大概思路为：研究区域覆盖多个轨道，也就包含多个HDF文件，需要进行按天的文件分类，然后进行每天的多个HDF文件的拼接、裁剪、重投影。

文件路径：

F:\GLASS\GLASS10E01.V40.A2008001.h27v06.2019363.hdf
F:\GLASS\GLASS10E01.V40.A2008001.h27v07.2019363.hdf
F:\GLASS\GLASS10E01.V40.A2008001.h28v07.2019363.hdf
F:\GLASS\GLASS10E01.V40.A2008001.h28v08.2019363.hdf

参数简介：

	GLASS：
		HDF文件的上一级文件夹
	A2008001： 
		目标文件日期
	h27v06：
		目标区域轨道号

合并以上四个轨道的HDF文件即可实现HDF文件的按天整理。

2. 代码

pro modfiles

  ;MRT Tools第一步
  ;输入文件夹下，所有“.hdf”文件路径
  ;输出的.hdf文件路径存放在.txt文件中，
  ;.txt文件名为天数
  ;.txt文件内容为该天的所有.hdf文件路径

  path='F:\GLASS\'  ;output path
  files=file_search(path,'*.hdf')  ;the files needed to be processed

  v06h27files=file_search(path,'GLASS10E01.V*h27v06*.hdf')
  v07h27files=file_search(path,'GLASS10E01.V*h27v07*.hdf')
  v07h28files=file_search(path,'GLASS10E01.V*h28v07*.hdf')
  v08h28files=file_search(path,'GLASS10E01.V*h28v08*.hdf')

  nf=n_elements(v06h27files)

  for i= 0 ,nf-1 do begin
    basename=file_basename(v06h27files[i])
    Text0='F:\GLASS\1\'+strmid(basename,16,7)+'.txt'   ;creat a *.txt file of every (h27+h28)

    openW,lun,Text0,/Get_lun

    t4=v06h27files[i]
    t5=v07h27files[i]
    
    t8=v07h28files[i]
    t9=v08h28files[i]

    printf,lun,t4
    printf,lun,t5
    
    printf,lun,t8
    printf,lun,t9

    free_lun,lun
  endfor

  print,'endi'

end

在 “ F:\GLASS\1 ” 文件夹下生成所有天数的txt文件，每个txt文件包含多个轨道的HDF文件。

二、HDF文件拼接

1. 原理简介

大概思路为：按天分类后的多个HDF文件保存在同一个txt文件中，获取每天的txt文件路径，结合MRT tools的mrtmosaic工具，构造出使用MRT进行MOSAIC的批处理命令。

批处理命令：

E:\Software\MRT\bin\mrtmosaic -i F:\GLASS\1\2008001.txt -o F:\GLASS\1\2008001_do_GLASS.hdf -s " 1 "

参数简介：

	E:\Software\MRT\bin\mrtmosaic：
		mrtmosaic工具安装路径
	-i： 
		文件读取选项
	F:\GLASS\1\2008001.txt：
		待处理HDF文件的集合
	-o： 
		文件输出选项
	-s " 1 "：
		标签选项。如果要提取的数据为分层存储，则按需更改标签数量

2. 代码

pro buildmosaic

  ;MRT Tools第二步
  ;调用MRT的mrtmosaic工具，生成批处理文件“mdo.bat”
  ;双击批处理命令后，仍然生成的是hdf文件

  ;Input file path
  inf_path='F:\GLASS\1\'
  ;Out flie path
  out_path='F:\GLASS\1\'
  
  files=file_search(inf_path,'*.txt')
  Text0=out_path+'mdo.bat'
  openW,lun,Text0,/Get_lun
  for i= 0 , n_elements(files)-1 do begin
    basename=file_basename(files[i])
    txt=strmid(basename,0,7)+'_do'
    Outfile = out_path + txt + '_GLASS.hdf'  
     
    ;GLASS FVC
    Text="E:\Software\MRT\bin\mrtmosaic -i " + files[i] + " -o " + Outfile + " -s " +'" 1 "' 

    ;change the 'bin' path and '1 1 0 0 0'
    print,text
    printf,lun,Text

  endfor
  free_lun,lun
  print,'endi'

end

在 “ F:\GLASS\1 ” 文件夹下生成mdo.bat文件，双击该批处理文件即可完成HDF文件拼接。

总结

MRT批处理文件的自动化实现。体会行政审批流程，抉择不同生活方式，面对不同人生挑战。

后记

写博客的初衷是分享我的一些经验，同时也方便自己在其他电脑上进行数据处理。帮了很多人，但评论区小伙伴也有遇到问题的，那么：知识付费，我的时间和经验正好可以解决你的问题。