#include <processor.h>
Inheritance diagram for EMAN::TestImageEllipse:
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 Processor * | NEW () |
Static Public Attributes | |
static const string | NAME = "testimage.ellipsoid" |
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.
virtual string EMAN::TestImageEllipse::get_desc | ( | ) | const [inline, virtual] |
Get the descrition of this specific processor.
This function must be overwritten by a subclass.
Implements EMAN::Processor.
Definition at line 6452 of file processor.h.
virtual string EMAN::TestImageEllipse::get_name | ( | ) | const [inline, virtual] |
Get the processor's name.
Each processor is identified by a unique name.
Implements EMAN::Processor.
Definition at line 6447 of file processor.h.
References NAME.
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.
Reimplemented from EMAN::Processor.
Definition at line 6462 of file processor.h.
References EMAN::EMObject::FLOAT, EMAN::TypeDict::put(), and EMAN::EMObject::TRANSFORM.
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 }
static Processor* EMAN::TestImageEllipse::NEW | ( | ) | [inline, static] |
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.
image | The image to be processed. |
Implements EMAN::Processor.
Definition at line 7928 of file processor.cpp.
References b, bi, 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.
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 }
const string TestImageEllipse::NAME = "testimage.ellipsoid" [static] |