RepositorySOMDataTerralib.h

#ifndef RepositorySOMDataTerralib_H
#define RepositorySOMDataTerralib_H

#include "TeAdoDB.h"
#include "TeDatabase.h"

#include "ISOMDataRepository.h"
#include "Data.h"
#include "Defs.h"

class RepositorySOMDataTerralib : public ISOMDataRepository {
public:
        void save( SOMData data );

        SOMData *  load( vector<string> params, int flag );
public:
        TMatrix getCentroidCoords( vector<string> params  );
};


TMatrix 
RepositorySOMDataTerralib::getCentroidCoords( vector<string> params  )
{
        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;
            exit(0);
    }
    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();
            exit(0);
    }

    // Retrieves the information about the polygons representation
    TeRepresentation* rep = layer1->getRepresentation(TePOLYGONS);
    if (!rep)
    {                                                                                               
            cout << "Layer has no polygons!" << endl;
            db_->close();
            exit(0);
    }
        
    string geomTableName = rep->tableName_;

    // Build the SQL to get all polygons that are inside a given box
    //TeBox bb(320000,7390000,338000,7400000);
    string q;
    q = "SELECT * FROM " + geomTableName + " WHERE ";
    //    q +=  db_->getSQLBoxWhere (bb, TePOLYGONS);
                q += where;

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

    // 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;
                exit(0);
    }

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

    // Shows the result
    for (int i=0; i<ps.size(); i++)
    {
                TePolygon pol = ps[i];
                cout << "Geometry id: " << pol[0].geomId() << " object id: " << pol.objectId() << endl;
                for (int j=1; j<pol.size(); j++)
                         cout << "\t Child id: " << pol[j].geomId() << endl;
     }

     // Close database
            db_->close();
        cout <<  "End! (Press enter)\n";
        cout.flush();
        getchar();
        return 0;
};


void 
RepositorySOMDataTerralib::save( SOMData data ) 
{
        cout << "Save Data\n";
};

SOMData *  
RepositorySOMDataTerralib::load( vector<string> params, int flag )
{
        
        TvecLabel lvec;
        MatrixInt mat;
        VecInt vec;
        MatrixInt::iterator im;
        VecInt::iterator        iv;

        string host = params[0];
    string dbname = params[1];
    string user = params[2];
    string pass = params[3];
        string table = params[4];
        string where = params[5];

        int dimension, i, j;

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

        //      4) Rettrieve some data from a table throught a SQL query 

    TeDatabasePortal* portal = db->getPortal();
    if (!portal)
    {
            cout << "Error getting a portal to dabase: " << db->errorMessage() << endl;
            db->close();
            return 0;
    }
        string variables = params[6];
        for (i=7; i<params.size(); i++) variables = variables + ", " + params[i];
    string sql = "SELECT " + variables + " FROM " + table + " WHERE " + where;

        dimension = params.size() - 7;

    if (portal->query(sql))
    {            
            while (portal->fetchRow())                      // retrieve the attributoes
            {                                   
                                        lvec.push_back( portal->getData(0) );

                                        vec.clear();
                                        for (int j=1; j<= dimension; j++)  vec.push_back( portal->getDouble(j) );                       

                                        mat.push_back( vec );
            }

                        TMatrix mdata = create_matrix( 1, mat.size(), 1, dimension );
                        i = j = 1;
                        
                        for (im=mat.begin(); im<mat.end(); im++) {
                                j = 1;                          
                                for (iv=(*im).begin(); iv < (*im).end(); iv++) 
                                {
                                        mdata[i][j] = *iv;                                                                      
                                        j++;
                                }
                                i++;
                        }

                        SOMData * som_data = new SOMData;
        
                        som_data->setDimension( dimension );
                        som_data->setDataSize( mat.size() );
                        som_data->setEntries( mdata );
                        som_data->setLabel( lvec );

                        delete portal;
                        db->close();

                        return som_data;

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

    delete portal;
    db->close();

        return 0;
}

#endif

---