code: def _get_distances_to_intersections(view_vector, sat_pos, ellipsoid):
      """Get distances to intersections with the earth's surface.
  
      Returns:
          Distances to two intersections with the surface.
      """
      a, b, c = _get_abc_helper(view_vector, sat_pos, ellipsoid)
      tmp = np.sqrt((b**2 - a * c))
      dist_1 = (-b + tmp) / a
      dist_2 = (-b - tmp) / a
      return dist_1, dist_2
  
imports: from collections import namedtuple
  import dask.array
  import numba
  import numpy
  from satpy.utils import get_legacy_chunk_size
label_desc: Use image node_id 25 for calling a global function with example usage: _get_distances_to_intersections(view_vector, sat_pos, ellipsoid) and returns: dist_1,  dist_2