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

06457                         {
06458                                 return "Insert a hollow ellipse into the image.";
06459                         }

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

References NAME.

06452                         {
06453                                 return NAME;
06454                         }

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

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

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

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

Definition at line 6461 of file processor.h.

06462                         {
06463                                 return new TestImageHollowEllipse();
06464                         }

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

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


Member Data Documentation

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

Definition at line 6481 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Thu May 3 10:10:55 2012 for EMAN2 by  doxygen 1.4.7