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

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

06415                         {
06416                                 return "Insert a hollow ellipse into the image.";
06417                         }

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

References NAME.

06410                         {
06411                                 return NAME;
06412                         }

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

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

06425                         {
06426                                 TypeDict d;
06427                                 d.put("xwidth", EMObject::FLOAT, "inner equatorial radii along x axes");
06428                                 d.put("ywidth", EMObject::FLOAT, "inner equatorial radii along y axes");
06429                                 d.put("zwidth", EMObject::FLOAT, "inner polar radius");
06430                                 d.put("a", EMObject::FLOAT, "outter equatorial radii along x axes");
06431                                 d.put("b", EMObject::FLOAT, "outter equatorial radii along y axes");
06432                                 d.put("c", EMObject::FLOAT, "outter polar radius");
06433                                 d.put("width",EMObject::FLOAT, "width - specify the width or specify each width explicitly - xwidth, ywidth, zwidth");
06434                                 d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06435                                 d.put("fill", EMObject::FLOAT, "value you want to fill in hollow ellipse, default to 1.0");
06436                                 return d;
06437                         }

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

Definition at line 6419 of file processor.h.

06420                         {
06421                                 return new TestImageHollowEllipse();
06422                         }

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

07614 {
07615         preprocess(image);
07616 
07617         float width = params.set_default("width",2.0f);
07618 
07619         float a2 = params.set_default("a",nx/2.0f-1.0f);
07620         float b2 = params.set_default("b",ny/2.0f-1.0f);
07621         float c2 = params.set_default("c",nz/2.0f-1.0f);
07622 
07623         float a1 = params.set_default("xwidth",a2-width);
07624         float b1 = params.set_default("ywidth",b2-width);
07625         float c1 = params.set_default("zwidth",c2-width);
07626 
07627         float fill = params.set_default("fill",1.0f);
07628         Transform* t;
07629         if (params.has_key("transform")) {
07630                 t = params["transform"];
07631         } else {
07632                 t = new Transform;
07633         }
07634 
07635 
07636         int mz = 2*(int)c2+1;
07637         if ( nz < mz ) mz = nz;
07638         int my = 2*(int)b2+1;
07639         if ( ny < my ) my = ny;
07640         int mx = 2*(int)a2+1;
07641         if ( nx < mx ) mx = nx;
07642 
07643         float ai1 = 1/(a1*a1);
07644         float bi1 = 1/(b1*b1);
07645         float ci1 = 1/(c1*c1);
07646 
07647         float ai2 = 1/(a2*a2);
07648         float bi2 = 1/(b2*b2);
07649         float ci2 = 1/(c2*c2);
07650 
07651         Vec3f origin(nx/2,ny/2,nz/2);
07652 
07653         float x2, y2, z2, r1,r2;
07654         int xl, yl, zl;
07655         for (int k = 0; k < mz; ++k) {
07656                 for (int j = 0; j < my; ++j) {
07657                         for (int i = 0; i < mx; ++i) {
07658                                 x2 = (float)(i - mx/2);
07659                                 y2 = (float)(j - my/2);
07660                                 z2 = (float)(k - mz/2);
07661                                 r1 = (x2*x2)*ai1 + (y2*y2)*bi1 + (z2*z2)*ci1;
07662                                 r2 = (x2*x2)*ai2 + (y2*y2)*bi2 + (z2*z2)*ci2;
07663                                 if (r2 <= 1 && r1 >= 1) {
07664 
07665                                         if (t != 0) {
07666                                                 Vec3f v(x2,y2,z2);
07667                                                 v = (*t)*v;
07668                                                 v += origin;
07669 
07670                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07671                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07672 
07673                                                 for( int kk = -1; kk <= 1; ++kk)
07674                                                         for( int jj = -1; jj <= 1; ++jj)
07675                                                                 for( int ii = -1; ii <= 1; ++ii) {
07676                                                                         xl = (int)v[0]+ii;
07677                                                                         yl = (int)v[1]+jj;
07678                                                                         zl = (int)v[2]+kk;
07679                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07680                                                                                 image->set_value_at(xl,yl,zl,1.0);
07681                                                                 }
07682                                         } else {
07683                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07684                                         }
07685                                 }
07686                         }
07687                 }
07688         }
07689 
07690         delete t;
07691 
07692         image->update();
07693 }


Member Data Documentation

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

Definition at line 6439 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Mon Jul 19 13:07:26 2010 for EMAN2 by  doxygen 1.4.4