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
|
/* ----------------------------- MNI Header -----------------------------------
@NAME :
@DESCRIPTION: Simplified whole file in a memory MINC access library, using minc_1_rw interface
@COPYRIGHT :
Copyright 2007 Vladimir Fonov, McConnell Brain Imaging Centre,
Montreal Neurological Institute, McGill University.
Permission to use, copy, modify, and distribute this
software and its documentation for any purpose and without
fee is hereby granted, provided that the above copyright
notice appear in all copies. The author and McGill University
make no representations about the suitability of this
software for any purpose. It is provided "as is" without
express or implied warranty.
---------------------------------------------------------------------------- */
#include <iostream>
#include <math.h>
#include "minc_1_simple_rw.h"
namespace minc
{
const double minc_eps=1e-5;
bool is_same(minc_1_reader& one,minc_1_reader& two,bool verbose)
{
if(one.dim_no()!=two.dim_no())
{
if(verbose)
std::cerr<<"Unequal number of dimensions !"<<std::endl;
return false;
}
for(int j=0;j<5;j++)
{
if(one.ndim(j)!=two.ndim(j))
{
if(verbose)
std::cerr<<"Unequal dimension sizes"<<std::endl;
return false;
}
if(fabs(one.nstart(j)-two.nstart(j))>minc_eps)
{
if(verbose)
std::cerr<<"Unequal dimension sarts"<<std::endl;
return false;
}
if(fabs(one.nspacing(j)-two.nspacing(j))>minc_eps)
{
if(verbose)
std::cerr<<"Unequal dimension steps"<<std::endl;
return false;
}
for(int i=0;i<3;i++)
if(fabs(one.ndir_cos(j,i)-two.ndir_cos(j,i))>minc_eps)
{
if(verbose)
std::cerr<<"Unequal direction cosines"<<std::endl;
return false;
}
}
return true;
}
}
|
