EMAN::TestImageHollowEllipse Class Reference

Generate an ellipse/ellipsoid image with an inner hollow ellipse/ellipsoid. More...

#include <processor.h>

Inheritance diagram for EMAN::TestImageHollowEllipse:

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

Collaboration graph
[legend]
List of all members.

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 ProcessorNEW ()

Static Public Attributes

static const string NAME = "testimage.ellipsoid.hollow"

Detailed Description

Generate an ellipse/ellipsoid image with an inner hollow ellipse/ellipsoid.

Parameters:
xwidth inner equatorial radii along x axes
ywidth inner equatorial radii along y axes
zwidth inner polar radius
a outter equatorial radii along x axes
b outter equatorial radii along y axes
c outter polar radius
width specify the width or specify each width explicitly - xwidth, ywidth, zwidth
transform Optionally transform the ellipse
fill value you want to fill in hollow ellipse, default to 1.0

Definition at line 6364 of file processor.h.


Member Function Documentation

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

Get the descrition of this specific processor.

This function must be overwritten by a subclass.

Returns:
The description of this processor.

Implements EMAN::Processor.

Definition at line 6374 of file processor.h.

06375                         {
06376                                 return "Insert a hollow ellipse into the image.";
06377                         }

virtual string EMAN::TestImageHollowEllipse::get_name (  )  const [inline, virtual]

Get the processor's name.

Each processor is identified by a unique name.

Returns:
The processor's name.

Implements EMAN::Processor.

Definition at line 6369 of file processor.h.

References NAME.

06370                         {
06371                                 return NAME;
06372                         }

virtual TypeDict EMAN::TestImageHollowEllipse::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.

Returns:
A dictionary containing the parameter info.

Reimplemented from EMAN::Processor.

Definition at line 6384 of file processor.h.

References EMAN::EMObject::FLOAT, EMAN::TypeDict::put(), and EMAN::EMObject::TRANSFORM.

06385                         {
06386                                 TypeDict d;
06387                                 d.put("xwidth", EMObject::FLOAT, "inner equatorial radii along x axes");
06388                                 d.put("ywidth", EMObject::FLOAT, "inner equatorial radii along y axes");
06389                                 d.put("zwidth", EMObject::FLOAT, "inner polar radius");
06390                                 d.put("a", EMObject::FLOAT, "outter equatorial radii along x axes");
06391                                 d.put("b", EMObject::FLOAT, "outter equatorial radii along y axes");
06392                                 d.put("c", EMObject::FLOAT, "outter polar radius");
06393                                 d.put("width",EMObject::FLOAT, "width - specify the width or specify each width explicitly - xwidth, ywidth, zwidth");
06394                                 d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06395                                 d.put("fill", EMObject::FLOAT, "value you want to fill in hollow ellipse, default to 1.0");
06396                                 return d;
06397                         }

static Processor* EMAN::TestImageHollowEllipse::NEW (  )  [inline, static]

Definition at line 6379 of file processor.h.

06380                         {
06381                                 return new TestImageHollowEllipse();
06382                         }

void TestImageHollowEllipse::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.

Parameters:
image The image to be processed.

Implements EMAN::Processor.

Definition at line 7506 of file processor.cpp.

References EMAN::Dict::has_key(), EMAN::Processor::params, EMAN::TestImageProcessor::preprocess(), EMAN::Dict::set_default(), EMAN::EMData::set_value_at(), t, EMAN::EMData::update(), and v.

07507 {
07508         preprocess(image);
07509 
07510         float width = params.set_default("width",2.0f);
07511 
07512         float a2 = params.set_default("a",nx/2.0f-1.0f);
07513         float b2 = params.set_default("b",ny/2.0f-1.0f);
07514         float c2 = params.set_default("c",nz/2.0f-1.0f);
07515 
07516         float a1 = params.set_default("xwidth",a2-width);
07517         float b1 = params.set_default("ywidth",b2-width);
07518         float c1 = params.set_default("zwidth",c2-width);
07519 
07520         float fill = params.set_default("fill",1.0f);
07521         Transform* t;
07522         if (params.has_key("transform")) {
07523                 t = params["transform"];
07524         } else {
07525                 t = new Transform;
07526         }
07527 
07528 
07529         int mz = 2*(int)c2+1;
07530         if ( nz < mz ) mz = nz;
07531         int my = 2*(int)b2+1;
07532         if ( ny < my ) my = ny;
07533         int mx = 2*(int)a2+1;
07534         if ( nx < mx ) mx = nx;
07535 
07536         float ai1 = 1/(a1*a1);
07537         float bi1 = 1/(b1*b1);
07538         float ci1 = 1/(c1*c1);
07539 
07540         float ai2 = 1/(a2*a2);
07541         float bi2 = 1/(b2*b2);
07542         float ci2 = 1/(c2*c2);
07543 
07544         Vec3f origin(nx/2,ny/2,nz/2);
07545 
07546         float x2, y2, z2, r1,r2;
07547         int xl, yl, zl;
07548         for (int k = 0; k < mz; ++k) {
07549                 for (int j = 0; j < my; ++j) {
07550                         for (int i = 0; i < mx; ++i) {
07551                                 x2 = (float)(i - mx/2);
07552                                 y2 = (float)(j - my/2);
07553                                 z2 = (float)(k - mz/2);
07554                                 r1 = (x2*x2)*ai1 + (y2*y2)*bi1 + (z2*z2)*ci1;
07555                                 r2 = (x2*x2)*ai2 + (y2*y2)*bi2 + (z2*z2)*ci2;
07556                                 if (r2 <= 1 && r1 >= 1) {
07557 
07558                                         if (t != 0) {
07559                                                 Vec3f v(x2,y2,z2);
07560                                                 v = (*t)*v;
07561                                                 v += origin;
07562 
07563                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07564                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07565 
07566                                                 for( int kk = -1; kk <= 1; ++kk)
07567                                                         for( int jj = -1; jj <= 1; ++jj)
07568                                                                 for( int ii = -1; ii <= 1; ++ii) {
07569                                                                         xl = (int)v[0]+ii;
07570                                                                         yl = (int)v[1]+jj;
07571                                                                         zl = (int)v[2]+kk;
07572                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07573                                                                                 image->set_value_at(xl,yl,zl,1.0);
07574                                                                 }
07575                                         } else {
07576                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07577                                         }
07578                                 }
07579                         }
07580                 }
07581         }
07582 
07583         delete t;
07584 
07585         image->update();
07586 }


Member Data Documentation

const string TestImageHollowEllipse::NAME = "testimage.ellipsoid.hollow" [static]

Definition at line 6399 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Tue May 25 17:17:37 2010 for EMAN2 by  doxygen 1.4.7