(phase flipping)

Collaboration diagram for (phase flipping):

Classes
class	ProgCTFShow
class	ProgCTFPhaseFlipping

Functions
void	actualPhaseFlip (MultidimArray< double > &I, CTFDescription ctf)

Detailed Description

Function Documentation

actualPhaseFlip()

void actualPhaseFlip	(	MultidimArray< double > &	I,
		CTFDescription	ctf
	)

Definition at line 88 of file ctf_phase_flip.cpp.

{
    // Perform the Fourier transform
    FourierTransformer transformer;
    MultidimArray< std::complex<double> > M_inFourier;
    transformer.FourierTransform(I,M_inFourier,false);

    Matrix1D<double> freq(2); // Frequencies for Fourier plane
    int yDim=YSIZE(I);
    int xDim=XSIZE(I);
    double iTm=1.0/ctf.Tm;
    ctf.phase_shift = (ctf.phase_shift*PI)/180;

    for (size_t i=0; i<YSIZE(M_inFourier); ++i)
    {
        FFT_IDX2DIGFREQ(i, yDim, YY(freq));
        YY(freq) *= iTm;
        for (size_t j=0; j<XSIZE(M_inFourier); ++j)
        {
            FFT_IDX2DIGFREQ(j, xDim, XX(freq));
            XX(freq) *= iTm;
            ctf.precomputeValues(XX(freq),YY(freq));
            if (ctf.getValuePureWithoutDampingAt()<0)
                DIRECT_A2D_ELEM(M_inFourier,i,j)*=-1;
        }
    }

    // Perform inverse Fourier transform and finish
    transformer.inverseFourierTransform();
}