CompareClassAveragesWithReprojections

ML2D with reprojections of your model (reference.vol) in Xmipp:

xmipp_header_extract -i classes.sel -o classes.doc

mkdir projections

xmipp_angular_project_library -i reference.vol -o projections/ref -sym c1h -sampling_rate 5 -experimental_images classes.doc -compute_neighbors -angular_distance -1

xmipp_angular_projection_matching -i classes.doc -ref projections/ref -o compare

xmipp_header_assign -i compare.doc -mirror

cat compare.doc | awk '{ if ($2==8) {print p0" 1"; printf "%s%06.0f%s\n","projections/ref",$8,".xmp 1"; } p0=$2}' >compare_classes_vs_projections.sel

xmipp_show -sel compare_classes_vs_projections.sel -w 2 &

The left-hand side of the show window will be the reference-free class averages, the right-hand side will be the most similar reprojection of your 3D model.

This sequence of commands was also implemented as a shell script: ClassAveragesVsReprojections.csh. Save this script in a file calledClassAveragesVsReprojections.csh, and execute it as follows:

chmod +x ClassAveragesVsReprojections.csh
./ClassAveragesVsReprojections.csh classes.sel reference.vol

-- Main.SjorsScheres - 04 Aug 2009