%matplotlib inline
import os
import numpy as np
from matplotlib import pyplot as plt

work_dir = os.path.expanduser('~/data/test/FernandinaSenDT128/mintpy/geo')
os.chdir(work_dir)
print('Go to directory', work_dir)

fname = os.path.join(work_dir, 'geo_velocity.grd')

# read
import gdal
ds = gdal.Open(fname, gdal.GA_ReadOnly)
data = np.array(ds.ReadAsArray(), dtype=np.float32)

# plot
plt.figure()
plt.imshow(data, cmap='jet')
plt.show()

# read
from scipy.io import netcdf_file as netcdf
with netcdf(fname, 'r') as f:
    width, length = f.variables['dimension'][:]
    data = np.array(f.variables['z'][:].reshape(length, width), dtype=np.float32)

# plot
plt.figure()
plt.imshow(data, cmap='jet')
plt.show()