Inertia Matrix computer Computes the Inertia Matrix for a 3-D volume. More...

#include <analyzer.h>

Inheritance diagram for EMAN::InertiaMatrixAnalyzer:

Collaboration diagram for EMAN::InertiaMatrixAnalyzer:

[legend]

Public Member Functions
	InertiaMatrixAnalyzer ()

virtual int	insert_image (EMData *image)
	insert a image to the list of input images More...

virtual vector< EMData * >	analyze ()
	main function for Analyzer, analyze input images and create output images More...

string	get_name () const
	Get the Analyzer's name. More...

string	get_desc () const
	Get the Analyzer's description. More...

TypeDict	get_param_types () const
	Get Analyzer parameter information in a dictionary. More...

Public Member Functions inherited from EMAN::Analyzer
	Analyzer ()

virtual	~Analyzer ()

virtual int	insert_images_list (vector< EMData * > image_list)
	insert a list of images to the list of input images More...

virtual void	set_params (const Dict &new_params)
	Set the Analyzer parameters using a key/value dictionary. More...

virtual Dict	get_params () const
	Get the Reconstructor's parameters in a key/value dictionary. More...

Static Public Member Functions
static Analyzer *	NEW ()

Static Public Attributes
static const string	NAME = "inertiamatrix"

Protected Attributes
int	verbose

vector< EMData * >	ret

Protected Attributes inherited from EMAN::Analyzer
Dict	params

vector< EMData * >	images

Detailed Description

Inertia Matrix computer Computes the Inertia Matrix for a 3-D volume.

Author: Steve Ludtke

Date: 12/29/2013

Parameters

verbose Display progress if set, more detail with larger numbers

Definition at line 127 of file analyzer.h.

Constructor & Destructor Documentation

◆ InertiaMatrixAnalyzer()

EMAN::InertiaMatrixAnalyzer::InertiaMatrixAnalyzer ( )

inline

Definition at line 130 of file analyzer.h.

130: verbose(0) {}

EMAN::InertiaMatrixAnalyzer::verbose

int verbose

Definition: analyzer.h:165

Referenced by NEW().

Member Function Documentation

◆ analyze()

vector< EMData * > InertiaMatrixAnalyzer::analyze ( )

virtual

main function for Analyzer, analyze input images and create output images

Returns: vector<EMData *> result os images analysis

Implements EMAN::Analyzer.

Definition at line 78 of file analyzer.cpp.

                                                {
        int verbose = params.set_default("verbose",0);
        EMData *mx = new EMData(3,3);   // result is a 3x3 matrix
        mx->to_zero();
        ret.push_back(mx);
 
        if (images.size()!=1) throw ImageDimensionException("Inertia matrix computation accepts only a single volume as input");
        int nx=images[0]->get_xsize();
        int ny=images[0]->get_ysize();
        int nz=images[0]->get_zsize();
        if (nz==1 || ny==1 || nz==1) throw ImageDimensionException("Map must be 3-D");
 
        if (verbose>0) printf("Inertia volume size: %d %d %d\n",nx,ny,nz);
 
        for (int z=0; z<nz; z++) {
                for (int y=0; y<ny; y++) {
                        for (int x=0; x<nx; x++) {
                                int xx=x-nx/2;
                                int yy=y-ny/2;
                                int zz=z-nz/2;
                                float v=images[0]->get_value_at(x,y,z);
                                mx->set_value_at(0,0,mx->get_value_at(0,0)+v*(yy*yy+zz*zz));
                                mx->set_value_at(0,1,mx->get_value_at(0,1)+v*(-xx*yy));
                                mx->set_value_at(0,2,mx->get_value_at(0,2)+v*(-xx*zz));
                                mx->set_value_at(1,0,mx->get_value_at(1,0)+v*(-xx*yy));
                                mx->set_value_at(1,1,mx->get_value_at(1,1)+v*(zz*zz+xx*xx));
                                mx->set_value_at(1,2,mx->get_value_at(1,2)+v*(-yy*zz));
                                mx->set_value_at(2,0,mx->get_value_at(2,0)+v*(-xx*zz));
                                mx->set_value_at(2,1,mx->get_value_at(2,1)+v*(-yy*zz));
                                mx->set_value_at(2,2,mx->get_value_at(2,2)+v*(xx*xx+yy*yy));
                        }
                }
        }
        mx->mult(1.0f/(nx*ny*nz));
 
        if (verbose>0) {
                printf("%1.3g\t%1.3g\t%1.3g\n",mx->get_value_at(0,0),mx->get_value_at(1,0),mx->get_value_at(2,0));
                printf("%1.3g\t%1.3g\t%1.3g\n",mx->get_value_at(0,1),mx->get_value_at(1,1),mx->get_value_at(2,1));
                printf("%1.3g\t%1.3g\t%1.3g\n",mx->get_value_at(0,2),mx->get_value_at(1,2),mx->get_value_at(2,2));
        }
 
        return ret;
}

