Rotational_Spectrum Class Reference

#include <rotational_spectrum.h>

Collaboration diagram for Rotational_Spectrum:

Public Member Functions
void	compute_rotational_spectrum (MultidimArray< double > &img, double xr1, double xr2, double xdr, double xr)
void	compute_rotational_spectrum (Cylindrical_Wave_Decomposition &cwd, double xr1, double xr2, double xdr, double xr)
	Compute rotational spectrum using CWD. More...

Public Attributes
int	ir
int	numin
	Minimum harmonics. More...
int	numax
	Maximum harmonics. More...
double	x0
	Center of symmetry (x). More...
double	y0
	Center of symmetry (y). More...
double	rl
	Minimum integration radius. More...
double	rh
	Maximum integration radius. More...
double	dr
	Integration increment. More...
MultidimArray< double >	rot_spectrum
	Rotational spectrum. More...

Friends
std::ostream &	operator<< (std::ostream &_out, const Rotational_Spectrum &_spt)
	Show. More...

Detailed Description

Rotational spectrum.

Example of use:

int main(int argc, char** argv)
{
    ImageXmipp I(argv[1]);

    int rl = 0;
    int rh = 22;
    int dr = 1;

    Rotational_Spectrum spt;
    spt.rl = rl;
    spt.rh = rh;
    spt.dr = dr;
    spt.numin = 1;
    spt.numax = 15;
    spt.x0 = (double) XSIZE(I()) / 2;
    spt.y0 = (double) YSIZE(I()) / 2;

    spt.compute_rotational_spectrum (I(), rl, rh, dr, rh - rl);

    std::cout << spt.rot_spectrum << std::endl;

    return 0;
}

Definition at line 114 of file rotational_spectrum.h.

Member Function Documentation

compute_rotational_spectrum() [1/2]

void Rotational_Spectrum::compute_rotational_spectrum	(	MultidimArray< double > &	img,
		double	xr1,
		double	xr2,
		double	xdr,
		double	xr
	)

Compute rotational spectrum of an image.

xr1 is the minimum integration radius. xr2 the maximum integration radius. xdr the increment, and xr the length of integration. Usually, xr=xr2-xr1.

Definition at line 403 of file rotational_spectrum.cpp.

{
    // Compute the cylindrical wave decomposition
    Cylindrical_Wave_Decomposition cwd;
    cwd.numin = numin;
    cwd.numax = numax;
    cwd.x0    = (x0 == -1) ? (double)XSIZE(img) / 2 : x0;
    cwd.y0    = (y0 == -1) ? (double)YSIZE(img) / 2 : y0;
    cwd.r1    = rl;
    cwd.r2    = rh;
    cwd.r3    = dr;
    cwd.compute_cwd(img);
    ir = cwd.ir;

    // Compute the rotational spectrum
    compute_rotational_spectrum(cwd, xr1, xr2, xdr, xr);
}

compute_rotational_spectrum() [2/2]

void Rotational_Spectrum::compute_rotational_spectrum	(	Cylindrical_Wave_Decomposition &	cwd,
		double	xr1,
		double	xr2,
		double	xdr,
		double	xr
	)

Compute rotational spectrum using CWD.

Definition at line 267 of file rotational_spectrum.cpp.

