#include <processor.h>
Inheritance diagram for EMAN::CCDNormProcessor:
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 string | get_desc () const |
Get the descrition of this specific processor. | |
virtual TypeDict | get_param_types () const |
Get processor parameter information in a dictionary. | |
Static Public Member Functions | |
Processor * | NEW () |
Static Public Attributes | |
const string | NAME = "filter.ccdnorm" |
width | number of pixels on either side of the seam to sample |
Definition at line 6558 of file processor.h.
|
Get the descrition of this specific processor. This function must be overwritten by a subclass.
Implements EMAN::Processor. Definition at line 6573 of file processor.h. 06574 { 06575 return "normalize the 4 quadrants of a CCD image"; 06576 }
|
|
Get the processor's name. Each processor is identified by a unique name.
Implements EMAN::Processor. Definition at line 6563 of file processor.h. 06564 {
06565 return NAME;
06566 }
|
|
Get processor parameter information in a dictionary. Each parameter has one record in the dictionary. Each record contains its name, data-type, and description.
Reimplemented from EMAN::Processor. Definition at line 6578 of file processor.h. References EMAN::TypeDict::put(). 06579 { 06580 TypeDict d; 06581 d.put("width", EMObject::INT, "number of pixels on either side of the seam to sample"); 06582 return d; 06583 }
|
|
Definition at line 6568 of file processor.h. 06569 { 06570 return new CCDNormProcessor(); 06571 }
|
|
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.
Implements EMAN::Processor. Definition at line 7793 of file processor.cpp. References EMAN::Log::error(), EMAN::EMData::get_col(), EMAN::EMData::get_data(), EMAN::EMData::get_row(), EMAN::EMData::get_value_at(), EMAN::EMData::get_xsize(), EMAN::EMData::get_ysize(), EMAN::EMData::get_zsize(), EMAN::Log::logger(), EMAN::Log::set_level(), EMAN::EMData::set_value_at_fast(), x, and y. 07794 { 07795 if (!image) { 07796 Log::logger()->set_level(Log::ERROR_LOG); 07797 Log::logger()->error("Null image during call to CCDNorm\n"); 07798 return; 07799 } 07800 if (image->get_zsize() > 1) { 07801 Log::logger()->set_level(Log::ERROR_LOG); 07802 Log::logger()->error("CCDNorm does not support 3d images\n"); 07803 return; 07804 } 07805 07806 int xs = image->get_xsize(); 07807 int ys = image->get_ysize(); 07808 07809 // width of sample area on either side of the seams 07810 int width = params["width"]; 07811 07812 width%=(xs > ys ? xs/2 : ys/2); // make sure width is a valid value 07813 if (width==0) { 07814 width=1; 07815 } 07816 07817 // get the 4 "seams" of the image 07818 float *left, *right, *top, *bottom; 07819 07820 double *temp; 07821 temp= (double*)malloc((xs > ys ? xs*width : ys*width)*sizeof(double)); 07822 if (temp==NULL) { 07823 Log::logger()->set_level(Log::ERROR_LOG); 07824 Log::logger()->error("Could not allocate enough memory during call to CCDNorm\n"); 07825 return; 07826 } 07827 07828 int x, y, z; 07829 07830 // the mean values of each seam and the average 07831 double mL,mR,mT,mB; 07832 07833 // how much to shift each seam 07834 double nl,nr,nt,nb; 07835 07836 // quad. shifting amount 07837 double q1,q2,q3,q4; 07838 07839 // calc. the mean for each quadrant 07840 for (z=0; z<width; z++) { 07841 left = image->get_col(xs/2 -1-z)->get_data(); 07842 for (y=0; y<ys; y++) 07843 temp[z*ys+y]=left[y]; 07844 } 07845 mL=gsl_stats_mean(temp,1,ys*width); 07846 07847 for (z=0; z<width; z++) { 07848 right = image->get_col(xs/2 +z)->get_data(); 07849 for (y=0; y<ys; y++) 07850 temp[z*ys+y]=right[y]; 07851 } 07852 mR=gsl_stats_mean(temp,1,ys*width); 07853 07854 for (z=0; z<width; z++) { 07855 top = image->get_row(ys/2 -1-z)->get_data(); 07856 for (x=0; x<xs; x++) 07857 temp[z*xs+x]=top[x]; 07858 } 07859 mT=gsl_stats_mean(temp,1,xs*width); 07860 07861 for (z=0; z<width; z++) { 07862 bottom = image->get_row(ys/2 +z)->get_data(); 07863 for (x=0; x<xs; x++) 07864 temp[z*xs+x]=bottom[x]; 07865 } 07866 mB=gsl_stats_mean(temp,1,xs*width); 07867 07868 free(temp); 07869 07870 nl=(mL+mR)/2-mL; 07871 nr=(mL+mR)/2-mR; 07872 nt=(mT+mB)/2-mT; 07873 nb=(mT+mB)/2-mB; 07874 07875 q1=nl+nt; 07876 q2=nr+nt; 07877 q3=nr+nb; 07878 q4=nl+nb; 07879 07880 // change the pixel values 07881 for (x = 0; x < xs / 2; x++) 07882 for (y = 0; y < ys / 2; y++) { 07883 image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q1)); 07884 } 07885 for (x = xs / 2; x < xs; x++) 07886 for (y = 0; y < ys / 2; y++) { 07887 image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q2)); 07888 } 07889 for (x = xs / 2; x < xs; x++) 07890 for (y = ys / 2; y < ys; y++) { 07891 image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q3)); 07892 } 07893 for (x = 0; x < xs / 2; x++) 07894 for (y = ys / 2; y < ys; y++) { 07895 image->set_value_at_fast(x, y, image->get_value_at(x, y) + static_cast<float>(q4)); 07896 } 07897 07898 }
|
|
Definition at line 210 of file processor.cpp. |