radon.h File Reference

#include <core/xmipp_image.h>

Include dependency graph for radon.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
void	Radon_Transform (const MultidimArray< double > &vol, double rot, double psi, MultidimArray< double > &RT)
void	Local_Radon_Transform (const MultidimArray< double > &vol, double rot, double tilt, int label, const MultidimArray< double > &vol_label, MultidimArray< double > &RT, MultidimArray< double > &RT_n)
void	Radon_Transform (const MultidimArray< double > &I, double rot_step, MultidimArray< double > &RT)

Function Documentation

Local_Radon_Transform()

void Local_Radon_Transform	(	const MultidimArray< double > &	vol,
		double	rot,
		double	tilt,
		int	label,
		const MultidimArray< double > &	vol_label,
		MultidimArray< double > &	RT,
		MultidimArray< double > &	RT_n
	)

Radon transform of a piece of volume along a direction

Definition at line 51 of file radon.cpp.

{
    MultidimArray<double>   Rvol;
    MultidimArray<double>   Lvol;

    // Align the Radon direction with the Z axis
    if (rot != 0 || tilt != 0)
    {
        Euler_rotate(vol, rot, tilt, 0.0F, Rvol);
        Euler_rotate(vol_label, rot, tilt, 0.0F, Lvol);
    }
    else
    {
        Rvol = vol;
        Lvol = vol_label;
    }

    // Project onto one line
    RT.initZeros(Rvol.zdim);
    STARTINGX(RT) = STARTINGZ(Rvol);
    RT_n = RT;

    for (int k = STARTINGZ(Rvol); k <= FINISHINGZ(Rvol); k++)
        for (int i = STARTINGY(Rvol); i <= FINISHINGY(Rvol); i++)
            for (int j = STARTINGX(Rvol); j <= FINISHINGX(Rvol); j++)
                if (A3D_ELEM(Lvol, k, i, j) == label)
                {
                    A1D_ELEM(RT, k) += A3D_ELEM(Rvol, k, i, j);
                    A1D_ELEM(RT_n, k)++;
                }
}

Radon_Transform() [1/2]

void Radon_Transform	(	const MultidimArray< double > &	vol,
		double	rot,
		double	psi,
		MultidimArray< double > &	RT
	)

Radon transform of a volume along a direction

Definition at line 30 of file radon.cpp.

{
    MultidimArray<double> Rvol;

    // Align the Radon direction with the Z axis
    if (rot != 0 || tilt != 0) Euler_rotate(vol, rot, tilt, 0.0F, Rvol);
    else Rvol = vol;

    // Project onto one line
    RT.initZeros(Rvol.zdim);
    STARTINGX(RT) = STARTINGZ(Rvol);

    for (int k = STARTINGZ(Rvol); k <= FINISHINGZ(Rvol); k++)
        for (int i = STARTINGY(Rvol); i <= FINISHINGY(Rvol); i++)
            for (int j = STARTINGX(Rvol); j <= FINISHINGX(Rvol); j++)
                A1D_ELEM(RT, k) += A3D_ELEM(Rvol, k, i, j);
}

Radon_Transform() [2/2]

void Radon_Transform	(	const MultidimArray< double > &	I,
		double	rot_step,
		MultidimArray< double > &	RT
	)

Radon transform of an image.

Definition at line 87 of file radon.cpp.

{
    MultidimArray<double> rot_I;
    RT.initZeros(CEIL(360.0 / rot_step), XSIZE(I));
    STARTINGX(RT) = STARTINGX(I);
    int iterRot = static_cast<int>(360.0 / rot_step);
    for (int l=0; l<iterRot; l++)
    {
        double rot = l*rot_step;
        // Rotate image
        rotate(xmipp_transformation::LINEAR, rot_I, I, rot);
        // Sum by columns
        FOR_ALL_ELEMENTS_IN_ARRAY2D(rot_I) RT(l, j) += rot_I(i, j);
    }
}