#include <geometry.h>

Collaboration diagram for Bspline_model:

Public Member Functions
double	evaluate (double x, double y) const
	Evaluate the model at the point (x,y) More...

Public Attributes
int	l0
	l0 More...

int	lF
	lF; More...

int	m0
	m0 More...

int	mF
	mF; More...

double	x0
	x0 More...

double	y0
	y0 More...

int	SplineDegree
	Order of the Bspline. More...

double	h_x
	Scale X. More...

double	h_y
	Scale Y. More...

MultidimArray< double >	c_ml

Detailed Description

Bspline model class

When you fit a Bspline model this is the type returned. You can use it to evaluate it anywhere.

The model is f(x,y)=sum_{l=l0}^{lF} {sum_{m=m0}^{mF} {c_{ml}Beta_n((x-x0)/h_x-l) Beta_n((y-y0)/h_y-m) } }.

The parameter n is the Bspline degree. l0, lF, m0 and mF are the Bspline indexes. hx and hy are related to the extent of the Bspline.

Definition at line 273 of file geometry.h.

Member Function Documentation

◆ evaluate()

double Bspline_model::evaluate	(	double	x,
		double	y
	)		const

Evaluate the model at the point (x,y)

Definition at line 1238 of file geometry.cpp.

 {
     int SplineDegree_1 = SplineDegree - 1;
     double x_arg = (x - x0) / h_x;
     double y_arg = (y - y0) / h_y;
 
     int l1 = CLIP(CEIL(x_arg - SplineDegree_1), l0, lF);
     int l2 = CLIP(l1 + SplineDegree, l0, lF);
     int m1 = CLIP(CEIL(y_arg - SplineDegree_1), m0, mF);
     int m2 = CLIP(m1 + SplineDegree, m0, mF);
     double columns = 0.0;
     for (int m = m1; m <= m2; m++)
     {
         double rows = 0.0;
         for (int l = l1; l <= l2; l++)
         {
             double xminusl = x_arg - (double)l;
             double Coeff = c_ml(m, l);
             switch (SplineDegree)
             {
             case 2:
                 rows += Coeff * Bspline02(xminusl);
                 break;
             case 3:
                 rows += Coeff * Bspline03(xminusl);
                 break;
             case 4:
                 rows += Coeff * Bspline04(xminusl);
                 break;
             case 5:
                 rows += Coeff * Bspline05(xminusl);
                 break;
             case 6:
                 rows += Coeff * Bspline06(xminusl);
                 break;
             case 7:
                 rows += Coeff * Bspline07(xminusl);
                 break;
             case 8:
                 rows += Coeff * Bspline08(xminusl);
                 break;
             case 9:
                 rows += Coeff * Bspline09(xminusl);
                 break;
             }
         }
 
         double yminusm = y_arg - (double)m;
         switch (SplineDegree)
         {
         case 2:
             columns += rows * Bspline02(yminusm);
             break;
         case 3:
             columns += rows * Bspline03(yminusm);
             break;
         case 4:
             columns += rows * Bspline04(yminusm);
             break;
         case 5:
             columns += rows * Bspline05(yminusm);
             break;
         case 6:
             columns += rows * Bspline06(yminusm);
             break;
         case 7:
             columns += rows * Bspline07(yminusm);
             break;
         case 8:
             columns += rows * Bspline08(yminusm);
             break;
         case 9:
             columns += rows * Bspline09(yminusm);
             break;
         }
     }
     return columns;
 }

Member Data Documentation

◆ c_ml

MultidimArray< double > Bspline_model::c_ml

Bspline coefficients, c_{ml}

The logical indexes of this matrix go from Y=[m0...mF] and X=[l0...lF]

Definition at line 302 of file geometry.h.

◆ h_x

double Bspline_model::h_x

Scale X.

Definition at line 294 of file geometry.h.

◆ h_y

double Bspline_model::h_y

Scale Y.

Definition at line 296 of file geometry.h.

◆ l0

int Bspline_model::l0

l0

Definition at line 277 of file geometry.h.

◆ lF

int Bspline_model::lF

lF;

Definition at line 279 of file geometry.h.

◆ m0

int Bspline_model::m0

m0

Definition at line 281 of file geometry.h.

◆ mF

int Bspline_model::mF

mF;

Definition at line 283 of file geometry.h.

◆ SplineDegree

int Bspline_model::SplineDegree

Order of the Bspline.

Definition at line 291 of file geometry.h.

◆ x0

double Bspline_model::x0

x0

Definition at line 286 of file geometry.h.

◆ y0

double Bspline_model::y0

y0

Definition at line 288 of file geometry.h.

The documentation for this class was generated from the following files:

xmippCore/core/geometry.h
xmippCore/core/geometry.cpp

Public Member Functions

Public Attributes