{
    double *e[MAX_HARMONIC];
    double *ep[MAX_HARMONIC];
    double *erp [MAX_HARMONIC];

    // Read the information from the Cylindrical Wave Decomposition .........
    auto c = (double *) calloc(5191, sizeof(double));
    auto s = (double *) calloc(5191, sizeof(double));
    if ((NULL == c) || (NULL == s))
        REPORT_ERROR(ERR_MEM_NOTENOUGH, "compute_rotational_spectrum::no memory");

    int numin  = cwd.numin;
    int numax  = cwd.numax;
    double rl  = cwd.r1;
    double rh  = cwd.r2;
    double dr  = cwd.r3;
    FOR_ALL_ELEMENTS_IN_ARRAY1D(cwd.out_ampcos)
    {
        c[i+1] = A1D_ELEM(cwd.out_ampcos, i);
        s[i+1] = A1D_ELEM(cwd.out_ampsin, i);
    }

    int n = numax - numin + 1;
    auto m = (int)((rh - rl) / dr + 1);
    for (int i = 1; i <= n; i++) {
        e[i] = (double *) calloc(m + 1, sizeof(double));
        if (NULL == e[i])
            REPORT_ERROR(ERR_MEM_NOTENOUGH, "compute_rotational_spectrum::no memory");
    }
    auto rv = (double *) calloc(m + 1, sizeof(double));
    auto st = (double *) calloc(m + 1, sizeof(double));
    if ((NULL == rv) || (NULL == st))
        REPORT_ERROR(ERR_MEM_NOTENOUGH, "compute_rotational_spectrum::no memory");

    // Computations .........................................................
    int j1 = 0;
    if (numin == 0)
        j1 = 2;
    else
        j1 = 1;
    int k = 0;
    for (int i = 1; i <= n; i++)
        for (int j = 1; j <= m; j++)
        {
            k++;
            e[i][j] = c[k] * c[k] + s[k] * s[k];
        }

    for (int i = 1; i <= m; i++)
    {
        rv[i] = rl + dr * (i - 1);
        st[i] = 0;
        for (int j = j1; j <= n; j++)
            st[i] += e[j][i];
    }

    auto ir1 = (int)((xr1 - rl) / dr + 1);
    if (ir1 < 1)
        ir1 = 1;
    auto ir2 = (int)((XMIPP_MIN(xr2, rh) - rl) / dr + 1);
    if (ir2 < ir1)
        ir2 = ir1;
    auto ndr = (int)(xdr / dr);
    if (ndr < 1)
        ndr = 1;
    int nr = 0;
    if (xr < 0)
        nr = m;
    else
        nr = (int)(xr / dr + 1);
    ir2 = XMIPP_MIN(ir2, m);
    int ncol = ir2 - nr + 1 - ir1;
    if (ncol < 0)
        ncol = 0;
    else
        ncol = 1 + ncol / ndr;
    if (ncol == 0)
        REPORT_ERROR(ERR_VALUE_INCORRECT, "compute_rotational_spectrum::Incorrect data");

    int nvez = (ncol - 1) / 13 + 1;
    for (int i = 1; i <= n; i++) {
        ep[i] = (double *) calloc(ncol + 1, sizeof(double));
        erp[i] = (double *) calloc(ncol + 1, sizeof(double));
        if ((NULL == ep[i]) || (NULL == erp[i]))
            REPORT_ERROR(ERR_MEM_NOTENOUGH, "compute_rotational_spectrum::no memory");
    }
    auto rp1 = (double *) calloc(ncol + 1, sizeof(double));
    auto rp2 = (double *) calloc(ncol + 1, sizeof(double));
    auto sp  = (double *) calloc(ncol + 1, sizeof(double));
    if ((NULL == rp1) || (NULL == rp2) || (NULL == sp))
        REPORT_ERROR(ERR_MEM_NOTENOUGH, "compute_rotational_spectrum::no memory");
    for (k = 1; k <= ncol; k++)
    {
        int irk = ir1 + (k - 1) * ndr - 1;
        rp1[k] = 10 * rv[irk+1];
        rp2[k] = 10 * rv[irk+nr];
        sp[k] = 0;
        for (int k1 = 1; k1 <= nr; k1++)
            sp[k] += st[irk+k1];
        for (int i = 1; i <= n; i++)
        {
            ep[i][k] = 0;
            for (int k1 = 1; k1 <= nr; k1 ++)
                ep [i][k] += e[i][irk+k1];
            erp[i][k] = 1000000. * ep[i][k] / sp[k];
        }
    }

    // Keep results .........................................................
    rot_spectrum.initZeros(n - j1 + 1);
    for (k = 1; k <= nvez; k++)
    {
        int k1 = 13 * (k - 1) + 1;
        for (int i = j1; i <= n; i++)
            rot_spectrum(i - j1) = erp[i][k1] / 10000;
    }

    // Free memory
    for (int i = 1; i <= n; i++)
    {
        free(ep[i]);
        free(erp [i]);
    }
    free(c);
    free(s);
    free(rv);
    free(st);
    free(rp1);
    free(rp2);
    free(sp);
}

Friends And Related Function Documentation

operator<<

std::ostream& operator<< ( std::ostream &  _out, const Rotational_Spectrum &  _spt  )

friend

Show.

Definition at line 253 of file rotational_spectrum.cpp.

{
    _out << "numin=" << _spt.numin << std::endl
         << "numax=" << _spt.numax << std::endl
         << "x0=" << _spt.x0 << std::endl
         << "y0=" << _spt.y0 << std::endl
         << "rl=" << _spt.rl << std::endl
         << "rh=" << _spt.rh << std::endl
         << "dr=" << _spt.dr << std::endl;
    return _out;
}

Member Data Documentation

dr

double Rotational_Spectrum::dr

Integration increment.

Definition at line 144 of file rotational_spectrum.h.

ir

int Rotational_Spectrum::ir

Ir.

This value is given by the Cylindrical Wave Decomposition. It is not necessary to fill this field if the rotational spectrum is computed on an image.

Definition at line 123 of file rotational_spectrum.h.

numax

int Rotational_Spectrum::numax

Maximum harmonics.

Definition at line 129 of file rotational_spectrum.h.

numin

int Rotational_Spectrum::numin

Minimum harmonics.

Definition at line 126 of file rotational_spectrum.h.

rh

double Rotational_Spectrum::rh

Maximum integration radius.

Definition at line 141 of file rotational_spectrum.h.

rl

double Rotational_Spectrum::rl

Minimum integration radius.

Definition at line 138 of file rotational_spectrum.h.

rot_spectrum

MultidimArray< double > Rotational_Spectrum::rot_spectrum

Rotational spectrum.

Definition at line 147 of file rotational_spectrum.h.

x0

double Rotational_Spectrum::x0

Center of symmetry (x).

Definition at line 132 of file rotational_spectrum.h.

y0

double Rotational_Spectrum::y0

Center of symmetry (y).

Definition at line 135 of file rotational_spectrum.h.

---

The documentation for this class was generated from the following files:

• xmipp/libraries/data/rotational_spectrum.h
• xmipp/libraries/data/rotational_spectrum.cpp