References ImageDimensionException, EMAN::Analyzer::images, EMAN::Analyzer::params, ret, EMAN::Dict::set_default(), verbose, x, and y.

◆ get_desc()

string EMAN::InertiaMatrixAnalyzer::get_desc ( ) const

inlinevirtual

Get the Analyzer's description.

Returns: The Analyzer's description.

Implements EMAN::Analyzer.

Definition at line 145 of file analyzer.h.

                {
                        return "Compute Inertia matrix for a volume";
                }

◆ get_name()

string EMAN::InertiaMatrixAnalyzer::get_name ( ) const

inlinevirtual

Get the Analyzer's name.

Each Analyzer is identified by a unique name.

Returns: The Analyzer's name.

Implements EMAN::Analyzer.

Definition at line 140 of file analyzer.h.

                {
                        return NAME;
                }

References NAME.

◆ get_param_types()

TypeDict EMAN::InertiaMatrixAnalyzer::get_param_types ( ) const

inlinevirtual

Get Analyzer 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.

Implements EMAN::Analyzer.

Definition at line 155 of file analyzer.h.

                {
                        TypeDict d;
                        d.put("verbose", EMObject::INT, "Display progress if set, more detail with larger numbers");
                        return d;
                }

References EMAN::EMObject::INT, and EMAN::TypeDict::put().

◆ insert_image()

virtual int EMAN::InertiaMatrixAnalyzer::insert_image ( EMData * image )

inlinevirtual

insert a image to the list of input images

Parameters

image

Returns: int 0 for success, <0 for fail

Implements EMAN::Analyzer.

Definition at line 132 of file analyzer.h.

                                                        {
                        images.push_back(image);
                        if (images.size()>1) { printf("InertiaMatrixAnalyzer takes only a single image\n"); return 1; }
                        return 0;
                }

References EMAN::Analyzer::images.

◆ NEW()

static Analyzer * EMAN::InertiaMatrixAnalyzer::NEW ( )

inlinestatic

Definition at line 150 of file analyzer.h.

                {
                        return new InertiaMatrixAnalyzer();
                }

References InertiaMatrixAnalyzer().

Member Data Documentation

◆ NAME

const string EMAN::InertiaMatrixAnalyzer::NAME = "inertiamatrix"

static

Definition at line 162 of file analyzer.h.

Referenced by get_name().

◆ ret

vector<EMData *> EMAN::InertiaMatrixAnalyzer::ret

protected

Definition at line 166 of file analyzer.h.

Referenced by analyze().

◆ verbose

int EMAN::InertiaMatrixAnalyzer::verbose

protected

Definition at line 165 of file analyzer.h.

Referenced by analyze().

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

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ InertiaMatrixAnalyzer()

Member Function Documentation

◆ analyze()

◆ get_desc()

◆ get_name()

◆ get_param_types()

◆ insert_image()

◆ NEW()

Member Data Documentation

◆ NAME

◆ ret

◆ verbose