X-Git-Url: http://git.kpe.io/?a=blobdiff_plain;f=src%2Fsample-ctsim.sh.in;fp=src%2Fsample-ctsim.sh.in;h=efc14e2fa724d4b0c9b67b0ebba74ab6d24c3a2b;hb=2f3d6e2580db607105bb072b13e4aff453ae4495;hp=0000000000000000000000000000000000000000;hpb=08f34bf3ba14d4f436f4d2ef0ee5af1d6eb266ac;p=ctsim.git

diff --git a/src/sample-ctsim.sh.in b/src/sample-ctsim.sh.in
new file mode 100755
index 0000000..efc14e2
--- /dev/null
+++ b/src/sample-ctsim.sh.in
@@ -0,0 +1,41 @@
+#!/bin/sh
+
+if test "$1" != "" ; then
+  bin=$1
+else
+  bin="@prefix@/bin/"
+fi
+
+if test "$1" = "clean" ; then
+  rm -f sample-phm.png sample-phm16.png sample-phm.if sample-pj.pj sample-pj.if sample-pj.png sample-pj16.png sample-rec.if sample-rec.png sample-rec16.png
+  exit
+fi
+
+# Generate phantom image
+
+${bin}phm2if sample-phm.if 256 256 --nsample 2 --phantom herman
+if [ -f sample-phm.if ] ; then
+  ${bin}if2img sample-phm.if sample-phm.png --format png
+  ${bin}if2img sample-phm.if sample-phm16.png --format png16
+fi
+
+# Simulate CT data collection and generate raysum sinugram for display
+${bin}phm2pj  sample-pj.pj 367 320 --nray 2  --phantom herman
+if [ -f sample-pj.pj ]; then
+  ${bin}pj2if  sample-pj.pj sample-pj.if
+fi
+if [ -f sample-pj.if ]; then
+  ${bin}if2img sample-pj.if sample-pj.png --format png
+  ${bin}if2img sample-pj.if sample-pj16.png --format png16
+fi
+
+# Reconstruct raysums and generate image for display
+${bin}pjrec   sample-pj.pj sample-rec.if 256 256 
+if [ -f sample-rec.if ]; then 
+  ${bin}if2img sample-rec.if sample-rec.png --format png
+  ${bin}if2img sample-rec.if sample-rec16.png --format png16
+
+  ${bin}ifinfo sample-phm.if sample-rec.if
+fi
+
+# Files sample-phm.png, sample-pj.png, and sample-rec.png are ready for display