code: def get_coefs(
      c_align,
      c_exp,
      tpz_size,
      nb_tpz,
      v_track,
      scans,
      scan_size,
      scan_offset,
  ):
      """Compute the coeffs in numpy domain."""
      nties = nb_tpz.item()
      tpz_size = tpz_size.item()
      v_scan = (
          np.arange(nties * tpz_size) % tpz_size + scan_offset
      ) / tpz_size
      s_scan, s_track = np.meshgrid(v_scan, v_track)
      s_track = s_track.reshape(scans, scan_size, nties, tpz_size)
      s_scan = s_scan.reshape(scans, scan_size, nties, tpz_size)
  
      c_align = c_align[np.newaxis, np.newaxis, :, np.newaxis]
      c_exp = c_exp[np.newaxis, np.newaxis, :, np.newaxis]
  
      a_scan = (
          s_scan
          + s_scan * (1 - s_scan) * c_exp
          + s_track * (1 - s_track) * c_align
      )
      a_track = s_track
      coef_a = (1 - a_track) * (1 - a_scan)
      coef_b = (1 - a_track) * a_scan
      coef_d = a_track * (1 - a_scan)
      coef_c = a_track * a_scan
      res = np.stack([coef_a, coef_b, coef_c, coef_d], axis=4).reshape(
          scans * scan_size, -1, 4
      )
      return res
  
imports: import logging
  from contextlib import suppress
  from datetime import datetime, timedelta
  import dask.array
  import h5py
  import numpy
  import xarray
  from satpy.readers.file_handlers import BaseFileHandler
  from satpy.readers.utils import np2str
  from satpy.utils import angle2xyz, get_legacy_chunk_size, lonlat2xyz, xyz2angle, xyz2lonlat
label_desc: Use image node_id 3 for calling a global function with example usage: get_coefs(c_align, c_exp, tpz_size, nb_tpz, v_track, scans, scan_size, scan_offset) and returns: res