diff options
Diffstat (limited to 'conversion/tor.cc')
-rw-r--r-- | conversion/tor.cc | 230 |
1 files changed, 230 insertions, 0 deletions
diff --git a/conversion/tor.cc b/conversion/tor.cc new file mode 100644 index 0000000..0663fd4 --- /dev/null +++ b/conversion/tor.cc @@ -0,0 +1,230 @@ +/* +* 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>. +*/ + +#include "conversion/tor.hh" +#include "raslib/endian.hh" +#include "raslib/rminit.hh" +#include "raslib/parseparams.hh" +#include "raslib/primitivetype.hh" + +void r_Conv_TOR::initTOR() +{ + if(params ==NULL) + { + params = new r_Parse_Params(3); + } +} + +r_Conv_TOR::r_Conv_TOR(const char* source, const r_Minterval& lengthordomain, const r_Type* tp) throw(r_Error) + : r_Convertor(source, lengthordomain, tp, true) + { + initTOR(); + } + +r_Conv_TOR::r_Conv_TOR(const char* source, const r_Minterval& lengthordomain, int tp) throw(r_Error) + : r_Convertor(source, lengthordomain, tp) + { + initTOR(); + } + +r_convDesc& +r_Conv_TOR::convertFrom(const char* options) throw (r_Error) + { + int swap = 0; + int rescale = 0; + char* domain = NULL; + params->add("swapendianness", &swap, r_Parse_Params::param_type_int); + params->add("rescale", &rescale, r_Parse_Params::param_type_int); + params->add("domain", &domain, r_Parse_Params::param_type_string); + params->process(options); + r_Minterval t; + RMInit::logOut << "r_Conv_TOR::convert swap " << swap << " rescale " << rescale << " domain " << (domain?domain:"NULL") << std::endl; + if (domain == NULL) + { + RMInit::logOut << "r_Conv_TOR::convertFrom no domain specified in options string" << std::endl; + throw r_Error(); + } + try { + t = r_Minterval(domain); + } + catch (r_Eno_interval& e) + { + RMInit::logOut << "r_Conv_TOR::convertFrom no correct domain specified in options string" << std::endl; + throw r_Error(); + } + if (t.dimension() != 2) + { + RMInit::logOut << "r_Conv_TOR::convertFrom domain in options string must have 2 dimensions" << std::endl; + throw r_Error(); + } + r_Range x = t[0].high() - t[0].low() + 1; + r_Range y = t[1].high() - t[1].low() + 1; + r_UShort* array = (r_UShort*)mystore.storage_alloc(x * y * sizeof(r_UShort)); + double max = 0; + double min = 0; + r_Char tc; + double maxu = 0; + double minu = 0; + r_UShort tu; + switch (desc.baseType) + { + case ctype_bool: + case ctype_char: + case ctype_uint8: + if (desc.srcInterv[0].high() - desc.srcInterv[0].low() + 1 != t.cell_count() * sizeof(r_Char)) + { + RMInit::logOut << "r_Conv_TOR::convertFrom the supplied data is " << desc.srcInterv[0].high() - desc.srcInterv[0].low() + 1 << " bytes, but " << t.cell_count() << " bytes were expected" << std::endl; + throw r_Error(); + } + get_limits(&tc, min, max); + get_limits(&tu, minu, maxu); + for(int i=0; i < x; i++) + { + for(int j=0; j < y; j++) + { + if (rescale) + array[i * x + j] = ((r_Octet*)desc.src)[j * x + i] / max * maxu; + else + array[i * x + j] = ((r_Octet*)desc.src)[j * x + i]; + } + } + break; + case ctype_uint16: + if (desc.srcInterv[0].high() - desc.srcInterv[0].low() + 1 != t.cell_count() * sizeof(r_UShort)) + { + RMInit::logOut << "r_Conv_TOR::convertFrom the supplied data is " << desc.srcInterv[0].high() - desc.srcInterv[0].low() + 1 << " bytes, but " << t.cell_count() << " bytes were expected" << std::endl; + throw r_Error(); + } + for(int i=0; i < x; i++) + { + for(int j=0; j < y; j++) + { + array[i * x + j] = ((r_UShort*)desc.src)[j * x + i]; + } + } + break; + default: + RMInit::logOut << "r_Conv_TOR unknown base type!" << std::endl; + throw r_Error(); + } + if (swap) + { + r_UShort* temp = (r_UShort*)mystore.storage_alloc(x * y * sizeof(r_UShort)); + r_Endian::swap_array(x * y * sizeof(r_UShort), array, temp); + mystore.storage_free(array); + array = temp; + } + + desc.dest = (char*)array; + desc.destInterv = t; + desc.destType = new r_Primitive_Type("UShort", r_Type::USHORT); + return desc; + } + +r_convDesc& +r_Conv_TOR::convertTo(const char* options) throw (r_Error) + { + int swap = 0; + int rescale = 0; + params->add("swapendianness", &swap, r_Parse_Params::param_type_int); + params->add("rescale", &rescale, r_Parse_Params::param_type_int); + params->process(options); + if (desc.srcInterv.dimension() != 2) + { + RMInit::logOut << "r_Conv_TOR::convertFrom domain in options string must have 2 dimensions" << std::endl; + throw r_Error(); + } + r_Range x = desc.srcInterv[0].high() - desc.srcInterv[0].low() + 1; + r_Range y = desc.srcInterv[1].high() - desc.srcInterv[1].low() + 1; + r_UShort* array = (r_UShort*)mystore.storage_alloc(x * y * sizeof(r_UShort)); + double max = 0; + double min = 0; + r_Char tc; + double maxu = 0; + double minu = 0; + r_UShort tu; + switch (desc.baseType) + { + case ctype_bool: + case ctype_char: + case ctype_uint8: + get_limits(&tc, min, max); + get_limits(&tu, minu, maxu); + for(int i=0; i < x; i++) + { + for(int j=0; j < y; j++) + { + if (rescale) + array[i * x + j] = ((r_Octet*)desc.src)[j * x + i] / max * maxu; + else + array[i * x + j] = ((r_Octet*)desc.src)[j * x + i]; + } + } + break; + case ctype_uint16: + for(int i=0; i < x; i++) + { + for(int j=0; j < y; j++) + { + array[i * x + j] = ((r_UShort*)desc.src)[j * x + i]; + } + } + break; + default: + RMInit::logOut << "r_Conv_TOR unknown base type!" << std::endl; + throw r_Error(); + } + if (swap) + { + r_UShort* temp =(r_UShort*)mystore.storage_alloc(x * y * sizeof(r_UShort)); + r_Endian::swap_array(x * y * sizeof(r_UShort), array, temp); + mystore.storage_free(array); + array = temp; + } + + desc.dest = (char*)array; + r_Minterval td(1); + td << r_Sinterval((r_Range)0, (r_Range)(desc.srcInterv.cell_count() * sizeof(r_UShort) - 1)); + desc.destInterv = td; + desc.destType = new r_Primitive_Type("Char", r_Type::CHAR); + return desc; + } + +const char* +r_Conv_TOR::get_name() const + { + return get_name_from_data_format(r_TOR); + } + +r_Data_Format +r_Conv_TOR::get_data_format() const + { + return r_TOR; + } + +r_Convertor* +r_Conv_TOR::clone() const + { + return new r_Conv_TOR(desc.src, desc.srcInterv, desc.srcType); + } + |