【IDL】重写HDF文件中数据集的方法

IDL中提供了一系列HDF_*函数对HDF文件进行读写、查询等操作。

遇到一个需求，使用IDL来重写HDF中的某个数据集。其实非常简单。只需要用到如下几个关键接口：

• HDF_SD_START
• HDF_SD_SELECT
• HDF_SD_GETDATA
• HDF_SD_ADDDATA

以一景MODIS MOD02数据为例，里边包含了经纬度信息。查看纬度（Latitude）数据如下所示，下面我们将纬度数据全部加1，重新写入到原始HDF文件中。

图：备份的原始HDF文件中的Latitude数据

运行如下代码：

file='D:\temp\MOD021KM.A2002248.0345.005.2007348121959.hdf'
;打开HDF文件，获取SD ID，注意使用关键字RDWR，可以读写HDF文件。
sd_id=HDF_SD_START(file, /RDWR)
;获取SD中Latitude的索引
latIndex = HDF_SD_NAMETOINDEX(sd_id,'Latitude')
;获得Latitude所在索引的ID
LatSD_ID = HDF_SD_SELECT(sd_id, latIndex)
;读取Latitude数据
HDF_SD_GETDATA, LatSD_ID, LatData
;处理Latitude数据
tempLatData = LatData+1
;将新数据写出
HDF_SD_ADDDATA, LatSD_ID, tempLatData
;关闭SD，关闭HDF文件
HDF_SD_ENDACCESS, LatSD_ID
HDF_SD_END, sd_id

运行之后，新的HDF数据如下所示：

图：更新后的数据

标签：数据,HDF,IDL,Latitude,sd,重写,ID,SD
来源： https://www.cnblogs.com/enviidl/p/16378471.html