]> Creatis software - clitk.git/blob - common/clitkListOfPair.cxx
new Airway filter (experimental)
[clitk.git] / common / clitkListOfPair.cxx
1 /*=========================================================================
2   Program:   vv                     http://www.creatis.insa-lyon.fr/rio/vv
3
4   Authors belong to:
5   - University of LYON              http://www.universite-lyon.fr/
6   - Léon Bérard cancer center       http://oncora1.lyon.fnclcc.fr
7   - CREATIS CNRS laboratory         http://www.creatis.insa-lyon.fr
8
9   This software is distributed WITHOUT ANY WARRANTY; without even
10   the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
11   PURPOSE.  See the copyright notices for more information.
12
13   It is distributed under dual licence
14
15   - BSD        See included LICENSE.txt file
16   - CeCILL-B   http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
17 ======================================================================-====*/
18 #ifndef CLITKLISTOFPAIR_CXX
19 #define CLITKLISTOFPAIR_CXX
20 /**
21    ------------------------------------------------=
22    * @file   clitkListOfPair.cxx
23    * @author David Sarrut <david.sarrut@creatis.insa-lyon.fr>
24    * @date   27 Feb 2007 09:55:56
25    *
26    * @brief
27    *
28    *
29    ------------------------------------------------=*/
30
31 #include "clitkListOfPair.h"
32
33 //--------------------------------------------------------------------
34 double clitk::convertValue(double v,
35                            const std::multimap<double, double> & conversionTable,
36                            bool linear)
37 {
38   std::multimap<double, double>::const_iterator i;
39   i = conversionTable.lower_bound(v);
40   if (i == conversionTable.end()) {
41     std::cerr << "The value " << v << " is out of the table" << std::endl;
42     exit(0);
43   }
44
45   double v2 = i->first;
46   double p2 = i->second;
47   if (i != conversionTable.begin()) i--;
48   double v1 = i->first;
49   double p1 = i->second;
50
51   // interpol
52   if (!linear) {
53     if ((v-v1) > (v2-v)) return p2;
54     else return p1;
55   } else {
56     double w = (v-v1)/(v2-v1);
57     return p1*(1.0-w)+w*p2;
58   }
59 }
60 //--------------------------------------------------------------------
61
62 #endif /* end #define CLITKLISTOFPAIR_CXX */
63