npsphinx prolog… I dont know what to write here….

get corners from the tiles in each bursts of the sub-swath

[5]:

import numpy as np
from matplotlib import pyplot as plt
from shapely.geometry import Polygon
all_poly = {}
all_poly_coords = {}
for bursti in range(ds.burst.size):
    for i_tile_sample in range(ds.tile_sample.size):
        for i_tile_line in range(ds.tile_line.size):
            #print('bursti',bursti,i_tile_sample,i_tile_line)
            lonco  = ds['corner_longitude'].isel({'burst':bursti,'tile_sample':i_tile_sample,'tile_line':i_tile_line}).values.ravel(order='C')
            latco  = ds['corner_latitude'].isel({'burst':bursti,'tile_sample':i_tile_sample,'tile_line':i_tile_line}).values.ravel(order='C')
            coords = np.stack([lonco,latco]).T
            #print('coords',coords)
            coords = coords[[0,1,3,2,0]]
            #print('coords2',coords)
            pp = Polygon(coords)
            #all_poly.append(pp)
            all_poly['burst%s_sample%s_line%s'%(bursti,i_tile_sample,i_tile_line)] = pp
            all_poly_coords['burst%s_sample%s_line%s'%(bursti,i_tile_sample,i_tile_line)] = coords
print('done',len(all_poly.keys()))

done 45

[6]:

import geoviews as gv
import holoviews as hv
gv.extension('bokeh')
from shapely.geometry import Polygon
import shapely.wkt

P = shapely.wkt.loads(ds.footprint)
#poly = Polygon(np.vstack([lons_tiles,lats_tiles]).T)
print(P)
foot = gv.Path(P).opts(color='magenta')
tmp =[]
for uu in  all_poly:
    tmp.append(gv.Path(all_poly[uu]))
(gv.tile_sources.Wikipedia * foot * hv.Overlay(tmp) ).opts(width=800,height=600,title='%s tiles SAR Sentinel-1 L1B IFREMER '%(cat_xspec))

POLYGON ((16.52569318605885 34.48393468420786, 15.47003602432591 34.63859437730419, 15.1537098076834 33.12316072296198, 16.19010348915945 32.96773923777913, 16.52569318605885 34.48393468420786))

/tmp/ipykernel_3173/1846620382.py:14: GeoviewsDeprecationWarning: 'Wikipedia' is deprecated and will be removed in version 1.11, use 'OSM' instead.
  (gv.tile_sources.Wikipedia * foot * hv.Overlay(tmp) ).opts(width=800,height=600,title='%s tiles SAR Sentinel-1 L1B IFREMER '%(cat_xspec))

[6]:

Example to display a IW L1B x-spectra

read L1B data

pick up intra or inter burst group

recombine real part and imaginary part (as recommanded by Fred)

get corners from the tiles in each bursts of the sub-swath

display pyplot figure of cross spectrum