packages/vtk/src/bbvtkRescaleSlopeIntercept.cxx

   1 //=====
   2 // Don't edit this file. This file is generated from xml description..
   3 //=====
   4 #include "bbvtkRescaleSlopeIntercept.h"
   5 #include "bbvtkPackage.h"
   6
   7
   8 namespace bbvtk
   9 {
  10
  11 BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,RescaleSlopeIntercept)
  12 BBTK_BLACK_BOX_IMPLEMENTATION(RescaleSlopeIntercept,bbtk::AtomicBlackBox);
  13 //=====
  14 // Don't edit this file. This file is generated from xml description..
  15 //=====
  16 void RescaleSlopeIntercept::Process()
  17 {
  18
  19 // THE MAIN PROCESSING METHOD BODY
  20 //   Here we simply set the input 'In' value to the output 'Out'
  21 //   And print out the output value
  22 // INPUT/OUTPUT ACCESSORS ARE OF THE FORM :
  23 //    void bbSet{Input|Output}NAME(const TYPE&)
  24 //    const TYPE& bbGet{Input|Output}NAME() const
  25 //    Where :
  26 //    * NAME is the name of the input/output
  27 //      (the one provided in the attribute 'name' of the tag 'input')
  28 //    * TYPE is the C++ type of the input/output
  29 //      (the one provided in the attribute 'type' of the tag 'input')
  30 //    bbSetOutputOut( bbGetInputIn() );
  31 //    std::cout << "Output value = " <<bbGetOutputOut() << std::endl;
  32
  33
  34         double range[2];
  35         double spc[3];
  36         int ext[6];
  37         double A=bbGetInputA();
  38         double B=bbGetInputB();
  39
  40         if (bbGetInputIn()!=0)
  41         {
  42                 if (bbGetInputOutputFormat()=="SAME")                                   outputformat = bbGetInputIn()->GetScalarType();
  43                 else if (bbGetInputOutputFormat()=="VTK_BIT")                   outputformat = VTK_BIT;
  44                 else if (bbGetInputOutputFormat()=="VTK_CHAR")                  outputformat = VTK_CHAR;
  45                 else if (bbGetInputOutputFormat()=="VTK_SIGNED_CHAR")   outputformat = VTK_SIGNED_CHAR;
  46                 else if (bbGetInputOutputFormat()=="VTK_UNSIGNED_CHAR") outputformat = VTK_UNSIGNED_CHAR;
  47                 else if (bbGetInputOutputFormat()=="VTK_UNSIGNED_SHORT") outputformat = VTK_UNSIGNED_SHORT;
  48                 else if (bbGetInputOutputFormat()=="VTK_INT")                   outputformat = VTK_INT;
  49                 else if (bbGetInputOutputFormat()=="VTK_UNSIGNED_INT")  outputformat = VTK_UNSIGNED_INT;
  50                 else if (bbGetInputOutputFormat()=="VTK_LONG")                  outputformat = VTK_LONG;
  51                 else if (bbGetInputOutputFormat()=="VTK_UNSIGNED_LONG") outputformat = VTK_UNSIGNED_LONG;
  52                 else if (bbGetInputOutputFormat()=="VTK_FLOAT")                 outputformat = VTK_FLOAT;
  53                 else if (bbGetInputOutputFormat()=="VTK_DOUBLE")                outputformat = VTK_DOUBLE;
  54
  55
  56                 bbGetInputIn()->GetScalarRange(range);
  57                 bbGetInputIn()->GetSpacing(spc);
  58                 bbGetInputIn()->GetExtent( ext );
  59
  60                 double range0plusrange1 = range[0]+range[1];
  61                 double difrange                 = range[1]-range[0];
  62                 double newdifrange              = B-A;
  63
  64                 double shift;
  65                 double scale;
  66
  67                 if (bbGetInputType()==0) // Rescale Slope/Intercept
  68                 {
  69                         //                      newValue = value*A+B;
  70                         scale = A;
  71                         shift = B;
  72                 } // if type 0
  73
  74                 if (bbGetInputType()==1) // Invert
  75                 {
  76                         //                      newValue = value*(-1) + (range0plusrange1);
  77                         scale=-1;
  78                         shift=range0plusrange1;
  79                 } // if type 1
  80
  81                 if (bbGetInputType()==2) // Redimension
  82                 {
  83                         //                      newValue = ((value-range[0])/difrange)*newdifrange + A;
  84                         scale=newdifrange/difrange;
  85                         shift= A - (range[0]*newdifrange/difrange);
  86                 } // if type 2
  87
  88                 if (bbGetInputType()==3) // Invert redimension
  89                 {
  90                         //                      newValue = ((value*(-1)+range[1])/difrange)*newdifrange + A;
  91                         scale=-newdifrange/difrange;
  92                         shift= A + (range[1]*newdifrange/difrange);
  93                 } // if type 3
  94
  95                 filter->SetInput( bbGetInputIn() );
  96                 filter->SetOutputScalarType( outputformat );
  97                 filter->SetScale( scale );
  98                 filter->SetShift( shift );
  99
 100                 mchange->SetInput( filter->GetOutput() );
 101                 mchange->SetExtentTranslation( -ext[0], -ext[2], -ext[4] );
 102                 mchange->SetOutputSpacing (spc);
 103                 mchange->Update();    //important
 104
 105                 bbSetOutputOut( mchange->GetOutput() );
 106
 107         } // if
 108 }
 109
 110 //=====
 111 // Don't edit this file. This file is generated from xml description..
 112 //=====
 113 void RescaleSlopeIntercept::bbUserSetDefaultValues()
 114 {
 115
 116 //  SET HERE THE DEFAULT INPUT/OUTPUT VALUES OF YOUR BOX
 117 //    Here we initialize the input 'In' to 0
 118         bbSetInputIn(NULL);
 119         bbSetInputType(0);
 120         bbSetInputA(1);
 121         bbSetInputB(0);
 122         bbSetInputOutputFormat("SAME");
 123
 124         filter  = NULL;
 125         mchange = NULL;
 126 }
 127 //=====
 128 // Don't edit this file. This file is generated from xml description..
 129 //=====
 130 void RescaleSlopeIntercept::bbUserInitializeProcessing()
 131 {
 132
 133 //  THE INITIALIZATION METHOD BODY :
 134 //    Here does nothing
 135 //    but this is where you should allocate the internal/output pointers
 136 //    if any
 137
 138         filter  = vtkImageShiftScale::New();
 139         mchange = vtkImageChangeInformation::New();
 140
 141
 142 }
 143 //=====
 144 // Don't edit this file. This file is generated from xml description..
 145 //=====
 146 void RescaleSlopeIntercept::bbUserFinalizeProcessing()
 147 {
 148
 149 //  THE FINALIZATION METHOD BODY :
 150 //    Here does nothing
 151 //    but this is where you should desallocate the internal/output pointers
 152 //    if any
 153         mchange->Delete();
 154         filter->Delete();
 155 }
 156
 157 }
 158 // EO namespace bbvtk
 159
 160