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
77
78
79
80
81
82
83
84
85
|
/* Copyright (c) 2008-2025 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 <cmath>
#include "command.h"
#include "image.h"
#include "image_helpers.h"
#include "algo/loop.h"
using namespace MR;
using namespace App;
void usage ()
{
AUTHOR = "Max Pietsch (maximilian.pietsch@kcl.ac.uk)";
SYNOPSIS = "Create bitwise checkerboard image";
ARGUMENTS
+ Argument ("input", "the input image to be used as a template.").type_image_in ()
+ Argument ("output", "the output binary image mask.").type_image_out ();
OPTIONS
+ Option ("tiles", "specify the number of tiles in any direction")
+ Argument ("value").type_integer()
+ Option ("invert", "invert output binary mask.")
+ Option ("nan", "use NaN as the output zero value.");
}
void run ()
{
size_t ntiles = 5;
auto opt = get_options ("tiles");
if (opt.size()) {
ntiles = opt[0][0];
}
bool invert = get_options ("invert").size();
const bool use_NaN = get_options ("nan").size();
auto in = Image<float>::open (argument[0]);
check_3D_nonunity (in);
size_t patchwidth_x = ceil((float) in.size(0) / ntiles);
size_t patchwidth_y = ceil((float) in.size(1) / ntiles);
size_t patchwidth_z = ceil((float) in.size(2) / ntiles);
Header header_out (in);
header_out.datatype() = use_NaN ? DataType::Float32 : DataType::Bit;
auto out = Image<float>::create (argument[1], header_out);
float zero = use_NaN ? NAN : 0.0;
float one = 1.0;
if (invert) std::swap (zero, one);
for (auto l = Loop(in) (in, out); l; ++l) {
const size_t x = in.index(0);
const size_t y = in.index(1);
const size_t z = in.index(2);
size_t xpatch = (x-(x%patchwidth_x))/patchwidth_x;
size_t ypatch = (y-(y%patchwidth_y))/patchwidth_y;
size_t zpatch = (z-(z%patchwidth_z))/patchwidth_z;
out.value() = (
(xpatch%2 + ypatch%2+ zpatch%2)%2==0) ? one : zero;
}
}
|
