code: def _mosaic_get_full_block(sorted_mosaics):
      
      # For each slice / mosaic create a data volume block
      data_blocks = []
      for index in range(0, len(sorted_mosaics)):
          data_blocks.append(_mosaic_to_block(sorted_mosaics[index]))
  
      # Add the data_blocks together to one 4d block
      size_x = numpy.shape(data_blocks[0])[0]
      size_y = numpy.shape(data_blocks[0])[1]
      size_z = numpy.shape(data_blocks[0])[2]
      size_t = len(data_blocks)
      full_block = numpy.zeros((size_x, size_y, size_z, size_t), dtype=data_blocks[0].dtype)
      for index in range(0, size_t):
          full_block[:, :, :, index] = data_blocks[index]
  
      # Apply the rescaling if needed
      common.apply_scaling(full_block, sorted_mosaics[0])
  
      return full_block
docs: Generate a full datablock containing all timepoints