FourierPixelInserter3DMode2 - encapsulates "method 2" for inserting a 2D Fourier slice into a 3D volume See comments in FourierPixelInserter3D for explanations. More...

#include <reconstructor_tools.h>

Inheritance diagram for EMAN::FourierInserter3DMode2:

[legend]

Collaboration diagram for EMAN::FourierInserter3DMode2:

[legend]

Public Member Functions
	FourierInserter3DMode2 ()

virtual	~FourierInserter3DMode2 ()

virtual bool	insert_pixel (const float &xx, const float &yy, const float &zz, const std::complex< float > dt, const float &weight=1.0)
	Insert a complex pixel [dt[0]+dt[1]i] at (float) coordinate [xx,yy,zz] with weighting into a discrete 3D volume. More...

virtual string	get_name () const
	Get the unique name of this class (especially for factory based instantiation access) More...

virtual string	get_desc () const
	Get a clear, concise description of this class. More...

	FourierInserter3DMode2 (const FourierInserter3DMode2 &)

FourierInserter3DMode2 &	operator= (const FourierInserter3DMode2 &)

Public Member Functions inherited from EMAN::FourierPixelInserter3D
	FourierPixelInserter3D ()
	Construct a FourierPixelInserter3D. More...

virtual	~FourierPixelInserter3D ()
	Desctruct a FourierPixelInserter3D. More...

TypeDict	get_param_types () const

virtual void	init ()

Public Member Functions inherited from EMAN::FactoryBase
	FactoryBase ()

virtual	~FactoryBase ()

Dict	get_params () const
	get a copy of the parameters of this class More...

void	set_params (const Dict &new_params)
	Set new parameters. More...

void	set_param (const string key, const EMObject val)

void	insert_params (const Dict &new_params)
	Insert parameters. More...

Dict	copy_relevant_params (const FactoryBase *const that) const

Static Public Member Functions
static FourierPixelInserter3D *	NEW ()

Static Public Attributes
static const string	NAME = "gauss_2"

Additional Inherited Members
Protected Attributes inherited from EMAN::FourierPixelInserter3D
float *	norm
	A pointer to the constructor argument normalize_values. More...

EMData *	data
	A pointer to the constructor argument real_data. More...

int	nx
	Image volume data sizes a convenience variable used here and there. More...

int	ny

int	nz

int	nxyz

int	nx2

int	ny2

int	nz2

int	subx0

int	suby0

int	subz0

int	fullnx

int	fullny

int	fullnz

Protected Attributes inherited from EMAN::FactoryBase
Dict	params
	This is the dictionary the stores the parameters of the object. More...

Detailed Description

FourierPixelInserter3DMode2 - encapsulates "method 2" for inserting a 2D Fourier slice into a 3D volume See comments in FourierPixelInserter3D for explanations.

Definition at line 193 of file reconstructor_tools.h.

Constructor & Destructor Documentation

◆ FourierInserter3DMode2() [1/2]

EMAN::FourierInserter3DMode2::FourierInserter3DMode2 ( )

inline

Definition at line 196 of file reconstructor_tools.h.

196{}

Referenced by NEW().

◆ ~FourierInserter3DMode2()

virtual EMAN::FourierInserter3DMode2::~FourierInserter3DMode2 ( )

inlinevirtual

Definition at line 197 of file reconstructor_tools.h.

197{}

◆ FourierInserter3DMode2() [2/2]

EMAN::FourierInserter3DMode2::FourierInserter3DMode2 ( const FourierInserter3DMode2 & )

Member Function Documentation

◆ get_desc()

virtual string EMAN::FourierInserter3DMode2::get_desc ( ) const

inlinevirtual

Get a clear, concise description of this class.

Returns: a clear, concise description of this class

Implements EMAN::FactoryBase.

Definition at line 211 of file reconstructor_tools.h.

                        {
                                return "Fourier pixel insertion using interpolation and the nearest 8 voxels";
                        }

◆ get_name()

virtual string EMAN::FourierInserter3DMode2::get_name ( ) const

inlinevirtual

Get the unique name of this class (especially for factory based instantiation access)

Returns: the unique name of this class

Implements EMAN::FactoryBase.

Definition at line 206 of file reconstructor_tools.h.

                        {
                                return NAME;
                        }

References NAME.

◆ insert_pixel()

