FourierPixelInserter3DMode8 - encapsulates "method 8" 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::FourierInserter3DMode8:

[legend]

Collaboration diagram for EMAN::FourierInserter3DMode8:

[legend]

Public Member Functions
	FourierInserter3DMode8 ()

virtual	~FourierInserter3DMode8 ()

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...

virtual void	init ()

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

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

TypeDict	get_param_types () const

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 = "experimental"

Private Member Functions
	FourierInserter3DMode8 (const FourierInserter3DMode8 &)

FourierInserter3DMode8 &	operator= (const FourierInserter3DMode8 &)

Private Attributes
int	mFreqCutoff

float	mDFreq

float *	W

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

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

Definition at line 414 of file reconstructor_tools.h.

Constructor & Destructor Documentation

◆ FourierInserter3DMode8() [1/2]

EMAN::FourierInserter3DMode8::FourierInserter3DMode8 ( )

inline

Definition at line 417 of file reconstructor_tools.h.

Referenced by NEW().

◆ ~FourierInserter3DMode8()

virtual EMAN::FourierInserter3DMode8::~FourierInserter3DMode8 ( )

inlinevirtual

Definition at line 421 of file reconstructor_tools.h.

                        {
                                if ( W != 0 )
                                        delete [] W;
                        }

References W.

◆ FourierInserter3DMode8() [2/2]

EMAN::FourierInserter3DMode8::FourierInserter3DMode8 ( const FourierInserter3DMode8 & )

private

Member Function Documentation

◆ get_desc()

virtual string EMAN::FourierInserter3DMode8::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 439 of file reconstructor_tools.h.

                        {
                                return "Fourier pixel insertion mode 8";
                        }

◆ get_name()

virtual string EMAN::FourierInserter3DMode8::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 434 of file reconstructor_tools.h.

                        {
                                return NAME;
                        }

References NAME.

◆ init()

void FourierInserter3DMode8::init ( )

virtual

Reimplemented from EMAN::FourierPixelInserter3D.

Definition at line 582 of file reconstructor_tools.cpp.

{
        FourierPixelInserter3D::init();
//      int P = (int)((1.0+0.25)*nx+1);
//      float r = (float)(nx+1)/(float)P;
//      mFreqCutoff = 2;
//      mDFreq = 0.2f;
//      if (W != 0) delete [] W;
//      W = Util::getBaldwinGridWeights(mFreqCutoff, (float)P, r,mDFreq,0.5f,0.2f);
 
}

References EMAN::FourierPixelInserter3D::init().

◆ insert_pixel()

bool FourierInserter3DMode8::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 594 of file reconstructor_tools.cpp.

{
        static FILE *out400=NULL,*out862=NULL,*out962,*out872,*out4093=NULL;
        
        int x0 = (int) floor(xx);
        int y0 = (int) floor(yy);
        int z0 = (int) floor(zz);
 
        if (out400==NULL) {
                out400=fopen("pxl4_0_0.txt","w");
                out862=fopen("pxl8_6_2.txt","w");
                out962=fopen("pxl9_6_2.txt","w");
                out872=fopen("pxl8_7_2.txt","w");
                out4093=fopen("pxl40_9_3.txt","w");
        }
                
        // 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=weight;
                                        gg = Util::fast_exp(-r *h)*weight;
//                                      gg = Util::fast_exp(-r / EMConsts::I2G)*weight;
//                                      gg = sqrt(Util::fast_exp(-r / EMConsts::I2G))*weight;
 
                                        size_t off;
                                        off=data->add_complex_at_fast(i,j,k,dt*gg);
//                                      off=data->add_complex_at(i,j,k,dt*gg);
                                        if (off!=nxyz) norm[off/2]+=gg;
                                        
                                        if (i==4&&j==0&&k==0) { fprintf(out400,"%1.4f\t%1.4f\t%1.4f\t%1.4f\t%1.4f\n",dt.real(),dt.imag(),gg,std::abs(dt),std::arg(dt)); fflush(out400); }
                                        if (i==8&&j==6&&k==2) { fprintf(out862,"%1.4f\t%1.4f\t%1.4f\t%1.4f\t%1.4f\n",dt.real(),dt.imag(),gg,std::abs(dt),std::arg(dt)); fflush(out862); }
                                        if (i==9&&j==6&&k==2) { fprintf(out962,"%1.4f\t%1.4f\t%1.4f\t%1.4f\t%1.4f\n",dt.real(),dt.imag(),gg,std::abs(dt),std::arg(dt)); fflush(out962); }
                                        if (i==8&&j==7&&k==2) { fprintf(out872,"%1.4f\t%1.4f\t%1.4f\t%1.4f\t%1.4f\n",dt.real(),dt.imag(),gg,std::abs(dt),std::arg(dt)); fflush(out872); }
                                        if (i==40&&j==9&&k==3) { fprintf(out4093,"%1.4f\t%1.4f\t%1.4f\t%1.4f\t%1.4f\n",dt.real(),dt.imag(),gg,std::abs(dt),std::arg(dt)); fflush(out4093); }
                                }
                        }
                }
                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;
        }
