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
|
#!/bin/bash
NUMPARAMS=$#
MAXITERATIONS=30x90x20
if [ $NUMPARAMS -lt 2 ]
then
echo " USAGE :: "
echo " sh $0 ImageDimension image.ext "
exit
fi
#initialization, here, is unbiased
DIM=$1
if [ ${#DIM} -gt 1 ]
then
echo " Problem with specified ImageDimension => User Specified Value = $DIM "
exit
fi
FIXED=$2
if [ ${#FIXED} -lt 1 -o ! -f $FIXED ]
then
echo " Problem with specified Fixed Image => User Specified Value = $FIXED "
exit
fi
OUTPUTNAME=` echo $FIXED | cut -d '.' -f 1 `
if [[ ! -s ${OUTPUTNAME}repaired.nii.gz ]] ; then
N3BiasFieldCorrection $DIM $FIXED ${OUTPUTNAME}repaired.nii.gz 4
fi
ThresholdImage $DIM ${OUTPUTNAME}repaired.nii.gz ${OUTPUTNAME}mask.nii.gz 0.1 99999
WM=${OUTPUTNAME}seg.nii.gz
if [[ ! -s $WM ]] ; then
Apocrita $DIM -x ${OUTPUTNAME}mask.nii.gz -m [ 0.5,1,0,0 ] -o [ ${OUTPUTNAME}seg.nii.gz ] -i Otsu[ ${OUTPUTNAME}repaired.nii.gz,3 ] -h 0
ThresholdImage $DIM $WM $WM 3 3
ImageMath $DIM $WM GetLargestComponent $WM
fi
OUT=${OUTPUTNAME}
if [[ ! -s ${OUT}max.nii.gz ]] ; then
# the closing operation - arachnoid surface
ImageMath 3 ${OUT}max.nii.gz MD $WM 15
ImageMath 3 ${OUT}max.nii.gz ME ${OUT}max.nii.gz 15
# unfolded or lissencephalic surface
# we would prefer to do this in a topology preserving way
ImageMath 3 ${OUT}min.nii.gz ME $WM 5
ImageMath 3 ${OUT}min.nii.gz MD ${OUT}min.nii.gz 5
# distance transform from arachnoid surface to inside of brain
# perhaps should be a geodesic distance transform
ImageMath 3 ${OUT}max.nii.gz Neg ${OUT}max.nii
ImageMath 3 ${OUT}dist.nii.gz D ${OUT}max.nii
fi
# multiply the distance transform by the "unfolded" surface
# this gives the "feature" surface/image that could be input
# to an image registration similarity metric
ImageMath 3 ${OUT}lohmann.nii.gz m ${OUT}dist.nii.gz ${OUT}min.nii.gz
# surface stuff for visualization
ImageMath 3 ${OUT}surf.nii.gz ME ${OUT}min.nii.gz 2
ImageMath 3 ${OUT}surf.nii.gz - ${OUT}min.nii.gz ${OUT}surf.nii.gz
ImageMath 3 ${OUT}maxvals.nii.gz m ${OUT}lohmann.nii.gz ${OUT}surf.nii.gz
# find sulci using curvature of the WM/lohmann surface
SurfaceCurvature ${OUT}lohmann.nii.gz ${OUT}lohmannk.nii.gz 1.5
ThresholdImage 3 ${OUT}lohmann.nii.gz ${OUT}lohmann.nii.gz 0.001 9999
ImageMath 3 ${OUT}lohmannk.nii.gz Normalize ${OUT}lohmannk.nii.gz
ImageMath 3 ${OUT}lohmannk.nii.gz m ${OUT}lohmannk.nii.gz ${OUT}lohmann.nii.gz
exit
|
