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

06459                         {
06460                                 return "Insert a hollow ellipse into the image.";
06461                         }

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

06454                         {
06455                                 return NAME;
06456                         }

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

References EMAN::TypeDict::put().

06469                         {
06470                                 TypeDict d;
06471                                 d.put("xwidth", EMObject::FLOAT, "inner equatorial radii along x axes");
06472                                 d.put("ywidth", EMObject::FLOAT, "inner equatorial radii along y axes");
06473                                 d.put("zwidth", EMObject::FLOAT, "inner polar radius");
06474                                 d.put("a", EMObject::FLOAT, "outter equatorial radii along x axes");
06475                                 d.put("b", EMObject::FLOAT, "outter equatorial radii along y axes");
06476                                 d.put("c", EMObject::FLOAT, "outter polar radius");
06477                                 d.put("width",EMObject::FLOAT, "width - specify the width or specify each width explicitly - xwidth, ywidth, zwidth");
06478                                 d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06479                                 d.put("fill", EMObject::FLOAT, "value you want to fill in hollow ellipse, default to 1.0");
06480                                 return d;
06481                         }

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

Definition at line 6463 of file processor.h.

06464                         {
06465                                 return new TestImageHollowEllipse();
06466                         }

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

07782 {
07783         preprocess(image);
07784 
07785         float width = params.set_default("width",2.0f);
07786 
07787         float a2 = params.set_default("a",nx/2.0f-1.0f);
07788         float b2 = params.set_default("b",ny/2.0f-1.0f);
07789         float c2 = params.set_default("c",nz/2.0f-1.0f);
07790 
07791         float a1 = params.set_default("xwidth",a2-width);
07792         float b1 = params.set_default("ywidth",b2-width);
07793         float c1 = params.set_default("zwidth",c2-width);
07794 
07795         float fill = params.set_default("fill",1.0f);
07796         Transform* t;
07797         if (params.has_key("transform")) {
07798                 t = params["transform"];
07799         } else {
07800                 t = new Transform;
07801         }
07802 
07803 
07804         int mz = 2*(int)c2+1;
07805         if ( nz < mz ) mz = nz;
07806         int my = 2*(int)b2+1;
07807         if ( ny < my ) my = ny;
07808         int mx = 2*(int)a2+1;
07809         if ( nx < mx ) mx = nx;
07810 
07811         float ai1 = 1/(a1*a1);
07812         float bi1 = 1/(b1*b1);
07813         float ci1 = 1/(c1*c1);
07814 
07815         float ai2 = 1/(a2*a2);
07816         float bi2 = 1/(b2*b2);
07817         float ci2 = 1/(c2*c2);
07818 
07819         Vec3f origin(nx/2,ny/2,nz/2);
07820 
07821         float x2, y2, z2, r1,r2;
07822         int xl, yl, zl;
07823         for (int k = 0; k < mz; ++k) {
07824                 for (int j = 0; j < my; ++j) {
07825                         for (int i = 0; i < mx; ++i) {
07826                                 x2 = (float)(i - mx/2);
07827                                 y2 = (float)(j - my/2);
07828                                 z2 = (float)(k - mz/2);
07829                                 r1 = (x2*x2)*ai1 + (y2*y2)*bi1 + (z2*z2)*ci1;
07830                                 r2 = (x2*x2)*ai2 + (y2*y2)*bi2 + (z2*z2)*ci2;
07831                                 if (r2 <= 1 && r1 >= 1) {
07832 
07833                                         if (t != 0) {
07834                                                 Vec3f v(x2,y2,z2);
07835                                                 v = (*t)*v;
07836                                                 v += origin;
07837 
07838                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07839                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07840 
07841                                                 for( int kk = -1; kk <= 1; ++kk)
07842                                                         for( int jj = -1; jj <= 1; ++jj)
07843                                                                 for( int ii = -1; ii <= 1; ++ii) {
07844                                                                         xl = (int)v[0]+ii;
07845                                                                         yl = (int)v[1]+jj;
07846                                                                         zl = (int)v[2]+kk;
07847                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07848                                                                                 image->set_value_at(xl,yl,zl,1.0);
07849                                                                 }
07850                                         } else {
07851                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07852                                         }
07853                                 }
07854                         }
07855                 }
07856         }
07857 
07858         delete t;
07859 
07860         image->update();
07861 }


Member Data Documentation

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

Definition at line 205 of file processor.cpp.


The documentation for this class was generated from the following files:
Generated on Mon Aug 13 13:42:29 2012 for EMAN2 by  doxygen 1.3.9.1