Чтение набора данных группы HDF5 с использованием Fortran - PullRequest
Новая
       19

Чтение набора данных группы HDF5 с использованием Fortran

0 голосов
/ 05 июля 2019

Я застрял при чтении набора данных членов группы файла HDF5 с использованием Fortran.

Я могу перечислить члена группы моего HDF5. Но я не могу получить доступ к данным участника в группе, 

program sds_info
use hdf5
implicit none


! Variables declaration
  CHARACTER*100 :: file_name
  CHARACTER*100 :: sds_name
  CHARACTER*100 :: gr_name
  INTEGER(HID_T):: file_id, gr_id, dset_id, attr_id
  INTEGER :: status, error, storage, nlinks,max_corder, attr_num

  REAL, DIMENSION(1) :: dset_data, data_out
  INTEGER, DIMENSION(1)   ::    buf

  INTEGER(HSIZE_T), DIMENSION(1)::  data_dims
  INTEGER(HSIZE_T), DIMENSION(1) ::dims
!
! varaibles to read a dataset in a group
  CHARACTER*100 :: ap_name
  integer(HID_T):: ap_id
  real, allocatable, dimension(:) :: ap
  integer(HSIZE_T), dimension(15624960) :: ap_dim
  integer :: nmembers ! Number of group members
  CHARACTER(LEN=20) :: name_buffer ! Buffer to hold object's name
  integer :: i
  integer :: type
! 
! Variables initalization
  file_name = "PVAR8.h5"
  sds_name = "time"
  gr_name = "part"
  attr_name = "attr1"
  ap_name="ap"

! Initialize the interface
  call h5open_f(status)
! Open an hdf5 file
  call h5fopen_f(file_name, H5F_ACC_RDWR_F, file_id, status)


! Open a group
  call h5gopen_f(file_id, gr_name, gr_id, status )
!  
! Open a dataset
  call h5dopen_f(file_id, sds_name, dset_id, error)
! Get the number of attributes
  call h5aget_num_attrs_f(dset_id, attr_num, error)
  print *, "attr_num ",attr_num
! Read the dataset
  call h5dread_f(dset_id, H5T_NATIVE_REAL, data_out, data_dims, error)
  print *, "data_out ",data_out


! Terminate access to the group
  call h5gclose_f(gr_id, error)
! Terminate access to the dataset
  call h5dclose_f(dset_id, error)
! Terminate access to the file
  call h5fclose_f(file_id, error)
! Close FORTRAN interface.
  call h5close_f(status)


end program sds_info

Я могу читать группу, но как получить доступ и прочитать данные члена группы в HDF5 с помощью Fortran?

1 Ответ

1 голос
/ 05 июля 2019

Если вы не обязаны использовать определенную технологию, посмотрите на HDFql , чтобы решить ваш вопрос.

Используя HDFql в Fortran, вы можете прочитать набор данных time (с типом данных real) хранится в группе part из файла PVAR8.h5, например: 

PROGRAM Test

    ! use HDFql module (make sure it can be found by the Fortran compiler)
    USE HDFql

    ! declare variables
    REAL(KIND = 8) :: value
    INTEGER :: state

    ! register variable "value" for subsequent use (by HDFql)
    state = hdfql_variable_transient_register(value)

    ! select (i.e. read) data from dataset "time" and populate variable "value" with it
    state = hdfql_execute("SELECT FROM PVAR8.h5 /part/time INTO MEMORY 0")

    ! display content of variable "value"
    WRITE(*, *) "Dataset value:", value

END PROGRAM

Дополнительные примеры HDFql можно найти здесь .

...