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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
|
tests/test-rtnlinv: traj scale phantom nufft resize rtnlinv fmac nrmse
set -e ; mkdir $(TESTS_TMP) ; cd $(TESTS_TMP) ;\
$(TOOLDIR)/traj -r -x128 -y21 -t5 traj.ra ;\
$(TOOLDIR)/scale 0.5 traj.ra traj2.ra ;\
$(TOOLDIR)/phantom -s8 -k -t traj2.ra ksp.ra ;\
$(TOOLDIR)/nufft -a traj.ra ksp.ra I.ra ;\
$(TOOLDIR)/rtnlinv -N -S -i9 -t traj2.ra ksp.ra r.ra c.ra ;\
$(TOOLDIR)/resize -c 0 64 1 64 c.ra c2.ra ;\
$(TOOLDIR)/fmac r.ra c2.ra x.ra ;\
$(TOOLDIR)/nufft traj2.ra x.ra k2.ra ;\
$(TOOLDIR)/nrmse -t 0.05 ksp.ra k2.ra ;\
rm *.ra ; cd .. ; rmdir $(TESTS_TMP)
touch $@
tests/test-rtnlinv-precomp: traj scale phantom ones repmat fft nufft rtnlinv nrmse
set -e ; mkdir $(TESTS_TMP) ; cd $(TESTS_TMP) ;\
$(TOOLDIR)/traj -r -x128 -y21 -t5 traj.ra ;\
$(TOOLDIR)/scale 0.5 traj.ra traj2.ra ;\
$(TOOLDIR)/phantom -s8 -k -t traj2.ra ksp.ra ;\
$(TOOLDIR)/ones 3 1 128 21 o.ra ;\
$(TOOLDIR)/repmat 10 5 o.ra o2.ra ;\
$(TOOLDIR)/nufft -a traj.ra o2.ra psf.ra ;\
$(TOOLDIR)/nufft -a traj.ra ksp.ra _adj.ra ;\
$(TOOLDIR)/scale 2 _adj.ra adj.ra ;\
$(TOOLDIR)/fft -u 7 psf.ra mtf.ra ;\
$(TOOLDIR)/fft -u 7 adj.ra ksp2.ra ;\
$(TOOLDIR)/scale 4. mtf.ra mtf2.ra ;\
$(TOOLDIR)/rtnlinv -w1. -N -i9 -p mtf2.ra ksp2.ra r1.ra c1.ra ;\
$(TOOLDIR)/rtnlinv -w1. -N -i9 -t traj2.ra ksp.ra r2.ra c2.ra ;\
$(TOOLDIR)/nrmse -t 0.000001 r2.ra r1.ra ;\
$(TOOLDIR)/nrmse -t 0.000001 c2.ra c1.ra ;\
rm *.ra ; cd .. ; rmdir $(TESTS_TMP)
touch $@
tests/test-rtnlinv-nlinv-noncart: traj scale phantom rtnlinv nlinv nrmse
set -e ; mkdir $(TESTS_TMP) ; cd $(TESTS_TMP) ;\
$(TOOLDIR)/traj -r -x128 -y21 traj.ra ;\
$(TOOLDIR)/scale 0.5 traj.ra traj2.ra ;\
$(TOOLDIR)/phantom -s8 -k -t traj2.ra ksp.ra ;\
$(TOOLDIR)/nlinv -w1. -N -i9 -t traj2.ra ksp.ra r1.ra c1.ra ;\
$(TOOLDIR)/rtnlinv -w1. -N -i9 -t traj2.ra ksp.ra r2.ra c2.ra ;\
$(TOOLDIR)/nrmse -t 0.000001 r2.ra r1.ra ;\
$(TOOLDIR)/nrmse -t 0.000001 c2.ra c1.ra ;\
rm *.ra ; cd .. ; rmdir $(TESTS_TMP)
touch $@
tests/test-rtnlinv-nlinv-pseudocart: scale phantom ones rtnlinv nlinv nrmse
set -e ; mkdir $(TESTS_TMP) ; cd $(TESTS_TMP) ;\
$(TOOLDIR)/phantom -s8 -k -x 128 ksp.ra ;\
$(TOOLDIR)/ones 2 128 128 psf.ra ;\
$(TOOLDIR)/nlinv -w1. -N -i9 -f1. -n ksp.ra r1.ra c1.ra ;\
$(TOOLDIR)/rtnlinv -w1. -N -i9 -f1. -p psf.ra ksp.ra r2.ra c2.ra ;\
$(TOOLDIR)/nrmse -t 0. r2.ra r1.ra ;\
$(TOOLDIR)/nrmse -t 0. c2.ra c1.ra ;\
rm *.ra ; cd .. ; rmdir $(TESTS_TMP)
touch $@
tests/test-rtnlinv-maps-dims: phantom ones rtnlinv show
set -e ; mkdir $(TESTS_TMP) ; cd $(TESTS_TMP) ;\
$(TOOLDIR)/phantom -s1 -k ksp.ra ;\
$(TOOLDIR)/ones 2 128 128 psf.ra ;\
$(TOOLDIR)/rtnlinv -m4 -S -i1 -ppsf.ra ksp.ra r.ra c.ra ;\
$(TOOLDIR)/rtnlinv -m4 -U -S -i1 -ppsf.ra ksp.ra rU.ra cU.ra ;\
$(TOOLDIR)/rtnlinv -m4 -N -S -i1 -ppsf.ra ksp.ra rN.ra cN.ra ;\
$(TOOLDIR)/rtnlinv -m4 -N -U -S -i1 -ppsf.ra ksp.ra rNU.ra cNU.ra ;\
[ 1 -eq `$(TOOLDIR)/show -d 4 r.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 c.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 rU.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 cU.ra` ] &&\
[ 1 -eq `$(TOOLDIR)/show -d 4 rN.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 cN.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 rNU.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 cNU.ra` ] &&\
true
touch $@
tests/test-rtnlinv-noncart-maps-dims: traj phantom rtnlinv show
set -e ; mkdir $(TESTS_TMP) ; cd $(TESTS_TMP) ;\
$(TOOLDIR)/traj -r -y21 traj.ra ;\
$(TOOLDIR)/phantom -s1 -k -ttraj.ra ksp.ra ;\
$(TOOLDIR)/rtnlinv -m4 -S -i1 -ttraj.ra ksp.ra r.ra c.ra ;\
$(TOOLDIR)/rtnlinv -m4 -U -S -i1 -ttraj.ra ksp.ra rU.ra cU.ra ;\
$(TOOLDIR)/rtnlinv -m4 -N -S -i1 -ttraj.ra ksp.ra rN.ra cN.ra ;\
$(TOOLDIR)/rtnlinv -m4 -N -U -S -i1 -ttraj.ra ksp.ra rNU.ra cNU.ra ;\
[ 1 -eq `$(TOOLDIR)/show -d 4 r.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 c.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 rU.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 cU.ra` ] &&\
[ 1 -eq `$(TOOLDIR)/show -d 4 rN.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 cN.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 rNU.ra` ] &&\
[ 4 -eq `$(TOOLDIR)/show -d 4 cNU.ra` ] &&\
true
touch $@
TESTS += tests/test-rtnlinv tests/test-rtnlinv-precomp tests/test-rtnlinv-nlinv-noncart tests/test-rtnlinv-nlinv-pseudocart
TESTS += tests/test-rtnlinv-maps-dims tests/test-rtnlinv-noncart-maps-dims
#TESTS += tests/test-rtnlinv-precomp
|
