Download this notebook from github.
illustration of the default Impulse Response (IR) files from xsarslc
[1]:
import xsarslc.get_config_infos
from importlib import reload
reload(xsarslc.get_config_infos)
import holoviews as hv
hv.extension('bokeh')
import xarray as xr
import os
[2]:
lst_nc = []
for unit in ['S1A','S1B']:
for sw in ['IW1','IW2','IW3']:
for pol in ['VV']:
fileir = xsarslc.get_config_infos.get_IR_file(unit=unit, subswath=sw, polarization=pol)
lst_nc.append(fileir)
print(fileir)
/home1/datahome/agrouaze/sources/git/xsar_slc/auxdata/S1A_IRs_IW1_VV.nc
/home1/datahome/agrouaze/sources/git/xsar_slc/auxdata/S1A_IRs_IW2_VV.nc
/home1/datahome/agrouaze/sources/git/xsar_slc/auxdata/S1A_IRs_IW3_VV.nc
/home1/datahome/agrouaze/sources/git/xsar_slc/auxdata/S1B_IRs_IW1_VV.nc
/home1/datahome/agrouaze/sources/git/xsar_slc/auxdata/S1B_IRs_IW2_VV.nc
/home1/datahome/agrouaze/sources/git/xsar_slc/auxdata/S1B_IRs_IW3_VV.nc
display IR in range dimension
[3]:
vari = 'range_IR'
all_cu = []
filter_out = ['HH','HV','VH']
for ff in lst_nc:
go = True
for fi in filter_out:
if fi in ff:
print(ff,'filtered out with',fi)
go = False
else:
pass
if go:
ds = xr.open_dataset(ff)
tmpcu = hv.Curve(ds[vari],label=os.path.basename(ff)).opts(tools=['hover'],show_legend=True)
all_cu.append(tmpcu)
print(len(all_cu))
hv.Overlay(all_cu).opts(width=1000,height=750,tools=['hover'],show_grid=True,show_legend=True)
6
[3]:
display IR in azimuth dimension
[4]:
vari = 'azimuth_IR'
all_cu = []
filter_out = ['HH','HV','VH']
for ff in lst_nc:
go = True
for fi in filter_out:
if fi in ff:
print(ff,'filtered out with',fi)
go = False
else:
pass
if go:
ds = xr.open_dataset(ff)
tmpcu = hv.Curve(ds[vari],label=os.path.basename(ff)).opts(tools=['hover'],show_legend=True)
all_cu.append(tmpcu)
print(len(all_cu))
hv.Overlay(all_cu).opts(width=1000,height=750,tools=['hover'],show_grid=True,show_legend=True)
6
[4]: