EMAN2
Public Member Functions | Static Public Member Functions | Static Public Attributes | Private Member Functions
EMAN::FourierInserter3DMode5 Class Reference

FourierPixelInserter3DMode5 - encapsulates "method 5" 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::FourierInserter3DMode5:
Inheritance graph
[legend]
Collaboration diagram for EMAN::FourierInserter3DMode5:
Collaboration graph
[legend]

List of all members.

Public Member Functions

 FourierInserter3DMode5 ()
virtual ~FourierInserter3DMode5 ()
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.
virtual string get_name () const
 Get the unique name of this class (especially for factory based instantiation access)
virtual string get_desc () const
 Get a clear, concise description of this class.

Static Public Member Functions

static FourierPixelInserter3DNEW ()

Static Public Attributes

static const string NAME = "gauss_5"

Private Member Functions

 FourierInserter3DMode5 (const FourierInserter3DMode5 &)
FourierInserter3DMode5operator= (const FourierInserter3DMode5 &)

Detailed Description

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

Definition at line 259 of file reconstructor_tools.h.


Constructor & Destructor Documentation

EMAN::FourierInserter3DMode5::FourierInserter3DMode5 ( ) [inline]

Definition at line 262 of file reconstructor_tools.h.

Referenced by NEW().

                        {
//                              gimx = EMAN::Interp::get_gimx();
                        }
virtual EMAN::FourierInserter3DMode5::~FourierInserter3DMode5 ( ) [inline, virtual]

Definition at line 267 of file reconstructor_tools.h.

                        {
                                // Don't delete gimx it causes a seg fault
//                              if ( gimx != 0 )
//                              {
//                                      delete gimx;
//                                      gimx = 0;
//                              }
                        }
EMAN::FourierInserter3DMode5::FourierInserter3DMode5 ( const FourierInserter3DMode5 ) [private]

Member Function Documentation

virtual string EMAN::FourierInserter3DMode5::get_desc ( ) const [inline, virtual]

Get a clear, concise description of this class.

Returns:
a clear, concise description of this class

Implements EMAN::FactoryBase.

Definition at line 289 of file reconstructor_tools.h.

                        {
                                return "Fourier pixel insertion mode 5";
                        }
virtual string EMAN::FourierInserter3DMode5::get_name ( ) const [inline, virtual]

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 284 of file reconstructor_tools.h.

References NAME.

                        {
                                return NAME;
                        }
bool FourierInserter3DMode5::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:
xxthe floating point x coordinate
yythe floating point y coordinate
zzthe floating point z coordinate
dtthe complex pixel value (dt[0] is real, dt[1] is imaginary)
weightthe 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 235 of file reconstructor_tools.cpp.

References EMAN::EMData::add_complex_at_fast(), EMAN::FourierPixelInserter3D::data, EMAN::Util::fast_exp(), EMAN::Util::hypot3(), EMAN::Util::hypot3sq(), EMAN::EMConsts::I3G, EMAN::FourierPixelInserter3D::norm, EMAN::FourierPixelInserter3D::nx2, EMAN::FourierPixelInserter3D::ny2, EMAN::FourierPixelInserter3D::nz2, and EMAN::FourierPixelInserter3D::subx0.

{
        int x0 = (int) floor(xx-1.5);
        int y0 = (int) floor(yy-1.5);
        int z0 = (int) floor(zz-1.5);
        
        if (subx0<0) {                  // normal full reconstruction
                if (x0<-nx2-4 || y0<-ny2-4 || z0<-nz2-4 || x0>nx2+3 || y0>ny2+3 || z0>nz2+3 ) return false;

                // no error checking on add_complex_fast, so we need to be careful here
                int x1=x0+4;
                int y1=y0+4;
                int z1=z0+4;
                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=2.0/EMConsts::I3G;
                float h=32.0f/((8.0f+Util::hypot3(xx,yy,zz))*EMConsts::I3G); 
                float w=weight/(1.0f+6.0f*Util::fast_exp(-h)+12*Util::fast_exp(-h*2.0f)+8*Util::fast_exp(-h*3.0f)+
                        6.0f*Util::fast_exp(-h*4.0f)+24.0f*Util::fast_exp(-h*5.0f)+24.0f*Util::fast_exp(-h*6.0f)+12.0f*Util::fast_exp(-h*8.0f)+
                        24.0f*Util::fast_exp(-h*9.0f)+8.0f*Util::fast_exp(-h*12.0f));   // approx normalization so higer radii aren't upweighted relative to lower due to wider Gaussian
                //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)*w;
//                                      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);
                                        norm[off/2]+=gg;
                                }
                        }
                }
                return true;
        }
        printf("region writing not supported in mode 3\n");
        return false;
}
static FourierPixelInserter3D* EMAN::FourierInserter3DMode5::NEW ( ) [inline, static]

Definition at line 279 of file reconstructor_tools.h.

References FourierInserter3DMode5().

                        {
                                return new FourierInserter3DMode5();
                        }
FourierInserter3DMode5& EMAN::FourierInserter3DMode5::operator= ( const FourierInserter3DMode5 ) [private]

Member Data Documentation

const string FourierInserter3DMode5::NAME = "gauss_5" [static]

Definition at line 294 of file reconstructor_tools.h.

Referenced by get_name().


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