EMAN2
Public Member Functions | Static Public Member Functions | Static Public Attributes
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:
xwidthinner equatorial radii along x axes
ywidthinner equatorial radii along y axes
zwidthinner polar radius
aoutter equatorial radii along x axes
boutter equatorial radii along y axes
coutter polar radius
widthspecify the width or specify each width explicitly - xwidth, ywidth, zwidth
transformOptionally transform the ellipse
fillvalue you want to fill in hollow ellipse, default to 1.0

Definition at line 6859 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 6869 of file processor.h.

                        {
                                return "Insert a hollow ellipse into the image.";
                        }
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 6864 of file processor.h.

References NAME.

                        {
                                return NAME;
                        }
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 6879 of file processor.h.

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

                        {
                                TypeDict d;
                                d.put("xwidth", EMObject::FLOAT, "inner equatorial radii along x axes");
                                d.put("ywidth", EMObject::FLOAT, "inner equatorial radii along y axes");
                                d.put("zwidth", EMObject::FLOAT, "inner polar radius");
                                d.put("a", EMObject::FLOAT, "outter equatorial radii along x axes");
                                d.put("b", EMObject::FLOAT, "outter equatorial radii along y axes");
                                d.put("c", EMObject::FLOAT, "outter polar radius");
                                d.put("width",EMObject::FLOAT, "width - specify the width or specify each width explicitly - xwidth, ywidth, zwidth");
                                d.put("transform", EMObject::TRANSFORM, "Optionally transform the ellipse");
                                d.put("fill", EMObject::FLOAT, "value you want to fill in hollow ellipse, default to 1.0");
                                return d;
                        }
static Processor* EMAN::TestImageHollowEllipse::NEW ( ) [inline, static]

Definition at line 6874 of file processor.h.

                        {
                                return new TestImageHollowEllipse();
                        }
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:
imageThe image to be processed.

Implements EMAN::Processor.

Definition at line 8284 of file processor.cpp.

References EMAN::Dict::has_key(), nx, ny, EMAN::Processor::params, EMAN::Dict::set_default(), EMAN::EMData::set_value_at(), t, EMAN::EMData::update(), and v.

{
        preprocess(image);

        float width = params.set_default("width",2.0f);

        float a2 = params.set_default("a",nx/2.0f-1.0f);
        float b2 = params.set_default("b",ny/2.0f-1.0f);
        float c2 = params.set_default("c",nz/2.0f-1.0f);

        float a1 = params.set_default("xwidth",a2-width);
        float b1 = params.set_default("ywidth",b2-width);
        float c1 = params.set_default("zwidth",c2-width);

        float fill = params.set_default("fill",1.0f);
        Transform* t;
        if (params.has_key("transform")) {
                t = params["transform"];
        } else {
                t = new Transform;
        }


        int mz = 2*(int)c2+1;
        if ( nz < mz ) mz = nz;
        int my = 2*(int)b2+1;
        if ( ny < my ) my = ny;
        int mx = 2*(int)a2+1;
        if ( nx < mx ) mx = nx;

        float ai1 = 1/(a1*a1);
        float bi1 = 1/(b1*b1);
        float ci1 = 1/(c1*c1);

        float ai2 = 1/(a2*a2);
        float bi2 = 1/(b2*b2);
        float ci2 = 1/(c2*c2);

        Vec3f origin(nx/2,ny/2,nz/2);

        float x2, y2, z2, r1,r2;
        int xl, yl, zl;
        for (int k = 0; k < mz; ++k) {
                for (int j = 0; j < my; ++j) {
                        for (int i = 0; i < mx; ++i) {
                                x2 = (float)(i - mx/2);
                                y2 = (float)(j - my/2);
                                z2 = (float)(k - mz/2);
                                r1 = (x2*x2)*ai1 + (y2*y2)*bi1 + (z2*z2)*ci1;
                                r2 = (x2*x2)*ai2 + (y2*y2)*bi2 + (z2*z2)*ci2;
                                if (r2 <= 1 && r1 >= 1) {

                                        if (t != 0) {
                                                Vec3f v(x2,y2,z2);
                                                v = (*t)*v;
                                                v += origin;

                                                // THIS ISN'T THE BEST STRATEGY BUT IT'S A STOP GAP. A FLOOD FILL IS PROBABLY BETTER
                                                // I fill in 3x3 cubes to make sure there are no gaps...

                                                for( int kk = -1; kk <= 1; ++kk)
                                                        for( int jj = -1; jj <= 1; ++jj)
                                                                for( int ii = -1; ii <= 1; ++ii) {
                                                                        xl = (int)v[0]+ii;
                                                                        yl = (int)v[1]+jj;
                                                                        zl = (int)v[2]+kk;
                                                                        if (xl >= 0 && xl < nx && yl >= 0 && yl < ny && zl >= 0 && zl < nz)
                                                                                image->set_value_at(xl,yl,zl,1.0);
                                                                }
                                        } else {
                                                image->set_value_at((int)x2+nx/2,(int)y2+ny/2,(int)z2+nz/2,fill);
                                        }
                                }
                        }
                }
        }

        delete t;

        image->update();
}

Member Data Documentation

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

Definition at line 6894 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files: