code: def _extract_two_byte_mask(
      data_a: np.ndarray,
      data_b: np.ndarray,
      bit_start: int,
      bit_count: int,
  ) -> np.ndarray:
      data_a = data_a.astype(np.uint16, copy=False)
      data_a = np.left_shift(data_a, 8)  # dataset_a << 8
      byte_dataset = np.bitwise_or(data_a, data_b).astype(np.uint16)
      shape = byte_dataset.shape
      # We replicate the concatenated byte with the right shape
      byte_dataset = np.repeat(
          np.repeat(byte_dataset, 4, axis=0), 4, axis=1
      )
      # All bits carry information, we update bit_start consequently
      bit_start = np.arange(16, dtype=np.uint16).reshape((4, 4))
      bit_start = np.tile(bit_start, (shape[0], shape[1]))
      return _bits_strip(bit_start, bit_count, byte_dataset)
  
imports: import logging
  import dask.array
  import numpy
  import xarray
  from satpy.readers.hdf4_utils import from_sds
  from satpy.readers.hdfeos_base import HDFEOSGeoReader
  from satpy.utils import get_legacy_chunk_size
label_desc: Use image node_id 2 for calling a global function with example usage: _extract_two_byte_mask(data_a, data_b, bit_start, bit_count) and returns: _bits_strip