#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 | |
Processor * | NEW () |
Static Public Attributes | |
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 5708 of file processor.h.
|
Implements EMAN::FourierAnlProcessor. Definition at line 6263 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::XYData::set_x(), EMAN::XYData::set_y(), and sqrt(). 06263 { 06264 // The radial mask comes in with the existing radial image profile 06265 // The radial mask runs from 0 to the 1-D Nyquist (it leaves out the corners in Fourier space) 06266 06267 EMData *to = params["to"]; 06268 XYData *sf = new XYData(); 06269 float apixto = to->get_attr("apix_x"); 06270 06271 06272 if (to->is_complex()) { 06273 vector<float> rd=to->calc_radial_dist(to->get_ysize()/2.0f,0,1.0f,1); 06274 for (size_t i=0; i<rd.size(); ++i) { 06275 sf->set_x(i,i/(apixto*2.0f*rd.size())); 06276 sf->set_y(i,rd[i]); 06277 } 06278 } 06279 else { 06280 EMData *tmp=to->do_fft(); 06281 vector<float> rd=tmp->calc_radial_dist(to->get_ysize()/2,0,1.0,1); 06282 for (size_t i=0; i<rd.size(); ++i) { 06283 sf->set_x(i,i/(apixto*2.0f*rd.size())); 06284 sf->set_y(i,rd[i]); 06285 } 06286 delete tmp; 06287 } 06288 06289 float apix=image->get_attr("apix_x"); 06290 06291 // sf->write_file("a.txt"); 06292 // Util::save_data(0,sf->get_x(1),radial_mask,"b.txt"); 06293 06294 int n = radial_mask.size(); 06295 for (int i=0; i<n; i++) { 06296 if (radial_mask[i]>0) radial_mask[i]= sqrt(sf->get_yatx(i/(apix*2.0f*n),false)/radial_mask[i]); 06297 else if (i>0) radial_mask[i]=radial_mask[i-1]; 06298 } 06299 06300 // Util::save_data(0,sf->get_x(1),radial_mask,"c.txt"); 06301 06302 delete sf; 06303 }
|
|
Get the descrition of this specific processor. This function must be overwritten by a subclass.
Implements EMAN::Processor. Definition at line 5717 of file processor.h. 05718 { 05719 return "Filters the image so its 1-D power spectrum matches a second image"; 05720 }
|
|
Get the processor's name. Each processor is identified by a unique name.
Implements EMAN::Processor. Definition at line 5712 of file processor.h. 05713 {
05714 return NAME;
05715 }
|
|
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 5727 of file processor.h. References EMAN::TypeDict::put(). 05728 { 05729 TypeDict d; 05730 d.put("to", EMObject::EMDATA, "The image to match with. Make sure apix values are correct."); 05731 return d; 05732 }
|
|
Definition at line 5722 of file processor.h. 05723 { 05724 return new MatchSFProcessor(); 05725 }
|
|
Definition at line 186 of file processor.cpp. |