bool FourierInserter3DMode2::insert_pixel	(	const float &	xx,
		const float &	yy,
		const float &	zz,
		const std::complex< float >	dt,
		const float &	weight = `1.0`
	)

virtual

Insert a complex pixel [dt[0]+dt[1]i] at (float) coordinate [xx,yy,zz] with weighting into a discrete 3D volume.

Parameters

xx	the floating point x coordinate
yy	the floating point y coordinate
zz	the floating point z coordinate
dt	the complex pixel value (dt[0] is real, dt[1] is imaginary)
weight	the weight to given to this complex pixel

Returns: A boolean that indicates the pixel has been inserted (or not)

Implements EMAN::FourierPixelInserter3D.

Definition at line 128 of file reconstructor_tools.cpp.

{
        int x0 = (int) floor(xx);
        int y0 = (int) floor(yy);
        int z0 = (int) floor(zz);
 
        // note that subnx differs in the inserters. In the reconstructors it subx0 is 0 for the full volume. Here it is -1
        if (subx0<0) {                  // normal full reconstruction
                if (x0<-nx2-1 || y0<-ny2-1 || z0<-nz2-1 || x0>nx2 || y0>ny2 || z0>nz2 ) return false;
 
                int x1=x0+1;
                int y1=y0+1;
                int z1=z0+1;
//              if (x0<-nx2) x0=-nx2;
//              if (x1>nx2) x1=nx2;
//              if (y0<-ny2) y0=-ny2;
//              if (y1>ny2) y1=ny2;
//              if (z0<-nz2) z0=-nz2;
//              if (z1>nz2) z1=nz2;
 
//              float h=2.0/((1.0+pow(Util::hypot3sq(xx,yy,zz),.5))*EMConsts::I2G);
                float h=1.0f/EMConsts::I2G;
                //size_t idx;
                float r, gg;
//              int pc=0;
                for (int k = z0 ; k <= z1; k++) {
                        for (int j = y0 ; j <= y1; j++) {
                                for (int i = x0; i <= x1; i ++) {
                                        r = Util::hypot3sq((float) i - xx, j - yy, k - zz);
                                        gg = Util::fast_exp(-r *h)*weight;
 
                                        size_t off;
                                        off=data->add_complex_at_fast(i,j,k,dt*gg);
                                        if (off!=nxyz) norm[off/2]+=gg;
//                                      if (off!=nxyz) norm[off/2]+=weight;     // experiment 6/1/20, use true Gaussian kernel, not just weight
                                }
                        }
                }
                return true;
        }
        else {                                  // for subvolumes, not optimized yet
                //size_t idx;
                float r, gg;
                int pc=0;
                for (int k = z0 ; k <= z0 + 1; k++) {
                        for (int j = y0 ; j <= y0 + 1; j++) {
                                for (int i = x0; i <= x0 + 1; i ++) {
                                        r = Util::hypot3sq((float) i - xx, j - yy, k - zz);
                                        gg = Util::fast_exp(-r / EMConsts::I2G)*weight;
 
                                        size_t off;
                                        if (subx0<0) off=data->add_complex_at(i,j,k,dt*gg);
                                        else off=data->add_complex_at(i,j,k,subx0,suby0,subz0,fullnx,fullny,fullnz,dt*gg);
                                        if (static_cast<int>(off)!=nxyz) { norm[off/2]+=gg; pc+=1; }
                                }
                        }
                }
 
                if (pc>0)  return true;
                return false;
        }
}

◆ NEW()

static FourierPixelInserter3D * EMAN::FourierInserter3DMode2::NEW ( )

inlinestatic

Definition at line 201 of file reconstructor_tools.h.

                        {
                                return new FourierInserter3DMode2();
                        }

References FourierInserter3DMode2().

◆ operator=()

FourierInserter3DMode2 & EMAN::FourierInserter3DMode2::operator= ( const FourierInserter3DMode2 & )

Member Data Documentation

◆ NAME

const string FourierInserter3DMode2::NAME = "gauss_2"

static

Definition at line 216 of file reconstructor_tools.h.

Referenced by get_name().

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

Public Member Functions

Static Public Member Functions

Static Public Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FourierInserter3DMode2() [1/2]

◆ ~FourierInserter3DMode2()

◆ FourierInserter3DMode2() [2/2]

Member Function Documentation

◆ get_desc()

◆ get_name()

◆ insert_pixel()

◆ NEW()

◆ operator=()

Member Data Documentation

◆ NAME