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) { }
}
}
|