【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数据如下所示:

图:更新后的数据

posted @   ENVI-IDL技术殿堂  阅读(301)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· DeepSeek 开源周回顾「GitHub 热点速览」
· 物流快递公司核心技术能力-地址解析分单基础技术分享
· .NET 10首个预览版发布:重大改进与新特性概览!
· AI与.NET技术实操系列(二):开始使用ML.NET
· 单线程的Redis速度为什么快?
点击右上角即可分享
微信分享提示