PRESTO pulsar 出版图编辑剪切

折叠，仅仅显示轮廓
prepfold -n 64 -npart 128 -accelcand 361 -accelfile NGC6517_20221231_DM182.50_ACCEL_20.cand -justprofs -noxwin -nosearch NGC6517_20221231_DM182.50.dat

prepfold -n 64 -npart 128 -par J1801-0857K_5.par -justprofs -noxwin -nosearch -topo -o NGC6517K 20200119_NGC6517K_DM182.50.dat

添加文字
convert NGC6517_20221231_DM182.50_ACCEL_Cand_361.pfd.png -gravity northwest -pointsize 50 -annotate +225+160 'NGC 6517R' NGC6517R.pfd.png
convert NGC6517R_PSR_1801-0857R.pfd.png -gravity northwest -pointsize 50 -annotate +225+150 'NGC 6517R' -pointsize 30 -annotate +130+180 'p = 6.585 ms; DM = 180.3 pc cm^{-3}' NGC6517RR.pfd.png

剪切

convert NGC6517R.pfd.png -gravity west -chop 3%x0 -gravity east -chop 75%x0 -gravity north -chop 0x6% -gravity south -chop 0x4% NGC6517R.pfd.png

并排放
左右并排
convert NGC6517R1.pfd.png NGC6517R2.pfd.png NGC6517R.png +append NGC6517R.png
上下并排
convert NGC6517R1.pfd.png NGC6517R2.pfd.png -append NGC6517R.png

average profile of pulsars from many observations
ls *.pfd > tmp.txt
sum_profiles.py -d 182.654793 -f -n 64 -t 20200123_NGC6517A_fitsfold_PSR_1801-0857A.pfd.bestprof tmp.txt
mv sum_profiles.bestprof A_sum_profiles.bestprof

obtain the average profile of pulsar form multiful observations

sum_profiles.py -d 182.45 -f -n 64 -t filename.pfd.bestprof tmp.txt (*.pfd file list of all observations)

sum_profiles.py -d 182.45 -f -n 64 -t filename.pfd.bestprof tmp.txt (*.pfd file list of all observations)

ls *.pfd > tmp.txt
sum_profiles.py -d 182.654793 -f -n 64 -t 20200123_NGC6517A_fitsfold_PSR_1801-0857A.pfd.bestprof tmp.txt
mv sum_profiles.bestprof A_sum_profiles.bestprof

c=colors[len(tempo_results.freqbands)][ii])
c=colors[len(tempo_results.freqbands)][ii], ms=15,capsize=5,elinewidth=3)

def create_plot():
# Set up the plot
#fig = plt.figure(figsize=(11,8.5))
fig = plt.figure(figsize=(20,7))

plt.savefig(savefn, orientation='landscape', papertype='letter')
plt.savefig(savefn, orientation='portrait', papertype='letter')

convert -append A_profile_res.png B_profile_res.png C_profile_res.png D_profile_res.png E_profile_res.png F_profile_res.png G_profile_res.png H_profile_res.png I_profile_res.png A_I.png

convert -append K_profile_res.png L_profile_res.png M_profile_res.png N_profile_res.png O_profile_res.png P_profile_res.png Q_profile_res.png R_profile_res.png K_R.png

convert A_I.png K_R.png +append A_R.png

convert A_R.png -background white -fill black -gravity south -splice 0x180 -pointsize 120 -draw "text +2200,-4140 'MJD'" -gravity east -splice 160x0 -gravity east -pointsize 150 -draw "rotate -90 text -1200,-4670 'Residuals (µs)'" A_R_res_profile.png