summaryrefslogtreecommitdiffstats
path: root/ds9_image.sh
blob: 76a03d56cef72eb1726001e397fff92c5ca7916d (plain)
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
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
#!/bin/sh
#
# Use DS9 to visualize the ACIS image, with useful/handy
# arguments/options passed.
# Also touch the output image filename for easier save.
#
# Aaron LI
# Created: 2016-04-16
# Updated: 2016-05-08
#

# Default parameters
FILE_PATTERN="img_c*_fill.fits"
REG_FILE=${REG_FILE:-r500.reg}
REG_FORMAT=${REG_FORMAT:-ds9}
REG_FOV=${REG_FOV:-skyfov.fits}
SMOOTH_RADIUS=${SMOOTH_RADIUS:-4}


case "$1" in
    -[hH]*)
        echo "Usage:"
        echo "    `basename $0` <img_dir1> ..."
        exit 1
        ;;
esac


analyze_path() {
    # extract `obs id' and `source name' from path
    echo "$@" | awk '
    # main part
    {
        if (NF==1) {
            ## oi & name
            input=($1 "/")
            if (input ~ /_oi/) {
                ## PATTERN: .../$name_oi$oi/...
                idx_oi = match(input, /oi[0-9]+/) + 2;    # "2" skip the "oi"
                len_oi = RLENGTH - 2;
                oi = substr(input, idx_oi, len_oi);
                idx_name = match(input, /\/[a-zA-Z0-9.+-]+_oi/) + 1;
                len_name = RLENGTH - 4;
                name = substr(input, idx_name, len_name);
                owner = "lwt";
            }
            else {
                ## PATTERN: .../$name/$oi/...
                idx_oi = match(input, /\/[0-9]+\//) + 1;
                len_oi = RLENGTH - 2;
                oi = substr(input, idx_oi, len_oi);
                idx_name1 = match(input, /\/[a-zA-Z0-9_.+-]+\/[0-9]+\//);
                len_name1 = RLENGTH;
                name1 = substr(input, idx_name1, len_name1);
                idx_name = match(name1, /\/[a-zA-Z0-9_.+-]+\//) + 1;
                len_name = RLENGTH - 2;
                name = substr(name1, idx_name, len_name);
                owner = "zzh";
            }
            ## output
            printf("input: %s\n", input)
            printf("oi: %s\nname: %s\nowner: %s\n", oi, name, owner)
        }
        else {
            printf("*** WARNING: invalid input: %s\n", $0)
        }
    }
    # END { }
    '
}


INIT_DIR=`pwd -P`
while [ ! -z "$1" ]; do
    imgdir="$1"
    shift
    cd ${imgdir}
    echo "====== ${PWD} ======"
    NAME=`analyze_path "${PWD}" | grep '^name:' | awk '{ print $2 }'`
    FILE=`\ls ${FILE_PATTERN}`
    [ -z "${FILE}" ] && continue
    echo "FITS file: ${FILE}"
    #
    PNG_SUFFIX=`echo "${FILE%.fits}" | sed -e 's/_c7//' -e 's/_c0-3//'`
    PNG_FILE="${NAME}_${PNG_SUFFIX}.png"
    [ ! -f "${PNG_FILE}" ] && touch ${PNG_FILE}
    echo "PNG file: ${PNG_FILE}"
    # FoV region
    if [ -f "${REG_FOV}" ]; then
        DS9_REG_FOV="-regions color green -regions ${REG_FOV}"
    fi
    #
    ds9 ${FILE} \
        -bin factor 1 \
        -smooth radius ${SMOOTH_RADIUS} -smooth yes \
        -scale asinh \
        -cmap sls \
        -geometry 1288x1026-0+0 \
        -regions format ${REG_FORMAT} \
        -regions color white -regions ${REG_FILE} ${DS9_REG_FOV}
    #
    cd ${INIT_DIR}
done