code: def _make_bti_dig_points(
      nasion,
      lpa,
      rpa,
      hpi,
      extra,
      convert=False,
      use_hpi=False,
      bti_dev_t=False,
      dev_ctf_t=False,
  ):
      from ..transforms import (
          Transform,
          combine_transforms,
          invert_transform,
      )
  
      _hpi = hpi if use_hpi else None
      info_dig, ctf_head_t = _call_make_dig_points(
          nasion, lpa, rpa, _hpi, extra, convert
      )
  
      if convert:
          t = combine_transforms(
              invert_transform(bti_dev_t), dev_ctf_t, "meg", "ctf_head"
          )
          dev_head_t = combine_transforms(t, ctf_head_t, "meg", "head")
      else:
          dev_head_t = Transform("meg", "head", trans=None)
  
      return info_dig, dev_head_t, ctf_head_t
  
imports: import heapq
  from collections import Counter
  import numpy
  from ..utils import Bunch, _check_fname, _validate_type, logger, verbose, warn
  from .constants import FIFF, _coord_frame_named
  from .tag import read_tag
  from .tree import dir_tree_find
  from .write import start_and_end_file, write_dig_points
label_desc: Use image node_id 13 for calling a global function with example usage: _make_bti_dig_points(nasion, lpa, rpa, hpi, extra, convert, use_hpi, bti_dev_t, dev_ctf_t) and returns: info_dig,  dev_head_t,  ctf_head_t