@@ -271,7 +271,7 @@ def specular_plot(self, cmap='YlOrRd', title='CYGNSS data', vmin=0,
271271 axis_label_flag = False , title_flag = True , indices = None ,
272272 save = None , lonrange = None , latrange = None ,
273273 truth_flag = False , return_flag = False , gpsid = None ,
274- gain = None , sat = None ):
274+ gain = None , sat = None , ** kwargs ):
275275 """
276276 Plots CYGNSS specular points on lat/lon axes using matplotlib's scatter
277277 object, which colors each point based on its wind speed value.
@@ -321,15 +321,16 @@ def specular_plot(self, cmap='YlOrRd', title='CYGNSS data', vmin=0,
321321 if basemap is None :
322322 sc = ax .scatter (ds .lon [ds .good ], ds .lat [ds .good ], c = ws [ds .good ],
323323 vmin = vmin , vmax = vmax , cmap = cmap , s = ms ,
324- marker = marker , edgecolors = ec )
324+ marker = marker , edgecolors = ec , ** kwargs )
325325 if lonrange is not None :
326326 ax .set_xlim (lonrange )
327327 if latrange is not None :
328328 ax .set_ylim (latrange )
329329 else :
330330 x , y = basemap (ds .lon [ds .good ], ds .lat [ds .good ])
331331 sc = basemap .scatter (x , y , c = ws [ds .good ], vmin = vmin , vmax = vmax ,
332- cmap = cmap , s = ms , marker = marker , edgecolors = ec )
332+ cmap = cmap , s = ms , marker = marker , edgecolors = ec ,
333+ ** kwargs )
333334 if colorbar_flag :
334335 plt .colorbar (sc , label = 'CYGNSS Wind Speed (m/s)' )
335336 if axis_label_flag :
@@ -534,8 +535,7 @@ class CygnssTrack(object):
534535 """
535536 Class to facilitate extraction of a single track of specular points
536537 from a CygnssSingleSat, CygnssMultiSat, or CygnssL2WindDisplay object.
537-
538-
538+
539539 Attributes
540540 ----------
541541 input = CygnssSubsection object
@@ -545,7 +545,7 @@ class CygnssTrack(object):
545545 lat = Latitudes of specular points
546546 rcg = Range-corrected gains of specular points
547547 datetimes = Datetime objects for specular points
548-
548+
549549 The following attributes are created by filter_track method:
550550 fws = Filtered wind speeds
551551 flon = Filtered longitudes
@@ -693,12 +693,12 @@ def get_tracks(data, indices=None, min_samples=10, verbose=False,
693693
694694 data = CygnssSingleSat, CygnssMultiSat, or CygnssL2WindDisplay object
695695 indices = Indices (2-element tuple) to use to limit the period of data
696- shown (i.e., limit by time). Not usually necessary unless
696+ shown (i.e., limit by time). Not usually necessary unless
697697 processing more than one day's worth of data.
698698 min_samples = Minimum allowable track size (number of specular points)
699699 verbose = Set to True for some text updates while running
700700 filter = Set to True to filter each track
701- window = Window length of filter, in number of specular points. Must be odd.
701+ window = Window length of filter, in # of specular points. Must be odd.
702702 """
703703 trl = []
704704 dts = get_datetime (data )
@@ -729,10 +729,12 @@ def get_tracks(data, indices=None, min_samples=10, verbose=False,
729729 dsc .tws = ds .tws [labels == element ]
730730 dsc .rcg = ds .lon [labels == element ]
731731 dsc .datetimes = ds .datetimes [labels == element ]
732+ dsc .sat = csat
733+ dsc .prn = gsat
732734 trl .append (dsc )
733735 if filter :
734736 for tr in trl :
735- tr .filter_track (window = window )
737+ tr .filter_track (window = window )
736738 return trl
737739
738740
0 commit comments