Read HDF5/MintPy file in Python¶

The file structure is described in https://mintpy.readthedocs.io/en/latest/api/data_structure/.

In [1]:

%matplotlib inline
import os
import h5py

fdir = os.path.expanduser('~/data/archives/Galapagos/GalapagosSenDT128/mintpy/geo')

1. Time series¶

Run mask.py geo_timeseries_ERA5_ramp_demErr.h5 -m geo_maskTempCoh.h5 to get the masked time series file.

In [2]:

fname = os.path.join(fdir, 'geo_timeseries_ERA5_ramp_demErr_msk.h5')
with h5py.File(fname,'r') as f:
    data = f['timeseries'][:]
    date_list = f['date'][:]
date_list = [x.decode('utf8') for x in date_list]

where:

date_list is a list of strings in size of num_date in YYYYMMDD for the SAR acquisition date,
data is a 3D np.ndarray in size of (num_date, length, width) for the cumulative displacement with respect to the 1st acquisition (which has all zero values in the file) in the unit of meters.

2. Velocity¶

Run mask.py geo_velocity.h5 -m geo_maskTempCoh.h5 to get the masked velocity file.

In [3]:

fname = os.path.join(fdir, 'geo_velocity_msk.h5')
with h5py.File(fname,'r') as f:
    data = f['velocity'][:]

where data is a 2D np.ndarray in size of (length, width) for the average linear velocity in the unit of meters per year.

In [ ]: