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

06287                         {
06288                                 return "Insert a hollow ellipse into the image.";
06289                         }

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

References NAME.

06282                         {
06283                                 return NAME;
06284                         }

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

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

06297                         {
06298                                 TypeDict d;
06299                                 d.put("xwidth", EMObject::FLOAT, "inner equatorial radii along x axes");
06300                                 d.put("ywidth", EMObject::FLOAT, "inner equatorial radii along y axes");
06301                                 d.put("zwidth", EMObject::FLOAT, "inner polar radius");
06302                                 d.put("a", EMObject::FLOAT, "outter equatorial radii along x axes");
06303                                 d.put("b", EMObject::FLOAT, "outter equatorial radii along y axes");
06304                                 d.put("c", EMObject::FLOAT, "outter polar radius");
06305                                 d.put("width",EMObject::FLOAT, "width - specify the width or specify each width explicitly - xwidth, ywidth, zwidth");
06306                                 d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06307                                 d.put("fill", EMObject::FLOAT, "value you want to fill in hollow ellipse, default to 1.0");
06308                                 return d;
06309                         }

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

Definition at line 6291 of file processor.h.

06292                         {
06293                                 return new TestImageHollowEllipse();
06294                         }

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

07579 {
07580         preprocess(image);
07581 
07582         float width = params.set_default("width",2.0f);
07583 
07584         float a2 = params.set_default("a",nx/2.0f-1.0f);
07585         float b2 = params.set_default("b",ny/2.0f-1.0f);
07586         float c2 = params.set_default("c",nz/2.0f-1.0f);
07587 
07588         float a1 = params.set_default("xwidth",a2-width);
07589         float b1 = params.set_default("ywidth",b2-width);
07590         float c1 = params.set_default("zwidth",c2-width);
07591 
07592         float fill = params.set_default("fill",1.0f);
07593         Transform* t;
07594         if (params.has_key("transform")) {
07595                 t = params["transform"];
07596         } else {
07597                 t = new Transform;
07598         }
07599 
07600 
07601         int mz = 2*(int)c2+1;
07602         if ( nz < mz ) mz = nz;
07603         int my = 2*(int)b2+1;
07604         if ( ny < my ) my = ny;
07605         int mx = 2*(int)a2+1;
07606         if ( nx < mx ) mx = nx;
07607 
07608         float ai1 = 1/(a1*a1);
07609         float bi1 = 1/(b1*b1);
07610         float ci1 = 1/(c1*c1);
07611 
07612         float ai2 = 1/(a2*a2);
07613         float bi2 = 1/(b2*b2);
07614         float ci2 = 1/(c2*c2);
07615 
07616         Vec3f origin(nx/2,ny/2,nz/2);
07617 
07618         float x2, y2, z2, r1,r2;
07619         int xl, yl, zl;
07620         for (int k = 0; k < mz; ++k) {
07621                 for (int j = 0; j < my; ++j) {
07622                         for (int i = 0; i < mx; ++i) {
07623                                 x2 = (float)(i - mx/2);
07624                                 y2 = (float)(j - my/2);
07625                                 z2 = (float)(k - mz/2);
07626                                 r1 = (x2*x2)*ai1 + (y2*y2)*bi1 + (z2*z2)*ci1;
07627                                 r2 = (x2*x2)*ai2 + (y2*y2)*bi2 + (z2*z2)*ci2;
07628                                 if (r2 <= 1 && r1 >= 1) {
07629 
07630                                         if (t != 0) {
07631                                                 Vec3f v(x2,y2,z2);
07632                                                 v = (*t)*v;
07633                                                 v += origin;
07634 
07635                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07636                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07637 
07638                                                 for( int kk = -1; kk <= 1; ++kk)
07639                                                         for( int jj = -1; jj <= 1; ++jj)
07640                                                                 for( int ii = -1; ii <= 1; ++ii) {
07641                                                                         xl = (int)v[0]+ii;
07642                                                                         yl = (int)v[1]+jj;
07643                                                                         zl = (int)v[2]+kk;
07644                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07645                                                                                 image->set_value_at(xl,yl,zl,1.0);
07646                                                                 }
07647                                         } else {
07648                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07649                                         }
07650                                 }
07651                         }
07652                 }
07653         }
07654 
07655         delete t;
07656 
07657         image->update();
07658 }


Member Data Documentation

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

Definition at line 6311 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Mon Mar 7 18:10:47 2011 for EMAN2 by  doxygen 1.4.7