Name

Fourier Ring/Shell Correlation with mask

Usage

frsc = fsc(img1, img2 ,[mask, w, filename])

Input

img1
first image
img2
second image
mask
mask file (does not have to be binary).
w
ring/shell width in Fourier space. Default: w=1.
filename
result of this function is stored in the filename provided.
Note
mask, w and filename are optional.

Output

frsc
a list with three columns:

Note:: This is the output when the result is not written to the file but just displayed on the screen when the script is executed.

filename
the name of the output text file where the results are stored in three columns:
  1. the absolute Fourier frequency
  2. FSC
  3. number of Fourier coefficients n within given ring/shell.

Note:: This is the output when the result is written to a file and when the contents of the file are displayed.

Description

The calculation is done in Fourier space as

where `hat(f)` and `hat(g)` are Fourier transforms of two input images and summation is performed over rings/shells in Fourier space at spatial frequency `|bbb{y}_n|` and within ring/shell width w.

Method

For each input image, the average pixel value is calculated at location where the mask value is 0.5 within one pixel wide band, this average is subtracted, both FFTs are calculated (without padding) and FSC is calculated for each ring. The format of input images is not changed.

Reference

W. O. Saxton and W. Baumeister, ‘‘The correlation averaging of a regularly arranged bacterial envelope protein,’’ J. Microsc. (Oxford) 127, 127–138 (1982).

P. A. Penczek, ‘‘Three-dimensional spectral signal-to-noise ratio for a class of reconstruction algorithms,’’ J. Struct. Biol. 138, 34–46 (2002).

Author / Maintainer

Pawel Penczek

Keywords

category 1
STATISTICS.
category 2
FOURIER.

Files

statistics.py

Maturity

stable
works.

Bugs

None. It is perfect.

fsc_mask (last edited 2013-07-01 13:13:02 by localhost)