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
|
/*
* Copyright (c) 2008-2018 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/
*
* MRtrix3 is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty
* of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
*
* For more details, see http://www.mrtrix.org/
*/
#include "command.h"
#include "image.h"
#include "algo/loop.h"
using namespace MR;
using namespace App;
void usage ()
{
AUTHOR = "Robert E. Smith (robert.smith@florey.edu.au)";
SYNOPSIS = "Print out the locations of all non-zero voxels in a mask image";
DESCRIPTION
+ "If no destination file is specified, the voxel locations will be "
"printed to stdout.";
ARGUMENTS
+ Argument ("input", "the input image.").type_image_in()
+ Argument ("output", "the (optional) output text file.").type_file_out().optional();
}
void run ()
{
auto H = Header::open (argument[0]);
if (H.datatype() != DataType::Bit)
WARN ("Input is not a genuine boolean mask image");
auto in = H.get_image<bool>();
vector< Eigen::ArrayXi > locations;
for (auto l = Loop(in) (in); l; ++l) {
if (in.value()) {
Eigen::ArrayXi this_voxel (in.ndim());
for (size_t axis = 0; axis != in.ndim(); ++axis)
this_voxel[axis] = in.index (axis);
locations.push_back (std::move (this_voxel));
}
}
Eigen::ArrayXXi prettyprint (locations.size(), in.ndim());
for (size_t row = 0; row != locations.size(); ++row)
prettyprint.row (row) = std::move (locations[row]);
INFO ("Printing locations of " + str(prettyprint.rows()) + " non-zero voxels");
if (argument.size() == 2)
save_matrix (prettyprint, argument[1]);
else
std::cout << prettyprint;
}
|
