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 6319 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 6329 of file processor.h.

06330                 {
06331                         return "Insert an ellipse into the image.";
06332                 }

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 6324 of file processor.h.

References NAME.

06325                 {
06326                         return NAME;
06327                 }

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 6339 of file processor.h.

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

06340                 {
06341                         TypeDict d;
06342                         d.put("a", EMObject::FLOAT, "equatorial radii along x axes");
06343                         d.put("b", EMObject::FLOAT, "equatorial radii along y axes");
06344                         d.put("c", EMObject::FLOAT, "polar radius");
06345                         d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06346                         d.put("fill", EMObject::FLOAT, "value you want to fill in ellipse, default to 1.0");
06347                         return d;
06348                 }

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

Definition at line 6334 of file processor.h.

06335                 {
06336                         return new TestImageEllipse();
06337                 }

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 7588 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.

07589 {
07590         preprocess(image);
07591 
07592 
07593         float a = params.set_default("a",nx/2.0f-1.0f);
07594         float b = params.set_default("b",ny/2.0f-1.0f);
07595         float c = params.set_default("c",nz/2.0f-1.0f);
07596         float fill = params.set_default("fill",1.0f);
07597         //bool hollow = params.set_default("hollow",false);
07598         Transform* t;
07599         if (params.has_key("transform")) {
07600                 t = params["transform"];
07601         } else {
07602                 t = new Transform;
07603         }
07604 
07605 
07606         int mz = 2*(int)c+1;
07607         if ( nz < mz ) mz = nz;
07608         int my = 2*(int)b+1;
07609         if ( ny < my ) my = ny;
07610         int mx = 2*(int)a+1;
07611         if ( nx < mx ) mx = nx;
07612 
07613 
07614         float ai = 1/(a*a);
07615         float bi = 1/(b*b);
07616         float ci = 1/(c*c);
07617 
07618         Vec3f origin(nx/2,ny/2,nz/2);
07619 
07620         float x2, y2, z2, r;
07621         int xl, yl, zl;
07622         for (int k = 0; k < mz; ++k) {
07623                 for (int j = 0; j < my; ++j) {
07624                         for (int i = 0; i < mx; ++i) {
07625                                 x2 = (float)(i - mx/2);
07626                                 y2 = (float)(j - my/2);
07627                                 z2 = (float)(k - mz/2);
07628                                 r = (x2*x2)*ai + (y2*y2)*bi + (z2*z2)*ci;
07629                                 if (r <= 1) {
07630 
07631                                         if (t != 0) {
07632                                                 Vec3f v(x2,y2,z2);
07633                                                 v = (*t)*v;
07634                                                 v += origin;
07635 
07636                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07637                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07638 
07639                                                 for( int kk = -1; kk <= 1; ++kk)
07640                                                         for( int jj = -1; jj <= 1; ++jj)
07641                                                                 for( int ii = -1; ii <= 1; ++ii) {
07642                                                                         xl = (int)v[0]+ii;
07643                                                                         yl = (int)v[1]+jj;
07644                                                                         zl = (int)v[2]+kk;
07645                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07646                                                                                 image->set_value_at(xl,yl,zl,fill);
07647                                                                 }
07648                                         } else {
07649                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07650                                         }
07651                                 }
07652                         }
07653                 }
07654         }
07655 
07656         delete t;
07657 
07658         image->update();
07659 }


Member Data Documentation

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

Definition at line 6350 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:36 2010 for EMAN2 by  doxygen 1.4.7