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

06498                         {
06499                                 return "Insert a hollow ellipse into the image.";
06500                         }

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

References NAME.

06493                         {
06494                                 return NAME;
06495                         }

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

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

06508                         {
06509                                 TypeDict d;
06510                                 d.put("xwidth", EMObject::FLOAT, "inner equatorial radii along x axes");
06511                                 d.put("ywidth", EMObject::FLOAT, "inner equatorial radii along y axes");
06512                                 d.put("zwidth", EMObject::FLOAT, "inner polar radius");
06513                                 d.put("a", EMObject::FLOAT, "outter equatorial radii along x axes");
06514                                 d.put("b", EMObject::FLOAT, "outter equatorial radii along y axes");
06515                                 d.put("c", EMObject::FLOAT, "outter polar radius");
06516                                 d.put("width",EMObject::FLOAT, "width - specify the width or specify each width explicitly - xwidth, ywidth, zwidth");
06517                                 d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06518                                 d.put("fill", EMObject::FLOAT, "value you want to fill in hollow ellipse, default to 1.0");
06519                                 return d;
06520                         }

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

Definition at line 6502 of file processor.h.

06503                         {
06504                                 return new TestImageHollowEllipse();
06505                         }

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

07847 {
07848         preprocess(image);
07849 
07850         float width = params.set_default("width",2.0f);
07851 
07852         float a2 = params.set_default("a",nx/2.0f-1.0f);
07853         float b2 = params.set_default("b",ny/2.0f-1.0f);
07854         float c2 = params.set_default("c",nz/2.0f-1.0f);
07855 
07856         float a1 = params.set_default("xwidth",a2-width);
07857         float b1 = params.set_default("ywidth",b2-width);
07858         float c1 = params.set_default("zwidth",c2-width);
07859 
07860         float fill = params.set_default("fill",1.0f);
07861         Transform* t;
07862         if (params.has_key("transform")) {
07863                 t = params["transform"];
07864         } else {
07865                 t = new Transform;
07866         }
07867 
07868 
07869         int mz = 2*(int)c2+1;
07870         if ( nz < mz ) mz = nz;
07871         int my = 2*(int)b2+1;
07872         if ( ny < my ) my = ny;
07873         int mx = 2*(int)a2+1;
07874         if ( nx < mx ) mx = nx;
07875 
07876         float ai1 = 1/(a1*a1);
07877         float bi1 = 1/(b1*b1);
07878         float ci1 = 1/(c1*c1);
07879 
07880         float ai2 = 1/(a2*a2);
07881         float bi2 = 1/(b2*b2);
07882         float ci2 = 1/(c2*c2);
07883 
07884         Vec3f origin(nx/2,ny/2,nz/2);
07885 
07886         float x2, y2, z2, r1,r2;
07887         int xl, yl, zl;
07888         for (int k = 0; k < mz; ++k) {
07889                 for (int j = 0; j < my; ++j) {
07890                         for (int i = 0; i < mx; ++i) {
07891                                 x2 = (float)(i - mx/2);
07892                                 y2 = (float)(j - my/2);
07893                                 z2 = (float)(k - mz/2);
07894                                 r1 = (x2*x2)*ai1 + (y2*y2)*bi1 + (z2*z2)*ci1;
07895                                 r2 = (x2*x2)*ai2 + (y2*y2)*bi2 + (z2*z2)*ci2;
07896                                 if (r2 <= 1 && r1 >= 1) {
07897 
07898                                         if (t != 0) {
07899                                                 Vec3f v(x2,y2,z2);
07900                                                 v = (*t)*v;
07901                                                 v += origin;
07902 
07903                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07904                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07905 
07906                                                 for( int kk = -1; kk <= 1; ++kk)
07907                                                         for( int jj = -1; jj <= 1; ++jj)
07908                                                                 for( int ii = -1; ii <= 1; ++ii) {
07909                                                                         xl = (int)v[0]+ii;
07910                                                                         yl = (int)v[1]+jj;
07911                                                                         zl = (int)v[2]+kk;
07912                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07913                                                                                 image->set_value_at(xl,yl,zl,1.0);
07914                                                                 }
07915                                         } else {
07916                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07917                                         }
07918                                 }
07919                         }
07920                 }
07921         }
07922 
07923         delete t;
07924 
07925         image->update();
07926 }


Member Data Documentation

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

Definition at line 6522 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Tue Jun 11 12:44:44 2013 for EMAN2 by  doxygen 1.4.7