Main Page | Modules | Namespace List | Class Hierarchy | Alphabetical List | Class List | Directories | File List | Namespace Members | Class Members | File Members

EMAN::TestImageEllipse Class Reference

Generate an ellipse or ellipsoid image. More...

#include <processor.h>

Inheritance diagram for EMAN::TestImageEllipse:

Inheritance graph
[legend]
Collaboration diagram for EMAN::TestImageEllipse:

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"

Detailed Description

Generate an ellipse or ellipsoid image.

Parameters:
a equatorial radii along x axes
b equatorial radii along y axes
c polar radius
transform Optionally transform the ellipse
fill value you want to fill in ellipse, default to 1.0

Definition at line 6442 of file processor.h.


Member Function Documentation

virtual string EMAN::TestImageEllipse::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 6452 of file processor.h.

06453                 {
06454                         return "Insert an ellipse into the image.";
06455                 }

virtual string EMAN::TestImageEllipse::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 6447 of file processor.h.

06448                 {
06449                         return NAME;
06450                 }

virtual TypeDict EMAN::TestImageEllipse::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 6462 of file processor.h.

References EMAN::TypeDict::put().

06463                 {
06464                         TypeDict d;
06465                         d.put("a", EMObject::FLOAT, "equatorial radius along x axes (major semiaxes)");
06466                         d.put("b", EMObject::FLOAT, "equatorial radius along y axes (minor semiaxes)");
06467                         d.put("c", EMObject::FLOAT, "polar radius for ellipsoid (x^2/a^2+y^2/b^2+z^2/c^2=1)");
06468                         d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
06469                         d.put("fill", EMObject::FLOAT, "value you want to fill in ellipse, default to 1.0");
06470                         return d;
06471                 }

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

Definition at line 6457 of file processor.h.

06458                 {
06459                         return new TestImageEllipse();
06460                 }

void TestImageEllipse::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 7928 of file processor.cpp.

References b, bi, 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.

07929 {
07930         preprocess(image);
07931 
07932 
07933         float a = params.set_default("a",nx/2.0f-1.0f);
07934         float b = params.set_default("b",ny/2.0f-1.0f);
07935         float c = params.set_default("c",nz/2.0f-1.0f);
07936         float fill = params.set_default("fill",1.0f);
07937         //bool hollow = params.set_default("hollow",false);
07938         Transform* t;
07939         if (params.has_key("transform")) {
07940                 t = params["transform"];
07941         } else {
07942                 t = new Transform;
07943         }
07944 
07945 
07946         int mz = 2*(int)c+1;
07947         if ( nz < mz ) mz = nz;
07948         int my = 2*(int)b+1;
07949         if ( ny < my ) my = ny;
07950         int mx = 2*(int)a+1;
07951         if ( nx < mx ) mx = nx;
07952 
07953 
07954         float ai = 1/(a*a);
07955         float bi = 1/(b*b);
07956         float ci = 1/(c*c);
07957 
07958         Vec3f origin(nx/2,ny/2,nz/2);
07959 
07960         float x2, y2, z2, r;
07961         int xl, yl, zl;
07962         for (int k = 0; k < mz; ++k) {
07963                 for (int j = 0; j < my; ++j) {
07964                         for (int i = 0; i < mx; ++i) {
07965                                 x2 = (float)(i - mx/2);
07966                                 y2 = (float)(j - my/2);
07967                                 z2 = (float)(k - mz/2);
07968                                 r = (x2*x2)*ai + (y2*y2)*bi + (z2*z2)*ci;
07969                                 if (r <= 1) {
07970 
07971                                         if (t != 0) {
07972                                                 Vec3f v(x2,y2,z2);
07973                                                 v = (*t)*v;
07974                                                 v += origin;
07975 
07976                                                 // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
07977                                                 // I fill in 3x3 cubes to make sure there are no gaps...
07978 
07979                                                 for( int kk = -1; kk <= 1; ++kk)
07980                                                         for( int jj = -1; jj <= 1; ++jj)
07981                                                                 for( int ii = -1; ii <= 1; ++ii) {
07982                                                                         xl = (int)v[0]+ii;
07983                                                                         yl = (int)v[1]+jj;
07984                                                                         zl = (int)v[2]+kk;
07985                                                                         if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
07986                                                                                 image->set_value_at(xl,yl,zl,fill);
07987                                                                 }
07988                                         } else {
07989                                                 image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
07990                                         }
07991                                 }
07992                         }
07993                 }
07994         }
07995 
07996         delete t;
07997 
07998         image->update();
07999 }


Member Data Documentation

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

Definition at line 205 of file processor.cpp.


The documentation for this class was generated from the following files:
Generated on Tue Jun 11 13:42:47 2013 for EMAN2 by  doxygen 1.3.9.1