code: def sen_seasonal_slopes(x):
      x = ma.array(x, subok=True, copy=False, ndmin=2)
      (n, _) = x.shape
      # Get list of slopes per season
      szn_slopes = ma.vstack(
          [
              (x[i + 1 :] - x[i]) / np.arange(1, n - i)[:, None]
              for i in range(n)
          ]
      )
      szn_medslopes = ma.median(szn_slopes, axis=0)
      medslope = ma.median(szn_slopes, axis=None)
      return szn_medslopes, medslope
  
imports: from __future__ import division, print_function, absolute_import
  import numpy
  from numpy import ndarray
  import numpy.ma
  from numpy.ma import masked, nomask
  from scipy._lib.six import iteritems
  import itertools
  import warnings
  from collections import namedtuple
  from .None import distributions
  import scipy.special
  from ._stats_mstats_common import _find_repeats, linregress, theilslopes, siegelslopes
label_desc: Use image node_id 19 for calling a global function with example usage: sen_seasonal_slopes(x) and returns: szn_medslopes,  medslope