IRVE.H

#ifndef IRVE_H
#define IRVE_H

#include <TeAdoDB.h>
#include <TeDatabase.h>
#include <TeProjection.h>
#include <TeGeneralizedProxMatrix.h>
#include <TeSTObject.h> 
#include "Som.h"


class IRVE {
public:
        IRVE( SOM * som, vector<string> params, TIntVector clusters, int c ): _som(som), _params(params), map_cluster(clusters), c(c) 
        {
                vector<string> vec;
                m = som->getMapcode()->getMapSize();
                
                num_data = 359;

                data_cluster    = create_intvector( 0, num_data );
                visited                 = create_intvector( 0, num_data );

                for (int i=0; i <= num_data; ++i ) visited[i] = data_cluster[i] = 0;
                                        
                p  = create_vector( 1, c );
                q  = create_vector( 1, c );

                for (int k = 1; k <= c; ++k ) p[k] = q[k] = 0;

                codebook_label = som->getMapcode()->getLabel();
                
                for (int j=1; j<=c; j++ ) clusters_label.push_back( vec );

                ReadProxMatrix();
        };      


        TMatrix Index();

private:

        void MarkNeighborsVisited( int i, int c );
        void Calculate_p_q();
        void AssociateCluster();

        bool ReadProxMatrix(); 
        
        vector< vector<string> >  clusters_label;

        TVector                 p;
        TVector                 q;
        TIntVector              map_cluster;
        TIntVector              data_cluster;
        TIntVector              visited;

        TMatrix                 neighbors;

        SOM                             *_som;
        vector<string>  _params;
        
        int                             m;      
        int                             num_data;
        int                             c;
        
        TLabel                  codebook_label;         
};


void 
IRVE::MarkNeighborsVisited( int i, int c )
{
        if ( ( data_cluster[i] == c ) && ( visited[ i ] == 0 ) ) {
                visited[ i ] = 1;               
                for (int j=0; j<=num_data; j++)
                        if ( neighbors[i][j] == 1 )
                                if ( data_cluster[j] == c )
                                        MarkNeighborsVisited( j, c );
        }
};


void 
IRVE::Calculate_p_q()
{
        for (int i=0; i<=num_data; i++ ) 
                if ( visited[i] == 0 ) {
                        ++q[data_cluster[i]];
                        MarkNeighborsVisited( i, data_cluster[i] );
                }
};

void 
IRVE::AssociateCluster()
{       
        for (int i = 0; i < m; ++i )
                for (TvecLabel::iterator ite = codebook_label[i].begin(); ite != codebook_label[i].end(); ++ite ) {                     
                        data_cluster[ strtoint(  *ite  ) ] = map_cluster[ i + 1 ];
                        ++p[ map_cluster[ i + 1 ] ];
                }
        
};


TMatrix 
IRVE::Index() 
{
        AssociateCluster();

        Calculate_p_q();

        TMatrix result = create_matrix( 0, c, 1, 3 );

        result[0][1] = 0.0;

        for (int i=1; i<=c; ++i ) {             
                result[ i ][ 1 ] = q[ i ];
                result[ i ][ 2 ] = p[ i ];

                if ( q[ i ] == p[ i ] )
                        result[ i ][ 3 ] = 0;
                else
                        result[ i ][ 3 ] = -( (q[i] - 1)/p[i]) + 1.0;
                result[ 0 ][ 1 ] += result[ i ][ 3 ];
        }

        result[ 0 ][ 1 ] /= c;

        return result;
};


bool 
IRVE::ReadProxMatrix() {
        
        string host = _params[0];
    string dbname = _params[1];
    string user = _params[2];
    string pass = _params[3];
        string layer = _params[4];
        string where = _params[5];

        // Opens a connection to a database accessible though ADO
    TeDatabase* db_ = new TeAdo();
    if (!db_->connect(host,user,pass,dbname,0))
    {
            cout << "Error: " << db_->errorMessage() << endl;
            return false;
    }
    //cout << "1) Opened connection to: " << dbname << " on server: " << host << endl << endl;


    // Loads a layer by its name 
    TeLayer *layer1 = new TeLayer( layer );

    if (!db_->loadLayer(layer1))                
    {                                                                                               
            cout << "Error: " << db_->errorMessage() << endl;
            db_->close();
            return false;
    }

    // Retrieves the information about the polygons representation
    TeRepresentation* rep = layer1->getRepresentation(TePOLYGONS);
    if (!rep)
    {                                                                                               
            cout << "Layer has no polygons!" << endl;
            db_->close();
            return false;
    }
        
    string geomTableName = rep->tableName_;
        
    string q;
        q = "SELECT * FROM " + geomTableName;
//      if ( !where.empty() )
//              q += " WHERE " + where;

    // Get a portal to the database
    TeDatabasePortal *portal = db_->getPortal();
    if (!portal)
    {                                                                                               
                cout << "Error trying to get a database portal" << endl;
                db_->close();
                return false;
    }

    // Submit the query and go to the first position of the portal 
    if (!portal->query(q) || !portal->fetchRow())
    {       
                // error trying to submit the query or no records returned
                delete portal;
                return false;
    }


    // Process all the polygons found
    TePolygonSet ps;
    bool flag = true;
    do {
                TePolygon poly; 
                flag = portal->fetchGeometry(poly); 
                ps.add(poly);}  
        while (flag);   // while there is a record to be read
    
        delete portal;

        
        TeSTObject<TePolygonSet> objects;               
        objects.setGeometry( ps );
        
        
        TeSTObjectSet<TePolygonSet> s( layer1, layer );
        s.insertSTObject( objects );
        
        TeGeneralizedProxMatrix g( s, TePOLYGONS, "Breymann" ); 
        
        int max = 359;
        int i, j, sum_tot = 0;
        TeNeighbours n;
        neighbors = create_matrix( 0, max, 0, max );
        num_data = max;

        for (i=0; i<=max; i++) for (j=0; j <= max; j++) neighbors[i][j] = 0;

        for (i=0; i<ps.size(); i++) 
        {
                TePolygon pol = ps[i]; 
                for (j=0; j<pol.size(); j++) { 
                        n = g.getNeighbours( Te2String( pol[j].geomId() ) ); 
                        for (int count = 0; count < n.size(); ++count )                 
                                neighbors[ pol[j].geomId() ][ strtoint( n.ObjectId(count) ) ] =  1;
                } 
        }
        db_->close();
        return true;
};



#endif

---