SingleSegmentation.h

#ifndef SingleSegmentation_H
#define SingleSegmentation_H


#include <TeAdoDB.h>
#include <TeDatabase.h>
#include "Som.h"


class SingleSegmentation  {
public:

        
        SingleSegmentation(){};

        
        void Segment( SOM* som, int );

        void createField( vector<string> params );

        


private:

        void HorizontalSegmentation();
        
        void VerticalSegmentation();

        void Diagonal01Segmentation();

        void Diagonal02Segmentation();
        
private:
        SOM * _som; 
        int m;
        TIntVector  _clusters;
        int                _num_clusters;
        
};


void 
SingleSegmentation::createField( vector<string> params ) {

        string host = params[0];
    string dbname = params[1];
    string user = params[2];
    string pass = params[3];
        string tableName = params[4];
        string fieldName = params[5];
        string labelField = params[6];

        TLabel label = _som->getMapcode()->getLabel();
        

        int i;

        // 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;
    }

        TeAttribute atr;
        atr.rep_.name_= fieldName;
    atr.rep_.type_ = TeREAL;
    atr.rep_.numChar_ = 15;
    atr.rep_.decimals_ = 0;
    
        if (!db->addColumn(tableName, atr.rep_))
        {
                cout << "Error adding column to table: " << db->errorMessage() << endl;
                db->close();
                return;
        }



    TeDatabasePortal* portal = db->getPortal();
    if (!portal)
    {
            cout << "Error getting a portal to dabase: " << db->errorMessage() << endl;
            db->close();
            return;
    }
                                
        string sql = "UPDATE " + tableName  + " SET " + fieldName + " = 0";
        

        if (portal->query(sql))
        {
                  cout << "Error submiting query: " << db->errorMessage() << endl;
          delete portal;
          db->close();
          return;
        }
        

        for (i=0; i< m; i++ ) {
                for (int j = 0; j < label[i].size(); j++ ) {

                        string value = inttostr(_clusters[i+1]);
                        string id    = label[i][j];
                        sql = "UPDATE " + tableName  + " SET " + fieldName + " = " + value + " WHERE " + labelField + " = " + "\""+ id + "\"";
                        //cout << sql << endl;
                        if (portal->query(sql))
                        {
                                cout << "Error submiting query: " << db->errorMessage() << endl;
                                delete portal;
                                db->close();
                                return;
                        }
                        
                }
        }

        delete portal;
    db->close();
        
};

void 
SingleSegmentation::HorizontalSegmentation() {
        
        int L = _som->getMapcode()->getDimensions(0);
        int C = _som->getMapcode()->getDimensions(1);
        
        _clusters = create_intvector( 1, m );

        int count = 1;
        for (int j=1; j<= C; j++ )
                for (int i=1; i<= L; i++ ) {
                        _clusters[ ((j - 1)*L + i) ] = count;
                        count++;
                }

        
};

void 
SingleSegmentation::VerticalSegmentation() {
        
        int L = _som->getMapcode()->getDimensions(0);
        int C = _som->getMapcode()->getDimensions(1);
        
        _clusters = create_intvector( 1, m );

        int count = 1;
        for (int i=1; i<= L; i++ )
                for (int j=1; j<= C; j++ ) {
                        _clusters[ ((j - 1)*L + i) ] = count;
                        ++count;
                }
        

};

void 
SingleSegmentation::Diagonal01Segmentation() {

        int L = _som->getMapcode()->getDimensions(0);
        int C = _som->getMapcode()->getDimensions(1);
        int i, j;
        
        _clusters = create_intvector( 1, m );
        for (i=1; i<=m; i++ ) _clusters[i] = 0;

        TMatrix mat = create_matrix( 1, L, 1, C );
        
        for (i=1; i<=L; i++) for (j=1; j<= C; j++ ) mat[i][j] = 1;
        
        int count = 1;
        bool flag = false;
        
        while ( count <= m ) {

                flag = false;
                for (i=1; i<=L; i++) {
                        for (j=1; j <= C; j++ )                 
                                if ( mat[i][j] == 1 ) {
                                        mat[i][j] = 0;
                                        flag = true;
                                        break;
                                }
                        if ( flag ) break;
                }

                _clusters[ ((j - 1)*L + i ) ] = count;
                ++count;
                while  (( j > 1 ) && ( i < L ) ) {
                                --j;
                                ++i;
                                mat[i][j] = 0;
                                _clusters[ ((j - 1)*L + i ) ] = count;
                                ++count;
                }
                
        }       
}

void 
SingleSegmentation::Diagonal02Segmentation() {

        int L = _som->getMapcode()->getDimensions(0);
        int C = _som->getMapcode()->getDimensions(1);
        int i, j;
        
        _clusters = create_intvector( 1, m );
        for (i=1; i<=m; i++ ) _clusters[i] = 0;

        TMatrix mat = create_matrix( 1, L, 1, C );
        
        for (i=1; i<=L; i++) for (j=1; j<= C; j++ ) mat[i][j] = 1;
        
        int count = 1;
        bool flag = false;
        
        while ( count <= m ) {

                flag = false;
                for (i=1; i<=L; i++) {
                        for (j=C; j >= 1; j-- )                 
                                if ( mat[i][j] == 1 ) {
                                        mat[i][j] = 0;
                                        flag = true;
                                        break;
                                }
                        if ( flag ) break;
                }

                _clusters[ ((j - 1)*L + i ) ] = count;
                ++count;
                while  (( j < C ) && ( i < L ) ) {
                                ++j;
                                ++i;
                                mat[i][j] = 0;
                                _clusters[ ((j - 1)*L + i ) ] = count;
                                ++count;
                }
                
        }               
};

void 
SingleSegmentation::Segment( SOM* som, int single_type ) {
        
        _som = som;
        m = som->getMapcode()->getMapSize();
        _num_clusters = m;
        
        switch( single_type ) {
                case 0 : HorizontalSegmentation(); break;
                case 1 : VerticalSegmentation(); break;
                case 2 : Diagonal01Segmentation();  break;
                case 3 : Diagonal02Segmentation(); break;
        };
        
};


#endif

---