#include <processor.h>
Inheritance diagram for EMAN::TestImageHollowEllipse:
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.hollow" |
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.
virtual string EMAN::TestImageHollowEllipse::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 6286 of file processor.h.
virtual string EMAN::TestImageHollowEllipse::get_name | ( | ) | const [inline, virtual] |
Get the processor's name.
Each processor is identified by a unique 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.
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.
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 }
const string TestImageHollowEllipse::NAME = "testimage.ellipsoid.hollow" [static] |