Milena (Olena)  User documentation 2.0a Id
 All Classes Namespaces Functions Variables Typedefs Enumerator Groups Pages
dicom/get_header.hh
1 // Copyright (C) 2010 EPITA Research and Development Laboratory (LRDE)
2 //
3 // This file is part of Olena.
4 //
5 // Olena is free software: you can redistribute it and/or modify it under
6 // the terms of the GNU General Public License as published by the Free
7 // Software Foundation, version 2 of the License.
8 //
9 // Olena is distributed in the hope that it will be useful,
10 // but WITHOUT ANY WARRANTY; without even the implied warranty of
11 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
12 // General Public License for more details.
13 //
14 // You should have received a copy of the GNU General Public License
15 // along with Olena. If not, see <http://www.gnu.org/licenses/>.
16 //
17 // As a special exception, you may use this file as part of a free
18 // software project without restriction. Specifically, if other files
19 // instantiate templates or use macros or inline functions from this
20 // file, or you compile this file and link it with other files to produce
21 // an executable, this file does not by itself cause the resulting
22 // executable to be covered by the GNU General Public License. This
23 // exception does not however invalidate any other reasons why the
24 // executable file might be covered by the GNU General Public License.
25 
26 #ifndef MLN_IO_DICOM_GET_HEADER_HH
27 # define MLN_IO_DICOM_GET_HEADER_HH
28 
32 
33 # include <iostream>
34 # include <fstream>
35 
36 # include <gdcm-2.0/gdcmReader.h>
37 # include <gdcm-2.0/gdcmImageReader.h>
38 # include <gdcm-2.0/gdcmWriter.h>
39 # include <gdcm-2.0/gdcmDataSet.h>
40 # include <gdcm-2.0/gdcmAttribute.h>
41 
42 # include <mln/core/concept/image.hh>
43 # include <mln/core/routine/initialize.hh>
44 # include <mln/core/box_runstart_piter.hh>
45 # include <mln/core/pixel.hh>
46 # include <mln/data/memcpy_.hh>
47 # include <mln/util/array.hh>
48 
49 namespace mln
50 {
51 
52  namespace io
53  {
54 
55  namespace dicom
56  {
57 
59  struct dicom_header
60  {
61  // The number of dimensions.
62  unsigned dim;
63 
64  // The size in each dimension.
66  };
67 
68 
70  dicom_header get_header(const std::string& filename);
71 
72 
73 # ifndef MLN_INCLUDE_ONLY
74 
75 
76  dicom_header get_header(const std::string& filename)
77  {
78  trace::entering("mln::io::dicom::get_header");
79 
80  dicom_header header;
81 
82  gdcm::ImageReader r;
83  r.SetFileName(filename.c_str());
84  if (!r.Read())
85  {
86  std::cerr << "error: cannot open file '" << filename << "'!";
87  abort();
88  }
89 
90  gdcm::Image& image = r.GetImage();
91 
92  header.dim = image.GetNumberOfDimensions();
93  const unsigned int* dims = image.GetDimensions();
94 
95  for (unsigned i = 2; i < header.dim; ++i)
96  header.size.append(dims[i]); // sli, ...
97  for (unsigned i = 0; i < 2; ++i)
98  header.size.append(dims[i]); // row, col
99 
100  trace::exiting("mln::io::dicom::get_header");
101  return header;
102  }
103 
104 
105 # endif // ! MLN_INCLUDE_ONLY
106 
107  } // end of namespace mln::io::dicom
108 
109  } // end of namespace mln::io
110 
111 } // end of namespace mln
112 
113 #endif // ! MLN_IO_DICOM_GET_HEADER_HH
114