BlobOperators.cpp

Go to the documentation of this file.

#include <limits.h>
#include "hrl_cvblobslib/BlobOperators.h"

/***************************************************************************
  Implementació de les classes per al càlcul de característiques sobre el blob

  Implementation of the helper classes to perform operations on blobs
/**************************************************************************/

double CBlobGetMoment::operator()(CBlob &blob)
{
        return blob.Moment(m_p, m_q);
}

double CBlobGetHullPerimeter::operator()(CBlob &blob)
{
        CvSeq *convexHull;
        double perimeter;

        convexHull = blob.GetConvexHull();

        if( convexHull )
                perimeter = fabs(cvArcLength(convexHull,CV_WHOLE_SEQ,1));
        else
                return 0;

        cvClearSeq(convexHull);

        return perimeter;
}

double CBlobGetHullArea::operator()(CBlob &blob)
{
        CvSeq *convexHull;
        double area;

        convexHull = blob.GetConvexHull();

        if( convexHull )
                area = fabs(cvContourArea(convexHull));
        else
                return 0;

        cvClearSeq(convexHull);

        return area;
}

double CBlobGetMinXatMinY::operator()(CBlob &blob)
{
        double result = LONG_MAX;
        
        CvSeqReader reader;
        CvPoint actualPoint;
        t_PointList externContour;
        
        externContour = blob.GetExternalContour()->GetContourPoints();
        if( !externContour ) return result;
        cvStartReadSeq( externContour, &reader);

        for( int i=0; i< externContour->total; i++)
        {
                CV_READ_SEQ_ELEM( actualPoint, reader);

                if( (actualPoint.y == blob.MinY()) && (actualPoint.x < result) )
                {
                        result = actualPoint.x;
                }       
        }

        return result;
}

double CBlobGetMinYatMaxX::operator()(CBlob &blob)
{
        double result = LONG_MAX;
        
        CvSeqReader reader;
        CvPoint actualPoint;
        t_PointList externContour;
        
        externContour = blob.GetExternalContour()->GetContourPoints();
        if( !externContour ) return result;
        cvStartReadSeq( externContour, &reader);

        for( int i=0; i< externContour->total; i++)
        {
                CV_READ_SEQ_ELEM( actualPoint, reader);

                if( (actualPoint.x == blob.MaxX()) && (actualPoint.y < result) )
                {
                        result = actualPoint.y;
                }       
        }

        return result;
}

double CBlobGetMaxXatMaxY::operator()(CBlob &blob)
{
        double result = LONG_MIN;
        
        CvSeqReader reader;
        CvPoint actualPoint;
        t_PointList externContour;
        
        externContour = blob.GetExternalContour()->GetContourPoints();
        if( !externContour ) return result;

        cvStartReadSeq( externContour, &reader);

        for( int i=0; i< externContour->total; i++)
        {
                CV_READ_SEQ_ELEM( actualPoint, reader);

                if( (actualPoint.y == blob.MaxY()) && (actualPoint.x > result) )
                {
                        result = actualPoint.x;
                }       
        }

        return result;
}

double CBlobGetMaxYatMinX::operator()(CBlob &blob)
{
        double result = LONG_MIN;
        
        CvSeqReader reader;
        CvPoint actualPoint;
        t_PointList externContour;
        
        externContour = blob.GetExternalContour()->GetContourPoints();
        if( !externContour ) return result;

        cvStartReadSeq( externContour, &reader);

        
        for( int i=0; i< externContour->total; i++)
        {
                CV_READ_SEQ_ELEM( actualPoint, reader);

                if( (actualPoint.x == blob.MinX()) && (actualPoint.y > result) )
                {
                        result = actualPoint.y;
                }       
        }

        return result;
}


double CBlobGetElongation::operator()(CBlob &blob)
{
        double ampladaC,longitudC,amplada,longitud;

        double tmp;

        tmp = blob.Perimeter()*blob.Perimeter() - 16*blob.Area();

        if( tmp > 0.0 )
                ampladaC = (double) (blob.Perimeter()+sqrt(tmp))/4;
        // error intrínsec en els càlculs de l'àrea i el perímetre 
        else
                ampladaC = (double) (blob.Perimeter())/4;

        if(ampladaC<=0.0) return 0;
        longitudC=(double) blob.Area()/ampladaC;

        longitud=MAX( longitudC , ampladaC );
        amplada=MIN( longitudC , ampladaC );

        return (double) longitud/amplada;
}

