diff options
Diffstat (limited to 'manuals_and_examples/examples/c++/avg-cell-red.cc')
-rw-r--r-- | manuals_and_examples/examples/c++/avg-cell-red.cc | 144 |
1 files changed, 144 insertions, 0 deletions
diff --git a/manuals_and_examples/examples/c++/avg-cell-red.cc b/manuals_and_examples/examples/c++/avg-cell-red.cc new file mode 100644 index 0000000..4896a83 --- /dev/null +++ b/manuals_and_examples/examples/c++/avg-cell-red.cc @@ -0,0 +1,144 @@ +/* +* This file is part of rasdaman community. +* +* Rasdaman community is free software: you can redistribute it and/or modify +* it under the terms of the GNU General Public License as published by +* the Free Software Foundation, either version 3 of the License, or +* (at your option) any later version. +* +* Rasdaman community is distributed in the hope that it will be useful, +* but WITHOUT ANY WARRANTY; without even the implied warranty of +* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the +* GNU General Public License for more details. +* +* You should have received a copy of the GNU General Public License +* along with rasdaman community. If not, see <http://www.gnu.org/licenses/>. +* +* Copyright 2003, 2004, 2005, 2006, 2007, 2008, 2009 Peter Baumann / +rasdaman GmbH. +* +* For more information please see <http://www.rasdaman.org> +* or contact Peter Baumann via <baumann@rasdaman.com>. +*/ +/************************************************************* + * + * + * PURPOSE: + * Example program for computing the avg cell value + * for each image in a given 2-D RGB collection. + * + * + ************************************************************/ + +#include <iostream> + +// Linux needs this for template instantiation +#ifdef EARLY_TEMPLATE +#define __EXECUTABLE__ +#ifdef __GNUG__ +#include "raslib/template_inst.hh" +#endif +#endif + + +// rasdl generated data type file: +// (contains predefined RGB pixel structure) +#include "basictypes.hh" + + +int main( int ac, char** av ) +{ + char rasmgrName[255]; + int rasmgrPort; + char baseName[255]; + char collName[255]; + char userName[255]; + char userPass[255]; + + double sum; + + if( ac != 7 ) + { + cout << "Usage: avg-cell-red rasmgr_name rasmgr_port base_name collection_name user_name user_passwd" << endl; + return -1; + } + + strcpy( rasmgrName, av[1] ); + rasmgrPort = strtoul( av[2], NULL, 0); + strcpy( baseName, av[3] ); + strcpy( collName, av[4] ); + strcpy( userName, av[5] ); + strcpy( userPass, av[6] ); + + r_Database database; + r_Transaction transaction; + r_Ref< r_Set< r_Ref< r_Marray< RGBPixel > > > > image_set; + r_Iterator< r_Ref< r_Marray< RGBPixel > > > iter; + r_Ref< r_Marray< RGBPixel > > image; + r_Range i,j; + long pixelcount; + + try + { + database.set_servername( rasmgrName, rasmgrPort ); + database.set_useridentification( userName, userPass ); + + cout << "Opening database " << baseName + << " on " << rasmgrName << "... " << flush; + database.open( baseName ); + cout << "OK" << endl; + + cout << "Starting read-only transaction ... " << flush; + transaction.begin( r_Transaction::read_only ); + cout << "OK" << endl; + + cout << "Looking up collection " << collName << " ..." << flush; + image_set = database.lookup_object( collName ); + cout << "OK" << endl; + + cout << "Collection contains " << image_set->cardinality() + << " entries" << endl; + + iter = image_set->create_iterator(); + for( iter.reset(); iter.not_done(); iter++ ) + { + RGBPixel pixel; + image = *iter; + sum = 0.0; + for ( i=image->spatial_domain()[0].low(); i<=image->spatial_domain()[0].high(); i++ ) + { + for ( j=image->spatial_domain()[1].low(); j<=image->spatial_domain()[1].high(); j++ ) + { + pixel = (*image)[ r_Point(i,j) ]; + sum += pixel.red; + } + } + pixelcount = + ( image->spatial_domain()[0].high() - image->spatial_domain()[0].low() + 1 ) + * ( image->spatial_domain()[1].high() - image->spatial_domain()[1].low() + 1 ); + cout << " avg over " << pixelcount + << " red pixels is " + << sum / pixelcount + << endl << flush; + } + + cout << "Committing transaction ... " << flush; + transaction.commit(); + cout << "OK" << endl; + + cout << "Closing database ... " << flush; + database.close(); + cout << "OK" << endl; + } + catch( r_Error& errorObj ) + { + cerr << errorObj.what() << endl; + return -1; + } + + return 0; +} + +/* + * end of avg-cell-red.cc + */ |