EMAN::TestImageEllipse Class Reference

Generate an ellipse or ellipsoid image. More...

#include <processor.h>

Inheritance diagram for EMAN::TestImageEllipse:

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

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"

Detailed Description

Generate an ellipse or ellipsoid image.

Parameters:
a equatorial radii along x axes
b equatorial radii along y axes
c polar radius
transform Optionally transform the ellipse
fill value you want to fill in ellipse, default to 1.0

Definition at line 6287 of file processor.h.


Member Function Documentation

virtual string EMAN::TestImageEllipse::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 6297 of file processor.h.

06298                 {
06299                         return "Insert an ellipse into the image.";
06300                 }

virtual string EMAN::TestImageEllipse::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 6292 of file processor.h.

References NAME.

06293                 {
06294                         return NAME;
06295                 }

virtual TypeDict EMAN::TestImageEllipse::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 6307 of file processor.h.

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

06308                 {
06309                         TypeDict d;
06310                         d.put("a", EMObject::FLOAT, "equatorial radius along x axes (major semiaxes)");
06311                         d.put("b", EMObject::FLOAT, "equatorial radius along y axes (minor semiaxes)");
06312                         d.put("c", EMObject::FLOAT, "polar radius for ellipsoid (x^2/a^2+y^2/b^2+z^2/c^2=1)");
06313                         d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06314                         d.put("fill", EMObject::FLOAT, "value you want to fill in ellipse, default to 1.0");
06315                         return d;
06316                 }

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

Definition at line 6302 of file processor.h.

06303                 {
06304                         return new TestImageEllipse();
06305                 }

void TestImageEllipse::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 7662 of file processor.cpp.

References b, bi, 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.

07663 {
07664         preprocess(image);
07665 
07666 
07667         float a = params.set_default("a",nx/2.0f-1.0f);
07668         float b = params.set_default("b",ny/2.0f-1.0f);
07669         float c = params.set_default("c",nz/2.0f-1.0f);
07670         float fill = params.set_default("fill",1.0f);
07671         //bool hollow = params.set_default("hollow",false);
07672         Transform* t;
07673         if (params.has_key("transform")) {
07674                 t = params["transform"];
07675         } else {
07676                 t = new Transform;
07677         }
07678 
07679 
07680         int mz = 2*(int)c+1;
07681         if ( nz < mz ) mz = nz;
07682         int my = 2*(int)b+1;
07683         if ( ny < my ) my = ny;
07684         int mx = 2*(int)a+1;
07685         if ( nx < mx ) mx = nx;
07686 
07687 
07688         float ai = 1/(a*a);
07689         float bi = 1/(b*b);
07690         float ci = 1/(c*c);
07691 
07692         Vec3f origin(nx/2,ny/2,nz/2);
07693 
07694         float x2, y2, z2, r;
07695         int xl, yl, zl;
07696         for (int k = 0; k < mz; ++k) {
07697                 for (int j = 0; j < my; ++j) {
07698                         for (int i = 0; i < mx; ++i) {
07699                                 x2 = (float)(i - mx/2);
07700                                 y2 = (float)(j - my/2);
07701                                 z2 = (float)(k - mz/2);
07702                                 r = (x2*x2)*ai + (y2*y2)*bi + (z2*z2)*ci;
07703                                 if (r <= 1) {
07704 
07705                                         if (t != 0) {
07706                                                 Vec3f v(x2,y2,z2);
07707                                                 v = (*t)*v;
07708                                                 v += origin;
07709 
07710                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07711                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07712 
07713                                                 for( int kk = -1; kk <= 1; ++kk)
07714                                                         for( int jj = -1; jj <= 1; ++jj)
07715                                                                 for( int ii = -1; ii <= 1; ++ii) {
07716                                                                         xl = (int)v[0]+ii;
07717                                                                         yl = (int)v[1]+jj;
07718                                                                         zl = (int)v[2]+kk;
07719                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07720                                                                                 image->set_value_at(xl,yl,zl,fill);
07721                                                                 }
07722                                         } else {
07723                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07724                                         }
07725                                 }
07726                         }
07727                 }
07728         }
07729 
07730         delete t;
07731 
07732         image->update();
07733 }


Member Data Documentation

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

Definition at line 6318 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Mon May 2 13:31:06 2011 for EMAN2 by  doxygen 1.4.7