#include <processor.h>
Inheritance diagram for EMAN::MatchSFProcessor:
Public Member Functions | |
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 | |
static Processor * | NEW () |
Static Public Attributes | |
static const string | NAME = "filter.matchto" |
Protected Member Functions | |
void | create_radial_func (vector< float > &radial_mask, EMData *image) const |
to | EMData object to match to. Make sure apix values are set properly |
Definition at line 5706 of file processor.h.
void MatchSFProcessor::create_radial_func | ( | vector< float > & | radial_mask, | |
EMData * | image | |||
) | const [protected, virtual] |
Implements EMAN::FourierAnlProcessor.
Definition at line 6259 of file processor.cpp.
References EMAN::EMData::calc_radial_dist(), EMAN::EMData::do_fft(), EMAN::EMData::get_attr(), EMAN::XYData::get_yatx(), EMAN::EMData::get_ysize(), EMAN::EMData::is_complex(), EMAN::Processor::params, EMAN::XYData::set_x(), EMAN::XYData::set_y(), and sqrt().
06259 { 06260 // The radial mask comes in with the existing radial image profile 06261 // The radial mask runs from 0 to the 1-D Nyquist (it leaves out the corners in Fourier space) 06262 06263 EMData *to = params["to"]; 06264 XYData *sf = new XYData(); 06265 float apixto = to->get_attr("apix_x"); 06266 06267 06268 if (to->is_complex()) { 06269 vector<float> rd=to->calc_radial_dist(to->get_ysize()/2.0f,0,1.0f,1); 06270 for (size_t i=0; i<rd.size(); ++i) { 06271 sf->set_x(i,i/(apixto*2.0f*rd.size())); 06272 sf->set_y(i,rd[i]); 06273 } 06274 } 06275 else { 06276 EMData *tmp=to->do_fft(); 06277 vector<float> rd=tmp->calc_radial_dist(to->get_ysize()/2,0,1.0,1); 06278 for (size_t i=0; i<rd.size(); ++i) { 06279 sf->set_x(i,i/(apixto*2.0f*rd.size())); 06280 sf->set_y(i,rd[i]); 06281 } 06282 delete tmp; 06283 } 06284 06285 float apix=image->get_attr("apix_x"); 06286 06287 // sf->write_file("a.txt"); 06288 // Util::save_data(0,sf->get_x(1),radial_mask,"b.txt"); 06289 06290 int n = radial_mask.size(); 06291 for (int i=0; i<n; i++) { 06292 if (radial_mask[i]>0) radial_mask[i]= sqrt(sf->get_yatx(i/(apix*2.0f*n),false)/radial_mask[i]); 06293 else if (i>0) radial_mask[i]=radial_mask[i-1]; 06294 } 06295 06296 // Util::save_data(0,sf->get_x(1),radial_mask,"c.txt"); 06297 06298 delete sf; 06299 }
virtual string EMAN::MatchSFProcessor::get_desc | ( | ) | const [inline, virtual] |
Get the descrition of this specific processor.
This function must be overwritten by a subclass.
Implements EMAN::Processor.
Definition at line 5715 of file processor.h.
virtual string EMAN::MatchSFProcessor::get_name | ( | ) | const [inline, virtual] |
Get the processor's name.
Each processor is identified by a unique name.
Implements EMAN::Processor.
Definition at line 5710 of file processor.h.
References NAME.
05711 { 05712 return NAME; 05713 }
virtual TypeDict EMAN::MatchSFProcessor::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.
Reimplemented from EMAN::FourierAnlProcessor.
Definition at line 5725 of file processor.h.
References EMAN::EMObject::EMDATA, and EMAN::TypeDict::put().
05726 { 05727 TypeDict d; 05728 d.put("to", EMObject::EMDATA, "The image to match with. Make sure apix values are correct."); 05729 return d; 05730 }
static Processor* EMAN::MatchSFProcessor::NEW | ( | ) | [inline, static] |
const string MatchSFProcessor::NAME = "filter.matchto" [static] |