Skip to content

Latest commit

 

History

History
91 lines (81 loc) · 3.09 KB

2021-09-29.md

File metadata and controls

91 lines (81 loc) · 3.09 KB

Plik wejściowy: a4_633M_3M0.raw

Próba konwersji do formatu .fil (filterbank) za pomocą bin2fil zakończyła się błędem, chociaż plik wyjściowy został utworzony:

Wartość parametru -f (center frequency) została dobrana eksperymentalnie. Przy podaniu wartości rzeczywistej (633MHz), bin2fil generował następujęce wpisy w nagłówku pliku .fil:

Central freq (MHz) = 637.5            
Low channel (MHz) = 636.06           
High channel (MHz) = 638.94  

Ostatecznie próba konwersji wyglądała tak:

$ time bin2fil -f 628.5 -sr 3000000 -o a4_633M_3M0.raw.fil a4_633M_3M0.raw 
conversion parameters:
conv_start     0
conv_len       7200
mjd_time       -1
filt_h         0.1
filt_l         60.0
ampli          3.0
channels       25
sample_rate    3000000
center_freq    628.5
source_name    B0329+54
source_ra      33259.37
source_de      543443.57
sr_error       3.0
input_file     a4_633M_3M0.raw
out_file       a4_633M_3M0.raw.fil
chunk_size     50000
Traceback (most recent call last):70%)  ( 99%)
  File "/home/thomson/radio/bin2fil-py/bin2fil.py", line 313, in <module>
    main()
  File "/home/thomson/radio/bin2fil-py/bin2fil.py", line 308, in main
    elaborate(chunk_size, debug=False)
  File "/home/thomson/radio/bin2fil-py/bin2fil.py", line 276, in elaborate
    d = d.reshape(-1, channels).transpose()
ValueError: cannot reshape array of size 1032112 into shape (25)

real	27m7,173s
user	7m40,673s
sys	0m24,698s

Inspekcja pliku za pomocą readfile z pakietu Presto wykazała:

$ readfile a4_633M_3M0.raw.fil 
Assuming the data is a SIGPROC filterbank file.


1: From the SIGPROC filterbank file 'a4_633M_3M0.raw.fil':
                  Telescope = Fake
                Source Name = B0329+54
                    Backend = Unknown
            Obs Date String = 2021-09-29T19:44:30.7096
             MJD start time = 59486.82257765707618
                   RA J2000 = 03:32:59.3700
             RA J2000 (deg) = 53.247375        
                  Dec J2000 = 54:34:43.5700
            Dec J2000 (deg) = 54.5787694444444 
                  Tracking? = True
              Azimuth (deg) = 7.410985e-323
           Zenith Ang (deg) = 6.947696e-310
            Number of polns = 2 (summed)
           Sample time (us) = 111111.444444444 
         Central freq (MHz) = 633              
          Low channel (MHz) = 631.56           
         High channel (MHz) = 634.44           
        Channel width (MHz) = 0.12             
         Number of channels = 25
      Total Bandwidth (MHz) = 3                
                       Beam = 1 of 1
            Beam FWHM (deg) = 1.000
         Spectra per subint = 2400
           Spectra per file = 1090850000
      Time per subint (sec) = 266.667466667
        Time per file (sec) = 121205919.172
            bits per sample = 8
          Are bytes signed? = False
          bytes per spectra = 25
        samples per spectra = 25
           bytes per subint = 60000
         samples per subint = 60000
                zero offset = 0                
           Invert the band? = False
       bytes in file header = 238