double CBlobGetCompactness::operator()(CBlob &blob)
{
        if( blob.Area() != 0.0 )
                return (double) pow(blob.Perimeter(),2)/(4*CV_PI*blob.Area());
        else
                return 0.0;
}

double CBlobGetRoughness::operator()(CBlob &blob)
{
        CBlobGetHullPerimeter getHullPerimeter = CBlobGetHullPerimeter();
        
        double hullPerimeter = getHullPerimeter(blob);

        if( hullPerimeter != 0.0 )
                return blob.Perimeter() / hullPerimeter;//HullPerimeter();

        return 0.0;
}

double CBlobGetLength::operator()(CBlob &blob)
{
        double ampladaC,longitudC;
        double tmp;

        tmp = blob.Perimeter()*blob.Perimeter() - 16*blob.Area();

        if( tmp > 0.0 )
                ampladaC = (double) (blob.Perimeter()+sqrt(tmp))/4;
        // error intrínsec en els càlculs de l'àrea i el perímetre 
        else
                ampladaC = (double) (blob.Perimeter())/4;

        if(ampladaC<=0.0) return 0;
        longitudC=(double) blob.Area()/ampladaC;

        return MAX( longitudC , ampladaC );
}

double CBlobGetBreadth::operator()(CBlob &blob)
{
        double ampladaC,longitudC;
        double tmp;

        tmp = blob.Perimeter()*blob.Perimeter() - 16*blob.Area();

        if( tmp > 0.0 )
                ampladaC = (double) (blob.Perimeter()+sqrt(tmp))/4;
        // error intrínsec en els càlculs de l'àrea i el perímetre 
        else
                ampladaC = (double) (blob.Perimeter())/4;

        if(ampladaC<=0.0) return 0;
        longitudC = (double) blob.Area()/ampladaC;

        return MIN( longitudC , ampladaC );
}

double CBlobGetDistanceFromPoint::operator()(CBlob &blob)
{
        double xmitjana, ymitjana;
        CBlobGetXCenter getXCenter;
        CBlobGetYCenter getYCenter;

        xmitjana = m_x - getXCenter( blob );
        ymitjana = m_y - getYCenter( blob );

        return sqrt((xmitjana*xmitjana)+(ymitjana*ymitjana));
}

double CBlobGetXYInside::operator()(CBlob &blob)
{
        if( blob.GetExternalContour()->GetContourPoints() )
        {
                return cvPointPolygonTest( blob.GetExternalContour()->GetContourPoints(), m_p,0) >= 0;
        }

        return 0;
}
#ifdef BLOB_OBJECT_FACTORY

