|
@@ -111,7 +111,7 @@ std::string Operation3D::writeInfos()
|
|
//############################## region feature ###############################
|
|
//############################## region feature ###############################
|
|
|
|
|
|
|
|
|
|
-double RegionFeat::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RegionFeat3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
return (*feats.rProbs)[(*feats.feats)(x,y,z,channel1)][channel2];
|
|
return (*feats.rProbs)[(*feats.feats)(x,y,z,channel1)][channel2];
|
|
}
|
|
}
|
|
@@ -122,7 +122,7 @@ double RegionFeat::getVal ( const Features &feats, const int &x, const int &y, c
|
|
//############################# two gray values ###############################
|
|
//############################# two gray values ###############################
|
|
|
|
|
|
|
|
|
|
-double Minus::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double Minus3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
@@ -133,7 +133,7 @@ double Minus::getVal ( const Features &feats, const int &x, const int &y, const
|
|
return v1-v2;
|
|
return v1-v2;
|
|
}
|
|
}
|
|
|
|
|
|
-double MinusAbs::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double MinusAbs3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
@@ -144,7 +144,7 @@ double MinusAbs::getVal ( const Features &feats, const int &x, const int &y, con
|
|
return abs(v1-v2);
|
|
return abs(v1-v2);
|
|
}
|
|
}
|
|
|
|
|
|
-double Addition::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double Addition3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
@@ -161,7 +161,7 @@ double Addition::getVal ( const Features &feats, const int &x, const int &y, con
|
|
//############################## one gray value ###############################
|
|
//############################## one gray value ###############################
|
|
|
|
|
|
|
|
|
|
-double Only1::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double Only13D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
@@ -175,21 +175,21 @@ double Only1::getVal ( const Features &feats, const int &x, const int &y, const
|
|
//############################ relative positions #############################
|
|
//############################ relative positions #############################
|
|
|
|
|
|
|
|
|
|
-double RelativeXPosition::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RelativeXPosition3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
return ( double ) x / ( double ) xsize;
|
|
return ( double ) x / ( double ) xsize;
|
|
}
|
|
}
|
|
|
|
|
|
-double RelativeYPosition::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RelativeYPosition3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
return ( double ) y / ( double ) ysize;
|
|
return ( double ) y / ( double ) ysize;
|
|
}
|
|
}
|
|
|
|
|
|
-double RelativeZPosition::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RelativeZPosition3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
@@ -204,7 +204,7 @@ double RelativeZPosition::getVal ( const Features &feats, const int &x, const in
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-double GlobalFeats::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double GlobalFeats3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int xsize, ysize, zsize;
|
|
int xsize, ysize, zsize;
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
getXYZ ( feats, xsize, ysize, zsize );
|
|
@@ -212,7 +212,7 @@ double GlobalFeats::getVal ( const Features &feats, const int &x, const int &y,
|
|
return feats.feats->getIntegralValue( 0, 0, 0, xsize - 1, ysize - 1, zsize - 1, channel1 );
|
|
return feats.feats->getIntegralValue( 0, 0, 0, xsize - 1, ysize - 1, zsize - 1, channel1 );
|
|
}
|
|
}
|
|
|
|
|
|
-void IntegralOps::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int _z2, int _channel1, int _channel2, int _featType )
|
|
|
|
|
|
+void IntegralOps3D::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int _z2, int _channel1, int _channel2, int _featType )
|
|
{
|
|
{
|
|
x1 = std::min ( _x1, _x2 );
|
|
x1 = std::min ( _x1, _x2 );
|
|
y1 = std::min ( _y1, _y2 );
|
|
y1 = std::min ( _y1, _y2 );
|
|
@@ -226,12 +226,12 @@ void IntegralOps::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int _z2, in
|
|
init = true;
|
|
init = true;
|
|
}
|
|
}
|
|
|
|
|
|
-double IntegralOps::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double IntegralOps3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
return feats.feats->getIntegralValue(x + x1, y + y1, z + z1, x + x2, y + y2, z + z2, channel1);
|
|
return feats.feats->getIntegralValue(x + x1, y + y1, z + z1, x + x2, y + y2, z + z2, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-void IntegralCenteredOps::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int _z2, int _channel1, int _channel2, int _featType )
|
|
|
|
|
|
+void IntegralCenteredOps3D::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int _z2, int _channel1, int _channel2, int _featType )
|
|
{
|
|
{
|
|
x1 = abs ( _x1 );
|
|
x1 = abs ( _x1 );
|
|
y1 = abs ( _y1 );
|
|
y1 = abs ( _y1 );
|
|
@@ -245,12 +245,12 @@ void IntegralCenteredOps::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int
|
|
init = true;
|
|
init = true;
|
|
}
|
|
}
|
|
|
|
|
|
-double IntegralCenteredOps::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double IntegralCenteredOps3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
return feats.feats->getIntegralValue(x - x1, y - y1, z - z1, x + x1, y + y1, z + z1, channel1);
|
|
return feats.feats->getIntegralValue(x - x1, y - y1, z - z1, x + x1, y + y1, z + z1, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-void BiIntegralCenteredOps::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int _z2, int _channel1, int _channel2, int _featType )
|
|
|
|
|
|
+void BiIntegralCenteredOps3D::set ( int _x1, int _y1, int _z1, int _x2, int _y2, int _z2, int _channel1, int _channel2, int _featType )
|
|
{
|
|
{
|
|
x1 = std::min ( abs ( _x1 ), abs ( _x2 ) );
|
|
x1 = std::min ( abs ( _x1 ), abs ( _x2 ) );
|
|
y1 = std::min ( abs ( _y1 ), abs ( _y2 ) );
|
|
y1 = std::min ( abs ( _y1 ), abs ( _y2 ) );
|
|
@@ -264,7 +264,7 @@ void BiIntegralCenteredOps::set ( int _x1, int _y1, int _z1, int _x2, int _y2, i
|
|
init = true;
|
|
init = true;
|
|
}
|
|
}
|
|
|
|
|
|
-double BiIntegralCenteredOps::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double BiIntegralCenteredOps3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
return feats.feats->getIntegralValue(x - x1, y - y1, z - z1, x + x1, y + y1, z + z1, channel1 ) - feats.feats->getIntegralValue(x - x2, y - y2, z - z2, x + x2, y + y2, z + z2, channel1);
|
|
return feats.feats->getIntegralValue(x - x1, y - y1, z - z1, x + x1, y + y1, z + z1, channel1 ) - feats.feats->getIntegralValue(x - x2, y - y2, z - z2, x + x2, y + y2, z + z2, channel1);
|
|
}
|
|
}
|
|
@@ -278,7 +278,7 @@ double BiIntegralCenteredOps::getVal ( const Features &feats, const int &x, cons
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-double HaarHorizontal::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double HaarHorizontal3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -290,7 +290,7 @@ double HaarHorizontal::getVal ( const Features &feats, const int &x, const int &
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, y, lrz, channel1 ) - feats.feats->getIntegralValue(tlx, y, tlz, lrx, lry, lrz, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, y, lrz, channel1 ) - feats.feats->getIntegralValue(tlx, y, tlz, lrx, lry, lrz, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double HaarVertical::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double HaarVertical3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -302,7 +302,7 @@ double HaarVertical::getVal ( const Features &feats, const int &x, const int &y,
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, x, lry, lrz, channel1) - feats.feats->getIntegralValue(x, tly, tlz, lrx, lry, lrz, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, x, lry, lrz, channel1) - feats.feats->getIntegralValue(x, tly, tlz, lrx, lry, lrz, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double HaarStacked::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double HaarStacked3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -314,7 +314,7 @@ double HaarStacked::getVal ( const Features &feats, const int &x, const int &y,
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, lry, z, channel1) - feats.feats->getIntegralValue(tlx, tly, z, lrx, lry, lrz, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, lry, z, channel1) - feats.feats->getIntegralValue(tlx, tly, z, lrx, lry, lrz, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double HaarDiagXY::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double HaarDiagXY3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -326,7 +326,7 @@ double HaarDiagXY::getVal ( const Features &feats, const int &x, const int &y, c
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, x, y, lrz, channel1) + feats.feats->getIntegralValue(x, y, tlz, lrx, lry, lrz, channel1) - feats.feats->getIntegralValue(tlx, y, tlz, x, lry, lrz, channel1) - feats.feats->getIntegralValue(x, tly, tlz, lrx, y, lrz, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, x, y, lrz, channel1) + feats.feats->getIntegralValue(x, y, tlz, lrx, lry, lrz, channel1) - feats.feats->getIntegralValue(tlx, y, tlz, x, lry, lrz, channel1) - feats.feats->getIntegralValue(x, tly, tlz, lrx, y, lrz, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double HaarDiagXZ::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double HaarDiagXZ3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -338,7 +338,7 @@ double HaarDiagXZ::getVal ( const Features &feats, const int &x, const int &y, c
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, x, lry, z, channel1) + feats.feats->getIntegralValue(x, tly, z, lrx, lry, lrz, channel1) - feats.feats->getIntegralValue(tlx, tly, z, x, lry, lrz, channel1) - feats.feats->getIntegralValue(x, tly, tlz, lrx, lry, z, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, x, lry, z, channel1) + feats.feats->getIntegralValue(x, tly, z, lrx, lry, lrz, channel1) - feats.feats->getIntegralValue(tlx, tly, z, x, lry, lrz, channel1) - feats.feats->getIntegralValue(x, tly, tlz, lrx, lry, z, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double HaarDiagYZ::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double HaarDiagYZ3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -350,7 +350,7 @@ double HaarDiagYZ::getVal ( const Features &feats, const int &x, const int &y, c
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, y, z, channel1) + feats.feats->getIntegralValue(tlx, y, z, lrx, lry, lrz, channel1) - feats.feats->getIntegralValue(tlx, tly, z, lrx, y, lrz, channel1) - feats.feats->getIntegralValue(tlx, y, tlz, lrx, lry, z, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, y, z, channel1) + feats.feats->getIntegralValue(tlx, y, z, lrx, lry, lrz, channel1) - feats.feats->getIntegralValue(tlx, tly, z, lrx, y, lrz, channel1) - feats.feats->getIntegralValue(tlx, y, tlz, lrx, lry, z, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double Haar3Horiz::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double Haar3Horiz3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -364,7 +364,7 @@ double Haar3Horiz::getVal ( const Features &feats, const int &x, const int &y, c
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, mtly, lrz, channel1) - feats.feats->getIntegralValue(tlx, mtly, tlz, lrx, mlry, lrz, channel1) + feats.feats->getIntegralValue(tlx, mlry, tlz, lrx, lry, lrz, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, lrx, mtly, lrz, channel1) - feats.feats->getIntegralValue(tlx, mtly, tlz, lrx, mlry, lrz, channel1) + feats.feats->getIntegralValue(tlx, mlry, tlz, lrx, lry, lrz, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double Haar3Vert::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double Haar3Vert3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -378,7 +378,7 @@ double Haar3Vert::getVal ( const Features &feats, const int &x, const int &y, co
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, mtlx, lry, lrz, channel1) - feats.feats->getIntegralValue(mtlx, tly, tlz, mlrx, lry, lrz, channel1) + feats.feats->getIntegralValue(mlrx, tly, tlz, lrx, lry, lrz, channel1);
|
|
return feats.feats->getIntegralValue(tlx, tly, tlz, mtlx, lry, lrz, channel1) - feats.feats->getIntegralValue(mtlx, tly, tlz, mlrx, lry, lrz, channel1) + feats.feats->getIntegralValue(mlrx, tly, tlz, lrx, lry, lrz, channel1);
|
|
}
|
|
}
|
|
|
|
|
|
-double Haar3Stack::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double Haar3Stack3D::getVal ( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
int tlx = x - abs(x1);
|
|
int tlx = x - abs(x1);
|
|
int tly = y - abs(y1);
|
|
int tly = y - abs(y1);
|
|
@@ -396,7 +396,7 @@ double Haar3Stack::getVal ( const Features &feats, const int &x, const int &y, c
|
|
|
|
|
|
|
|
|
|
//############################### Ray Features ################################
|
|
//############################### Ray Features ################################
|
|
-double RayDiff::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RayDiff3D::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
|
|
|
|
double dist1 = feats.feats->get( x, y, z, feats.feats->channels()-24+z1 );
|
|
double dist1 = feats.feats->get( x, y, z, feats.feats->channels()-24+z1 );
|
|
@@ -408,17 +408,17 @@ double RayDiff::getVal( const Features &feats, const int &x, const int &y, const
|
|
return 0.0;
|
|
return 0.0;
|
|
}
|
|
}
|
|
|
|
|
|
-double RayDist::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RayDist3D::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
return feats.feats->get( x, y, z, feats.feats->channels()-24+z1 );
|
|
return feats.feats->get( x, y, z, feats.feats->channels()-24+z1 );
|
|
}
|
|
}
|
|
|
|
|
|
-double RayNorm::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RayNorm3D::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
return feats.feats->get( x, y, z, feats.feats->channels()-16+z1 );
|
|
return feats.feats->get( x, y, z, feats.feats->channels()-16+z1 );
|
|
}
|
|
}
|
|
|
|
|
|
-double RayOrient::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
|
|
|
|
+double RayOrient3D::getVal( const Features &feats, const int &x, const int &y, const int &z )
|
|
{
|
|
{
|
|
return feats.feats->get( x, y, z, feats.feats->channels()-8+z1 );
|
|
return feats.feats->get( x, y, z, feats.feats->channels()-8+z1 );
|
|
}
|
|
}
|