ClassAveragesVsReprojections

#!/usr/bin/env csh
if ( $# == 2 ) then

 set selfile=$1
 set volume=$2
 xmipp_header_extract -i $selfile -o classes.doc
 mkdir projections
 xmipp_angular_project_library -i $volume -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 &
else
 echo "Usage:  ClassAveragesVsReprojections selfile_with_class_averages reference_volume"
endif

-- Main.SjorsScheres - 04 Aug 2009