void RegistraTotsOperadors( t_OperadorBlobFactory &fabricaOperadorsBlob )
{
        // blob shape
        fabricaOperadorsBlob.Register( CBlobGetArea().GetNom(), Type2Type<CBlobGetArea>());
        fabricaOperadorsBlob.Register( CBlobGetBreadth().GetNom(), Type2Type<CBlobGetBreadth>());
        fabricaOperadorsBlob.Register( CBlobGetCompactness().GetNom(), Type2Type<CBlobGetCompactness>());
        fabricaOperadorsBlob.Register( CBlobGetElongation().GetNom(), Type2Type<CBlobGetElongation>());
        fabricaOperadorsBlob.Register( CBlobGetExterior().GetNom(), Type2Type<CBlobGetExterior>());
        fabricaOperadorsBlob.Register( CBlobGetLength().GetNom(), Type2Type<CBlobGetLength>());
        fabricaOperadorsBlob.Register( CBlobGetPerimeter().GetNom(), Type2Type<CBlobGetPerimeter>());
        fabricaOperadorsBlob.Register( CBlobGetRoughness().GetNom(), Type2Type<CBlobGetRoughness>());

        // blob color
        fabricaOperadorsBlob.Register( CBlobGetMean(NULL).GetNom(), Type2Type<CBlobGetMean>());
        fabricaOperadorsBlob.Register( CBlobGetStdDev(NULL).GetNom(), Type2Type<CBlobGetStdDev>());

        // extern pixels
        fabricaOperadorsBlob.Register( CBlobGetExternPerimeterRatio().GetNom(), Type2Type<CBlobGetExternPerimeterRatio>());
        fabricaOperadorsBlob.Register( CBlobGetExternHullPerimeterRatio().GetNom(), Type2Type<CBlobGetExternHullPerimeterRatio>());
        fabricaOperadorsBlob.Register( CBlobGetExternPerimeter().GetNom(), Type2Type<CBlobGetExternPerimeter>());
        

        // hull 
        fabricaOperadorsBlob.Register( CBlobGetHullPerimeter().GetNom(), Type2Type<CBlobGetHullPerimeter>());
        fabricaOperadorsBlob.Register( CBlobGetHullArea().GetNom(), Type2Type<CBlobGetHullArea>());
        

        // elipse info
        fabricaOperadorsBlob.Register( CBlobGetMajorAxisLength().GetNom(), Type2Type<CBlobGetMajorAxisLength>());
        fabricaOperadorsBlob.Register( CBlobGetMinorAxisLength().GetNom(), Type2Type<CBlobGetMinorAxisLength>());
        fabricaOperadorsBlob.Register( CBlobGetAxisRatio().GetNom(), Type2Type<CBlobGetAxisRatio>());
        fabricaOperadorsBlob.Register( CBlobGetOrientation().GetNom(), Type2Type<CBlobGetOrientation>());
        fabricaOperadorsBlob.Register( CBlobGetOrientationCos().GetNom(), Type2Type<CBlobGetOrientationCos>());
        fabricaOperadorsBlob.Register( CBlobGetAreaElipseRatio().GetNom(), Type2Type<CBlobGetAreaElipseRatio>());

        // min an max
        fabricaOperadorsBlob.Register( CBlobGetMaxX().GetNom(), Type2Type<CBlobGetMaxX>());
        fabricaOperadorsBlob.Register( CBlobGetMaxY().GetNom(), Type2Type<CBlobGetMaxY>());
        fabricaOperadorsBlob.Register( CBlobGetMinX().GetNom(), Type2Type<CBlobGetMinX>());
        fabricaOperadorsBlob.Register( CBlobGetMinY().GetNom(), Type2Type<CBlobGetMinY>());

        fabricaOperadorsBlob.Register( CBlobGetMaxXatMaxY().GetNom(), Type2Type<CBlobGetMaxXatMaxY>());
        fabricaOperadorsBlob.Register( CBlobGetMaxYatMinX().GetNom(), Type2Type<CBlobGetMaxYatMinX>());
        fabricaOperadorsBlob.Register( CBlobGetMinXatMinY().GetNom(), Type2Type<CBlobGetMinXatMinY>());
        fabricaOperadorsBlob.Register( CBlobGetMinYatMaxX().GetNom(), Type2Type<CBlobGetMinYatMaxX>());

        // coordinate info
        fabricaOperadorsBlob.Register( CBlobGetXYInside().GetNom(), Type2Type<CBlobGetXYInside>());
        fabricaOperadorsBlob.Register( CBlobGetDiffY().GetNom(), Type2Type<CBlobGetDiffY>());
        fabricaOperadorsBlob.Register( CBlobGetDiffX().GetNom(), Type2Type<CBlobGetDiffX>());
        fabricaOperadorsBlob.Register( CBlobGetXCenter().GetNom(), Type2Type<CBlobGetXCenter>());
        fabricaOperadorsBlob.Register( CBlobGetYCenter().GetNom(), Type2Type<CBlobGetYCenter>());
        fabricaOperadorsBlob.Register( CBlobGetDistanceFromPoint().GetNom(), Type2Type<CBlobGetDistanceFromPoint>());

        // moments
        fabricaOperadorsBlob.Register( CBlobGetMoment().GetNom(), Type2Type<CBlobGetMoment>());

}

#endif  //BLOB_OBJECT_FACTORY