code: def expand(data, coefs, scans, scan_size):
      """Perform the expansion in numpy domain."""
      data = data.reshape(data.shape[:-1])
  
      coefs = coefs.reshape(scans, scan_size, data.shape[1] - 1, -1, 4)
  
      coef_a = coefs[:, :, :, :, 0]
      coef_b = coefs[:, :, :, :, 1]
      coef_c = coefs[:, :, :, :, 2]
      coef_d = coefs[:, :, :, :, 3]
  
      corner_coefficients = (coef_a, coef_b, coef_c, coef_d)
      fdata = _interpolate_data(data, corner_coefficients, scans)
      return fdata.reshape(scans * scan_size, -1)
  
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 4 for calling a global function with example usage: expand(data, coefs, scans, scan_size) and returns: fdata