IDL读写hdf5数据

转载于:http://rssthome.com/forum.php?mod=viewthread&tid=50&extra=page=1

pro hdf5
    filter=['*.HDF','*.*']
    files=DIALOG_PICKFILE(/MULTIPLE_FILES,FILTER=filter,title='打开hdf5数据',GET_PATH=path)
    CD,path
    if (files[0] EQ '') then return
    storage_path=dialog_pickfile(/directory,title='选择存储路径')
    n=n_elements(files)

    for i=0,n-1 do begin
    ;-----------------------------读取hdf5--------------------------------
    basename=file_basename(files[i],'.HDF')
    ;H5F_OPEN打开hdf5
    file_id = H5F_OPEN(files[i])
    ;H5G_OPEN打开group
    group_id=H5G_OPEN(file_id,'Data')
    ;H5D_OPEN打开数据集
    dataset_id=H5D_OPEN(group_id,"EV_1KM_RefSB")
    ;H5D_READ读取数据集数据,存入SdsData中
    SdsData=H5D_READ(dataset_id)
    weishu=size(SdsData,/dimensions)
    ;读取属性数据
    Attribute_num= H5A_GET_NUM_ATTRS(dataset_id)
    for j=0,Attribute_num-1 do begin
        ;H5A_OPEN_IDX按照属性索引打开属性
        Attribute_id= H5A_OPEN_IDX(dataset_id,j)
        ;H5A_OPEN_Name按照属性名字打开属性
        ;Attribute_id= H5A_OPEN_Name(dataset_id,'units')
        ;得到属性的名字
        Attribute_name= H5A_GET_NAME(Attribute_id)
        ;读取属性
        AttData=H5A_READ(Attribute_id)
    endfor

    ;-----------------------------输出hdf5--------------------------------
    outname=storage_path+basename+'_out2.HDF'
    ;创建hdf5文件
    outfid=H5F_CREATE(outname)
    ;创建group
    OutGroupId = H5G_CREATE(outfid,"FrameSampleMeasurements")     
    ;创建数据集
     datatype_id = H5T_IDL_CREATE(SdsData)
    dataspace_id = H5S_CREATE_SIMPLE(weishu)
    Outdataset_id = H5D_CREATE(OutGroupId,"sample_measurements_o2",datatype_id,dataspace_id)
    OutData=SdsData
    H5D_WRITE,Outdataset_id,OutData
    ;写数据集的属性
    Attribute_num=H5A_GET_NUM_ATTRS(dataset_id)
    for k=0,Attribute_num-1 do begin
        ;获取属性
        sds_att_id=H5A_OPEN_IDX(dataset_id,k)
        sds_att_name=H5A_GET_NAME(sds_att_id)
        sds_att_value=H5A_READ(sds_att_id)
        sds_att_type=H5A_GET_TYPE(sds_att_id)
        ;写出属性
        ;属性dataspace的大小必须和属性包含的内容的个数一致
        dataspace_id = H5S_CREATE_SIMPLE(n_elements(sds_att_value))
        attr_id = H5A_CREATE(Outdataset_id,sds_att_name,sds_att_type,dataspace_id)
        H5A_WRITE,attr_id,sds_att_value
        H5A_CLOSE,attr_id
    endfor
  • 2
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
针对您的问题,我可以提供以下步骤: 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进行几何校正的步骤。

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

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

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值