polyTEM.image.fft.frc

polyTEM.image.fft.frc(fft1, fft2, dx, processor)

Fourier Ring Coefficient based on Baxter, et al. https://doi.org/10.1016/j.jsb.2009.02.012

FRC(k) =

rac{Sigma_k fft1(k) * conj(fft2(k)}

{np.sqrt(Sigma_k magnitude(fft2(k))**2) * np.sqrt(Sigma_k magnitude(fft1(k))**2)} where k is the spectral frequency k is equivalent to u in other modules, = 1/d

BUG: ASSUMES SQUARE FFT of same shapes

Arguments:

fft1: [real, imag] components of fourier transform of image 1 fft2: [real, imag] components of fourier transform of image 2

Returns:

FRC_u: np.ndarray of reciprocal space 1/d units FRC_counts: the FRC coefficient FRC(q)