-j 截一段数据,跳过哪一部分
-l 取多长的数据,单位s 秒
--drop 搜索最大限度,一般设置比真实信号脉冲宽度要大
--overlap 重叠率,10%;避免边缘信噪比低的被删除,重叠率低有些脉冲会搜不到
--zapthred 取2或1 差不多,IQR四分位距,统计量
--rfi kadaneF 4 4 加速,对数据没有影响!数字越大越快,但是太大可能会跳过些干扰。
-i 对于单波束,设置为1
--rfi zap 频带消干扰
--snrloss 滑动窗口大小决定损失率,默认,不用改
--maxcand 最大搜寻候选体数量,比如100,50
--minpts 聚类算法中一类最小数量
--baseline 消基线,不改,前均值,后幅度 (0.0;0.1)
--z (--rfi)-zdot 消掉0dm的干扰,效果可能不太好
--repeater 还没写好
后面为输出的各项设置
--cont 数据连起来,否则会停掉搜寻,需要加上
--ddplan 一般不需要加;
--td 2 时间分辨率乘2
guppi_59696_10845_J0534+2200_search_0005_0017.fits
transientx_fil --psrfits -f guppi_59696_10845_J0534+2200_search_0005_0017.fits --dms 40 --ndm 1000 --thre 10 --ddm 0.3 -t 4 -v -j 0 0 -l 1 --drop --overlap 0.1 --zapthre 2
判断是干扰还是真实信号
1.看脉冲宽度,周期 (脉冲宽度占周期10%)
2.最高频率,最低频率--色散量--来算时间延迟(斜线)
图上倾斜造成的延迟量与色散造成的延迟量如果相等,可能是假信号。--消干扰后,干扰的低频先到达,高频后到达。
ps:
1. 周期性的窄带干扰,很明显。FRB也不会很窄
2. 图像的背景色都是归一化后呈现的结果,并无什么区别
guppi_59696_06233_J0332+5434_search_test_0004_0001.fits
斜线一般都是干扰
!!!检查一下终端:陡峭--脉冲频率变化太快,阶跃怎么出现的,可用vela数据试一下
频率通道带宽多少???1Mhz左右
验证是否是真实信号,DM不能太小了
算延迟,检查是否是干扰
可以先消干扰(比如presto),然后用transientx搜寻!
防止信号别漏掉!
20221020 10_20_SGR J1935+2154
guppi_59872_33313_SGRJ1935+2154_search_0007_0007.fits
transientx_fil -t 3 -v -j 0 0 -l 2 --fd 2 --td 2 -s J1935 -o J1935 --maxw 0.1 --drop --overlap 0.1 --zapthre 2 --dms 300 --ddm 0.1 --ndm 1000 --rfi kadaneF 4 4 --thre 7 -i 1 --psrfits guppi_59872_33313_SGRJ1935+2154_search_0007_0027.fits
取得时间可能太小,导致的!
transientx_fil --psrfits -f guppi_59872_33313_SGRJ1935+2154_search_0007_0026.fits --dms 200 --ndm 1000 --thre 10 --ddm 0.3 -t 4 -v -j 0 0 -l 1 --drop --overlap 0.1 --zapthre 2
避免运算量太大,下采样(--td)取大一些,色散步长宽一些!!
主要看脉冲来判断!
找出搜索最优的图像--一个脉冲聚类算出两个值(一个高字节,一个低字节)-比较两个脉冲的mjd,判断是否是同一个,接近可能为同一个。
exemple:
### transientx_fil: search for single pulse
```
transientx_fil -v -o test -t 4 --zapthre 3.0 --fd 1 --overlap 0.1 --ddplan ddplan.txt --thre 7 --maxw 0.1 --snrloss 0.1 -l 2.0 --drop -z kadaneF 8 4 zdot -f *.fil
```
### replot_fil: filter out the bad candidates
```
replot_fil -v -t 4 --zapthre 3.0 --td 1 --fd 1 --dmcutoff 3 --widthcutoff 0.1 --snrcutoff 7 --snrloss 0.1 --zap --zdot --kadane 8 4 7 --candfile test.cands --clean -f *.fil
```
testing:
transientx_fil -t 4 -v -j 0 0 -l 2 --fd 2 --td 2 -s J1935+2154 -o J1935+2154 --maxw 0.1 --drop --overlap 0.1 --zapthre 2 --dms 300 --ddm 1 --ndm 30 --rfi kadaneF 4 4 --thre 7 -i 1 --rfi zap 2174 2190 -f J1935+2154_59873.417742448495119_6.fil
transientx_fil -t 4 -v -j 0 0 -l 2 --fd 2 --td 2 -s J1935+2154 -o J1935+2154 --maxw 0.1 --drop --overlap 0.1 --zapthre 2 --dms 320 --ddm 1 --ndm 20 --rfi kadaneF 4 4 --thre 7 -i 1 --rfi zap 2174 2190 -f J1935+2154_59903.475594585237559_17_n17.fil