ENVI_IDL:如何读取HDF4数据?

我将以MODIS气溶胶产品进行一些简单的分析.如果你觉得稍有繁琐, 你可以直接阅读下方关于HDF4的文件读取的一些常见函数说明.

file_id = hdf_sd_start(file_path, /read)

功能: 打开HDF4文件,返回该文件的id(内存地址,以后访问该文件可以通过其)

解释: 第一个参数传入hdf4文件的绝对路径(如果该hdf4文件在目前pro文件所在文件夹中, 那么相对路径也是可以的),第二个参数传入打开hdf4文件的方式,除了目前传入的/read(只读),你还可以传入/rdwr(可读可写),/create(创建新的hdf4文件)。函数返回该文件的内存地址

hdf_sd_info, file_id, dataset_num(接收数据的变量), attr_num(接收数据的变量)

功能:获取hdf4文件内的数据集个数,全局属性个数

解释:第一个参数传入hdf4文件的id;第二个参数传入一个变量名,用于接收该文件的数据集个数;第三个参数传入一个变量名,用于接收该文件的全局属性个数。</

  • 0
    点赞
  • 6
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 5
    评论
针对您的问题,我可以提供以下步骤: 1. 打开ENVI,选择File -> Open External File,选择HDF5数据文件。 2. 打开IDL环境,在命令行中输入以下脚本: ``` envi, /restore_base_save_files file = FILEPATH('input_file.h5', SUBDIRECTORY=['data'], $ ROOT_DIR=['/']) hdf5 = HDF(/read, FILENAME=file) data = hdf5->GetData() metadata = hdf5->GetMetadata() lat = metadata[0].Data.Latitude lon = metadata[0].Data.Longitude ``` 这个脚本将读取HDF5数据文件并提取数据和元数据。其中,lat和lon是数据的经纬度信息。 3. 创建GLT文件。在IDL中输入以下脚本: ``` gltpath = FILEPATH('output_glt.glt', SUBDIRECTORY=['data'], $ ROOT_DIR=['/']) glt = ENVIGeometricCorrectionRaster(gltpath, $ DATASET_NAME='input_file', $ SPATIAL_REF_OBJ=hdf5->GetSpatialRef(), $ SPATIAL_REF_CODE=hdf5->GetSpatialRefCode(), $ PIXEL_SIZE=hdf5->GetPixelSize(), $ UL_LATITUDE=max(lat), $ UL_LONGITUDE=min(lon), $ ROTATION_ANGLE=0.0, $ OUTPUT_PIXEL_SIZE=hdf5->GetPixelSize(), $ OUTPUT_COORDINATE_SYSTEM_CODE=hdf5->GetSpatialRefCode(), $ OUTPUT_INTERLEAVE=hdf5->GetInterleave(), $ OUTPUT_DATA_TYPE=hdf5->GetDataType()) ``` 这个脚本将创建一个GLT文件并设置几何校正参数。其中,glt是一个ENVIGeometricCorrectionRaster对象,表示GLT文件;gltpath是GLT文件的保存路径;DATASET_NAME是输入数据集的名称;SPATIAL_REF_OBJ和SPATIAL_REF_CODE是输入数据的空间参考对象和空间参考代码;PIXEL_SIZE是输入数据的像素大小;UL_LATITUDE和UL_LONGITUDE是输入数据的左上角经纬度坐标;ROTATION_ANGLE是输入数据的旋转角度;OUTPUT_PIXEL_SIZE、OUTPUT_COORDINATE_SYSTEM_CODE、OUTPUT_INTERLEAVE和OUTPUT_DATA_TYPE是输出数据的相关参数。 4. 运行几何校正。在IDL中输入以下脚本: ``` outputpath = FILEPATH('output_file.h5', SUBDIRECTORY=['data'], $ ROOT_DIR=['/']) geocorrect = ENVIGeometricCorrection(input_file, $ OUTPUT_URI=outputpath, $ GLT=glt, $ SPATIAL_REF_OBJ=hdf5->GetSpatialRef(), $ SPATIAL_REF_CODE=hdf5->GetSpatialRefCode(), $ PIXEL_SIZE=hdf5->GetPixelSize(), $ OUTPUT_INTERLEAVE=hdf5->GetInterleave(), $ OUTPUT_DATA_TYPE=hdf5->GetDataType()) ``` 这个脚本将执行几何校正操作,并将结果保存到指定路径的HDF5文件中。其中,geocorrect是一个ENVIGeometricCorrection对象,表示几何校正操作;input_file是输入数据文件;OUTPUT_URI是输出数据文件的保存路径;GLT、SPATIAL_REF_OBJ、SPATIAL_REF_CODE、PIXEL_SIZE、OUTPUT_INTERLEAVE和OUTPUT_DATA_TYPE是与GLT文件相同的参数。 以上就是使用ENVI IDLHDF5数据建立GLT进行几何校正的步骤。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 5
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

炒茄子

不装逼我浑身难受aaa

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值