1. 课堂内容
OMI数据的读取与输出，这里实际考察如何取读取HDF5文件以及HDF5文件输出(这里输出是以tiff格式为例)
2. 知识储备
HDF5文件打开的相关函数
3. 编程
; 本函数用于打开HDF5文件并获取数据集 function get_hdf5_ds, file_path, ds_path ; 输入数据集所在文件的路径以及 --数据集的路径(关于这个看主程序的注释) ; 获取文件的id————》通过h5f_open()函数获取 file_id = h5f_open(file_path) ; 传入文件的路径 ; 获取数据集的id————》通过h5d_open()函数获取(另外需要说一下， h5d表示hdf(简称h)5(简称5)文件中的data set(简称d)的相关操作) ds_id = h5d_open(file_id, ds_path) ; 传入数据集所在文件的id，传入数据集的名称 ; 获取数据集的数据————》通过h5d_read()函数获取 data = h5d_read(ds_id) ; 传入数据集的id ; 打开的文件以及数据集需要全部关闭当不用的时候，这是习惯和态度问题(其实是为了后续你再使用这些文件避免一些错误和问题) h5d_close, ds_id ; 这里是h5d, 下面是h5f，两个函数是不一样的，需要注意一下 h5f_close, file_id ; 返回数据集的数据 return, data ; 销毁数据变量(避免空占内存——但是我是存疑的，return之后函数内部的所有变量不都会被销毁吗？还需要自己操作吗) data = !null end pro week_three_study1 ; 本程序解决如何读取OMI数据以及它的输出(OMI数据这里是HDF5格式,实际上是学习如何读取和输出HDF5文件) ; 获取文件路径以及数据集的路径(通过hdf explorer软件轻松获取) file_path = 'D:\IDL_program\experiment_data\chapter_2\NO2\2017\OMI-Aura_L3-OMNO2d_2017m0101_v003-2018m0627t042221.he5' ds_path = '/HDFEOS/GRIDS/ColumnAmountNO2/Data Fields/ColumnAmountNO2TropCloudScreened' ; 这里说明一下，hdf4不需要数据集在文件里面的目录，但是hdf5需要输入数据集在文件里面的目录(原因不知道，我有猜测但是不确定) ; 获取输出的文件夹路径 out_dir_path = 'D:\IDL_program\experiment_data\chapter_2\output' ; 这里假定你不知道路径中output文件夹是否存在，且你并不愿意亲自查看，可以进行如下操作 if file_test(out_dir_path, /directory) eq 0 then begin ; 这里有一个函数需要认识一下————》file_test()——》检查文件是否存在(默认检查文件，需要检查文件夹需要加上参数/directory) ; 该函数会返回0或者1，0表示不存在，1表示存在 ; 创建文件夹 file_mkdir, out_dir_path endif ; 获取输出的文件的名称 out_name = file_basename(file_path, '.he5') + '.tiff' ; 输出为tiff格式可以进行查看 ; 这里file_basename()传入file_path已经可以将文件名称输出，但是我们不需要后缀名.he5 ; 所以传入字符串'.he5'函数会将已经获取的文件名称中这个字符串删除并返回 ; 输出的路径 = 输出文件夹的路径 + 输出文件的名称 out_path = out_dir_path + '\' + out_name ; 获取数据集的数据 data = get_hdf5_ds(file_path, ds_path) ; 输出 write_tiff, out_path, data, /float ;传入输出的路径、输出的数据、另外这里默认是输出int型，但是我们通过hdf explorer发现该数据集的type是float，所需要指定一下 end
编译+运行之后输出的结果：
查看输出的TIFF文件：
4. 题外话(HDF5下的中文路径问题)
这里我遇到一个小小的不兼容，原来我们使用hdf_sd_start()函数打开HDF4文件时
即使传入的路径是中文，那么也是可以打开这个文件的。
但是，在h5f_open()函数中打开HDF5文件传入中文路径就会运行报错，显示无法打开文件，只需要将中文路径修改为英文路径即可(所以，养成写英文路径的习惯吧，也算学了学英语)
报错显示如下：
