ctf_estimate_psd_with_arma.cpp File Reference

#include "ctf_estimate_psd_with_arma.h"
#include <core/args.h>
#include <core/xmipp_fftw.h>
#include <data/filters.h>

Include dependency graph for ctf_estimate_psd_with_arma.cpp:

Go to the source code of this file.

Functions
void	First_Quadrant_Neighbors (int N, int M, MultidimArray< double > &Neighbors)
void	Second_Quadrant_Neighbors (int N, int M, MultidimArray< double > &Neighbors)
void	CausalARMA (MultidimArray< double > &Img, ARMA_parameters &prm)
void	ARMAFilter (MultidimArray< double > &Img, MultidimArray< double > &Filter, ARMA_parameters &prm)

Function Documentation

First_Quadrant_Neighbors()

void First_Quadrant_Neighbors	(	int	N,
		int	M,
		MultidimArray< double > &	Neighbors
	)

Definition at line 52 of file ctf_estimate_psd_with_arma.cpp.

{
    long NumberOfPoints = (N + 1) * M;
    int n;
    Neighbors.resize(NumberOfPoints, 3);

    n = 0; // Number of neighbors found so far.

    for (int p = N; p >= 0; p--)
    {
        for (int q = M; q > 0; q--)
        {
            A2D_ELEM(Neighbors,n, 0) = (double)p;
            A2D_ELEM(Neighbors,n, 1) = (double)q;
            n++;
        }
    }
}

Second_Quadrant_Neighbors()

void Second_Quadrant_Neighbors	(	int	N,
		int	M,
		MultidimArray< double > &	Neighbors
	)

Definition at line 72 of file ctf_estimate_psd_with_arma.cpp.

{
    long NumberOfPoints = N * (M + 1);
    int n;
    Neighbors.resize(NumberOfPoints, 3);

    n = 0; // Number of neighbors found so far.

    for (int p = N; p >= 1; p--)
    {
        for (int q = 0; q >= -M; q--)
        {
            A2D_ELEM(Neighbors,n, 0) = (double)p;
            A2D_ELEM(Neighbors,n, 1) = (double)q;
            n++;
        }
    }
}