Basis function

Collaboration diagram for Basis function:

Modules
	Blobs
	Splines

Classes
class	Basis

Functions
void	createZernike3DBasis (const MultidimArray< double > &Vin, MultidimArray< double > &Vbasis, int l1, int n, int l2, int m, int Rmax=-1)

Detailed Description

This class defines the basis function to use for the reconstruction. Currently, valid basis functions are blobs and voxels.

Function Documentation

createZernike3DBasis()

void createZernike3DBasis	(	const MultidimArray< double > &	Vin,
		MultidimArray< double > &	Vbasis,
		int	l1,
		int	n,
		int	l2,
		int	m,
		int	Rmax = -1
	)

Create a Zernike Spherical Basis volume. With the same size as an input volume. Z(l,n,m)=R_n^l(r)Y_l^m(x,y,z) for l=0,1,2,..., n=0,1,2,..., m=-l,...,l The basis volume is created up to Rmax, by default half the size of the input volume

Definition at line 406 of file basis.cpp.

{
    Vbasis.initZeros(Vin);
    Vbasis.setXmippOrigin();
    if (Rmax<0)
        Rmax=XSIZE(Vin)/2;
    double Rmax2=Rmax*Rmax;
    double iRmax=1.0/Rmax;
    for (int k=STARTINGZ(Vbasis); k<=FINISHINGZ(Vbasis); k++)
    {
        double k2=k*k;
        if (k2>Rmax2)
            continue;
        double kr=k*iRmax;
        for (int i=STARTINGY(Vbasis); i<=FINISHINGY(Vbasis); i++)
        {
            double k2i2=k2+i*i;
            if (k2i2>Rmax2)
                continue;
            double ir=i*iRmax;
            for (int j=STARTINGX(Vbasis); j<=FINISHINGX(Vbasis); j++)
            {
                double r2=k2i2+j*j;
                if (r2>Rmax2)
                    continue;
                double jr=j*iRmax;
                A3D_ELEM(Vbasis,k,i,j)=ZernikeSphericalHarmonics(l1,n,l2,m,jr,ir,kr,sqrt(r2)*iRmax);
            }
        }

    }
}