EMAN2
Public Member Functions | Static Public Member Functions | Protected Member Functions
EMAN::FourierAnlProcessor Class Reference

Same as FourierProcessor, except first computes the current image radial power spectrum and passes it to the processor (no radial oversampling, number of elements = radius) More...

#include <processor.h>

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

List of all members.

Public Member Functions

void process_inplace (EMData *image)
 To process an image in-place.
TypeDict get_param_types () const
 Get processor parameter information in a dictionary.

Static Public Member Functions

static string get_group_desc ()
 Get the description of this group of processors.

Protected Member Functions

virtual void preprocess (EMData *)
virtual void create_radial_func (vector< float > &radial_mask, EMData *image) const =0

Detailed Description

Same as FourierProcessor, except first computes the current image radial power spectrum and passes it to the processor (no radial oversampling, number of elements = radius)

Parameters:
cutoff_absProcessor radius in terms of Nyquist (0-.5).
cutoff_pixelsWidth in Fourier pixels (0 - size()/2).
cutoff_freqResolution in 1/A (0 - 1 / size*apix).
apixOverride A/pix in the image header (changes x,y and z).

Definition at line 394 of file processor.h.


Member Function Documentation

virtual void EMAN::FourierAnlProcessor::create_radial_func ( vector< float > &  radial_mask,
EMData image 
) const [protected, pure virtual]
static string EMAN::FourierAnlProcessor::get_group_desc ( ) [inline, static]

Get the description of this group of processors.

This function is defined in a parent class. It gives a introduction to a group of processors.

Returns:
The description of this group of processors.

Reimplemented from EMAN::Processor.

Definition at line 399 of file processor.h.

                {
                        return "Fourier Filter processors are a group of processor in the frequency domain. Before using such processors on an image, the image must be transformed from real space to the fourier space. FourierProcessor class is the base class of fourier space processors. Each specific processor is either a lowpass filter processor, or a highpass filter processor, or neighter. The unit of lowpass and highpass parameters are in terms of Nyquist, valid range is [0,0.5]. ";
                }
TypeDict EMAN::FourierAnlProcessor::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.

Reimplemented in EMAN::LowpassAutoBProcessor, EMAN::MatchSFProcessor, and EMAN::SetSFProcessor.

Definition at line 404 of file processor.h.

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

                {
                        TypeDict d;
                        d.put("cutoff_abs", EMObject::FLOAT, "Processor radius in terms of Nyquist (0-.5)");
                        d.put("cutoff_pixels", EMObject::FLOAT, " Width in Fourier pixels (0 - size()/2)");
                        d.put("cutoff_freq", EMObject::FLOAT, "1/Resolution in 1/A (0 - 1 / 2*apix). eg - a 20 A filter is cutoff_freq=0.05");
                        d.put("apix", EMObject::FLOAT, " Override A/pix in the image header (changes x,y and z)");
                        d.put("return_radial", EMObject::BOOL, "Return the radial filter function as an attribute (filter_curve)");
                        d.put("interpolate", EMObject::BOOL, "Whether or not to interpolate the radial scaling function. Default=false. Prb should be true.");
                        return d;
                }
virtual void EMAN::FourierAnlProcessor::preprocess ( EMData ) [inline, protected, virtual]

Reimplemented in EMAN::LowpassAutoBProcessor, and EMAN::HighpassAutoPeakProcessor.

Definition at line 417 of file processor.h.

Referenced by process_inplace().

{}
void FourierAnlProcessor::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 616 of file processor.cpp.

References EMAN::EMData::apply_radial_func(), EMAN::EMData::calc_radial_dist(), create_radial_func(), EMAN::EMData::do_fft(), EMAN::EMData::do_ift(), EMAN::EMData::get_data(), EMAN::EMData::get_xsize(), EMAN::EMData::get_ysize(), EMAN::EMData::get_zsize(), EMAN::EMData::is_complex(), LOGWARN, EMAN::Processor::params, preprocess(), EMAN::EMData::set_attr(), EMAN::Dict::set_default(), and EMAN::EMData::update().

{
        if (!image) {
                LOGWARN("NULL Image");
                return;
        }

        preprocess(image);

//      int array_size = FFTRADIALOVERSAMPLE * image->get_ysize();
//      float step=0.5f/array_size;
//
//      vector < float >yarray(array_size);

        bool return_radial=(bool)params.set_default("return_radial",0);
        bool interpolate=(bool)params.set_default("interpolate",0);

        if (image->is_complex()) {
                vector <float>yarray = image->calc_radial_dist(image->get_ysize()/2,0,1.0,1);
                create_radial_func(yarray,image);
                image->apply_radial_func(0, 0.5f/yarray.size(), yarray);
                if (return_radial) image->set_attr("filter_curve",yarray);
        }
        else {
                EMData *fft = image->do_fft();
                vector <float>yarray = fft->calc_radial_dist(fft->get_ysize()/2,0,1.0,1);
                create_radial_func(yarray,image);
                // 4/30/10 stevel turned off interpolation to fix problem with matched filter
                // 9/12/14 stevel, not sure why I turned off interp. Seems to cause rather than fix problems. Adding option to enable. Afraid to turn it on
                fft->apply_radial_func(0, 0.5f/yarray.size(), yarray,0);
                EMData *ift = fft->do_ift();

                memcpy(image->get_data(),ift->get_data(),ift->get_xsize()*ift->get_ysize()*ift->get_zsize()*sizeof(float));
                if (return_radial) image->set_attr("filter_curve",yarray);

                //ift->update(); Unecessary

                delete fft;
                delete ift;

        }

        image->update();
}

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