Download this notebook from github.

SAR Sentinel-1 cross spectrum computation from large image Interferometric Wide-swath SLC

[1]:

import os
import xsarslc.get_test_files
import logging
import xsar
from importlib import reload
reload(logging)
logging.basicConfig(level=logging.INFO)
logging.info('start notebook')
import xsarslc.processing.xspectra as proc
import warnings
warnings.simplefilter(action='ignore', category=FutureWarning)
warnings.filterwarnings("ignore")

INFO:root:start notebook

get datatree with digital numbers and annotations from IW SLC product

[2]:

import xsarslc
import xsar
import os
#one_tiff = '/home/datawork-cersat-public/cache/project/mpc-sentinel1/data/esa/sentinel-1a/L1/IW/S1A_IW_SLC__1S/2022/127/S1A_IW_SLC__1SDV_20220507T162411_20220507T162439_043107_0525DE_734D.SAFE/measurement/s1a-iw1-slc-vv-20220507t162411-20220507t162439-043107-0525de-004.tiff'
#tmp = xsar.get_test_file('S1A_IW_SLC__1SDV_20170907T103019_20170907T103047_018268_01EB76_5F55.SAFE')
fullpathsafeSLC = xsar.get_test_file('S1A_IW_SLC__1SDV_20181008T134439_20181008T134506_024045_02A0A6_65C8.SAFE') # only one subswath IW1 zipped

#fullpathsafeSLC = os.path.dirname(os.path.dirname(one_tiff))
#imagette_number = os.path.basename(one_tiff).split('-')[1].replace('iw','')
imagette_number = 1
str_gdal = 'SENTINEL1_DS:%s:IW%s'%(fullpathsafeSLC,imagette_number)
print('str_gdal',str_gdal)
xsarobj = xsar.Sentinel1Dataset(str_gdal) #,resolution='30m'
dt = xsarobj.datatree
dt.load()
dt

str_gdal SENTINEL1_DS:/home1/scratch/agrouaze/xsardatasync/xsardata/S1A_IW_SLC__1SDV_20181008T134439_20181008T134506_024045_02A0A6_65C8.SAFE:IW1

[2]:

<xarray.DatasetView>
Dimensions:  ()
Data variables:
    *empty*
Attributes: (12/15)
    name:              SENTINEL1_DS:/home1/scratch/agrouaze/xsardatasync/xsar...
    short_name:        SENTINEL1_DS:S1A_IW_SLC__1SDV_20181008T134439_20181008...
    product:           SLC
    safe:              S1A_IW_SLC__1SDV_20181008T134439_20181008T134506_02404...
    swath:             IW
    multidataset:      False
    ...                ...
    start_date:        2018-10-08 13:44:40.984312
    stop_date:         2018-10-08 13:45:06.117599
    footprint:         POLYGON ((-114.6305305568416 32.60521729566573, -115.5...
    coverage:          171km * 88km (line * sample )
    orbit_pass:        Descending
    platform_heading:  -167.2011573025802

cross spectrum in inter bursts parts of the IW product

[3]:

inter_xs = proc.compute_IW_subswath_interburst_xspectra(dt,
                                        tile_width={'sample': 20.e3, 'line': 20.e3},
                                        tile_overlap={'sample': 0, 'line': 0},
                                        periodo_width = {'sample': 2000, 'line': 1200},
                                        periodo_overlap = {'sample': 1000, 'line': 600},
                                        dev=True,polarization='VV')
inter_xs

INFO:root:reduce number of burst -> 1
0it [00:00, ?it/s]

[3]: