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

EMAN::WaveletProcessor Class Reference

Perform a Wavelet transform using GSL. More...

#include <processor.h>

Inheritance diagram for EMAN::WaveletProcessor:

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

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 TypeDict get_param_types () const
 Get processor parameter information in a dictionary.
virtual string get_desc () const
 Get the descrition of this specific processor.

Static Public Member Functions

ProcessorNEW ()

Static Public Attributes

const string NAME = "basis.wavelet"

Detailed Description

Perform a Wavelet transform using GSL.

Author:
Steve Ludtke <sludtke@bcm.edu>
Date:
10/15/2006
Parameters:
type "daub", "harr", or "bspl"
dir 1 for forward transform, -1 for inverse transform
ord for Daubechies (4,6,8,...,20), for Harr (2), for B-Splines (103, 105, 202, 204, 206, 208, 301, 303, 305 307, 309)

Definition at line 6719 of file processor.h.


Member Function Documentation

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

06744                 {
06745                         return "Computes the DWT (discrete wavelet transform) of an image in one of 3 possible bases";
06746                 }

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

Referenced by process_inplace().

06725                 {
06726                         return NAME;
06727                 }

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

References EMAN::TypeDict::put().

06735                 {
06736                         TypeDict d;
06737                         d.put("type", EMObject::STRING, "'daub', 'harr' or 'bspl'");
06738                         d.put("dir", EMObject::INT, "1 for forward transform, -1 for inverse transform");
06739                         d.put("ord", EMObject::INT, "Daubechies (4,6,8,...,20), for Harr (2), for B-Splines (103, 105, 202, 204, 206, 208, 301, 303, 305 307, 309)");
06740                         return d;
06741                 }

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

Definition at line 6729 of file processor.h.

06730                 {
06731                         return new WaveletProcessor();
06732                 }

void WaveletProcessor::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 8271 of file processor.cpp.

References get_name(), EMAN::EMData::get_value_at(), EMAN::EMData::get_xsize(), EMAN::EMData::get_ysize(), EMAN::EMData::get_zsize(), ImageDimensionException, InvalidStringException, EMAN::Util::IsPower2(), LOGERR, nx, ny, and EMAN::EMData::set_value_at_fast().

08272 {
08273         if (image->get_zsize() != 1) {
08274                         LOGERR("%s Processor doesn't support 3D", get_name().c_str());
08275                         throw ImageDimensionException("3D model not supported");
08276         }
08277 
08278         int i,nx,ny;
08279         const gsl_wavelet_type * T;
08280         nx=image->get_xsize();
08281         ny=image->get_ysize();
08282 
08283         if (nx != ny && ny!=1) throw ImageDimensionException("Wavelet transform only supports square images");
08284 //      float l=log((float)nx)/log(2.0f);
08285 //      if (l!=floor(l)) throw ImageDimensionException("Wavelet transform size must be power of 2");
08286         if( !Util::IsPower2(nx) )  throw ImageDimensionException("Wavelet transform size must be power of 2");
08287 
08288         // Unfortunately GSL works only on double() arrays
08289         // eventually we should put our own wavelet code in here
08290         // but this will work for now
08291         double *cpy = (double *)malloc(nx*ny*sizeof(double));
08292 
08293         for (i=0; i<nx*ny; i++) cpy[i]=image->get_value_at(i,0,0);
08294 
08295         string tp = (const char*)params["type"];
08296         if (tp=="daub") T=gsl_wavelet_daubechies;
08297         else if (tp=="harr") T=gsl_wavelet_haar;
08298         else if (tp=="bspl") T=gsl_wavelet_bspline;
08299         else throw InvalidStringException(tp,"Invalid wavelet name, 'daub', 'harr' or 'bspl'");
08300 
08301         int K=(int)params["ord"];
08302         gsl_wavelet_direction dir;
08303         if ((int)params["dir"]==1) dir=forward;
08304         else dir=backward;
08305 
08306         gsl_wavelet *w = gsl_wavelet_alloc(T, K);
08307         gsl_wavelet_workspace *work = gsl_wavelet_workspace_alloc(nx);
08308 
08309         if (ny==1) gsl_wavelet_transform (w,cpy, 1, nx, dir, work);
08310         else gsl_wavelet2d_transform (w, cpy, nx,nx,ny, dir, work);
08311 
08312         gsl_wavelet_workspace_free (work);
08313         gsl_wavelet_free (w);
08314 
08315         for (i=0; i<nx*ny; i++) image->set_value_at_fast(i,0,0,static_cast<float>(cpy[i]));
08316 
08317         free(cpy);
08318 }


Member Data Documentation

const string WaveletProcessor::NAME = "basis.wavelet" [static]
 

Definition at line 212 of file processor.cpp.


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