docs: def orientnii(imfile):
      '''Get the orientation from NIfTI sform.  Not fully functional yet.'''
      strorient = ['L-R', 'S-I', 'A-P']
      niiorient = []
      niixyz = np.zeros(3,dtype=np.int8)
      if os.path.isfile(imfile):
          nim = nib.load(imfile)
          pct = nim.get_data()
          A = nim.get_sform()
          for i in range(3):
              niixyz[i] = np.argmax(abs(A[i,:-1]))
              niiorient.append( strorient[ niixyz[i] ] )
          print niiorient