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

06367                 {
06368                         return "Insert an ellipse into the image.";
06369                 }

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

References NAME.

06362                 {
06363                         return NAME;
06364                 }

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

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

06377                 {
06378                         TypeDict d;
06379                         d.put("a", EMObject::FLOAT, "equatorial radius along x axes (major semiaxes)");
06380                         d.put("b", EMObject::FLOAT, "equatorial radius along y axes (minor semiaxes)");
06381                         d.put("c", EMObject::FLOAT, "polar radius for ellipsoid (x^2/a^2+y^2/b^2+z^2/c^2=1)");
06382                         d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06383                         d.put("fill", EMObject::FLOAT, "value you want to fill in ellipse, default to 1.0");
06384                         return d;
06385                 }

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

Definition at line 6371 of file processor.h.

06372                 {
06373                         return new TestImageEllipse();
06374                 }

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

07777 {
07778         preprocess(image);
07779 
07780 
07781         float a = params.set_default("a",nx/2.0f-1.0f);
07782         float b = params.set_default("b",ny/2.0f-1.0f);
07783         float c = params.set_default("c",nz/2.0f-1.0f);
07784         float fill = params.set_default("fill",1.0f);
07785         //bool hollow = params.set_default("hollow",false);
07786         Transform* t;
07787         if (params.has_key("transform")) {
07788                 t = params["transform"];
07789         } else {
07790                 t = new Transform;
07791         }
07792 
07793 
07794         int mz = 2*(int)c+1;
07795         if ( nz < mz ) mz = nz;
07796         int my = 2*(int)b+1;
07797         if ( ny < my ) my = ny;
07798         int mx = 2*(int)a+1;
07799         if ( nx < mx ) mx = nx;
07800 
07801 
07802         float ai = 1/(a*a);
07803         float bi = 1/(b*b);
07804         float ci = 1/(c*c);
07805 
07806         Vec3f origin(nx/2,ny/2,nz/2);
07807 
07808         float x2, y2, z2, r;
07809         int xl, yl, zl;
07810         for (int k = 0; k < mz; ++k) {
07811                 for (int j = 0; j < my; ++j) {
07812                         for (int i = 0; i < mx; ++i) {
07813                                 x2 = (float)(i - mx/2);
07814                                 y2 = (float)(j - my/2);
07815                                 z2 = (float)(k - mz/2);
07816                                 r = (x2*x2)*ai + (y2*y2)*bi + (z2*z2)*ci;
07817                                 if (r <= 1) {
07818 
07819                                         if (t != 0) {
07820                                                 Vec3f v(x2,y2,z2);
07821                                                 v = (*t)*v;
07822                                                 v += origin;
07823 
07824                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07825                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07826 
07827                                                 for( int kk = -1; kk <= 1; ++kk)
07828                                                         for( int jj = -1; jj <= 1; ++jj)
07829                                                                 for( int ii = -1; ii <= 1; ++ii) {
07830                                                                         xl = (int)v[0]+ii;
07831                                                                         yl = (int)v[1]+jj;
07832                                                                         zl = (int)v[2]+kk;
07833                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07834                                                                                 image->set_value_at(xl,yl,zl,fill);
07835                                                                 }
07836                                         } else {
07837                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07838                                         }
07839                                 }
07840                         }
07841                 }
07842         }
07843 
07844         delete t;
07845 
07846         image->update();
07847 }


Member Data Documentation

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

Definition at line 6387 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Thu Nov 17 12:46:59 2011 for EMAN2 by  doxygen 1.4.7