code: def dct_ii_cut_basis(volume_times, cut_period):
      """DCT-II regressors with periods >= `cut_period`
  
      See: http://en.wikipedia.org/wiki/Discrete_cosine_transform#DCT-II
  
      Parameters
      ----------
      volume_times : array-like
          Times of acquisition of each volume.  Must be regular and continuous
          otherwise we raise an error.
      cut_period: float
          Cut period (wavelength) of the low-pass filter (in time units).
  
      Returns
      -------
      cdrift: array shape (n_scans, n_drifts)
          DCT-II drifts plus a constant regressor in the final column.  Constant
          regressor always present, regardless of `cut_period`.
      """
      N = len(volume_times)
      hfcut = 1.0 / cut_period
      dt = volume_times[1] - volume_times[0]
      # Such that hfcut = 1/(2*dt) yields N
      order = int(np.floor(2 * N * hfcut * dt))
      # Always return constant column
      if order == 0:
          return np.ones((N, 1))
      basis = np.ones((N, order))
      basis[:, :-1] = dct_ii_basis(volume_times, order, normcols=True)[
          :, 1:
      ]
      return basis
  
label_desc: Use image node_id 2 for calling a global function with example usage: dct_ii_cut_basis(volume_times, cut_period) and returns: basis, np