SMAP L4 dataΒΆ

This example code illustrates how to access and visualize a SMAP L4 data.

#Add data file
print 'Read data...'
fn = 'D:/Temp/hdf/SMAP_L4_SM_gph_20150414T013000_Vb1010_001.h5'
f = addfile(fn)
lon = f['cell_lon'][:,:]
lat = f['cell_lat'][::-1,:]
#vname = 'surface_temp'
#vname = 'sm_surface_wetness'
#vname = 'sm_surface'
vname = 'sm_profile'
v = f[vname]
data = v[::-1,:]
longname = v.attrvalue('long_name')[0]
units = v.attrvalue('units')[0]

#Interpolate data to grid
print 'Interpolate data to grid...'
lon1 = linspace(lon.min(), lon.max(), lon.dimlen(1))
lat1 = linspace(lat.min(), lat.max(), lat.dimlen(0))
data1 = griddata((lon, lat), data, xi=(lon1, lat1), method='surface')[0]

#Plot
print 'Plot...'
axesm()
geoshow('country', edgecolor='k')
layer = imshow(lon1, lat1, data1, 20, cmap='wcgyr_1000', interpolation='bilinear')
colorbar(layer, label=units)
title('{0}\n {1}'.format(fn, longname))
print 'Finish!'
../../../_images/smap_l4.png