subtom_parallel_eigenvolumes
Computes projections of data onto Eigenvectors.
subtom_parallel_eigenvolumes(
'all_motl_fn_prefix', all_motl_fn_prefix ('combinedmotl/allmotl'),
'ptcl_fn_prefix', ptcl_fn_prefix ('subtomograms/subtomo'),
'eig_vec_fn_prefix', eig_vec_fn_prefix ('class/eigvec_pca'),
'eig_val_fn_prefix', eig_val_fn_prefix ('class/eigval_pca'),
'xmatrix_fn_prefix', xmatrix_fn_prefix ('class/xmatrix_pca'),
'eig_vol_fn_prefix', eig_vol_fn_prefix ('class/eigvol_pca'),
'mask_fn', mask_fn ('none'),
'iteration', iteration (1),
'num_xmatrix_batch', num_xmatrix_batch (1),
'process_idx', process_idx (1))
Calculates the summed projections of particles onto previously determined Eigen
(or left singular) vectors, by means of an also previously calculated X-matrix
to produce Eigenvolumes which can then be used to determine which vectors can
best influence classification. The Eigenvectors are named based on
eig_vec_fn_prefix
and iteration
and the X-Matrix is named based on
xmatrix_fn_prefix
, iteration
, and process_idx
. The Eigenvolumes are
also masked by the file specified by mask_fn
. The Eigenvolumes are split
into num_xmatrix_batch
sums, which is the same number of batches that the
X-Matrix was broken into in its computation. process_idx
is a counter that
designates the current batch being determined. The output sum Eigenvolume will
be written out as specified by eig_vol_fn_prefix
, iteration
,
process_idx
and # where the # is the particular Eigenvolume being written
out.
Example
subtom_parallel_eigenvolumes(
'all_motl_fn_prefix', 'combinedmotl/allmotl', ...
'ptcl_fn_prefix', 'subtomograms/subtomo', ...
'eig_vec_fn_prefix', 'class/eigvec', ...
'eig_val_fn_prefix', 'class/eigval', ...
'xmatrix_fn_prefix', 'class/xmatrix', ...
'eig_vol_fn_prefix', 'class/eigvol', ...
'mask_fn', 'otherinputs/classification_mask.em', ...
'iteration', 1, ...
'num_xmatrix_batch', 100, ...
'process_idx', 1)