#!/bin/sh

if test "$1" != "" ; then
  bin=$1
else
  bin="@prefix@/bin/"
fi

# Generate phantom image

${bin}phm2sdf sample-phm.sdf 256 256 --nsample 2 --phantom herman
if [ -f sample-phm.sdf ] ; then
  ${bin}sdf2img sample-phm.sdf sample-phm.png --format png
  ${bin}sdf2img sample-phm.sdf sample-phm16.png --format png16
fi

# Simulate CT data collection and generate raysum sinugram for display
${bin}phm2rs  sample-rs.rs 367 320 --nray 2  --phantom herman
if [ -f sample-rs.rs ]; then
  ${bin}rs2sdf  sample-rs.rs sample-rs.sdf
fi
if [ -f sample-rs.sdf ]; then
  ${bin}sdf2img sample-rs.sdf sample-rs.png --format png
  ${bin}sdf2img sample-rs.sdf sample-rs16.png --format png16
fi

# Reconstruct raysums and generate image for display
${bin}ctrec   sample-rs.rs sample-rec.sdf 256 256 
if [ -f sample-rec.sdf ]; then 
  ${bin}sdf2img sample-rec.sdf sample-rec.png --format png
  ${bin}sdf2img sample-rec.sdf sample-rec16.png --format png16
fi

# Files sample-phm.png, sample-rs.png, and sample-rec.png are ready for display