Add scripts "ciao_img_rotcrop.sh" and "ciao_mkimg_det.sh".

Writen for ZHANG Chenhao to generate rotated and cropped image for
power spectrum analysis.

WARNING: the image generated by "ciao_mkimg_det.sh" in the DET
coordinate is very different than that in SKY coordinate. Watch out for
scientific usage!!!

1 files changed, 46 insertions, 0 deletions

diff --git a/scripts/ciao_mkimg_det.sh b/scripts/ciao_mkimg_det.sh
new file mode 100755
index 0000000..3e2d881
--- /dev/null
+++ b/scripts/ciao_mkimg_det.sh
@@ -0,0 +1,46 @@
+#!/bin/sh
+#
+# Make an image from evt2 file by binning the 'TDET' coordinates.
+# If the 'TDET' coordinate information is not available (e.g., blanksky),
+# the 'DET' coordinate is used for binning.
+#
+# Aaron LI
+# 2015/08/21
+#
+
+if [ $# -ne 2 ]; then
+    printf "Usage:\n"
+    printf "    `basename $0` <input:evt2> <output:img>\n"
+    exit 1
+fi
+
+EVT2="$1"
+OUTIMG="$2"
+
+# Test whether 'TEDT' coordinate exists?
+if dmlist ${EVT2} cols | grep -q tdet; then
+    COORDX="tdetx"
+    COORDY="tdety"
+else
+    printf "WARNING: tdet coordinate NOT exist! use det instead.\n"
+    COORDX="detx"
+    COORDY="dety"
+fi
+
+# Get TDET/DET coordinate min & max values
+# COORDX:
+punlearn dmstat
+dmstat "${EVT2}[cols ${COORDX}]" >/dev/null 2>&1
+XMIN=`pget dmstat out_min`
+XMAX=`pget dmstat out_max`
+dmstat "${EVT2}[cols ${COORDY}]" >/dev/null 2>&1
+YMIN=`pget dmstat out_min`
+YMAX=`pget dmstat out_max`
+
+BINSPEC="[bin ${COORDX}=${XMIN}:${XMAX}:1,${COORDY}=${YMIN}:${YMAX}:1]"
+
+punlearn dmcopy
+CMD="dmcopy \"${EVT2}${BINSPEC}\" ${OUTIMG}"
+printf "CMD: ${CMD}\n"
+eval ${CMD}
+