EMAN::Wiener2DFourierProcessor Class Reference

Wiener filter based on a Ctf object either in the image header. More...

#include <processor.h>

Inheritance diagram for EMAN::Wiener2DFourierProcessor:

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

Collaboration graph
[legend]
List of all members.

Public Member Functions

string get_name () const
 Get the processor's name.
virtual EMDataprocess (const EMData *const image)
 To proccess an image out-of-place.
void process_inplace (EMData *image)
 To process an image in-place.
void set_params (const Dict &new_params)
 Set the processor parameters using a key/value dictionary.
TypeDict get_param_types () const
 Get processor parameter information in a dictionary.
string get_desc () const
 Get the descrition of this specific processor.

Static Public Member Functions

static ProcessorNEW ()

Static Public Attributes

static const string NAME = "filter.wiener2d"

Protected Attributes

Ctfctf

Detailed Description

Wiener filter based on a Ctf object either in the image header.

Parameters:
ctf[in] A Ctf object to use
Author:
Steve Ludtke
Date:
2008/11/03

Definition at line 820 of file processor.h.


Member Function Documentation

string EMAN::Wiener2DFourierProcessor::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 851 of file processor.h.

00852                 {
00853                         return "Applies a 2-D Wiener filter to an image based on its Ctf parameters";
00854                 }

string EMAN::Wiener2DFourierProcessor::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 823 of file processor.h.

References NAME.

00824                 {
00825                         return NAME;
00826                 }

TypeDict EMAN::Wiener2DFourierProcessor::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 839 of file processor.h.

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

00840                 {
00841                         TypeDict d;
00842                         d.put("ctf", EMObject::EMDATA, "Ctf object to use for Wiener filter parameters");
00843                         return d;
00844                 }

static Processor* EMAN::Wiener2DFourierProcessor::NEW (  )  [inline, static]

Definition at line 846 of file processor.h.

00847                 {
00848                         return new Wiener2DFourierProcessor();
00849                 }

EMData * Wiener2DFourierProcessor::process ( const EMData *const   image  )  [virtual]

To proccess an image out-of-place.

For those processors which can only be processed out-of-place, override this function to give the right behavior.

Parameters:
image The image will be copied, actual process happen on copy of image.
Returns:
the image processing result, may or may not be the same size of the input image

Reimplemented from EMAN::Processor.

Definition at line 1078 of file processor.cpp.

References EMAN::Ctf::compute_2d_complex(), ctf, EMAN::Ctf::CTF_WIENER_FILTER, EMAN::EMData::do_ift(), in, and EMAN::EMData::mult().

Referenced by process_inplace().

01079 {
01080         const EMData *in2 = 0;
01081         if (in->is_complex()) in2=in;
01082         else in=in->do_fft();
01083 
01084         EMData *filt = in->copy_head();
01085         Ctf *ictf = ctf;
01086 
01087         if (!ictf) ctf=(Ctf *)in->get_attr("ctf");
01088 
01089         ictf->compute_2d_complex(filt,Ctf::CTF_WIENER_FILTER);
01090         filt->mult(*in2);
01091         EMData *ret=filt->do_ift();
01092 
01093         delete filt;
01094         if (!in->is_complex()) delete in2;
01095 
01096         if(!ictf && ctf) {delete ctf; ctf=0;}
01097         return(ret);
01098 /*      const EMData *fft;
01099         float *fftd;
01100         int f=0;
01101 
01102         if (!image) {
01103                 LOGWARN("NULL Image");
01104                 return ret;
01105         }
01106 
01107         if (!image->is_complex()) {
01108                 fft = image->do_fft();
01109                 fftd = fft->get_data();
01110                 f=1;
01111         }
01112         else {
01113                 fft=image;
01114                 fftd=image->get_data();
01115         }
01116         powd=image->get_data();
01117 
01118         int bad=0;
01119         for (int i=0; i<image->get_xsize()*image->get_ysize(); i+=2) {
01120                 snr=(fftd[i]*fftd[i]+fftd[i+1]*fftd[i+1]-powd[i])/powd[i];
01121                 if (snr<0) { bad++; snr=0; }
01122 
01123         }
01124 
01125         print("%d bad pixels\n",snr);
01126 */      return ret;
01127 
01128 }

void Wiener2DFourierProcessor::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 1130 of file processor.cpp.

References EMAN::EMData::get_data(), EMAN::EMData::get_xsize(), EMAN::EMData::get_ysize(), EMAN::EMData::get_zsize(), process(), and EMAN::EMData::update().

01130                                                             {
01131         EMData *tmp=process(image);
01132         memcpy(image->get_data(),tmp->get_data(),image->get_xsize()*image->get_ysize()*image->get_zsize()*sizeof(float));
01133         delete tmp;
01134         image->update();
01135         return;
01136 }

void EMAN::Wiener2DFourierProcessor::set_params ( const Dict new_params  )  [inline, virtual]

Set the processor parameters using a key/value dictionary.

Parameters:
new_params A dictionary containing the new parameters.

Reimplemented from EMAN::Processor.

Definition at line 832 of file processor.h.

References ctf, and EMAN::Processor::params.

00833                 {
00834                         params = new_params;
00835                         ctf = params["ctf"];
00836 //                      printf("%s %f\n",params.keys()[0].c_str(),lowpass);
00837                 }


Member Data Documentation

Ctf* EMAN::Wiener2DFourierProcessor::ctf [protected]

Definition at line 859 of file processor.h.

Referenced by process(), and set_params().

const string Wiener2DFourierProcessor::NAME = "filter.wiener2d" [static]

Definition at line 856 of file processor.h.

Referenced by get_name().


The documentation for this class was generated from the following files:
Generated on Thu Nov 17 12:45:53 2011 for EMAN2 by  doxygen 1.4.7