//      int x0 = (int) floor(qx);
//      int y0 = (int) floor(qy);
//      int z0 = (int) floor(qz);
 
//      int sizeW = (int)(1+2*mFreqCutoff/mDFreq);
//      int sizeWmid = sizeW/2;
 
//      for (int z = z0-mFreqCutoff; z < z0+mFreqCutoff; ++z){
//              for (int y = y0-mFreqCutoff; y < y0+mFreqCutoff; ++y){
//                      for (int x = x0-mFreqCutoff; x < x0+mFreqCutoff; ++x){
//                              if ( x < 0 || x >= nx ) continue;
//                              if ( y < 0 || y >= ny ) continue;
//                              if ( z < 0 || z >= nz ) continue;
//                              float dist = (float)((x-x0)*(x-x0)+(y-y0)*(y-y0)+(z-z0)*(z-z0));
//                              dist = sqrtf(dist);
//                              // We enforce a spherical kernel
//                              if ( dist > mFreqCutoff ) continue;
//                              int idx = (int)(sizeWmid + dist/mDFreq);
//                              if (idx >= sizeW) throw;
//                              float residual = dist/mDFreq - (int)(dist/mDFreq);
//                              if ( fabs(residual) > 1) throw;
//
//                              float factor = W[idx]*(1.0f-residual)+W[idx+1]*residual*weight;
//
//                              size_t k = z*nx*ny + y*nx + 2*x;
//
// //                           float c = Util::agauss(1, x-x0,y-y0,z-z0, EMConsts::I2G);
//                              rdata[k] += fq[0]*factor;
//                              rdata[k+1] += fq[1]*factor;
//
//
//                              norm[k/2] += weight;
//
//                      }
//              }
//      }
 
        return true;
}

◆ NEW()

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

inlinestatic

Definition at line 429 of file reconstructor_tools.h.

                        {
                                return new FourierInserter3DMode8();
                        }

References FourierInserter3DMode8().

◆ operator=()

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

private

Member Data Documentation

◆ mDFreq

float EMAN::FourierInserter3DMode8::mDFreq

private

Definition at line 450 of file reconstructor_tools.h.

◆ mFreqCutoff

int EMAN::FourierInserter3DMode8::mFreqCutoff

private

Definition at line 449 of file reconstructor_tools.h.

◆ NAME

const string FourierInserter3DMode8::NAME = "experimental"

static

Definition at line 446 of file reconstructor_tools.h.

Referenced by get_name().

◆ W

float* EMAN::FourierInserter3DMode8::W

private

Definition at line 455 of file reconstructor_tools.h.

Referenced by ~FourierInserter3DMode8().

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

Public Member Functions

Static Public Member Functions

Static Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FourierInserter3DMode8() [1/2]

◆ ~FourierInserter3DMode8()

◆ FourierInserter3DMode8() [2/2]

Member Function Documentation

◆ get_desc()

◆ get_name()

◆ init()

◆ insert_pixel()

◆ NEW()

◆ operator=()

Member Data Documentation

◆ mDFreq

◆ mFreqCutoff

◆ NAME

◆ W