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
|
/*
* 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 "datatype.h"
#include "image.h"
#include "fixel/helpers.h"
#include "diff_images.h"
using namespace MR;
using namespace App;
void usage ()
{
AUTHOR = "David Raffelt (david.raffelt@florey.edu.au) and Robert E. Smith (robert.smith@florey.edu.au)";
SYNOPSIS = "Compare two fixel images for differences, within specified tolerance";
ARGUMENTS
+ Argument ("fixel1", "fixel directory.").type_directory_in()
+ Argument ("fixel2", "another fixel directory.").type_directory_in();
OPTIONS
+ Testing::Diff_Image_Options;
}
void run ()
{
std::string fixel_directory1 = argument[0];
Fixel::check_fixel_directory (fixel_directory1);
std::string fixel_directory2 = argument[1];
Fixel::check_fixel_directory (fixel_directory2);
if (fixel_directory1 == fixel_directory2)
throw Exception ("Input fixel directorys are the same");
auto dir_walker1 = Path::Dir (fixel_directory1);
std::string fname;
while ((fname = dir_walker1.read_name()).size()) {
auto in1 = Image<cdouble>::open (Path::join (fixel_directory1, fname));
std::string filename2 = Path::join (fixel_directory2, fname);
if (!Path::exists (filename2))
throw Exception ("File (" + fname + ") exists in fixel directory (" + fixel_directory1 + ") but not in fixel directory (" + fixel_directory2 + ") ");
auto in2 = Image<cdouble>::open (filename2);
Testing::diff_images (in1, in2);
}
auto dir_walker2 = Path::Dir (fixel_directory2);
while ((fname = dir_walker2.read_name()).size()) {
std::string filename1 = Path::join (fixel_directory1, fname);
if (!Path::exists (filename1))
throw Exception ("File (" + fname + ") exists in fixel directory (" + fixel_directory2 + ") but not in fixel directory (" + fixel_directory1 + ") ");
}
CONSOLE ("data checked OK");
}
|
