docs: def is_orthogonal(dicoms, log_details=False):
      """
      Validate that volume is orthonormal
  
      :param dicoms: check that we have a volume without skewing
      """
      first_image_orient1 = numpy.array(dicoms[0].ImageOrientationPatient)[0:3]
      first_image_orient2 = numpy.array(dicoms[0].ImageOrientationPatient)[3:6]
      first_image_pos = numpy.array(dicoms[0].ImagePositionPatient)
  
      last_image_pos = numpy.array(dicoms[-1].ImagePositionPatient)
  
      first_image_dir = numpy.cross(first_image_orient1, first_image_orient2)
      first_image_dir /= numpy.linalg.norm(first_image_dir)
  
      combined_dir = last_image_pos - first_image_pos
      combined_dir /= numpy.linalg.norm(combined_dir)
  
      if not numpy.allclose(first_image_dir, combined_dir, rtol=0.05, atol=0.05) \
              and not numpy.allclose(first_image_dir, -combined_dir, rtol=0.05, atol=0.05):
          if log_details:
              logger.warning('Orthogonality check failed: non cubical image')
              logger.warning('---------------------------------------------------------')
              logger.warning(first_image_dir)
              logger.warning(combined_dir)
              logger.warning('---------------------------------------------------------')
          return False
      return True