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

06412                         {
06413                                 return "Insert a hollow ellipse into the image.";
06414                         }

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

References NAME.

06407                         {
06408                                 return NAME;
06409                         }

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

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

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

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

Definition at line 6416 of file processor.h.

06417                         {
06418                                 return new TestImageHollowEllipse();
06419                         }

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

07695 {
07696         preprocess(image);
07697 
07698         float width = params.set_default("width",2.0f);
07699 
07700         float a2 = params.set_default("a",nx/2.0f-1.0f);
07701         float b2 = params.set_default("b",ny/2.0f-1.0f);
07702         float c2 = params.set_default("c",nz/2.0f-1.0f);
07703 
07704         float a1 = params.set_default("xwidth",a2-width);
07705         float b1 = params.set_default("ywidth",b2-width);
07706         float c1 = params.set_default("zwidth",c2-width);
07707 
07708         float fill = params.set_default("fill",1.0f);
07709         Transform* t;
07710         if (params.has_key("transform")) {
07711                 t = params["transform"];
07712         } else {
07713                 t = new Transform;
07714         }
07715 
07716 
07717         int mz = 2*(int)c2+1;
07718         if ( nz < mz ) mz = nz;
07719         int my = 2*(int)b2+1;
07720         if ( ny < my ) my = ny;
07721         int mx = 2*(int)a2+1;
07722         if ( nx < mx ) mx = nx;
07723 
07724         float ai1 = 1/(a1*a1);
07725         float bi1 = 1/(b1*b1);
07726         float ci1 = 1/(c1*c1);
07727 
07728         float ai2 = 1/(a2*a2);
07729         float bi2 = 1/(b2*b2);
07730         float ci2 = 1/(c2*c2);
07731 
07732         Vec3f origin(nx/2,ny/2,nz/2);
07733 
07734         float x2, y2, z2, r1,r2;
07735         int xl, yl, zl;
07736         for (int k = 0; k < mz; ++k) {
07737                 for (int j = 0; j < my; ++j) {
07738                         for (int i = 0; i < mx; ++i) {
07739                                 x2 = (float)(i - mx/2);
07740                                 y2 = (float)(j - my/2);
07741                                 z2 = (float)(k - mz/2);
07742                                 r1 = (x2*x2)*ai1 + (y2*y2)*bi1 + (z2*z2)*ci1;
07743                                 r2 = (x2*x2)*ai2 + (y2*y2)*bi2 + (z2*z2)*ci2;
07744                                 if (r2 <= 1 && r1 >= 1) {
07745 
07746                                         if (t != 0) {
07747                                                 Vec3f v(x2,y2,z2);
07748                                                 v = (*t)*v;
07749                                                 v += origin;
07750 
07751                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07752                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07753 
07754                                                 for( int kk = -1; kk <= 1; ++kk)
07755                                                         for( int jj = -1; jj <= 1; ++jj)
07756                                                                 for( int ii = -1; ii <= 1; ++ii) {
07757                                                                         xl = (int)v[0]+ii;
07758                                                                         yl = (int)v[1]+jj;
07759                                                                         zl = (int)v[2]+kk;
07760                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07761                                                                                 image->set_value_at(xl,yl,zl,1.0);
07762                                                                 }
07763                                         } else {
07764                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07765                                         }
07766                                 }
07767                         }
07768                 }
07769         }
07770 
07771         delete t;
07772 
07773         image->update();
07774 }


Member Data Documentation

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

Definition at line 6436 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:47:00 2011 for EMAN2 by  doxygen 1.4.7