【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 @   地理遥感生态网平台  阅读(17)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· winform 绘制太阳,地球,月球 运作规律
· AI与.NET技术实操系列(五):向量存储与相似性搜索在 .NET 中的实现
· 超详细:普通电脑也行Windows部署deepseek R1训练数据并当服务器共享给他人
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 上周热点回顾(3.3-3.9)
点击右上角即可分享
微信分享提示