Milena (Olena)
User documentation 2.0a Id
|
00001 // Copyright (C) 2010 EPITA Research and Development Laboratory (LRDE) 00002 // 00003 // This file is part of Olena. 00004 // 00005 // Olena is free software: you can redistribute it and/or modify it under 00006 // the terms of the GNU General Public License as published by the Free 00007 // Software Foundation, version 2 of the License. 00008 // 00009 // Olena is distributed in the hope that it will be useful, 00010 // but WITHOUT ANY WARRANTY; without even the implied warranty of 00011 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU 00012 // General Public License for more details. 00013 // 00014 // You should have received a copy of the GNU General Public License 00015 // along with Olena. If not, see <http://www.gnu.org/licenses/>. 00016 // 00017 // As a special exception, you may use this file as part of a free 00018 // software project without restriction. Specifically, if other files 00019 // instantiate templates or use macros or inline functions from this 00020 // file, or you compile this file and link it with other files to produce 00021 // an executable, this file does not by itself cause the resulting 00022 // executable to be covered by the GNU General Public License. This 00023 // exception does not however invalidate any other reasons why the 00024 // executable file might be covered by the GNU General Public License. 00025 00026 #ifndef MLN_IO_DICOM_GET_HEADER_HH 00027 # define MLN_IO_DICOM_GET_HEADER_HH 00028 00032 00033 # include <iostream> 00034 # include <fstream> 00035 00036 # include <gdcm-2.0/gdcmReader.h> 00037 # include <gdcm-2.0/gdcmImageReader.h> 00038 # include <gdcm-2.0/gdcmWriter.h> 00039 # include <gdcm-2.0/gdcmDataSet.h> 00040 # include <gdcm-2.0/gdcmAttribute.h> 00041 00042 # include <mln/core/concept/image.hh> 00043 # include <mln/core/routine/initialize.hh> 00044 # include <mln/core/box_runstart_piter.hh> 00045 # include <mln/core/pixel.hh> 00046 # include <mln/data/memcpy_.hh> 00047 # include <mln/util/array.hh> 00048 00049 namespace mln 00050 { 00051 00052 namespace io 00053 { 00054 00055 namespace dicom 00056 { 00057 00059 struct dicom_header 00060 { 00061 // The number of dimensions. 00062 unsigned dim; 00063 00064 // The size in each dimension. 00065 util::array<unsigned> size; 00066 }; 00067 00068 00070 dicom_header get_header(const std::string& filename); 00071 00072 00073 # ifndef MLN_INCLUDE_ONLY 00074 00075 00076 dicom_header get_header(const std::string& filename) 00077 { 00078 trace::entering("mln::io::dicom::get_header"); 00079 00080 dicom_header header; 00081 00082 gdcm::ImageReader r; 00083 r.SetFileName(filename.c_str()); 00084 if (!r.Read()) 00085 { 00086 std::cerr << "error: cannot open file '" << filename << "'!"; 00087 abort(); 00088 } 00089 00090 gdcm::Image& image = r.GetImage(); 00091 00092 header.dim = image.GetNumberOfDimensions(); 00093 const unsigned int* dims = image.GetDimensions(); 00094 00095 for (unsigned i = 2; i < header.dim; ++i) 00096 header.size.append(dims[i]); // sli, ... 00097 for (unsigned i = 0; i < 2; ++i) 00098 header.size.append(dims[i]); // row, col 00099 00100 trace::exiting("mln::io::dicom::get_header"); 00101 return header; 00102 } 00103 00104 00105 # endif // ! MLN_INCLUDE_ONLY 00106 00107 } // end of namespace mln::io::dicom 00108 00109 } // end of namespace mln::io 00110 00111 } // end of namespace mln 00112 00113 #endif // ! MLN_IO_DICOM_GET_HEADER_HH 00114