DGtal  1.4.beta
testGenericReader.cpp
Go to the documentation of this file.
1 
31 #include <iostream>
32 #include "DGtal/base/Common.h"
33 #include "DGtal/images/ImageContainerBySTLVector.h"
34 #include "DGtal/io/readers/GenericReader.h"
35 #include "DGtal/helpers/StdDefs.h"
36 #include "ConfigTest.h"
37 
38 
39 
40 
41 
43 
45 // Functions for testing class GenericReader.
47 
51 bool testGenericReader()
52 {
53  unsigned int nbok = 0;
54  unsigned int nb = 0;
55 
56  DGtal::trace.beginBlock ( "Testing 2D/3D Image Reader" );
57 
58  typedef DGtal::ImageContainerBySTLVector<DGtal::Z3i::Domain, unsigned char> Image3D;
59  typedef DGtal::ImageContainerBySTLVector<DGtal::Z2i::Domain, unsigned char> Image2D;
60 
61 #ifdef WITH_ITK
62  typedef DGtal::ImageContainerBySTLVector<DGtal::Z3i::Domain, DGtal::uint16_t>
63  Image3D16bits;
64 #endif // WITH_ITK
65 
66  typedef DGtal::ImageContainerBySTLVector<DGtal::Z3i::Domain, DGtal::uint32_t> Image3D32bits;
67  typedef DGtal::ImageContainerBySTLVector<DGtal::Z3i::Domain, DGtal::uint64_t> Image3D64bits;
68  typedef DGtal::ImageContainerBySTLVector<DGtal::Z2i::Domain, unsigned int> Image2D32bits;
69 
70 
71  std::string filenameImage1 = testPath + "samples/cat10.vol";
72  Image3D anImportedImage1= DGtal::GenericReader<Image3D>::import(filenameImage1);
73  DGtal::Z3i::Domain domain1 = anImportedImage1.domain();
74  unsigned int size0Img1= domain1.upperBound()[0]-domain1.lowerBound()[0]+1;
75  unsigned int size1Img1= domain1.upperBound()[1]-domain1.lowerBound()[1]+1;
76  unsigned int size2Img1= domain1.upperBound()[2]-domain1.lowerBound()[2]+1;
77  DGtal::trace.info()<<"Vol image read: size[0]:" << size0Img1 ;
78  DGtal::trace.info()<<"size[1]: " << size1Img1;
79  DGtal::trace.info()<<"size[2]: " << size2Img1 << std::endl;
80  nbok += (size0Img1==40 && size1Img1==40 && size2Img1==40) ? 1 : 0;
81  nb++;
82  std::string filenameImage0 = testPath + "samples/test.longvol";
83  Image3D64bits anImportedImage0= DGtal::GenericReader<Image3D64bits>::import(filenameImage0);
84  DGtal::Z3i::Domain domain0 = anImportedImage0.domain();
85  unsigned int size0Img0= domain0.upperBound()[0]-domain0.lowerBound()[0]+1;
86  unsigned int size1Img0= domain0.upperBound()[1]-domain0.lowerBound()[1]+1;
87  unsigned int size2Img0= domain0.upperBound()[2]-domain0.lowerBound()[2]+1;
88  DGtal::trace.info()<<"Longvol image read: size[0]:" << size0Img0 ;
89  DGtal::trace.info()<<"size[1]: " << size1Img0;
90  DGtal::trace.info()<<"size[2]: " << size2Img0 << std::endl;
91  nbok += (size0Img0==16 && size1Img0==16 && size2Img0==16) ? 1 : 0;
92  nb++;
93  std::string filenameImage2 = testPath + "samples/cat10.pgm3d";
94  Image3D anImportedImage2= DGtal::GenericReader<Image3D>::import(filenameImage2);
95  DGtal::Z3i::Domain domain2 = anImportedImage2.domain();
96  unsigned int size0Img2= domain2.upperBound()[0]-domain2.lowerBound()[0]+1;
97  unsigned int size1Img2= domain2.upperBound()[1]-domain2.lowerBound()[1]+1;
98  unsigned int size2Img2= domain2.upperBound()[2]-domain2.lowerBound()[2]+1;
99  DGtal::trace.info()<<"Pgm3D image read: size[0]:" << size0Img2 ;
100  DGtal::trace.info()<<"size[1]: " << size1Img2;
101  DGtal::trace.info()<<"size[2]: " << size2Img2 << std::endl;;
102  nbok += (size0Img2==40 && size1Img2==40 && size2Img2==40) ? 1 : 0;
103  nb++;
104 
105  typedef DGtal::functors::Rescaling<unsigned char ,unsigned char > RescalFCT;
106  Image3D imageanImportedImage2 = DGtal::GenericReader<Image3D>::importWithValueFunctor( filenameImage2, RescalFCT(0, 120,
107  0, 255) );
108  domain2 = anImportedImage2.domain();
109  size0Img2= domain2.upperBound()[0]-domain2.lowerBound()[0]+1;
110  size1Img2= domain2.upperBound()[1]-domain2.lowerBound()[1]+1;
111  size2Img2= domain2.upperBound()[2]-domain2.lowerBound()[2]+1;
112  DGtal::trace.info()<<"Pgm3D image read (with scale functor) : size[0]:" << size0Img2 ;
113  DGtal::trace.info()<<"size[1]: " << size1Img2;
114  DGtal::trace.info()<<"size[2]: " << size2Img2 << std::endl;;
115  nbok += (size0Img2==40 && size1Img2==40 && size2Img2==40) ? 1 : 0;
116  nb++;
117 
118 
119 #ifdef WITH_HDF5
120  std::string filenameImageh5 = testPath + "samples/cat10.h5";
121  Image3D anImportedImageh5= DGtal::GenericReader<Image3D>::import(filenameImageh5);
122  DGtal::Z3i::Domain domainh5 = anImportedImageh5.domain();
123  unsigned int size0Imgh5= domainh5.upperBound()[0]-domainh5.lowerBound()[0]+1;
124  unsigned int size1Imgh5= domainh5.upperBound()[1]-domainh5.lowerBound()[1]+1;
125  unsigned int size2Imgh5= domainh5.upperBound()[2]-domainh5.lowerBound()[2]+1;
126  DGtal::trace.info()<<"HDF5 3D image read: size[0]:" << size0Imgh5;
127  DGtal::trace.info()<<"size[1]: " << size1Imgh5;
128  DGtal::trace.info()<<"size[2]: " << size2Imgh5 << std::endl;;
129  nbok += (size0Imgh5==40 && size1Imgh5==40 && size2Imgh5==40) ? 1 : 0;
130  nb++;
131 #endif
132  std::string filenameImage3 = testPath + "samples/contourS.pgm";
133  Image2D anImportedImage3= DGtal::GenericReader<Image2D>::import(filenameImage3);
134  DGtal::Z2i::Domain domain3 = anImportedImage3.domain();
135  unsigned int size0Img3= domain3.upperBound()[0]-domain3.lowerBound()[0]+1;
136  unsigned int size1Img3= domain3.upperBound()[1]-domain3.lowerBound()[1]+1;
137 
138  DGtal::trace.info()<<"Pgm image read: size[0]:" << size0Img3 ;
139  DGtal::trace.info()<<"size[1]: " << size1Img3 << std::endl;;
140  nbok += (size0Img3==185 && size1Img3==85 ) ? 1 : 0;
141  nb++;
142 
143  std::string filenameImage4 = testPath + "samples/raw32bits5x5x5.raw";
144  Image3D32bits anImportedImage4= DGtal::GenericReader<Image3D32bits>::import(filenameImage4, 5, 5, 5);
145  DGtal::Z3i::Domain domain4 = anImportedImage4.domain();
146  unsigned int size0Img4= domain4.upperBound()[0]-domain4.lowerBound()[0]+1;
147  unsigned int size1Img4= domain4.upperBound()[1]-domain4.lowerBound()[1]+1;
148  unsigned int size2Img4= domain4.upperBound()[2]-domain4.lowerBound()[2]+1;
149  DGtal::trace.info()<<"Raw32 bits 3D image read: size[0]:" << size0Img4;
150  DGtal::trace.info()<<"size[1]: " << size1Img4;
151  DGtal::trace.info()<<"size[2]: " << size2Img4 << std::endl;
152  DGtal::trace.info()<<"Image value of Point (2,3,4): " << anImportedImage4(DGtal::Z3i::Point(2,3,4)) << " (should be 250000*2*3*4) "<<std::endl;
153  nbok += (size0Img4==5 && size1Img4==5 && size2Img4==5 && anImportedImage4(DGtal::Z3i::Point(2,3,4))==250000*2*3*4) ? 1 : 0;
154  nb++;
155 
156  std::string filenameImage5 = testPath + "samples/raw32bits5x5.raw";
157  Image2D32bits anImportedImage5= DGtal::GenericReader<Image2D32bits>::import(filenameImage5, 5, 5);
158  DGtal::Z2i::Domain domain5 = anImportedImage5.domain();
159  unsigned int size0Img5= domain5.upperBound()[0]-domain5.lowerBound()[0]+1;
160  unsigned int size1Img5= domain5.upperBound()[1]-domain5.lowerBound()[1]+1;
161  DGtal::trace.info()<<"Raw32 bits 3D image read: size[0]:" << size0Img4;
162  DGtal::trace.info()<<"size[1]: " << size1Img5 << std::endl;
163  DGtal::trace.info()<<"Image value of Point (2,3): " << anImportedImage5(DGtal::Z2i::Point(2,3)) << " (should be 250000*2*3*4) "<<std::endl;
164  nbok += (size0Img5==5 && size1Img5==5 && anImportedImage5(DGtal::Z2i::Point(2,3))==250000*2*3*4) ? 1 : 0;
165  nb++;
166 
167 #ifdef WITH_ITK
168  std::string filenameImage6 = testPath + "samples/lobsterCroped16b.mhd";
169  typedef DGtal::functors::Rescaling<uint16_t, unsigned char> RescalFCT2;
170  RescalFCT2 resc2 = RescalFCT2( 0, 65535, 0, 255 );
171  Image3D anImportedImage6 =
172  DGtal::GenericReader<Image3D>::importWithValueFunctor( filenameImage6,
173  resc2 );
174  DGtal::Z3i::Domain domain6 = anImportedImage6.domain();
175  unsigned int size0Img6 =
176  domain6.upperBound()[ 0 ] - domain6.lowerBound()[ 0 ] + 1;
177  unsigned int size1Img6 =
178  domain6.upperBound()[ 1 ] - domain6.lowerBound()[ 1 ] + 1;
179  unsigned int size2Img6 =
180  domain6.upperBound()[ 2 ] - domain6.lowerBound()[ 2 ] + 1;
181 
182  DGtal::trace.info() << "mhd 16 bits 3D image read: size[0]:" << size0Img6;
183  DGtal::trace.info() << "size[1]: " << size1Img6 << std::endl;
184  DGtal::trace.info() << "Image value of Point (35,29, 3): "
185  << anImportedImage6( DGtal::Z3i::Point( 35, 29, 3 ) )
186  << " (should be" << resc2( 60400 ) << std::endl;
187  nbok +=
188  ( size0Img6 == 51 && size1Img6 == 91 && size2Img6 == 31 &&
189  anImportedImage6( DGtal::Z3i::Point( 35, 29, 3 ) ) == resc2( 60400 ) )
190  ? 1
191  : 0;
192  nb++;
193 
194  Image3D16bits anImportedImage7 =
195  DGtal::GenericReader<Image3D16bits>::import( filenameImage6 );
196  DGtal::trace.info() << "Image value of Point (35,29, 3): "
197  << anImportedImage7( DGtal::Z3i::Point( 35, 29, 3 ) )
198  << " (should be" << 60400 << " )" << std::endl;
199  nbok +=
200  ( anImportedImage7( DGtal::Z3i::Point( 35, 29, 3 ) ) == 60400 ) ? 1 : 0;
201  nb++;
202 #endif // WITH_ITK
203 
204  DGtal::trace.info() << "(" << nbok << "/" << nb << ") " << std::endl;
205  DGtal::trace.endBlock();
206 
207  return nbok == nb;
208 }
209 
211 // Standard services - public :
212 
213 int main( int argc, char** argv )
214 {
215  DGtal::trace.beginBlock ( "Testing class GenericReader" );
216  DGtal::trace.info() << "Args:";
217  for ( int i = 0; i < argc; ++i )
218  DGtal::trace.info() << " " << argv[ i ];
219  DGtal::trace.info() << std::endl;
220 
221  bool res = testGenericReader(); // && ... other tests
222  DGtal::trace.emphase() << ( res ? "Passed." : "Error." ) << std::endl;
223  DGtal::trace.endBlock();
224  return res ? 0 : 1;
225 }
226 // //
DGtal::HyperRectDomain::lowerBound
const Point & lowerBound() const
DGtal::HyperRectDomain::upperBound
const Point & upperBound() const
DGtal::Trace::beginBlock
void beginBlock(const std::string &keyword="")
DGtal::Trace::emphase
std::ostream & emphase()
DGtal::Trace::info
std::ostream & info()
DGtal::Trace::endBlock
double endBlock()
DGtal::trace
Trace trace
Definition: Common.h:153
DGtal::GenericReader::import
static TContainer import(const std::string &filename, std::vector< unsigned int > dimSpace=std::vector< unsigned int >())
DGtal::GenericReader::importWithValueFunctor
static TContainer importWithValueFunctor(const std::string &filename, const TFunctor &aFunctor, std::vector< unsigned int > dimSpace=std::vector< unsigned int >())
Definition: GenericReader.h:129
DGtal::functors::Rescaling
Aim: Functor allowing to rescale a value. Values of the initial scale [initMin,initMax] are rescaled ...
Definition: BasicFunctors.h:906
main
int main(int argc, char **argv)
Definition: testGenericReader.cpp:213
testGenericReader
bool testGenericReader()
Definition: testGenericReader.cpp:51