bool NX::ViewFrustum::Visible(const NX::Sphere &sphere, const unsigned int mask){
const NX::vector<float, 3> center = sphere.GetCenter();
const float r = sphere.GetRadius();
if(mask & NX::VF_VT_LEFT){
if(NX::Dot(GetLeftPlane().GetNormal(), center) + GetLeftPlane().GetDistFromOriginal() < -r){
return false;
}
}
if(mask & NX::VF_VT_RIGHT){
if(NX::Dot(GetRightPlane().GetNormal(), center) + GetRightPlane().GetDistFromOriginal() < -r){
return false;
}
}
if(mask & NX::VF_VT_TOP){
if(NX::Dot(GetTopPlane().GetNormal(), center) + GetTopPlane().GetDistFromOriginal() < -r){
return false;
}
}
if(mask & NX::VF_VT_BOTTOM){
if(NX::Dot(GetBottomPlane().GetNormal(), center) + GetBottomPlane().GetDistFromOriginal() < -r){
return false;
}
}
if(mask & NX::VF_VT_FRONT){
if(NX::Dot(GetFrontPlane().GetNormal(), center) + GetFrontPlane().GetDistFromOriginal() < -r){
return false;
}
}
if(mask & NX::VF_VT_BACK){
if(NX::Dot(GetBackPlane().GetNormal(), center) + GetBackPlane().GetDistFromOriginal() < -r){
return false;
}
}
return true;
}
bool NX::ViewFrustum::Visible(const NX::Circle &circle, const unsigned int mask){
const NX::vector<float, 3> normal = circle.GetNormal();
const float r = circle.GetRadius();
#undef NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST
#define NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(plane) \
const float RRef = r * NX::NXAbs(NX::Dot(plane.GetNormal(), normal));\
if(NX::Dot(plane.GetNormal(), circle.GetCenter()) + plane.GetDistFromOriginal() < - RRef)\
return false;
if(mask & NX::VF_VT_LEFT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetLeftPlane());
}
if(mask & NX::VF_VT_RIGHT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetRightPlane());
}
if(mask & NX::VF_VT_TOP){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetTopPlane());
}
if(mask & NX::VF_VT_BOTTOM){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetBottomPlane());
}
if(mask & NX::VF_VT_FRONT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetFrontPlane());
}
if(mask & NX::VF_VT_BACK){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetBackPlane());
}
return true;
#undef NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST
}
std::vector<vgl_plane_3d<double> > GetAllPlanes(const vgl_box_3d<double> &B)
{
std::vector<vgl_plane_3d<double> > Planes(6);
Planes[0] = GetFrontPlane(B);
Planes[1] = GetBackPlane(B);
Planes[2] = GetLeftPlane(B);
Planes[3] = GetRightPlane(B);
Planes[4] = GetTopPlane(B);
Planes[5] = GetBottomPlane(B);
return Planes;
}
bool NX::ViewFrustum::Visible(const NX::Ellipsoid &ellipsoid, const unsigned int mask){
const NX::vector<float, 3> R = ellipsoid.GetAxisX() * ellipsoid.GetAxisXLength();
const NX::vector<float, 3> S = ellipsoid.GetAxisY() * ellipsoid.GetAxisYLength();
const NX::vector<float, 3> T = ellipsoid.GetAxisZ() * ellipsoid.GetAxisZLength();
#undef NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST
#define NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(plane) \
const NX::vector<float, 3> N = plane.GetNormal();\
const float r = NX::Dot(R, N);\
const float s = NX::Dot(S, N);\
const float t = NX::Dot(T, N);\
const float RRef = std::sqrt(r * r + s * s + t * t);\
if(NX::Dot(N, ellipsoid.GetCenter()) + plane.GetDistFromOriginal() <= -RRef){\
return false;\
}
if(mask & NX::VF_VT_LEFT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetLeftPlane());
}
if(mask & NX::VF_VT_RIGHT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetRightPlane());
}
if(mask & NX::VF_VT_FRONT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetFrontPlane());
}
if(mask & NX::VF_VT_BACK){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetBackPlane());
}
if(mask & NX::VF_VT_TOP){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetTopPlane());
}
if(mask & NX::VF_VT_BOTTOM){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetBottomPlane());
}
return true;
#undef NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST
}
bool NX::ViewFrustum::Visible(const NX::Ellipse &ellipse, const unsigned int mask){
const NX::vector<float, 3> R = ellipse.GetLongAxis() * ellipse.GetLongAxisLength();
const NX::vector<float, 3> S = ellipse.GetShortAxis() * ellipse.GetShortAxisLength();
const NX::vector<float, 3> C = ellipse.GetCenter();
#undef NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST
#define NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(plane) \
const NX::vector<float, 3> N = plane.GetNormal();\
const float rn = NX::Dot(R, N);\
const float rs = NX::Dot(S, N);\
const float RRef = std::sqrt(rn * rn + rs * rs);\
if(NX::Dot(N, C) + plane.GetDistFromOriginal() < -RRef){\
return false;\
}
if(mask & NX::VF_VT_LEFT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetLeftPlane());
}
if(mask & NX::VF_VT_RIGHT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetRightPlane());
}
if(mask & NX::VF_VT_FRONT){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetFrontPlane());
}
if(mask & NX::VF_VT_BACK){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetBackPlane());
}
if(mask & NX::VF_VT_TOP){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetTopPlane());
}
if(mask & NX::VF_VT_BOTTOM){
NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST(GetBottomPlane());
}
return true;
#undef NX_VIEWFRUSTUM_POSITIVE_SIDE_TEST
}