Public Types |
Public Member Functions |
Static Public Member Functions |
Public Attributes |
Private Types |
Static Private Attributes
CubicMACFieldInterp< Data_T > Class Template Reference
#include <FieldInterp.h>
Inheritance diagram for CubicMACFieldInterp< Data_T >:
Public Types | |
| typedef CubicMACFieldInterp | class_type |
| typedef boost::intrusive_ptr < CubicMACFieldInterp > | Ptr |
| typedef Data_T | value_type |
Public Member Functions | |
| Data_T | sample (const MACField< Data_T > &data, const V3d &vsP) const |
Static Public Member Functions | |
| static const char * | classType () |
| static const char * | staticClassName () |
Public Attributes | |
| DEFINE_FIELD_RTTI_CONCRETE_CLASS | |
Private Types | |
| typedef RefBase | base |
| Convenience typedef for referring to base class. | |
Static Private Attributes | |
| static TemplatedFieldType < CubicMACFieldInterp< Data_T > > | ms_classType |
Detailed Description
template<class Data_T>
class CubicMACFieldInterp< Data_T >
Definition at line 434 of file FieldInterp.h.
Member Typedef Documentation
template<class Data_T >
| typedef Data_T CubicMACFieldInterp< Data_T >::value_type |
Definition at line 440 of file FieldInterp.h.
template<class Data_T >
| typedef boost::intrusive_ptr<CubicMACFieldInterp> CubicMACFieldInterp< Data_T >::Ptr |
Reimplemented from RefBase.
Definition at line 441 of file FieldInterp.h.
template<class Data_T >
| typedef CubicMACFieldInterp CubicMACFieldInterp< Data_T >::class_type |
Definition at line 445 of file FieldInterp.h.
template<class Data_T >
typedef RefBase CubicMACFieldInterp< Data_T >::base [private] |
Convenience typedef for referring to base class.
Definition at line 471 of file FieldInterp.h.
Member Function Documentation
template<class Data_T >
| static const char* CubicMACFieldInterp< Data_T >::staticClassName | ( | ) | [inline, static] |
Definition at line 448 of file FieldInterp.h.
{
return "CubicMACFieldInterp";
}
template<class Data_T >
| static const char* CubicMACFieldInterp< Data_T >::classType | ( | ) | [inline, static] |
Reimplemented from RefBase.
Definition at line 453 of file FieldInterp.h.
{
return CubicMACFieldInterp<Data_T>::ms_classType.name();
}
template<class Data_T >
| Data_T CubicMACFieldInterp< Data_T >::sample | ( | const MACField< Data_T > & | data, |
| const V3d & | vsP | ||
| ) | const |
Definition at line 1205 of file FieldInterp.h.
References FieldRes::dataWindow(), monotonicCubicInterpolant(), MACField< Data_T >::u(), MACField< Data_T >::v(), and MACField< Data_T >::w().
{
typedef typename Data_T::BaseType T;
const Box3i &dataWindow = data.dataWindow();
// Pixel centers are at .5 coordinates
// NOTE: Don't use contToDisc for this, we're looking for sample
// point locations, not coordinate shifts.
Data_T ret;
// X component ---
V3d clampedVsP(std::max(0.5, vsP.x),
std::max(0.5, vsP.y),
std::max(0.5, vsP.z));
FIELD3D_VEC3_T<double> p(vsP.x,
clampedVsP.y - 0.5,
clampedVsP.z - 0.5);
// Lower left corner
V3i c(static_cast<int>(floor(p.x)),
static_cast<int>(floor(p.y)),
static_cast<int>(floor(p.z)));
FIELD3D_VEC3_T<double> t(p - static_cast<FIELD3D_VEC3_T<double> >(c));
{
// Clamp the coordinates
int im, jm, km;
im = std::max(dataWindow.min.x, std::min(c.x, dataWindow.max.x + 1));
jm = std::max(dataWindow.min.y, std::min(c.y, dataWindow.max.y));
km = std::max(dataWindow.min.z, std::min(c.z, dataWindow.max.z));
int im_1, jm_1, km_1;
im_1 = std::max(dataWindow.min.x, std::min(im - 1, dataWindow.max.x + 1));
jm_1 = std::max(dataWindow.min.y, std::min(jm - 1, dataWindow.max.y));
km_1 = std::max(dataWindow.min.z, std::min(km - 1, dataWindow.max.z));
int im1, jm1, km1;
im1 = std::max(dataWindow.min.x, std::min(im + 1, dataWindow.max.x + 1));
jm1 = std::max(dataWindow.min.y, std::min(jm + 1, dataWindow.max.y));
km1 = std::max(dataWindow.min.z, std::min(km + 1, dataWindow.max.z));
int im2, jm2, km2;
im2 = std::max(dataWindow.min.x, std::min(im + 2, dataWindow.max.x + 1));
jm2 = std::max(dataWindow.min.y, std::min(jm + 2, dataWindow.max.y));
km2 = std::max(dataWindow.min.z, std::min(km + 2, dataWindow.max.z));
T z11 = monotonicCubicInterpolant(data.u(im_1, jm_1, km_1),
data.u(im_1, jm_1, km),
data.u(im_1, jm_1, km1),
data.u(im_1, jm_1, km2), t.z);
T z12 = monotonicCubicInterpolant(data.u(im_1, jm, km_1),
data.u(im_1, jm, km),
data.u(im_1, jm, km1),
data.u(im_1, jm, km2), t.z);
T z13 = monotonicCubicInterpolant(data.u(im_1, jm1, km_1),
data.u(im_1, jm1, km),
data.u(im_1, jm1, km1),
data.u(im_1, jm1, km2), t.z);
T z14 = monotonicCubicInterpolant(data.u(im_1, jm2, km_1),
data.u(im_1, jm2, km),
data.u(im_1, jm2, km1),
data.u(im_1, jm2, km2), t.z);
T z21 = monotonicCubicInterpolant(data.u(im, jm_1, km_1),
data.u(im, jm_1, km),
data.u(im, jm_1, km1),
data.u(im, jm_1, km2), t.z);
T z22 = monotonicCubicInterpolant(data.u(im, jm, km_1),
data.u(im, jm, km),
data.u(im, jm, km1),
data.u(im, jm, km2), t.z);
T z23 = monotonicCubicInterpolant(data.u(im, jm1, km_1),
data.u(im, jm1, km),
data.u(im, jm1, km1),
data.u(im, jm1, km2), t.z);
T z24 = monotonicCubicInterpolant(data.u(im, jm2, km_1),
data.u(im, jm2, km),
data.u(im, jm2, km1),
data.u(im, jm2, km2), t.z);
T z31 = monotonicCubicInterpolant(data.u(im1, jm_1, km_1),
data.u(im1, jm_1, km),
data.u(im1, jm_1, km1),
data.u(im1, jm_1, km2), t.z);
T z32 = monotonicCubicInterpolant(data.u(im1, jm, km_1),
data.u(im1, jm, km),
data.u(im1, jm, km1),
data.u(im1, jm, km2), t.z);
T z33 = monotonicCubicInterpolant(data.u(im1, jm1, km_1),
data.u(im1, jm1, km),
data.u(im1, jm1, km1),
data.u(im1, jm1, km2), t.z);
T z34 = monotonicCubicInterpolant(data.u(im1, jm2, km_1),
data.u(im1, jm2, km),
data.u(im1, jm2, km1),
data.u(im1, jm2, km2), t.z);
T z41 = monotonicCubicInterpolant(data.u(im2, jm_1, km_1),
data.u(im2, jm_1, km),
data.u(im2, jm_1, km1),
data.u(im2, jm_1, km2), t.z);
T z42 = monotonicCubicInterpolant(data.u(im2, jm, km_1),
data.u(im2, jm, km),
data.u(im2, jm, km1),
data.u(im2, jm, km2), t.z);
T z43 = monotonicCubicInterpolant(data.u(im2, jm1, km_1),
data.u(im2, jm1, km),
data.u(im2, jm1, km1),
data.u(im2, jm1, km2), t.z);
T z44 = monotonicCubicInterpolant(data.u(im2, jm2, km_1),
data.u(im2, jm2, km),
data.u(im2, jm2, km1),
data.u(im2, jm2, km2), t.z);
T y1 = monotonicCubicInterpolant(z11, z12, z13, z14, t.y);
T y2 = monotonicCubicInterpolant(z21, z22, z23, z24, t.y);
T y3 = monotonicCubicInterpolant(z31, z32, z33, z34, t.y);
T y4 = monotonicCubicInterpolant(z41, z42, z43, z44, t.y);
ret.x = monotonicCubicInterpolant(y1, y2, y3, y4, t.x);
}
// Y component ---
p.setValue(clampedVsP.x - 0.5, vsP.y , clampedVsP.z - 0.5);
// Lower left corner
c.x = static_cast<int>(floor(p.x));
c.y = static_cast<int>(floor(p.y));
c.z = static_cast<int>(floor(p.z));
t.setValue(p - static_cast<FIELD3D_VEC3_T<double> >(c));
{
// Clamp the coordinates
int im, jm, km;
im = std::max(dataWindow.min.x, std::min(c.x, dataWindow.max.x));
jm = std::max(dataWindow.min.y, std::min(c.y, dataWindow.max.y + 1));
km = std::max(dataWindow.min.z, std::min(c.z, dataWindow.max.z));
int im_1, jm_1, km_1;
im_1 = std::max(dataWindow.min.x, std::min(im - 1, dataWindow.max.x));
jm_1 = std::max(dataWindow.min.y, std::min(jm - 1, dataWindow.max.y + 1));
km_1 = std::max(dataWindow.min.z, std::min(km - 1, dataWindow.max.z));
int im1, jm1, km1;
im1 = std::max(dataWindow.min.x, std::min(im + 1, dataWindow.max.x));
jm1 = std::max(dataWindow.min.y, std::min(jm + 1, dataWindow.max.y + 1));
km1 = std::max(dataWindow.min.z, std::min(km + 1, dataWindow.max.z));
int im2, jm2, km2;
im2 = std::max(dataWindow.min.x, std::min(im + 2, dataWindow.max.x));
jm2 = std::max(dataWindow.min.y, std::min(jm + 2, dataWindow.max.y + 1));
km2 = std::max(dataWindow.min.z, std::min(km + 2, dataWindow.max.z));
T z11 = monotonicCubicInterpolant(data.v(im_1, jm_1, km_1),
data.v(im_1, jm_1, km),
data.v(im_1, jm_1, km1),
data.v(im_1, jm_1, km2), t.z);
T z12 = monotonicCubicInterpolant(data.v(im_1, jm, km_1),
data.v(im_1, jm, km),
data.v(im_1, jm, km1),
data.v(im_1, jm, km2), t.z);
T z13 = monotonicCubicInterpolant(data.v(im_1, jm1, km_1),
data.v(im_1, jm1, km),
data.v(im_1, jm1, km1),
data.v(im_1, jm1, km2), t.z);
T z14 = monotonicCubicInterpolant(data.v(im_1, jm2, km_1),
data.v(im_1, jm2, km),
data.v(im_1, jm2, km1),
data.v(im_1, jm2, km2), t.z);
T z21 = monotonicCubicInterpolant(data.v(im, jm_1, km_1),
data.v(im, jm_1, km),
data.v(im, jm_1, km1),
data.v(im, jm_1, km2), t.z);
T z22 = monotonicCubicInterpolant(data.v(im, jm, km_1),
data.v(im, jm, km),
data.v(im, jm, km1),
data.v(im, jm, km2), t.z);
T z23 = monotonicCubicInterpolant(data.v(im, jm1, km_1),
data.v(im, jm1, km),
data.v(im, jm1, km1),
data.v(im, jm1, km2), t.z);
T z24 = monotonicCubicInterpolant(data.v(im, jm2, km_1),
data.v(im, jm2, km),
data.v(im, jm2, km1),
data.v(im, jm2, km2), t.z);
T z31 = monotonicCubicInterpolant(data.v(im1, jm_1, km_1),
data.v(im1, jm_1, km),
data.v(im1, jm_1, km1),
data.v(im1, jm_1, km2), t.z);
T z32 = monotonicCubicInterpolant(data.v(im1, jm, km_1),
data.v(im1, jm, km),
data.v(im1, jm, km1),
data.v(im1, jm, km2), t.z);
T z33 = monotonicCubicInterpolant(data.v(im1, jm1, km_1),
data.v(im1, jm1, km),
data.v(im1, jm1, km1),
data.v(im1, jm1, km2), t.z);
T z34 = monotonicCubicInterpolant(data.v(im1, jm2, km_1),
data.v(im1, jm2, km),
data.v(im1, jm2, km1),
data.v(im1, jm2, km2), t.z);
T z41 = monotonicCubicInterpolant(data.v(im2, jm_1, km_1),
data.v(im2, jm_1, km),
data.v(im2, jm_1, km1),
data.v(im2, jm_1, km2), t.z);
T z42 = monotonicCubicInterpolant(data.v(im2, jm, km_1),
data.v(im2, jm, km),
data.v(im2, jm, km1),
data.v(im2, jm, km2), t.z);
T z43 = monotonicCubicInterpolant(data.v(im2, jm1, km_1),
data.v(im2, jm1, km),
data.v(im2, jm1, km1),
data.v(im2, jm1, km2), t.z);
T z44 = monotonicCubicInterpolant(data.v(im2, jm2, km_1),
data.v(im2, jm2, km),
data.v(im2, jm2, km1),
data.v(im2, jm2, km2), t.z);
T y1 = monotonicCubicInterpolant(z11, z12, z13, z14, t.y);
T y2 = monotonicCubicInterpolant(z21, z22, z23, z24, t.y);
T y3 = monotonicCubicInterpolant(z31, z32, z33, z34, t.y);
T y4 = monotonicCubicInterpolant(z41, z42, z43, z44, t.y);
ret.y = monotonicCubicInterpolant(y1, y2, y3, y4, t.x);
}
// Z component ---
p.setValue(clampedVsP.x - 0.5 , clampedVsP.y - 0.5, vsP.z);
// Lower left corner
c.x = static_cast<int>(floor(p.x));
c.y = static_cast<int>(floor(p.y));
c.z = static_cast<int>(floor(p.z));
t.setValue(p - static_cast<FIELD3D_VEC3_T<double> >(c));
{
// Clamp the coordinates
int im, jm, km;
im = std::max(dataWindow.min.x, std::min(c.x, dataWindow.max.x));
jm = std::max(dataWindow.min.y, std::min(c.y, dataWindow.max.y));
km = std::max(dataWindow.min.z, std::min(c.z, dataWindow.max.z + 1));
int im_1, jm_1, km_1;
im_1 = std::max(dataWindow.min.x, std::min(im - 1, dataWindow.max.x));
jm_1 = std::max(dataWindow.min.y, std::min(jm - 1, dataWindow.max.y));
km_1 = std::max(dataWindow.min.z, std::min(km - 1, dataWindow.max.z + 1));
int im1, jm1, km1;
im1 = std::max(dataWindow.min.x, std::min(im + 1, dataWindow.max.x));
jm1 = std::max(dataWindow.min.y, std::min(jm + 1, dataWindow.max.y));
km1 = std::max(dataWindow.min.z, std::min(km + 1, dataWindow.max.z + 1));
int im2, jm2, km2;
im2 = std::max(dataWindow.min.x, std::min(im + 2, dataWindow.max.x));
jm2 = std::max(dataWindow.min.y, std::min(jm + 2, dataWindow.max.y));
km2 = std::max(dataWindow.min.z, std::min(km + 2, dataWindow.max.z + 1));
T z11 = monotonicCubicInterpolant(data.w(im_1, jm_1, km_1),
data.w(im_1, jm_1, km),
data.w(im_1, jm_1, km1),
data.w(im_1, jm_1, km2), t.z);
T z12 = monotonicCubicInterpolant(data.w(im_1, jm, km_1),
data.w(im_1, jm, km),
data.w(im_1, jm, km1),
data.w(im_1, jm, km2), t.z);
T z13 = monotonicCubicInterpolant(data.w(im_1, jm1, km_1),
data.w(im_1, jm1, km),
data.w(im_1, jm1, km1),
data.w(im_1, jm1, km2), t.z);
T z14 = monotonicCubicInterpolant(data.w(im_1, jm2, km_1),
data.w(im_1, jm2, km),
data.w(im_1, jm2, km1),
data.w(im_1, jm2, km2), t.z);
T z21 = monotonicCubicInterpolant(data.w(im, jm_1, km_1),
data.w(im, jm_1, km),
data.w(im, jm_1, km1),
data.w(im, jm_1, km2), t.z);
T z22 = monotonicCubicInterpolant(data.w(im, jm, km_1),
data.w(im, jm, km),
data.w(im, jm, km1),
data.w(im, jm, km2), t.z);
T z23 = monotonicCubicInterpolant(data.w(im, jm1, km_1),
data.w(im, jm1, km),
data.w(im, jm1, km1),
data.w(im, jm1, km2), t.z);
T z24 = monotonicCubicInterpolant(data.w(im, jm2, km_1),
data.w(im, jm2, km),
data.w(im, jm2, km1),
data.w(im, jm2, km2), t.z);
T z31 = monotonicCubicInterpolant(data.w(im1, jm_1, km_1),
data.w(im1, jm_1, km),
data.w(im1, jm_1, km1),
data.w(im1, jm_1, km2), t.z);
T z32 = monotonicCubicInterpolant(data.w(im1, jm, km_1),
data.w(im1, jm, km),
data.w(im1, jm, km1),
data.w(im1, jm, km2), t.z);
T z33 = monotonicCubicInterpolant(data.w(im1, jm1, km_1),
data.w(im1, jm1, km),
data.w(im1, jm1, km1),
data.w(im1, jm1, km2), t.z);
T z34 = monotonicCubicInterpolant(data.w(im1, jm2, km_1),
data.w(im1, jm2, km),
data.w(im1, jm2, km1),
data.w(im1, jm2, km2), t.z);
T z41 = monotonicCubicInterpolant(data.w(im2, jm_1, km_1),
data.w(im2, jm_1, km),
data.w(im2, jm_1, km1),
data.w(im2, jm_1, km2), t.z);
T z42 = monotonicCubicInterpolant(data.w(im2, jm, km_1),
data.w(im2, jm, km),
data.w(im2, jm, km1),
data.w(im2, jm, km2), t.z);
T z43 = monotonicCubicInterpolant(data.w(im2, jm1, km_1),
data.w(im2, jm1, km),
data.w(im2, jm1, km1),
data.w(im2, jm1, km2), t.z);
T z44 = monotonicCubicInterpolant(data.w(im2, jm2, km_1),
data.w(im2, jm2, km),
data.w(im2, jm2, km1),
data.w(im2, jm2, km2), t.z);
T y1 = monotonicCubicInterpolant(z11, z12, z13, z14, t.y);
T y2 = monotonicCubicInterpolant(z21, z22, z23, z24, t.y);
T y3 = monotonicCubicInterpolant(z31, z32, z33, z34, t.y);
T y4 = monotonicCubicInterpolant(z41, z42, z43, z44, t.y);
ret.z = monotonicCubicInterpolant(y1, y2, y3, y4, t.x);
}
return ret;
}
Member Data Documentation
template<class Data_T >
| CubicMACFieldInterp< Data_T >::DEFINE_FIELD_RTTI_CONCRETE_CLASS |
Definition at line 446 of file FieldInterp.h.
template<class Data_T >
TemplatedFieldType<CubicMACFieldInterp<Data_T> > CubicMACFieldInterp< Data_T >::ms_classType [static, private] |
Definition at line 466 of file FieldInterp.h.
The documentation for this class was generated from the following file:
- export/FieldInterp.h
