xmipp_resolution_fsc (v3.0)


Calculate the resolution of one or more volumes or images with respect to a single reference.

Three methods are employed:

* Differential Phase Residual (DPR)
In this method the the resolution is defined as the spatial frequency at which the average phase discrepancy between the two transforms exceeds 45 degrees.
* Fourier Ring Correlation (FRC)
In this method the resolution is defined as the spatial frequency at which annular samplings of the two Fourier transforms register negligible cross-correlation.
* Spectral Signal-to-Noise Ratio (SSNR)
This method is based on the measurement of the signal-to-noise ratio as a function of spatial frequency. The SSNR is determined by comparing the Fourier transform of individual images with that of the global average image. (this option is only available for 2D images not for volumes, see the program ssnr if you are interested in computing the signal to noise ratio in 3D).

To calculate "the resolution of a reconstruction", use --set_of_images, where the program divides the corresponding set of projection images into two subsets and reconstructs two volumes from these subsets. This program may then be used to calculate the DPR and FRC between these two volumes. The resulting plots are commonly used to assess the high-resolution limit of the initial reconstruction (see Frank's book for more details).

The program writes out filename.frc files, for each input volume or image, or selfilename.frc, the set_of_images mode. These ACSII files contain the DPR, FRC and SSNR as a function of resolution (in 1/Angstrom). The .frc files also contain a column for the FRC expected for pure noise.

See also


-i <input_file> ( requires --ref )
either an image/volume or a selection file
or --set_of_images <selfile>
selfile containing a set of 2D-images
--oroot <root_file=>
Root of the output metadata. If not set, input file rootname is taken.
-o <output_file=>
Output file name.
--ref <input_file>
filename for reference image/volume
--sampling_rate, -s <Ts=1>
Pixel size (Angstrom)
for 2D-images: do not apply transformation stored in the header
compute dpr, by default only frc is computed
--max_sam <max_sr=-1>
set fsc to 0 for frequencies above this one (Angstrom), -1 -> all fequencies

Examples and notes

Resolution of subset2.vol volume with respect to subset1.vol reference volume using 5.6 pixel size (in Angstrom):
xmipp_resolution_fsc --ref subset1.vol  -i subset2.vol --sam 5.6 
Resolution of a set of images using 5.6 pixel size (in Angstrom):
xmipp_resolution_fsc --set_of_images selfile.sel --sam 5.6

User's comments