packages/vtk/src/bbvtkImageBoundaries.cxx

   1 /*
   2  # ---------------------------------------------------------------------
   3  #
   4  # Copyright (c) CREATIS (Centre de Recherche en Acquisition et Traitement de l'Image
   5  #                        pour la SantÈ)
   6  # Authors : Eduardo Davila, Frederic Cervenansky, Claire Mouton
   7  # Previous Authors : Laurent Guigues, Jean-Pierre Roux
   8  # CreaTools website : www.creatis.insa-lyon.fr/site/fr/creatools_accueil
   9  #
  10  #  This software is governed by the CeCILL-B license under French law and
  11  #  abiding by the rules of distribution of free software. You can  use,
  12  #  modify and/ or redistribute the software under the terms of the CeCILL-B
  13  #  license as circulated by CEA, CNRS and INRIA at the following URL
  14  #  http://www.cecill.info/licences/Licence_CeCILL-B_V1-en.html
  15  #  or in the file LICENSE.txt.
  16  #
  17  #  As a counterpart to the access to the source code and  rights to copy,
  18  #  modify and redistribute granted by the license, users are provided only
  19  #  with a limited warranty  and the software's author,  the holder of the
  20  #  economic rights,  and the successive licensors  have only  limited
  21  #  liability.
  22  #
  23  #  The fact that you are presently reading this means that you have had
  24  #  knowledge of the CeCILL-B license and that you accept its terms.
  25  # ------------------------------------------------------------------------ */
  26
  27
  28 //=====
  29 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
  30 //=====
  31 #include "bbvtkImageBoundaries.h"
  32 #include "bbvtkPackage.h"
  33 #include "creaVtk_MACROS.h"
  34
  35 #include <omp.h>
  36 namespace bbvtk
  37 {
  38
  39 BBTK_ADD_BLACK_BOX_TO_PACKAGE(vtk,ImageBoundaries)
  40 BBTK_BLACK_BOX_IMPLEMENTATION(ImageBoundaries,bbtk::AtomicBlackBox);
  41 //=====
  42 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
  43 //=====
  44 void ImageBoundaries::Process()
  45 {
  46         printf("bbtk:vtk::ImageBoundaries::Process start\n");
  47         if (bbGetInputIn()!=NULL){
  48                 int i,j,k;
  49                 int ext[6];
  50                 double value;
  51
  52 //EED 2017-01-01 Migration VTK7
  53 #if (VTK_MAJOR_VERSION <= 5)
  54                 bbGetInputIn()->GetWholeExtent(ext);
  55 #endif
  56 #if (VTK_MAJOR_VERSION >= 6)
  57                 bbGetInputIn()->GetExtent(ext);
  58 #endif
  59
  60                 int maxX        = ext[1]-ext[0]+1;
  61                 int maxY        = ext[3]-ext[2]+1;
  62                 int maxZ        = ext[5]-ext[4]+1;
  63                 int bXMax       = 0;
  64                 int bYMax       = 0;
  65                 int bZMax       = 0;
  66                 int bXMin       = maxX-1;
  67                 int bYMin       = maxY-1;
  68                 int bZMin       = maxZ-1;
  69
  70 /*
  71                 #pragma omp parallel for
  72                 for (int i=0 ; i<maxX ; i++)
  73                 {
  74                         int j;
  75                         int k;
  76                         for (j=0 ; j<maxY ; j++)
  77                         {
  78                                 for (k=0 ; k<maxZ ; k++)
  79                                 {
  80                                         double value = bbGetInputIn()->GetScalarComponentAsDouble(i, j, k,0);
  81                                         if ( (value>=bbGetInputLowerValue()) && (value<=bbGetInputUpperValue())  )
  82                                         {
  83                                                 if (i<bXMin) bXMin=i;
  84                                                 if (j<bYMin) bYMin=j;
  85                                                 if (k<bZMin) bZMin=k;
  86                                                 if (i>bXMax) bXMax=i;
  87                                                 if (j>bYMax) bYMax=j;
  88                                                 if (k>bZMax) bZMax=k;
  89                                         }
  90                                 } // for k
  91                         } // for j
  92                 } // i
  93 */
  94
  95                 DEF_POINTER_IMAGE_VTK_CREA(vI,ssI,pI,stI,bbGetInputIn())
  96                 long int index=0;
  97                 for ( k=0 ; k<maxZ ; k++)
  98                 {
  99                         for (j=0 ; j<maxY ; j++)
 100                         {
 101                                 for (i=0 ; i<maxX ; i++)
 102                                 {
 103 //                                      double value = bbGetInputIn()->GetScalarComponentAsDouble(i, j, k,0);
 104                                         GETVALUE2_VTK_CREA(vI,pI,stI,index)
 105                                         index++;
 106                                         if ( (vI>=bbGetInputLowerValue()) && (vI<=bbGetInputUpperValue())  )
 107                                         {
 108                                                 if (i<bXMin) bXMin=i;
 109                                                 if (j<bYMin) bYMin=j;
 110                                                 if (k<bZMin) bZMin=k;
 111                                                 if (i>bXMax) bXMax=i;
 112                                                 if (j>bYMax) bYMax=j;
 113                                                 if (k>bZMax) bZMax=k;
 114                                         }
 115                                 } // for k
 116                         } // for j
 117                 } // i
 118                 std::vector<int> tmpIndex;
 119                 std::vector<int> tmpSize;
 120                 if (bXMin>bXMax){
 121                         tmpIndex.push_back(0);
 122                         tmpIndex.push_back(0);
 123                         tmpIndex.push_back(0);
 124                         tmpSize.push_back(1);
 125                         tmpSize.push_back(1);
 126                         tmpSize.push_back(1);
 127                 } else {
 128                         tmpIndex.push_back(bXMin);
 129                         tmpIndex.push_back(bYMin);
 130                         tmpIndex.push_back(bZMin);
 131                         tmpSize.push_back(bXMax-bXMin+1);
 132                         tmpSize.push_back(bYMax-bYMin+1);
 133                         tmpSize.push_back(bZMax-bZMin+1);
 134                 }
 135                 bbSetOutputIndex(tmpIndex);
 136                 bbSetOutputSize(tmpSize);
 137         } else {
 138                 printf("bbtk warnning: bbtk:vtk:ImageBoundaries box, Input image (In) not initialized..\n");
 139         }
 140 }
 141 //=====
 142 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
 143 //=====
 144 void ImageBoundaries::bbUserSetDefaultValues()
 145 {
 146
 147
 148 //  SET HERE THE DEFAULT INPUT/OUTPUT VALUES OF YOUR BOX
 149    bbSetInputIn(NULL);
 150         bbSetInputLowerValue(0);
 151         bbSetInputUpperValue(5000);
 152
 153         std::vector<int> tmpIndex;
 154         tmpIndex.push_back(0);
 155         tmpIndex.push_back(0);
 156         tmpIndex.push_back(0);
 157
 158         std::vector<int> tmpSize;
 159         tmpSize.push_back(1);
 160         tmpSize.push_back(1);
 161         tmpSize.push_back(1);
 162
 163         bbSetOutputIndex(tmpIndex);
 164         bbSetOutputSize(tmpSize);
 165 }
 166 //=====
 167 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
 168 //=====
 169 void ImageBoundaries::bbUserInitializeProcessing()
 170 {
 171
 172 //  THE INITIALIZATION METHOD BODY :
 173 //    Here does nothing
 174 //    but this is where you should allocate the internal/output pointers
 175 //    if any
 176
 177
 178 }
 179 //=====
 180 // Before editing this file, make sure it's a file of your own (i.e.: it wasn't generated from xml description; if so : your modifications will be lost)
 181 //=====
 182 void ImageBoundaries::bbUserFinalizeProcessing()
 183 {
 184
 185 //  THE FINALIZATION METHOD BODY :
 186 //    Here does nothing
 187 //    but this is where you should desallocate the internal/output pointers
 188 //    if any
 189
 190 }
 191 }
 192 // EO namespace bbvtk
 193
 194