Adding phantoms

bkg.out

@@ -0,0 +1 @@
+1 0 0 25 1

fisher1D.py

@@ -0,0 +1,48 @@
+import numpy
+import matplotlib.pyplot
+import scipy.special
+
+def generateAlgebra(w,o,bw,FWHM,R):
+    sigma=FWHM/2.35
+    ctr=0.5*o.size
+    #w=numpy.zeros([n,n])
+    q=bw/numpy.sqrt(2)
+    for i in numpy.arange(o.size):
+        for j in numpy.arange(o.size):
+            w[i,j]=0.5*(scipy.special.erf(q*(j-i+0.5)/sigma)-
+                scipy.special.erf(q*(j-i-0.5)/sigma))
+
+    #o=numpy.zeros(n)
+    p=0
+    for i in numpy.arange(o.size):
+        if numpy.abs(i-ctr)<R/bw:
+            o[i]=1
+            p+=1;
+    o/=numpy.sqrt(p)
+
+
+def SNR(o,w):
+    return numpy.dot(numpy.dot(o,w),o)/\
+            numpy.sqrt(numpy.dot(numpy.dot(numpy.dot(o,w),w),o))
+
+def Noise(o,w):
+    return numpy.sqrt(numpy.dot(numpy.dot(numpy.dot(o,w),w),o))
+
+def Signal(o,w):
+    return numpy.dot(numpy.dot(o,w),o)
+
+def graphSNR(n,bw,x,R,var):
+    s=numpy.zeros(x.size)
+    w=numpy.zeros([n,n])
+    o=numpy.zeros(n)
+    for i in numpy.arange(x.size):
+        generateAlgebra(w,o,bw,x[i],R)
+        if var=="SNR":
+            s[i]=SNR(o,w)
+        if var=="Signal":
+            s[i]=Signal(o,w)
+        if var=="Noise":
+            s[i]=Noise(o,w)
+    return s
+
+

phantom_0_rw_0p6.out

@@ -0,0 +1,74 @@
+1 5.06906 4.20893 0.6 1
+1 3.22525 5.74528 0.6 1
+1 7.32148 5.03754 0.6 1
+1 5.47767 6.57389 0.6 1
+1 3.63386 8.11024 0.6 1
+1 9.5739 5.86615 0.6 1
+1 7.73009 7.4025 0.6 1
+1 5.88628 8.93885 0.6 1
+1 4.04248 10.4752 0.6 1
+1 11.8263 6.69476 0.6 1
+1 9.98252 8.23111 0.6 1
+1 8.13871 9.76746 0.6 1
+1 6.2949 11.3038 0.6 1
+1 4.45109 12.8402 0.6 1
+1 14.0787 7.52338 0.6 1
+1 12.2349 9.05972 0.6 1
+1 10.3911 10.5961 0.6 1
+1 8.54732 12.1324 0.6 1
+1 6.70351 13.6688 0.6 1
+1 4.8597 15.2051 0.6 1
+1 16.3312 8.35199 0.6 1
+1 14.4874 9.88834 0.6 1
+1 12.6435 11.4247 0.6 1
+1 10.7997 12.961 0.6 1
+1 8.95593 14.4974 0.6 1
+1 7.11212 16.0337 0.6 1
+1 5.26831 17.5701 0.6 1
+1 16.7398 10.7169 0.6 1
+1 14.896 12.2533 0.6 1
+1 13.0522 13.7896 0.6 1
+1 11.2084 15.326 0.6 1
+1 9.36454 16.8623 0.6 1
+1 7.52073 18.3987 0.6 1
+1 6.79604 -2.79306 0.8 1
+1 7.33875 0.360586 0.8 1
+1 9.25581 -4.83988 0.8 1
+1 9.79853 -1.68624 0.8 1
+1 10.3412 1.4674 0.8 1
+1 11.7156 -6.88671 0.8 1
+1 12.2583 -3.73306 0.8 1
+1 12.801 -0.579422 0.8 1
+1 13.3437 2.57422 0.8 1
+1 14.1754 -8.93353 0.8 1
+1 14.7181 -5.77989 0.8 1
+1 15.2608 -2.62625 0.8 1
+1 15.8035 0.527395 0.8 1
+1 16.3462 3.68104 0.8 1
+1 17.1779 -7.82672 0.8 1
+1 17.7206 -4.67307 0.8 1
+1 18.2633 -1.51943 0.8 1
+1 18.806 1.63421 0.8 1
+1 1.12952 -9.98407 1.2 1
+1 5.63215 -8.32085 1.2 1
+1 0.318599 -14.7151 1.2 1
+1 4.82123 -13.0519 1.2 1
+1 9.32386 -11.3886 1.2 1
+1 -0.492318 -19.4461 1.2 1
+1 4.01031 -17.7829 1.2 1
+1 8.51294 -16.1196 1.2 1
+1 13.0156 -14.4564 1.2 1
+1 -3.83164 -4.6172 1.6 1
+1 -9.83367 -6.83884 1.6 1
+1 -4.90866 -10.9259 1.6 1
+1 -15.8357 -9.06047 1.6 1
+1 -10.9107 -13.1476 1.6 1
+1 -5.98568 -17.2346 1.6 1
+1 -7.09324 1.2353 2 1
+1 -13.2324 6.36465 2 1
+1 -14.605 -1.51673 2 1
+1 -2.61918 7.13442 2.4 1
+1 -0.978426 16.5932 2.4 1
+1 -9.99032 13.2847 2.4 1
+2 21.4 22.6 1
+1 0 0 25 0.01