#include <util.h>
Inheritance diagram for EMAN::Util::KaiserBessel:
Public Member Functions | |
KaiserBessel (float alpha_, int K, float r_, float v_, int N_, float vtable_=0.f, int ntable_=5999) | |
virtual | ~KaiserBessel () |
float | I0table_maxerror () |
Compute the maximum error in the table. | |
vector< float > | dump_table () |
virtual float | sinhwin (float x) const |
Kaiser-Bessel Sinh window function. | |
virtual float | i0win (float x) const |
Kaiser-Bessel I0 window function. | |
float | i0win_tab (float x) const |
Kaiser-Bessel I0 window function (uses table lookup). | |
int | get_window_size () const |
Return the size of the I0 window. | |
kbsinh_win | get_kbsinh_win () |
Sinh window function object factory. | |
kbi0_win | get_kbi0_win () |
I0 window function object factory. | |
Protected Member Functions | |
virtual void | build_I0table () |
2*pi*alpha*r*vadjust | |
Protected Attributes | |
float | alpha |
float | v |
float | r |
int | N |
Kaiser-Bessel parameters. | |
int | K |
size in Ix-space | |
float | vtable |
I0 window size. | |
int | ntable |
table I0 non-zero domain maximum | |
vector< float > | i0table |
float | dtable |
float | alphar |
table spacing | |
float | fac |
alpha*r | |
float | vadjust |
2*pi*alpha*r*v | |
float | facadj |
float | fltb |
Tabulate I0 window for speed. | |
Classes | |
class | kbi0_win |
I0 window function object. More... | |
class | kbsinh_win |
Sinh window function object. More... |
(It's a class so that the windowing parameters may be instantiated and held in the instance object.)
The I0 version can be tabulated and interpolated upon demand, but the max error needs to be checked. The "vtable" parameter corresponds to the maximum value of x for which the I0 window is non-zero. Setting "vtable" different from "v" corresponds to a change in units of x. In practice, it is often handy to replace x in some sort of absolute units with x described in terms of grid intervals.
The get_kbsinh_win and get_kbi0_win functions return single-argument function objects, which is what a generic routine is likely to want.
Definition at line 305 of file util.h.
Util::KaiserBessel::KaiserBessel | ( | float | alpha_, | |
int | K, | |||
float | r_, | |||
float | v_, | |||
int | N_, | |||
float | vtable_ = 0.f , |
|||
int | ntable_ = 5999 | |||
) |
Definition at line 1981 of file util_sparx.cpp.
References alpha, alphar, build_I0table(), fac, facadj, K, r, twopi, v, vadjust, and vtable.
01983 : alpha(alpha_), v(v_), r(r_), N(N_), K(K_), vtable(vtable_), 01984 ntable(ntable_) { 01985 // Default values are alpha=1.25, K=6, r=0.5, v = K/2 01986 if (0.f == v) v = float(K)/2; 01987 if (0.f == vtable) vtable = v; 01988 alphar = alpha*r; 01989 fac = static_cast<float>(twopi)*alphar*v; 01990 vadjust = 1.0f*v; 01991 facadj = static_cast<float>(twopi)*alphar*vadjust; 01992 build_I0table(); 01993 }
virtual EMAN::Util::KaiserBessel::~KaiserBessel | ( | ) | [inline, virtual] |
void Util::KaiserBessel::build_I0table | ( | ) | [protected, virtual] |
2*pi*alpha*r*vadjust
Reimplemented in EMAN::Util::FakeKaiserBessel.
Definition at line 2004 of file util_sparx.cpp.
References facadj, fltb, i0table, K, N, ntable, EMAN::Util::round(), sqrt(), and vadjust.
Referenced by KaiserBessel().
02004 { 02005 i0table.resize(ntable+1); // i0table[0:ntable] 02006 int ltab = int(round(float(ntable)/1.25f)); 02007 fltb = float(ltab)/(K/2); 02008 float val0 = static_cast<float>(gsl_sf_bessel_I0(facadj)); 02009 for (int i=ltab+1; i <= ntable; i++) i0table[i] = 0.f; 02010 for (int i=0; i <= ltab; i++) { 02011 float s = float(i)/fltb/N; 02012 if (s < vadjust) { 02013 float rt = sqrt(1.f - pow(s/vadjust, 2)); 02014 i0table[i] = static_cast<float>(gsl_sf_bessel_I0(facadj*rt))/val0; 02015 } else { 02016 i0table[i] = 0.f; 02017 } 02018 // cout << " "<<s*N<<" "<<i0table[i] <<endl; 02019 } 02020 }
vector<float> EMAN::Util::KaiserBessel::dump_table | ( | ) | [inline] |
kbi0_win EMAN::Util::KaiserBessel::get_kbi0_win | ( | ) | [inline] |
kbsinh_win EMAN::Util::KaiserBessel::get_kbsinh_win | ( | ) | [inline] |
int EMAN::Util::KaiserBessel::get_window_size | ( | ) | const [inline] |
Return the size of the I0 window.
Definition at line 349 of file util.h.
Referenced by EMAN::EMData::extract_plane(), EMAN::EMData::extract_plane_rect(), EMAN::EMData::extract_plane_rect_fast(), EMAN::EMData::extractline(), EMAN::EMData::extractpoint(), EMAN::EMData::get_pixel_conv(), EMAN::Util::get_pixel_conv_new(), EMAN::Util::get_pixel_conv_new_background(), EMAN::EMData::rot_scale_conv(), and EMAN::EMData::rot_scale_conv7().
float Util::KaiserBessel::I0table_maxerror | ( | ) |
Compute the maximum error in the table.
Definition at line 2022 of file util_sparx.cpp.
References i0table, and ntable.
02022 { 02023 float maxdiff = 0.f; 02024 for (int i = 1; i <= ntable; i++) { 02025 float diff = fabs(i0table[i] - i0table[i-1]); 02026 if (diff > maxdiff) maxdiff = diff; 02027 } 02028 return maxdiff; 02029 }
float Util::KaiserBessel::i0win | ( | float | x | ) | const [virtual] |
Kaiser-Bessel I0 window function.
Reimplemented in EMAN::Util::FakeKaiserBessel.
Definition at line 1995 of file util_sparx.cpp.
References facadj, sqrt(), and vadjust.
Referenced by EMAN::Processor::EMFourierFilterFunc().
01995 { 01996 float val0 = float(gsl_sf_bessel_I0(facadj)); 01997 float absx = fabs(x); 01998 if (absx > vadjust) return 0.f; 01999 float rt = sqrt(1.f - pow(absx/vadjust, 2)); 02000 float res = static_cast<float>(gsl_sf_bessel_I0(facadj*rt))/val0; 02001 return res; 02002 }
float EMAN::Util::KaiserBessel::i0win_tab | ( | float | x | ) | const [inline] |
Kaiser-Bessel I0 window function (uses table lookup).
Definition at line 336 of file util.h.
Referenced by EMAN::EMData::extract_plane(), EMAN::EMData::extract_plane_rect(), EMAN::EMData::extract_plane_rect_fast(), EMAN::EMData::extractline(), EMAN::EMData::extractpoint(), EMAN::Util::extractpoint2(), EMAN::EMData::get_pixel_conv(), EMAN::Util::get_pixel_conv_new(), EMAN::Util::get_pixel_conv_new_background(), EMAN::EMData::rot_scale_conv(), and EMAN::EMData::rot_scale_conv7().
float Util::KaiserBessel::sinhwin | ( | float | x | ) | const [virtual] |
Kaiser-Bessel Sinh window function.
Reimplemented in EMAN::Util::FakeKaiserBessel.
Definition at line 2031 of file util_sparx.cpp.
References alphar, fac, and sqrt().
Referenced by EMAN::EMData::divkbsinh(), EMAN::EMData::divkbsinh_rect(), and EMAN::Processor::EMFourierFilterFunc().
02031 { 02032 float val0 = sinh(fac)/fac; 02033 float absx = fabs(x); 02034 if (0.0 == x) { 02035 float res = 1.0f; 02036 return res; 02037 } else if (absx == alphar) { 02038 return 1.0f/val0; 02039 } else if (absx < alphar) { 02040 float rt = sqrt(1.0f - pow((x/alphar), 2)); 02041 float facrt = fac*rt; 02042 float res = (sinh(facrt)/facrt)/val0; 02043 return res; 02044 } else { 02045 float rt = sqrt(pow((x/alphar),2) - 1.f); 02046 float facrt = fac*rt; 02047 float res = (sin(facrt)/facrt)/val0; 02048 return res; 02049 } 02050 }
float EMAN::Util::KaiserBessel::alpha [protected] |
float EMAN::Util::KaiserBessel::alphar [protected] |
table spacing
Definition at line 315 of file util.h.
Referenced by KaiserBessel(), EMAN::Util::FakeKaiserBessel::sinhwin(), and sinhwin().
float EMAN::Util::KaiserBessel::dtable [protected] |
float EMAN::Util::KaiserBessel::fac [protected] |
alpha*r
Definition at line 316 of file util.h.
Referenced by KaiserBessel(), EMAN::Util::FakeKaiserBessel::sinhwin(), and sinhwin().
float EMAN::Util::KaiserBessel::facadj [protected] |
Definition at line 318 of file util.h.
Referenced by EMAN::Util::FakeKaiserBessel::build_I0table(), build_I0table(), EMAN::Util::FakeKaiserBessel::i0win(), i0win(), and KaiserBessel().
float EMAN::Util::KaiserBessel::fltb [protected] |
Tabulate I0 window for speed.
Definition at line 320 of file util.h.
Referenced by EMAN::Util::FakeKaiserBessel::build_I0table(), and build_I0table().
vector<float> EMAN::Util::KaiserBessel::i0table [protected] |
Definition at line 313 of file util.h.
Referenced by EMAN::Util::FakeKaiserBessel::build_I0table(), build_I0table(), and I0table_maxerror().
int EMAN::Util::KaiserBessel::K [protected] |
size in Ix-space
Definition at line 310 of file util.h.
Referenced by EMAN::Util::FakeKaiserBessel::build_I0table(), build_I0table(), and KaiserBessel().
int EMAN::Util::KaiserBessel::N [protected] |
Kaiser-Bessel parameters.
Definition at line 309 of file util.h.
Referenced by EMAN::Util::FakeKaiserBessel::build_I0table(), and build_I0table().
int EMAN::Util::KaiserBessel::ntable [protected] |
table I0 non-zero domain maximum
Definition at line 312 of file util.h.
Referenced by EMAN::Util::FakeKaiserBessel::build_I0table(), build_I0table(), and I0table_maxerror().
float EMAN::Util::KaiserBessel::r [protected] |
float EMAN::Util::KaiserBessel::v [protected] |
float EMAN::Util::KaiserBessel::vadjust [protected] |
2*pi*alpha*r*v
Definition at line 317 of file util.h.
Referenced by EMAN::Util::FakeKaiserBessel::build_I0table(), build_I0table(), EMAN::Util::FakeKaiserBessel::i0win(), i0win(), and KaiserBessel().
float EMAN::Util::KaiserBessel::vtable [protected] |