#include <processor.h>
Inheritance diagram for EMAN::TestImagePureGaussian:
Public Member Functions | |
virtual void | process_inplace (EMData *image) |
To process an image in-place. | |
virtual string | get_name () const |
Get the processor's name. | |
virtual string | get_desc () const |
Get the descrition of this specific processor. | |
virtual TypeDict | get_param_types () const |
Get processor parameter information in a dictionary. | |
Static Public Member Functions | |
static Processor * | NEW () |
Static Public Attributes | |
static const string | NAME = "testimage.puregaussian" |
x_sigma | sigma value for this Gaussian blob on x direction | |
y_sigma | sigma value for this Gaussian blob on y direction | |
z_sigma | sigma value for this Gaussian blob on z direction | |
x_center | center for this Gaussian blob on x direction | |
y_center | center for this Gaussian blob on y direction | |
z_center | center for this Gaussian blob on z direction |
Definition at line 5753 of file processor.h.
virtual string EMAN::TestImagePureGaussian::get_desc | ( | ) | const [inline, virtual] |
Get the descrition of this specific processor.
This function must be overwritten by a subclass.
Implements EMAN::Processor.
Definition at line 5763 of file processor.h.
virtual string EMAN::TestImagePureGaussian::get_name | ( | ) | const [inline, virtual] |
Get the processor's name.
Each processor is identified by a unique name.
Implements EMAN::Processor.
Definition at line 5758 of file processor.h.
References NAME.
05759 { 05760 return NAME; 05761 }
virtual TypeDict EMAN::TestImagePureGaussian::get_param_types | ( | ) | const [inline, virtual] |
Get processor parameter information in a dictionary.
Each parameter has one record in the dictionary. Each record contains its name, data-type, and description.
Reimplemented from EMAN::Processor.
Definition at line 5773 of file processor.h.
References EMAN::EMObject::FLOAT, and EMAN::TypeDict::put().
05774 { 05775 TypeDict d; 05776 d.put("x_sigma", EMObject::FLOAT, "sigma value for this Gaussian blob on x direction"); 05777 d.put("y_sigma", EMObject::FLOAT, "sigma value for this Gaussian blob on y direction"); 05778 d.put("z_sigma", EMObject::FLOAT, "sigma value for this Gaussian blob on z direction"); 05779 d.put("x_center", EMObject::FLOAT, "center for this Gaussian blob on x direction" ); 05780 d.put("y_center", EMObject::FLOAT, "center for this Gaussian blob on y direction" ); 05781 d.put("z_center", EMObject::FLOAT, "center for this Gaussian blob on z direction" ); 05782 return d; 05783 }
static Processor* EMAN::TestImagePureGaussian::NEW | ( | ) | [inline, static] |
Definition at line 5768 of file processor.h.
05769 { 05770 return new TestImagePureGaussian(); 05771 }
void TestImagePureGaussian::process_inplace | ( | EMData * | image | ) | [virtual] |
To process an image in-place.
For those processors which can only be processed out-of-place, override this function to just print out some error message to remind user call the out-of-place version.
image | The image to be processed. |
Implements EMAN::Processor.
Definition at line 7204 of file processor.cpp.
References EMAN::EMData::get_xsize(), EMAN::EMData::get_ysize(), EMAN::EMData::get_zsize(), norm(), EMAN::Processor::params, pi, EMAN::TestImageProcessor::preprocess(), sqrt(), EMAN::EMData::update(), and x.
07205 { 07206 preprocess(image); 07207 07208 float x_sigma = params["x_sigma"]; 07209 float y_sigma = params["y_sigma"]; 07210 float z_sigma = params["z_sigma"]; 07211 07212 float x_center = params["x_center"]; 07213 float y_center = params["y_center"]; 07214 float z_center = params["z_center"]; 07215 07216 int nx = image->get_xsize(); 07217 int ny = image->get_ysize(); 07218 int nz = image->get_zsize(); 07219 07220 float x_twosig2 = 2*x_sigma*x_sigma; 07221 float y_twosig2 = 2*y_sigma*y_sigma; 07222 float z_twosig2 = 2*z_sigma*z_sigma; 07223 07224 float sr2pi = sqrt( 2.0f*(float)pi ); 07225 float norm = 1.0f/ ( x_sigma*sr2pi ); 07226 if (ny > 1) { 07227 norm *= 1.0f/ ( y_sigma*sr2pi ); 07228 if (nz > 1) norm *= 1.0f/ ( z_sigma*sr2pi ); 07229 } 07230 07231 float z, y, x, sum, val; 07232 for (int iz=0; iz < nz; ++iz) { 07233 z = static_cast<float>(iz) - z_center; 07234 for (int iy=0; iy < ny; ++iy) { 07235 y = static_cast<float>(iy) - y_center; 07236 for (int ix=0; ix < nx; ++ix) { 07237 x = static_cast<float>(ix) - x_center; 07238 sum = x*x/x_twosig2 + y*y/y_twosig2 + z*z/z_twosig2; 07239 val = norm*exp(-sum); 07240 (*image)(ix,iy,iz) = val; 07241 } 07242 } 07243 } 07244 image->update(); 07245 }
const string TestImagePureGaussian::NAME = "testimage.puregaussian" [static] |