Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members

EMAN::FourierInserter3DMode3 Class Reference

FourierPixelInserter3DMode3 - encapsulates "method 3" 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::FourierInserter3DMode3:

Inheritance graph
[legend]
Collaboration diagram for EMAN::FourierInserter3DMode3:

Collaboration graph
[legend]
List of all members.

Public Member Functions

 FourierInserter3DMode3 ()
virtual ~FourierInserter3DMode3 ()
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

FourierPixelInserter3DNEW ()

Static Public Attributes

const string NAME = "gauss_3"

Private Member Functions

 FourierInserter3DMode3 (const FourierInserter3DMode3 &)
FourierInserter3DMode3operator= (const FourierInserter3DMode3 &)

Detailed Description

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

Definition at line 225 of file reconstructor_tools.h.


Constructor & Destructor Documentation

EMAN::FourierInserter3DMode3::FourierInserter3DMode3  )  [inline]
 

Definition at line 228 of file reconstructor_tools.h.

00228 {}

virtual EMAN::FourierInserter3DMode3::~FourierInserter3DMode3  )  [inline, virtual]
 

Definition at line 229 of file reconstructor_tools.h.

00229 {}

EMAN::FourierInserter3DMode3::FourierInserter3DMode3 const FourierInserter3DMode3  )  [private]
 


Member Function Documentation

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

00244                         {
00245                                 return "Fourier pixel insertion using a 3x3x3 Gaussian kernel";
00246                         }

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

00239                         {
00240                                 return NAME;
00241                         }

bool FourierInserter3DMode3::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 185 of file reconstructor_tools.cpp.

References EMAN::EMData::add_complex_at_fast(), dt, EMAN::Util::fast_exp(), EMAN::Util::hypot3(), EMAN::Util::hypot3sq(), and weight.

00186 {
00187         int x0 = (int) floor(xx-.5);
00188         int y0 = (int) floor(yy-.5);
00189         int z0 = (int) floor(zz-.5);
00190         
00191         if (subx0<0) {                  // normal full reconstruction
00192                 if (x0<-nx2-2 || y0<-ny2-2 || z0<-nz2-2 || x0>nx2+1 || y0>ny2+1 || z0>nz2+1 ) return false;
00193 
00194                 // no error checking on add_complex_fast, so we need to be careful here
00195                 int x1=x0+2;
00196                 int y1=y0+2;
00197                 int z1=z0+2;
00198                 if (x0<-nx2) x0=-nx2;
00199                 if (x1>nx2) x1=nx2;
00200                 if (y0<-ny2) y0=-ny2;
00201                 if (y1>ny2) y1=ny2;
00202                 if (z0<-nz2) z0=-nz2;
00203                 if (z1>nz2) z1=nz2;
00204                 
00205 //              float h=2.0/((1.0+pow(Util::hypot3sq(xx,yy,zz),.5))*EMConsts::I2G);
00206 //              float h=2.0/EMConsts::I3G;
00207                 float h=32.0f/((8.0f+Util::hypot3(xx,yy,zz))*EMConsts::I3G);
00208 //              float w=weight;
00209                 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));     // approx normalization so higer radii aren't upweighted relative to lower due to wider Gaussian
00210                 //size_t idx;
00211                 float r, gg;
00212 //              int pc=0;
00213                 for (int k = z0 ; k <= z1; k++) {
00214                         for (int j = y0 ; j <= y1; j++) {
00215                                 for (int i = x0; i <= x1; i ++) {
00216                                         r = Util::hypot3sq((float) i - xx, (float)j - yy, (float)k - zz);
00217 //                                      gg=weight;
00218                                         gg = Util::fast_exp(-r *h)*w;
00219 //                                      gg = Util::fast_exp(-r / EMConsts::I2G)*weight;
00220 //                                      gg = sqrt(Util::fast_exp(-r / EMConsts::I2G))*weight;
00221                                         
00222                                         size_t off;
00223                                         off=data->add_complex_at_fast(i,j,k,dt*gg);
00224                                         norm[off/2]+=gg;
00225                                 }
00226                         }
00227                 }
00228                 return true;
00229         }
00230         printf("region writing not supported in mode 3\n");
00231         return false;
00232 }

FourierPixelInserter3D* EMAN::FourierInserter3DMode3::NEW  )  [inline, static]
 

Definition at line 233 of file reconstructor_tools.h.

00234                         {
00235                                 return new FourierInserter3DMode3();
00236                         }

FourierInserter3DMode3& EMAN::FourierInserter3DMode3::operator= const FourierInserter3DMode3  )  [private]
 


Member Data Documentation

const string FourierInserter3DMode3::NAME = "gauss_3" [static]
 

Definition at line 43 of file reconstructor_tools.cpp.


The documentation for this class was generated from the following files:
Generated on Tue Jun 11 13:49:30 2013 for EMAN2 by  doxygen 1.3.9.1