Mapping

下载示例数据

链接：https://pan.cstcloud.cn/s/u6u3iyfdRqo

$ ls RAW_data/*
RAW_data/KY:
M33_OTF_2021_07_31_05_14_00_000.xlsx

RAW_data/M33_OTF:
20210731

噪音管定标、RA-DEC

  • 方法一

 1#!/bin/bash
 2
 3# enable **/
 4shopt -s globstar
 5# 1. Temperature calibration
 6echo "Temperature calibration:"
 7# only need input chunk *0001.fits
 8files="$(ls RAW_data/M33_OTF/20210731/M33_OTF_1_*W_0001.fits)"
 9# check files
10printf "${files}"
11echo
12# commands
13# python -m hifast sep -p 3 同时执行5个 python -m hifast.sep
14python -m hifast sep -p 1 \
15    $files \
16   -d 4 -m 4 -n 596 --step 1 \
17   --frange 1350 1450 \
18   --smooth gaussian --s_sigma 2 \
19   --check_cal A --pcal_vary_lim_bin 0.02 \
20   --merge_pcals True --method_merge median --method_interp linear \
21   --save_pcals True \
22   --outdir 'output/%(project)s/%(date)s'
23
24# RA-DEC
25# # only need beam 01 specs_T.hdf5 file
26python -m hifast radec **/*-M01*-specs_T.hdf5

 1#!/bin/bash
 2
 3# enable **/
 4shopt -s globstar
 5# 1. Temperature calibration
 6echo "Temperature calibration:"
 7# only need input chunk *0001.fits
 8files="$(ls RAW_data/M33_OTF/20210731/M33_OTF_1_*W_0001.fits)"
 9# check files
10printf "${files}"
11echo
12# commands
13# python -m hifast sep -p 3 同时执行5个 python -m hifast.sep
14python -m hifast sep -p 3 \
15    $files \
16   -d 4 -m 4 -n 596 --step 1 \
17   --frange 1350 1450 \
18   --smooth gaussian --s_sigma 2 \
19   --check_cal A --pcal_vary_lim_bin 0.02 \
20   --merge_pcals False --cal_dis_lim 2 \
21   --save_pcals True \
22   --outdir 'output/%(project)s/%(date)s'
23
24# RA-DEC
25# # only need beam 01 specs_T.hdf5 file
26python -m hifast radec **/*-M01*-specs_T.hdf5

基线、驻波、RFI ...

  • 方法一

    用去基线的方法扣除系统温度

run-s2-v1.sh S2-sw.ini S2-rfi.ini

 1#!/bin/bash
 2
 3# enable **/
 4shopt -s globstar
 5
 6files="$(ls output/M33_OTF_1_MultiBeamOTF/20210731/M33_OTF_1_MultiBeamOTF-M*_W-*-specs_T.hdf5)"
 7# check fpatten
 8printf "${files}"
 9echo
10
11# using hifast.sh to run multi files
12# commands
13commands=$(cat <<EOF
14python -m hifast.flux  | 
15python -m hifast.bld  | --nproc 5 --frange 1400 1440 \
16                        --method PLS-asym2 --lam 1e9 \
17                        --s_method_freq gaussian --s_sigma_freq 3 \
18                        --s_method_t gaussian --s_sigma_t 11 \
19                        --exclude_type auto2
20# The `hifast.bld` process below can be replaced by using the `--post_method` parameter in the previous `hifast.bld` command
21python -m hifast.bld  | --nproc 5 \
22                        --method poly-asym2 --deg 2 \
23                        --s_method_freq gaussian --s_sigma_freq 1 \
24                        --exclude_type auto2
25python -m hifast.rfi  | -c conf/S2-rfi.ini
26python -m hifast.sw   | -c conf/S2-sw.ini
27python -m hifast.multi | --vtype optical --frame LSRK --merge_polar True --replace_rfi True
28EOF
29)
30# check commands
31echo "commands: ""$commands"
32# run multiple files
33hifast.sh "$files" -c "$commands" -n 5

  • 方法二

    减去Ref来扣除系统温度。S2-sw.ini 和 S2-rfi.ini 这两个参数文件与方法一的一致即可。

    run-s2-v1.sh

 1#!/bin/bash
 2
 3# enable **/
 4shopt -s globstar
 5
 6files="$(ls output/M33_OTF_1_MultiBeamOTF/20210731/M33_OTF_1_MultiBeamOTF-M*_W-*-specs_T.hdf5)"
 7# check fpatten
 8printf "${files}"
 9echo
10
11# using hifast.sh to run multi files
12# commands
13commands=$(cat <<EOF
14python -m hifast.flux  | 
15python -m hifast.bld  | --nproc 5 --frange 1400 1440 \
16                        --method PLS-asym2 --lam 1e9 \
17                        --s_method_freq gaussian --s_sigma_freq 3 \
18                        --njoin 50 \
19                        --exclude_type auto2
20# The `hifast.bld` process below can be replaced by using the `--post_method` parameter in the previous `hifast.bld` command
21python -m hifast.bld  | --nproc 5 \
22                        --method poly-asym2 --deg 2 \
23                        --s_method_freq gaussian --s_sigma_freq 1 \
24                        --exclude_type auto2
25python -m hifast.rfi  | -c conf/S2-rfi.ini
26# `--no_bld True` to keep baseline and only remove standing wave
27python -m hifast.sw   | -c conf/S2-sw.ini --no_bld True
28# sustract "ref" observation
29python -m hifast.ref | --method MedMed --nsection 11 --npart 8 \
30                       --post_method poly-asym2 --post_deg 3 \
31                       --post_s_method_freq gaussian --post_s_sigma_freq 1 \
32                       --post_exclude_type auto2
33python -m hifast.multi | --vtype optical --frame LSRK --merge_polar True --replace_rfi True
34EOF
35)
36# check commands
37echo "commands: ""$commands"
38# run multiple files
39hifast.sh "$files" -c "$commands" -n 5

Imaging

1#!/bin/bash
2
3python -m hifast.cube \
4     --outname cube.fits  \
5     -m bessel_gaussian  --type3 vopt \
6     --range3 -1000 2000 \
7     --share-mem False --step 38 --nproc 40 \
8     output/M33_OTF_1_MultiBeamOTF/*/*fc.hdf5