summaryrefslogtreecommitdiffstats log msg author committer range
path: root/rasodmg/interesttiling.hh
blob: 37d6eca7a358ac5099a63a881cc72d6c600019bc (plain)
 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186  /* * 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 . * * Copyright 2003, 2004, 2005, 2006, 2007, 2008, 2009 Peter Baumann / rasdaman GmbH. * * For more information please see * or contact Peter Baumann via . / /** * INCLUDE: dirtiling.hh * * MODULE: rasodmg * CLASS: r_Interest_Tiling * * COMMENTS: * None */ #ifndef _R_INTERESTTILING_HH_ #define _R_INTERESTTILING_HH_ class r_Dir_Decompose; class r_Interest_Tiling; #include "rasodmg/tiling.hh" //@ManMemo: Module: {\bf rasodmg} /*@Doc: This class implements the "Interesting Areas" tiling algorithm. The user specifies which areas are of interest to him and tiling is performed accordingly. Example of usage: The user says that [10:20, 50:60] and [18:50, 65:70] are of interest to him in the [0:1000, 0:1000] domain. Tiling should be done according to this setting. To specify this with this class the following code would work: ... r_Minterval domain(2); domain << r_Sinterval(0L, 1000L) << r_Sinterval(0L, 1000L); r_Minteval interest1(2); interest1 << r_Sinterval(10L, 20L) << r_Sinterval(50L, 60L); r_Minterval interest2(2); interest2 << r_Sinterval(18L, 50L) << r_Sinterval(65L, 70L); std::vector interest_areas; interest_areas.insert_element(interest1); interest_areas.insert_element(interest2); r_Interest_Tiling(interest_areas); ... */ /** * \ingroup Rasodmgs */ class r_Interest_Tiling : public r_Dimension_Tiling { public: /// Possible strategies to limit the tilesize enum Tilesize_Limit {NO_LIMIT, REGROUP, SUB_TILING, REGROUP_AND_SUBTILING, NUMBER}; /// read everything from an encoded string /// e.g. "2;[0:9,0:9];[100:109,0:9];100;REGROUPSUBTILING" r_Interest_Tiling(const char* encoded) throw (r_Error); r_Interest_Tiling(r_Dimension dim, const std::vector& interest_areas, r_Bytes ts = RMInit::clientTileSize, Tilesize_Limit strat = SUB_TILING) throw (r_Error); /** It takes as parameter a list containing the areas of interest to the user and also the tilesize to be used. The constructor also takes as parameter the tilesize limitation strategy, that can be the following: NO_LIMIT: The generated blocks can have any size. REGROUP: Only when performing grouping/merging of blocks, the size of the resulting block of two merges is checked against tilesize. If it's bigger, they are not merged. Blocks larger than tilesize may exist (for instance, if the user specifies an interest area larger then tilesize). SUB_TILING: In this strategie, regrouping is done regardless of the size of the generated blocks. After all the blocks are created, sub-tiling is performed on those whose size is larger than tilesize. REGROUP_AND_SUBTILING: This combines the last two strategies. When merging blocks, blocks larger than tilesize are never created and, when the final blocks are all created, sub-tiling is performed on those whose size is larger then tilesize. An exception is thrown when the dimensions of the domains does not match the specified dimension. */ virtual ~r_Interest_Tiling(); std::vector* compute_tiles(const r_Minterval& obj_domain, r_Bytes cell_size) const throw (r_Error); virtual void print_status(std::ostream& os) const; virtual r_Tiling* clone() const; virtual r_Tiling_Scheme get_tiling_scheme() const; //@ManMemo: Module: {\bf raslib} /** Get a tilesize limit for a tilisize limit name */ static r_Interest_Tiling::Tilesize_Limit get_tilesize_limit_from_name(const char* name); //@ManMemo: Module: {\bf raslib} /** Get a tilisize limit name for a tilesize limit */ static const char* get_name_from_tilesize_limit(Tilesize_Limit tsl); static const char* description; protected: // methods. /// Given a domain and a set of interest areas (internal) gener. partition std::vector* make_partition(const r_Minterval& domain) const; /// The Block types (A, B or C) enum Blocks_Type {BLOCKS_A, BLOCKS_B, BLOCKS_C }; /** A Blocks ---> Blocks that belong to two or more interest zones (non-groupable) B Blocks ---> Blocks within the same interest zone (groupable by interest zone) C Blocks ---> Blocks outside interest zones (groupable) */ /// Merge as many blocks together in a list as possible std::vector* group(std::vector& blocks, r_Bytes typelen, Blocks_Type btype) const; /** Parameters: the block list, the lenght of the base cells and the type of the block (A, B or C) */ /// Tilesize limitation strategie Tilesize_Limit ts_strat; /// Interest areas std::vector iareas; //@ManMemo: Module: {\bf raslib} /** The names of all tilesizelimit types, to avoid redundant storage and inconsistencies. The variable name convention is the prefix tilisizelimit_name_ followed by the name of the data format in lower case, i.e. for NOLIMIT tilisizelimit_name_nolimit. In addition there's an array of names all_tilesizelimit_names where the tilesize limit can be used as index to get the name. */ static const char* tilesizelimit_name_nolimit; static const char* tilesizelimit_name_regroup; static const char* tilesizelimit_name_subtiling; static const char* tilesizelimit_name_regroupandsubtiling; static const char* all_tilesizelimit_names[r_Interest_Tiling::NUMBER]; }; #endif