Skip to content

Latest commit

 

History

History
152 lines (121 loc) · 5.42 KB

2021-09-26.md

File metadata and controls

152 lines (121 loc) · 5.42 KB
Raw
  • Plik wejściowy: p4_626M_2M4.raw
  • Rozmiar pliku: 17206607872
  • Częstoliwość: 626MHz
  • Próbkowanie: 2.4MSPS

Próba interpretacji formatu

Nagrałem w GQRXie 10 sekund danych, co zajęło 206MB przy próbkowaniu 2.4MSPS, czyli ok. 20MB/sek. Najsensowniesze wydaje się założenie, że każda próbka I/Q zajmuje 8 bajtów, co przy 2.4MSPS daje długość 19.2MB/sek. Wynika z tego, że plik p4_626M_2M4.raw był nagrywany przez 896.177 sekund, czyli niecałe 15 minut.

Ok, z tego wynika kilka ważnych wniosków:

  • każda próbka IQ jest zapisana na 8 bajtach, po 4 bajty na I/Q
  • niestety, nie ma dokumentacji do GQRX, ale jest taki wątek: gqrx-sdr/gqrx#250 Wynika z niego, że to jest typ complex.

Poprawiona interpretacja formatu

  • Zakładając, że plik trwał godzinę: jedna próbka ma 2 bajty (RTL próbkuje po 8 bitów na każdy sygnał I oraz Q).

Eksperyment 1: konwersja za pomocą mmsdr2fil

Użyta komenda: ./mmsdr2fil ../data/2021-09-26/p4_626M_2M4.raw 626000000 2400000 2

Plik wynikowy: a.fil

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

1: From the SIGPROC filterbank file 'a.fil':
                  Telescope = Fake
                Source Name = B0329+54
                    Backend = Unknown
            Obs Date String = 2021-10-27T19:48:59.29
             MJD start time = 59514.82568622715917
                   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) = 4.940656e-324
           Zenith Ang (deg) = 6.949681e-310
            Number of polns = 2 (summed)
           Sample time (us) = 250              
         Central freq (MHz) = 627.2            
          Low channel (MHz) = 626.15           
         High channel (MHz) = 628.25           
        Channel width (MHz) = 0.3              
         Number of channels = 8
      Total Bandwidth (MHz) = 2.4              
                       Beam = 1 of 1
            Beam FWHM (deg) = 1.000
         Spectra per subint = 2400
           Spectra per file = 14338839
      Time per subint (sec) = 0.6
        Time per file (sec) = 3584.70975
            bits per sample = 8
          Are bytes signed? = False
          bytes per spectra = 8
        samples per spectra = 8
           bytes per subint = 19200
         samples per subint = 19200
                zero offset = 0                
           Invert the band? = False
       bytes in file header = 238

Wynik analizy rfifind:

$ rfifind -time 10 -o rtl10 a.fil 


               Pulsar Data RFI Finder
                 by Scott M. Ransom

Assuming the data are SIGPROC filterbank format...
Reading SIGPROC filterbank data from 1 file:
  'a.fil'

    Number of files = 1
       Num of polns = 2 (summed)
  Center freq (MHz) = 627.2
    Num of channels = 8
    Sample time (s) = 0.00025       
     Spectra/subint = 2400
   Total points (N) = 14338839
     Total time (s) = 3584.70975    
     Clipping sigma = 6.000
   Invert the band? = False
          Byteswap? = False
     Remove zeroDM? = False

File  Start Spec   Samples     Padding        Start MJD
----  ----------  ----------  ----------  --------------------
1              0    14338839           0  59514.82568622715917

Analyzing data sections of length 40800 points (10.2 sec).
  Prime factors are:  2 2 2 2 2 3 5 5 17 
  WARNING:  The largest prime factor is pretty big!  This will
            cause the FFTs to take a long time to compute.  I
            recommend choosing a different -time value.

Writing mask data  to 'rtl10_rfifind.mask'.
Writing  RFI data  to 'rtl10_rfifind.rfi'.
Writing statistics to 'rtl10_rfifind.stats'.

Massaging the data ...

Amount Complete = 100%
There are 1960 RFI instances.


Total number of intervals in the data:  2816

  Number of padded intervals:        8  ( 0.284%)
  Number of  good  intervals:        0  ( 0.000%)
  Number of  bad   intervals:     2808  (99.716%)

  Ten most significant birdies:
#  Sigma     Period(ms)      Freq(Hz)       Number 
----------------------------------------------------
1  30.92     10.0024(85)     99.976(85)     5384    
2  26.51     6.6672(54)      149.99(12)     4526    
3  23.93     0.727260(45)    1375.024(85)   41566   
4  22.73     20.002(17)      49.994(42)     2528    
5  21.47     5.0007(34)      199.97(14)     6906    
6  19.67     3.3336(17)      299.98(16)     4927    
7  16.14     2.2223(12)      449.98(24)     2355    
8  15.79     31.2404         32.0098        2       
9  15.47     2.00017(77)     499.96(19)     2625    
10 14.66     2.5003(14)      399.94(22)     2925    

  Ten most numerous birdies:
#  Number    Period(ms)      Freq(Hz)       Sigma 
----------------------------------------------------
1  47448     0.820967(46)    1218.075(68)   10.41   
2  41566     0.727260(45)    1375.024(85)   23.93   
3  36704     0.710774(21)    1406.916(41)   8.81    
4  8717      0.65108(15)     1535.90(35)    5.23    
5  8064      0.65267(23)     1532.18(55)    5.29    
6  6906      5.0007(34)      199.97(14)     21.47   
7  6152      0.843051(65)    1186.168(92)   5.32    
8  5384      10.0024(85)     99.976(85)     30.92   
9  4927      3.3336(17)      299.98(16)     19.67   
10 4526      6.6672(54)      149.99(12)     26.51   

Done.

Żaden interwał nie został zaraportowany jako good. Wszystkie są albo bad (99.99%) albo użyte jako padding.