#include <processor.h>
Inheritance diagram for EMAN::RotationalSubstractProcessor:
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. | |
Static Public Member Functions | |
static Processor * | NEW () |
Static Public Attributes | |
static const string | NAME = "math.rotationalsubtract" |
Definition at line 4462 of file processor.h.
virtual string EMAN::RotationalSubstractProcessor::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 4477 of file processor.h.
virtual string EMAN::RotationalSubstractProcessor::get_name | ( | ) | const [inline, virtual] |
Get the processor's name.
Each processor is identified by a unique name.
Implements EMAN::Processor.
Definition at line 4467 of file processor.h.
References NAME.
04468 { 04469 return NAME; 04470 }
static Processor* EMAN::RotationalSubstractProcessor::NEW | ( | ) | [inline, static] |
Definition at line 4472 of file processor.h.
04473 { 04474 return new RotationalSubstractProcessor(); 04475 }
void RotationalSubstractProcessor::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.
image | The image to be processed. |
Implements EMAN::Processor.
Definition at line 4244 of file processor.cpp.
References EMAN::EMData::calc_radial_dist(), dist(), EMAN::EMData::get_data(), EMAN::EMData::get_ndim(), EMAN::EMData::get_xsize(), EMAN::EMData::get_ysize(), ImageDimensionException, EMAN::EMData::is_complex(), LOGWARN, rdata, EMAN::EMData::update(), and x.
04245 { 04246 if (!image || image->is_complex()) { 04247 LOGWARN("only works on real image. do nothing."); 04248 return; 04249 } 04250 04251 if (image->get_ndim() != 2) throw ImageDimensionException("This processor works only for 2D images"); 04252 04253 float *rdata = image->get_data(); 04254 int nx = image->get_xsize(); 04255 int ny = image->get_ysize(); 04256 04257 vector < float >dist = image->calc_radial_dist(nx / 2, 0, 1,0); 04258 04259 int c = 0; 04260 for (int y = 0; y < ny; y++) { 04261 for (int x = 0; x < nx; x++, c++) { 04262 #ifdef _WIN32 04263 float r = (float) _hypot(x - nx / 2, y - ny / 2); 04264 #else 04265 float r = (float) hypot(x - nx / 2, y - ny / 2); 04266 #endif 04267 int i = (int) floor(r); 04268 r -= i; 04269 if (i >= 0 && i < nx / 2 - 1) { 04270 rdata[c] -= dist[i] * (1.0f - r) + dist[i + 1] * r; 04271 } 04272 else { 04273 rdata[c] = 0; 04274 } 04275 } 04276 } 04277 04278 image->update(); 04279 }
const string RotationalSubstractProcessor::NAME = "math.rotationalsubtract" [static] |