File: mosaic.cpp

package info (click to toggle)
mrtrix3 3.0.4-1
  • links: PTS, VCS
  • area: main
  • in suites: forky, sid, trixie
  • size: 13,712 kB
  • sloc: cpp: 129,776; python: 9,494; sh: 593; makefile: 234; xml: 47
file content (76 lines) | stat: -rw-r--r-- 2,516 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
/* Copyright (c) 2008-2022 the MRtrix3 contributors.
 *
 * This Source Code Form is subject to the terms of the Mozilla Public
 * License, v. 2.0. If a copy of the MPL was not distributed with this
 * file, You can obtain one at http://mozilla.org/MPL/2.0/.
 *
 * Covered Software is provided under this License on an "as is"
 * basis, without warranty of any kind, either expressed, implied, or
 * statutory, including, without limitation, warranties that the
 * Covered Software is free of defects, merchantable, fit for a
 * particular purpose or non-infringing.
 * See the Mozilla Public License v. 2.0 for more details.
 *
 * For more details, see http://www.mrtrix.org/.
 */

#include <limits>

#include "app.h"
#include "progressbar.h"
#include "header.h"
#include "image_io/mosaic.h"

namespace MR
{
  namespace ImageIO
  {


    void Mosaic::load (const Header& header, size_t)
    {
      if (files.empty())
        throw Exception ("no files specified in header for image \"" + header.name() + "\"");

      assert (header.datatype().bits() > 1);

      size_t bytes_per_segment = header.datatype().bytes() * segsize;
      if (files.size() * bytes_per_segment > std::numeric_limits<size_t>::max())
        throw Exception ("image \"" + header.name() + "\" is larger than maximum accessible memory");

      DEBUG ("loading mosaic image \"" + header.name() + "\"...");
      addresses.resize (1);
      addresses[0].reset (new uint8_t [files.size() * bytes_per_segment]);
      if (!addresses[0])
        throw Exception ("failed to allocate memory for image \"" + header.name() + "\"");

      ProgressBar progress ("reformatting DICOM mosaic images", slices*files.size());
      uint8_t* data = addresses[0].get();
      for (size_t n = 0; n < files.size(); n++) {
        File::MMap file (files[n], false, false, m_xdim * m_ydim * header.datatype().bytes());
        size_t nx = 0, ny = 0;
        for (size_t z = 0; z < slices; z++) {
          size_t ox = nx*xdim;
          size_t oy = ny*ydim;
          for (size_t y = 0; y < ydim; y++) {
            memcpy (data, file.address() + header.datatype().bytes() * (ox + m_xdim* (y+oy)), xdim * header.datatype().bytes());
            data += xdim * header.datatype().bytes();
          }
          nx++;
          if (nx >= m_xdim / xdim) {
            nx = 0;
            ny++;
          }
          ++progress;
        }
      }

      segsize = std::numeric_limits<size_t>::max();
    }

    void Mosaic::unload (const Header& header) { }

  }
}