EMAN::CCDNormProcessor Class Reference

Try to normalize the 4 quadrants of a CCD image. More...

#include <processor.h>

Inheritance diagram for EMAN::CCDNormProcessor:

Collaboration diagram for EMAN::CCDNormProcessor:

Public Member Functions
virtual void	process_inplace (EMData *image)
	To process an image in-place. More...
virtual string	get_name () const
	Get the processor's name. More...
virtual string	get_desc () const
	Get the descrition of this specific processor. More...
virtual TypeDict	get_param_types () const
	Get processor parameter information in a dictionary. More...
Public Member Functions inherited from EMAN::Processor
virtual	~Processor ()
virtual EMData *	process (const EMData *const image)
	To proccess an image out-of-place. More...
virtual void	process_list_inplace (vector< EMData * > &images)
	To process multiple images using the same algorithm. More...
virtual Dict	get_params () const
	Get the processor parameters in a key/value dictionary. More...
virtual void	set_params (const Dict &new_params)
	Set the processor parameters using a key/value dictionary. More...

Static Public Member Functions
static Processor *	NEW ()
Static Public Member Functions inherited from EMAN::Processor
static string	get_group_desc ()
	Get the description of this group of processors. More...
static void	EMFourierFilterInPlace (EMData *fimage, Dict params)
	Compute a Fourier-filter processed image in place. More...
static EMData *	EMFourierFilter (EMData *fimage, Dict params)
	Compute a Fourier-processor processed image without altering the original image. More...

Static Public Attributes
static const string	NAME = "filter.ccdnorm"

Additional Inherited Members
Public Types inherited from EMAN::Processor
enum	fourier_filter_types {   TOP_HAT_LOW_PASS , TOP_HAT_HIGH_PASS , TOP_HAT_BAND_PASS , TOP_HOMOMORPHIC ,   GAUSS_LOW_PASS , GAUSS_HIGH_PASS , GAUSS_BAND_PASS , GAUSS_INVERSE ,   GAUSS_HOMOMORPHIC , BUTTERWORTH_LOW_PASS , BUTTERWORTH_HIGH_PASS , BUTTERWORTH_HOMOMORPHIC ,   KAISER_I0 , KAISER_SINH , KAISER_I0_INVERSE , KAISER_SINH_INVERSE ,   SHIFT , TANH_LOW_PASS , TANH_HIGH_PASS , TANH_HOMOMORPHIC ,   TANH_BAND_PASS , RADIAL_TABLE , CTF_ }
	Fourier filter Processor type enum. More...
Protected Attributes inherited from EMAN::Processor
Dict	params

Detailed Description

Try to normalize the 4 quadrants of a CCD image.

Author

Deepy Mann dsman.nosp@m.n@bc.nosp@m.m.tmc.nosp@m..edu

Date

9-2005

Parameters

width

number of pixels on either side of the seam to sample

Definition at line 8921 of file processor.h.

Member Function Documentation

get_desc()

virtual string EMAN::CCDNormProcessor::get_desc ( ) const

inlinevirtual

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 8936 of file processor.h.

                {
                        return "normalize the 4 quadrants of a CCD image";
                }

get_name()

virtual string EMAN::CCDNormProcessor::get_name ( ) const

inlinevirtual

Get the processor's name.

Each processor is identified by a unique name.

Returns

The processor's name.

Implements EMAN::Processor.

Definition at line 8926 of file processor.h.

                {
                        return NAME;
                }

References NAME.

get_param_types()

virtual TypeDict EMAN::CCDNormProcessor::get_param_types ( ) const

inlinevirtual

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 8941 of file processor.h.

                {
                        TypeDict d;
                        d.put("width", EMObject::INT, "number of pixels on either side of the seam to sample");
                        return d;
                }

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

NEW()

static Processor * EMAN::CCDNormProcessor::NEW ( )

inlinestatic

Definition at line 8931 of file processor.h.

                {
                        return new CCDNormProcessor();
                }

process_inplace()

void CCDNormProcessor::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 10255 of file processor.cpp.

{
        if (!image) {
          Log::logger()->set_level(Log::ERROR_LOG);
          Log::logger()->error("Null image during call to CCDNorm\n");
          return;
        }
        if (image->get_zsize() > 1) {
          Log::logger()->set_level(Log::ERROR_LOG);
          Log::logger()->error("CCDNorm does not support 3d images\n");
          return;
        }
 
        int xs = image->get_xsize();
        int ys = image->get_ysize();
 
        // width of sample area on either side of the seams
        int width = params["width"];
 
        width%=(xs > ys ? xs/2 : ys/2);  // make sure width is a valid value
        if (width==0) {
          width=1;
        }
 
        // get the 4 "seams" of the image
        float *left, *right, *top, *bottom;
 
        double *temp;
        temp= (double*)malloc((xs > ys ? xs*width : ys*width)*sizeof(double));
        if (temp==NULL) {
          Log::logger()->set_level(Log::ERROR_LOG);
          Log::logger()->error("Could not allocate enough memory during call to CCDNorm\n");
          return;
        }
 
        int x, y, z;
 
        // the mean values of each seam and the average
        double mL,mR,mT,mB;
 
        // how much to shift each seam
        double nl,nr,nt,nb;
 
        // quad. shifting amount
        double q1,q2,q3,q4;
 
        // calc. the mean for each quadrant
        for (z=0; z<width; z++) {
          left = image->get_col(xs/2 -1-z)->get_data();
          for (y=0; y<ys; y++)
                temp[z*ys+y]=left[y];
        }
        mL=gsl_stats_mean(temp,1,ys*width);
 
        for (z=0; z<width; z++) {
          right = image->get_col(xs/2 +z)->get_data();
          for (y=0; y<ys; y++)
                temp[z*ys+y]=right[y];
        }
        mR=gsl_stats_mean(temp,1,ys*width);
 
        for (z=0; z<width; z++) {
          top = image->get_row(ys/2 -1-z)->get_data();
          for (x=0; x<xs; x++)
                temp[z*xs+x]=top[x];
        }
        mT=gsl_stats_mean(temp,1,xs*width);
 
        for (z=0; z<width; z++) {
          bottom = image->get_row(ys/2 +z)->get_data();
          for (x=0; x<xs; x++)
                temp[z*xs+x]=bottom[x];
        }
        mB=gsl_stats_mean(temp,1,xs*width);
 
        free(temp);
 
        nl=(mL+mR)/2-mL;
        nr=(mL+mR)/2-mR;
        nt=(mT+mB)/2-mT;
        nb=(mT+mB)/2-mB;
 
        q1=nl+nt;
        q2=nr+nt;
        q3=nr+nb;
        q4=nl+nb;
 
        // change the pixel values
        for (x = 0; x < xs / 2; x++)
          for (y = 0; y < ys / 2; y++) {
                image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q1));
          }
        for (x = xs / 2; x < xs; x++)
          for (y = 0; y < ys / 2; y++) {
                image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q2));
          }
        for (x = xs / 2; x < xs; x++)
          for (y = ys / 2; y < ys; y++) {
                image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q3));
          }
        for (x = 0; x < xs / 2; x++)
          for (y = ys / 2; y < ys; y++) {
                image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q4));
          }
 
}

References EMAN::Log::error(), EMAN::Log::ERROR_LOG, EMAN::Log::logger(), EMAN::Processor::params, EMAN::Log::set_level(), x, and y.

Member Data Documentation

NAME

const string CCDNormProcessor::NAME = "filter.ccdnorm"

static

Definition at line 8948 of file processor.h.

Referenced by get_name().

---

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

• processor.h
• processor.cpp