code: def _create_affine(x_axis, y_axis, z_axis, image_pos, voxel_sizes):
      
  
      # Create affine matrix (http://nipy.sourceforge.net/nibabel/dicom/dicom_orientation.html#dicom-slice-affine)
  
      affine = numpy.array(
          [[x_axis[0] * voxel_sizes[0], y_axis[0] * voxel_sizes[1], z_axis[0] * voxel_sizes[2], image_pos[0]],
           [x_axis[1] * voxel_sizes[0], y_axis[1] * voxel_sizes[1], z_axis[1] * voxel_sizes[2], image_pos[1]],
           [x_axis[2] * voxel_sizes[0], y_axis[2] * voxel_sizes[1], z_axis[2] * voxel_sizes[2], image_pos[2]],
           [0, 0, 0, 1]])
      return affine
docs: Function to generate the affine matrix for a dicom series
      This method was based on (http://nipy.org/nibabel/dicom/dicom_orientation.html)
  
      :param sorted_dicoms: list with sorted dicom files