Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members

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

ProcessorNEW ()

Static Public Attributes

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

06202                         {
06203                                 return "Insert a hollow ellipse into the image.";
06204                         }

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

06197                         {
06198                                 return NAME;
06199                         }

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

References EMAN::TypeDict::put().

06212                         {
06213                                 TypeDict d;
06214                                 d.put("xwidth", EMObject::FLOAT, "inner equatorial radii along x axes");
06215                                 d.put("ywidth", EMObject::FLOAT, "inner equatorial radii along y axes");
06216                                 d.put("zwidth", EMObject::FLOAT, "inner polar radius");
06217                                 d.put("a", EMObject::FLOAT, "outter equatorial radii along x axes");
06218                                 d.put("b", EMObject::FLOAT, "outter equatorial radii along y axes");
06219                                 d.put("c", EMObject::FLOAT, "outter polar radius");
06220                                 d.put("width",EMObject::FLOAT, "width - specify the width or specify each width explicitly - xwidth, ywidth, zwidth");
06221                                 d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06222                                 d.put("fill", EMObject::FLOAT, "value you want to fill in hollow ellipse, default to 1.0");
06223                                 return d;
06224                         }

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

Definition at line 6206 of file processor.h.

06207                         {
06208                                 return new TestImageHollowEllipse();
06209                         }

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 7511 of file processor.cpp.

References EMAN::Dict::has_key(), nx, ny, EMAN::TestImageProcessor::preprocess(), EMAN::Dict::set_default(), EMAN::EMData::set_value_at(), t, EMAN::EMData::update(), v, and EMAN::Vec3f.

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


Member Data Documentation

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

Definition at line 200 of file processor.cpp.


The documentation for this class was generated from the following files:
Generated on Thu Dec 9 13:48:24 2010 for EMAN2 by  doxygen 1.3.9.1