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
|
/* 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/.
*/
#ifndef __image_transform_h__
#define __image_transform_h__
#include "types.h"
namespace MR
{
class Transform { MEMALIGN(Transform)
public:
//! An object for transforming between voxel, scanner and image coordinate spaces
template <class HeaderType>
Transform (const HeaderType& header) :
voxelsize (header.spacing(0), header.spacing(1), header.spacing(2)),
voxel2scanner (header.transform() * voxelsize),
scanner2voxel (voxel2scanner.inverse()),
image2scanner (header.transform()),
scanner2image (image2scanner.inverse()) { }
Transform (const Transform&) = default;
Transform (Transform&&) = default;
Transform& operator= (const Transform&) = delete;
Transform& operator= (Transform&&) = delete;
const Eigen::DiagonalMatrix<default_type, 3> voxelsize;
const transform_type voxel2scanner, scanner2voxel, image2scanner, scanner2image;
template <class HeaderType>
static inline transform_type get_default (const HeaderType& header) {
transform_type M;
M.setIdentity();
M.translation() = Eigen::Vector3d (
-0.5 * (header.size (0)-1) * header.spacing (0),
-0.5 * (header.size (1)-1) * header.spacing (1),
-0.5 * (header.size (2)-1) * header.spacing (2)
);
return M;
}
};
}
#endif
|