code: def read_taiwan_ntu_dsi():
      """Load Taiwan NTU dataset.
  
      Returns
      -------
      img : obj,
          Nifti1Image
      gtab : obj,
          GradientTable
  
      """
      fraw, fbval, fbvec = get_fnames("taiwan_ntu_dsi")
      bvals, bvecs = read_bvals_bvecs(fbval, fbvec)
      bvecs[1:] = (
          bvecs[1:]
          / np.sqrt(np.sum(bvecs[1:] * bvecs[1:], axis=1))[:, None]
      )
      gtab = gradient_table(bvals, bvecs)
      img = nib.load(fraw)
      return img, gtab
  
imports: import os
  import contextlib
  import logging
  import json
  import os.path
  from os.path import join
  from hashlib import md5
  from shutil import copyfileobj
  from tqdm.auto import tqdm
  import numpy
  import nibabel
  import tarfile
  import tempfile
  import zipfile
  from dipy.core.gradients import gradient_table, gradient_table_from_gradient_strength_bvecs
  from dipy.io.gradients import read_bvals_bvecs
  from dipy.io.image import load_nifti, load_nifti_data, save_nifti
  from dipy.io.streamline import load_trk
  from dipy.utils.optpkg import optional_package, TripWire
  from urllib.request import urlopen
label_desc: Use image node_id 19 for calling a global function with example usage: read_taiwan_ntu_dsi() and returns: img,  gtab