rds demo

import%20marimo%0A%0A__generated_with%20%3D%20%220.25.1%22%0Aapp%20%3D%20marimo.App()%0A%0A%0A%40app.cell%0Adef%20_()%3A%0A%20%20%20%20import%20os%0A%20%20%20%20import%20urllib.request%0A%20%20%20%20import%20matplotlib.pyplot%20as%20plt%0A%20%20%20%20import%20numpy%20as%20np%0A%20%20%20%20%23from%20gnuradio.filter%20import%20firdes%0A%20%20%20%20from%20scipy.signal%20import%20resample_poly%2C%20firwin%0A%20%20%20%20from%20matplotlib.animation%20import%20FuncAnimation%0A%0A%20%20%20%20%23recording%20%3D%20np.fromfile('%2Fhome%2Fmarc%2FDownloads%2Ffm_clip_for_rds.iq'%2C%20dtype%3Dnp.complex64)%20%23%20med%20SNR%0A%20%20%20%20%23recording%20%3D%20np.fromfile('%2Fhome%2Fmarc%2FDownloads%2Ffm_rds_250k.iq'%2C%20dtype%3Dnp.complex64)%20%23%20high%20SNR%0A%20%20%20%20%23%20high%20SNR%2C%20shorter.%20%20If%20it's%20not%20in%20the%20current%20directory%20it%20gets%20downloaded%20from%20GitHub%2C%20which%20is%20also%20how%20the%20web-based%20marimo%20notebook%20gets%20it%0A%20%20%20%20filename%20%3D%20'fm_rds_250k_1Msamples.iq'%0A%20%20%20%20if%20os.path.exists(filename)%3A%0A%20%20%20%20%20%20%20%20recording%20%3D%20np.fromfile(filename%2C%20dtype%3Dnp.complex64)%0A%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20recording%20%3D%20np.frombuffer(urllib.request.urlopen('https%3A%2F%2Fraw.githubusercontent.com%2F777arc%2F498x%2Fmaster%2F'%20%2B%20filename).read()%2C%20dtype%3Dnp.complex64)%0A%0A%20%20%20%20%23%20MAKE%20MY%20OWN%20RECORDING%20OF%20A%20COOLER%20STATION%2C%20JUST%20MAKE%20SURE%20ITS%20HIGH%20SNR%20AND%20SCALE%20IT%20TO%20THE%20SAME%20SIGNAL%20LEVEL%20AS%20THIS%20WORKING%20ONE%0A%20%20%20%20return%20FuncAnimation%2C%20firwin%2C%20np%2C%20plt%2C%20recording%2C%20resample_poly%0A%0A%0A%40app.cell%0Adef%20_(FuncAnimation%2C%20firwin%2C%20freq_log%2C%20np%2C%20plt%2C%20recording%2C%20resample_poly)%3A%0A%20%20%20%20%23%20Everything%20from%20the%20FM%20demod%20through%20the%20BPSK%20demod%20is%20one%20cell%2C%20because%20x%20gets%20reused%20for%20each%20step%20of%20the%20signal%20chain%0A%20%20%20%20sample_rate%20%3D%20250e3%0A%20%20%20%20center_freq%20%3D%2099.5e6%0A%0A%20%20%20%20%23%20Plot%20time%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20plt.plot(recording)%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%20%20%20%20%23%20PSD%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20PSD%20%3D%2010*np.log10(np.abs(np.fft.fftshift(np.fft.fft(recording)))**2)%0A%20%20%20%20%20%20%20%20PSD%20%3D%20PSD%5B%3A%3A100%5D%0A%20%20%20%20%20%20%20%20PSD%20%3D%20PSD%20-%20np.max(PSD)%0A%20%20%20%20%20%20%20%20f%20%3D%20np.linspace(sample_rate%2F-2%2C%20sample_rate%2F2%2C%20len(PSD))%2F1e3%0A%20%20%20%20%20%20%20%20plt.plot(f%2C%20PSD)%0A%20%20%20%20%20%20%20%20plt.axis(%5B-125%2C%20125%2C%20-55%2C%201%5D)%0A%20%20%20%20%20%20%20%20plt.xlabel(%22Frequency%20%5BkHz%5D%22)%0A%20%20%20%20%20%20%20%20plt.ylabel(%22PSD%20Before%20FM%20Demod%20%5BdB%5D%22)%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%0A%20%20%20%20%23%20Quadrature%20Demod%0A%20%20%20%20x%20%3D%200.5%20*%20np.angle(recording%5B0%3A-1%5D%20*%20np.conj(recording%5B1%3A%5D))%20%23%20see%20https%3A%2F%2Fwiki.gnuradio.org%2Findex.php%2FQuadrature_Demod%0A%0A%20%20%20%20%23%20PSD%0A%20%20%20%20if%20True%3A%0A%20%20%20%20%20%20%20%20print(len(x))%0A%20%20%20%20%20%20%20%20PSD%20%3D%2010*np.log10(np.abs(np.fft.fftshift(np.fft.fft(x)))**2%20%2F%20len(x)%20%2F%20sample_rate)%0A%20%20%20%20%20%20%20%20PSD%20%3D%20PSD%5B%3A%3A10%5D%0A%20%20%20%20%20%20%20%20PSD%20%3D%20PSD%5Blen(PSD)%2F%2F2%3A%5D%0A%20%20%20%20%20%20%20%20PSD%20%3D%20PSD%20-%20np.max(PSD)%0A%20%20%20%20%20%20%20%20f%20%3D%20np.linspace(0%2C%20sample_rate%2F2%2C%20len(PSD))%2F1e3%0A%20%20%20%20%20%20%20%20plt.plot(f%2C%20PSD)%0A%20%20%20%20%20%20%20%20plt.axis(%5B0%2C%2080%2C%20-70%2C%200%5D)%0A%20%20%20%20%20%20%20%20plt.xlabel(%22Frequency%20%5BkHz%5D%22)%0A%20%20%20%20%20%20%20%20plt.ylabel(%22PSD%20After%20FM%20Demod%20%5BdB%5D%22)%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%20%20%20%20%23%20Spectrogram%20(once%20i%20get%20a%20higher%20SNR%20better%20looking%20recording%20I%20can%20include%20spectrogram%20towards%20the%20start%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20fft_size%20%3D%201024%0A%20%20%20%20%20%20%20%20num_rows%20%3D%20int(np.floor(len(x)%2Ffft_size))%0A%20%20%20%20%20%20%20%20spectrogram%20%3D%20np.zeros((num_rows%2C%20fft_size))%0A%20%20%20%20%20%20%20%20for%20i%20in%20range(num_rows)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20spectrogram%5Bi%2C%3A%5D%20%3D%20np.log10(np.abs(np.fft.fftshift(np.fft.fft(x%5Bi*fft_size%3A(i%2B1)*fft_size%5D)))**2)%0A%0A%20%20%20%20%20%20%20%20plt.imshow(spectrogram%2C%20aspect%3D'auto'%2C%20extent%20%3D%20%5B-sample_rate%2F2%2F1e6%2C%20sample_rate%2F2%2F1e6%2C%200%2C%20len(x)%2Fsample_rate%5D)%0A%20%20%20%20%20%20%20%20plt.xlabel(%22Frequency%20%5BMHz%5D%22)%0A%20%20%20%20%20%20%20%20plt.ylabel(%22Time%20%5Bs%5D%22)%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%0A%20%20%20%20%23%20Freq%20shift%0A%20%20%20%20N%20%3D%20len(x)%0A%20%20%20%20f_o%20%3D%20-57e3%20%23%20amount%20we%20need%20to%20shift%20by%2C%20because%20it's%20a%20real%20signal%20being%20fed%20it%2C%20in%20doesn't%20matter%20if%20this%20is%20-%20or%20%2B%0A%20%20%20%20t%20%3D%20np.arange(N)%2Fsample_rate%0A%20%20%20%20x%20%3D%20x%20*%20np.exp(2j*np.pi*f_o*t)%20%23%20down%20shift%0A%0A%20%20%20%20%23%20Low%20Pass%20Filter%0A%20%20%20%20%23taps%20%3D%20firdes.low_pass(1.0%2C%20sample_rate%2C%207.5e3%2C%205e3)%20%23%20same%20as%20original%20flowgraph%0A%20%20%20%20taps%20%3D%20firwin(numtaps%3D101%2C%20cutoff%3D7.5e3%2C%20fs%3Dsample_rate)%0A%20%20%20%20x%20%3D%20np.convolve(x%2C%20taps%2C%20'valid')%0A%0A%20%20%20%20%23%20Decimate%20by%2010%2C%20now%20that%20we%20filtered%20and%20there%20wont%20be%20aliasing%0A%20%20%20%20x%20%3D%20x%5B%3A%3A10%5D%0A%20%20%20%20sample_rate%20%3D%2025e3%0A%0A%20%20%20%20%23%20At%20this%20point%20we're%20at%2025kHz%20sample%20rate%2C%20due%20to%20the%20decimate%20by%2010%0A%0A%20%20%20%20%23%20Resample%20to%2019kHz%0A%20%20%20%20x%20%3D%20resample_poly(x%2C%2019%2C%2025)%20%0A%20%20%20%20sample_rate%20%3D%2019e3%0A%0A%20%20%20%20%23%20Filters%20time%2Ffreq%20response%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20plt.plot(taps%2C'.-')%0A%20%20%20%20%20%20%20%20plt.show()%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20print(taps)%0A%20%20%20%20%20%20%20%20fft%20%3D%20np.abs(np.fft.fftshift(np.fft.fft(taps)))**2%0A%20%20%20%20%20%20%20%20fft%20%3D%20fft%5Blen(fft)%2F%2F2%3A%5D%0A%20%20%20%20%20%20%20%20f%20%3D%20np.linspace(0%2C%20sample_rate%2F2%2C%20len(fft))%2F1e3%0A%20%20%20%20%20%20%20%20plt.plot(f%2Cfft)%0A%20%20%20%20%20%20%20%20plt.grid()%0A%20%20%20%20%20%20%20%20plt.xlabel('Frequency%20%5BkHz%5D')%0A%20%20%20%20%20%20%20%20plt.ylabel('Frequency%20Response%20of%20Filter')%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%0A%0A%20%20%20%20%23%20TODO%20use%20https%3A%2F%2Fgithub.com%2Fwavewalkerdsp%2FblogDownloads%2Fblob%2Fmain%2FsrrcDesign.py%20or%20another%20to%20make%20RRC%20filter%20and%20freq%20shift%20it%20like%20in%20filters%20chapter%0A%0A%0A%20%20%20%20%23%20Normalize%20signal%2C%20to%20emulate%20AGC%2C%20this%20isn't%20actually%20required%20for%20any%20future%20step%2C%20they%20all%20scale%20to%20amplitude%20just%20fine%0A%20%20%20%20%23x%20%2F%3D%20np.mean(np.abs(x))%20%23%20take%20magnitude%2C%20find%20mean%2C%20normalize%20to%20that%0A%0A%20%20%20%20%23%20Plot%20constellation%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20plt.plot(np.real(x)%2C%20np.imag(x)%2C%20'.')%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%20%20%20%20%23%20For%20now%20just%20read%20in%20samples%20saved%20after%20the%20AGC%20that%20def%20worked%20in%20GR%0A%20%20%20%20%23x%20%3D%20np.fromfile('%2Fhome%2Fmarc%2FDownloads%2Ftest.iq'%2C%20dtype%3Dnp.complex64)%0A%0A%20%20%20%20%23%20Next%20part%20is%20symbol%20sync%2C%20using%20what%20we%20did%20in%20sync%20chapter%0A%20%20%20%20samples%20%3D%20x%20%23%20for%20the%20sake%20of%20matching%20the%20sync%20chapter%0A%20%20%20%20samples_interpolated%20%3D%20resample_poly(samples%2C%2032%2C%201)%20%23%20we'll%20use%2032%20as%20the%20interpolation%20factor%2C%20arbitrarily%20chosen%0A%20%20%20%20sps%20%3D%2016%0A%20%20%20%20mu%20%3D%200.01%20%23%20initial%20estimate%20of%20phase%20of%20sample%0A%20%20%20%20out%20%3D%20np.zeros(len(samples)%20%2B%2010%2C%20dtype%3Dnp.complex64)%0A%20%20%20%20out_rail%20%3D%20np.zeros(len(samples)%20%2B%2010%2C%20dtype%3Dnp.complex64)%20%23%20stores%20values%2C%20each%20iteration%20we%20need%20the%20previous%202%20values%20plus%20current%20value%0A%20%20%20%20i_in%20%3D%200%20%23%20input%20samples%20index%0A%20%20%20%20i_out%20%3D%202%20%23%20output%20index%20(let%20first%20two%20outputs%20be%200)%0A%20%20%20%20while%20i_out%20%3C%20len(samples)%20and%20i_in%2B32%20%3C%20len(samples)%3A%0A%20%20%20%20%20%20%20%20out%5Bi_out%5D%20%3D%20samples_interpolated%5Bi_in*32%20%2B%20int(mu*32)%5D%20%23%20grab%20what%20we%20think%20is%20the%20%22best%22%20sample%0A%20%20%20%20%20%20%20%20out_rail%5Bi_out%5D%20%3D%20int(np.real(out%5Bi_out%5D)%20%3E%200)%20%2B%201j*int(np.imag(out%5Bi_out%5D)%20%3E%200)%0A%20%20%20%20%20%20%20%20x%20%3D%20(out_rail%5Bi_out%5D%20-%20out_rail%5Bi_out-2%5D)%20*%20np.conj(out%5Bi_out-1%5D)%0A%20%20%20%20%20%20%20%20y%20%3D%20(out%5Bi_out%5D%20-%20out%5Bi_out-2%5D)%20*%20np.conj(out_rail%5Bi_out-1%5D)%0A%20%20%20%20%20%20%20%20mm_val%20%3D%20np.real(y%20-%20x)%0A%20%20%20%20%20%20%20%20mu%20%2B%3D%20sps%20%2B%200.01*mm_val%0A%20%20%20%20%20%20%20%20i_in%20%2B%3D%20int(np.floor(mu))%20%23%20round%20down%20to%20nearest%20int%20since%20we%20are%20using%20it%20as%20an%20index%0A%20%20%20%20%20%20%20%20mu%20%3D%20mu%20-%20np.floor(mu)%20%23%20remove%20the%20integer%20part%20of%20mu%0A%20%20%20%20%20%20%20%20i_out%20%2B%3D%201%20%23%20increment%20output%20index%0A%20%20%20%20x%20%3D%20out%5B2%3Ai_out%5D%20%23%20remove%20the%20first%20two%2C%20and%20anything%20after%20i_out%20(that%20was%20never%20filled%20out)%0A%0A%20%20%20%20%23%20Plot%20constellation%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20plt.plot(np.real(x%5B1000%3A1200%5D)%2C%20np.imag(x%5B1000%3A1200%5D)%2C%20'.')%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%20%20%20%20%23%20Animate%20constellation%0A%20%20%20%20if%20False%3A%20%0A%20%20%20%20%20%20%20%20subset%20%3D%20x%5B0%3A2000%5D%0A%20%20%20%20%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots()%0A%20%20%20%20%20%20%20%20fig.set_tight_layout(True)%0A%20%20%20%20%20%20%20%20line%2C%20%3D%20ax.plot(%5B0%2C%200%2C%200%2C%200%2C%200%5D%2C%20%5B0%2C%200%2C%200%2C%200%2C%200%5D%2C%20'.')%0A%20%20%20%20%20%20%20%20ax.axis(%5B-0.02%2C%200.02%2C%20-0.02%2C%200.02%5D)%0A%20%20%20%20%20%20%20%20subset%20%3D%20np.concatenate((np.zeros(100)%2C%20subset))%20%23%20Add%20zeros%20at%20the%20beginning%20so%20that%20when%20gif%20loops%20it%20has%20a%20transition%20period%0A%20%20%20%20%20%20%20%20def%20update(i)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20i%20%3D%20int(i)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(i)%0A%20%20%20%20%20%20%20%20%20%20%20%20line.set_xdata(%5Bnp.real(subset%5Bi*10%3A(i%2B2)*10%5D)%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20line.set_ydata(%5Bnp.imag(subset%5Bi*10%3A(i%2B2)*10%5D)%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20line%2C%20ax%0A%20%20%20%20%20%20%20%20anim%20%3D%20FuncAnimation(fig%2C%20update%2C%20frames%3Dnp.arange(0%2C%20len(subset)%2F10%20-%202)%2C%20interval%3D20)%0A%20%20%20%20%20%20%20%20anim.save('%2Ftmp%2Fconstellation-animated.gif'%2C%20dpi%3D80%2C%20writer%3D'imagemagick')%0A%0A%0A%0A%20%20%20%20%23%20Plot%20time%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20plt.plot(np.real(x)%2C%20'.-')%0A%20%20%20%20%20%20%20%20plt.plot(np.imag(x)%2C%20'.-')%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%0A%20%20%20%20%23%20Fine%20freq%20sync%0A%20%20%20%20samples%20%3D%20x%0A%20%20%20%20N%20%3D%20len(samples)%0A%20%20%20%20phase%20%3D%200%0A%20%20%20%20freq%20%3D%200%0A%20%20%20%20%23%20These%20next%20two%20params%20is%20what%20to%20adjust%2C%20to%20make%20the%20feedback%20loop%20faster%20or%20slower%20(which%20impacts%20stability)%0A%20%20%20%20alpha%20%3D%20100.0%0A%20%20%20%20beta%20%3D%200.5%0A%20%20%20%20out%20%3D%20np.zeros(N%2C%20dtype%3Dnp.complex64)%0A%20%20%20%20%23freq_log%20%3D%20%5B%5D%0A%20%20%20%20for%20i%20in%20range(N)%3A%0A%20%20%20%20%20%20%20%20out%5Bi%5D%20%3D%20samples%5Bi%5D%20*%20np.exp(-1j*phase)%20%23%20adjust%20the%20input%20sample%20by%20the%20inverse%20of%20the%20estimated%20phase%20offset%0A%20%20%20%20%20%20%20%20error%20%3D%20np.real(out%5Bi%5D)%20*%20np.imag(out%5Bi%5D)%20%23%20This%20is%20the%20error%20formula%20for%202nd%20order%20Costas%20Loop%20(e.g.%20for%20BPSK)%0A%0A%20%20%20%20%20%20%20%20%23%20Advance%20the%20loop%20(recalc%20phase%20and%20freq%20offset)%0A%20%20%20%20%20%20%20%20freq%20%2B%3D%20(beta%20*%20error)%0A%20%20%20%20%20%20%20%20%23freq_log.append(freq%20*%20sample_rate%20%2F%20(2*np.pi))%20%23%20convert%20from%20angular%20velocity%20to%20Hz%20for%20logging%0A%20%20%20%20%20%20%20%20phase%20%2B%3D%20freq%20%2B%20(alpha%20*%20error)%0A%0A%20%20%20%20%20%20%20%20%23%20Optional%3A%20Adjust%20phase%20so%20its%20always%20between%200%20and%202pi%2C%20recall%20that%20phase%20wraps%20around%20every%202pi%0A%20%20%20%20%20%20%20%20while%20phase%20%3E%3D%202*np.pi%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20phase%20-%3D%202*np.pi%0A%20%20%20%20%20%20%20%20while%20phase%20%3C%200%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20phase%20%2B%3D%202*np.pi%0A%20%20%20%20x%20%3D%20out%0A%0A%0A%20%20%20%20%23%20Animate%20constellation%0A%20%20%20%20if%20False%3A%20%0A%20%20%20%20%20%20%20%20subset%20%3D%20x%5B0%3A2000%5D%0A%20%20%20%20%20%20%20%20fig%2C%20ax%20%3D%20plt.subplots()%0A%20%20%20%20%20%20%20%20fig.set_tight_layout(True)%0A%20%20%20%20%20%20%20%20line%2C%20%3D%20ax.plot(%5B0%2C%200%2C%200%2C%200%2C%200%5D%2C%20%5B0%2C%200%2C%200%2C%200%2C%200%5D%2C%20'.')%0A%20%20%20%20%20%20%20%20ax.axis(%5B-0.02%2C%200.02%2C%20-0.02%2C%200.02%5D)%0A%20%20%20%20%20%20%20%20subset%20%3D%20np.concatenate((np.zeros(100)%2C%20subset))%20%23%20Add%20zeros%20at%20the%20beginning%20so%20that%20when%20gif%20loops%20it%20has%20a%20transition%20period%0A%20%20%20%20%20%20%20%20def%20update(i)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20i%20%3D%20int(i)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(i)%0A%20%20%20%20%20%20%20%20%20%20%20%20line.set_xdata(%5Bnp.real(subset%5Bi*10%3A(i%2B2)*10%5D)%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20line.set_ydata(%5Bnp.imag(subset%5Bi*10%3A(i%2B2)*10%5D)%5D)%0A%20%20%20%20%20%20%20%20%20%20%20%20return%20line%2C%20ax%0A%20%20%20%20%20%20%20%20anim%20%3D%20FuncAnimation(fig%2C%20update%2C%20frames%3Dnp.arange(0%2C%20len(subset)%2F10%20-%202)%2C%20interval%3D20)%0A%20%20%20%20%20%20%20%20anim.save('%2Ftmp%2Fconstellation-animated-postcostas.gif'%2C%20dpi%3D80%2C%20writer%3D'imagemagick')%0A%0A%0A%20%20%20%20%23%20Plot%20freq%20over%20time%20to%20see%20how%20long%20it%20takes%20to%20hit%20the%20right%20offset%2C%20this%20is%20what%20we%20look%20at%20to%20tweak%20alpha%20and%20beta%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20plt.plot(freq_log%5B%3A%3A100%5D%2C'.-')%0A%20%20%20%20%20%20%20%20plt.xlabel('Sample')%0A%20%20%20%20%20%20%20%20plt.ylabel('Frequency%20%5BHz%5D')%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%20%20%20%20%23%20Plot%20time%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20plt.plot(np.real(x)%2C%20'.-')%0A%20%20%20%20%20%20%20%20plt.plot(np.imag(x)%2C%20'.-')%0A%20%20%20%20%20%20%20%20plt.show()%0A%0A%20%20%20%20%23%20Demod%20BPSK%0A%20%20%20%20bits%20%3D%20(np.real(x)%20%3E%200).astype(int)%20%23%201's%20and%200's%0A%0A%20%20%20%20%23%20Differential%20decoding%2C%20so%20that%20it%20doesn't%20matter%20whether%20our%20BPSK%20was%20180%20degrees%20rotated%20without%20us%20realizing%20it%0A%20%20%20%20print(bits%5B0%3A10%5D)%0A%20%20%20%20bits%20%3D%20(bits%5B1%3A%5D%20-%20bits%5B0%3A-1%5D)%20%25%202%0A%20%20%20%20bits%20%3D%20bits.astype(np.uint8)%20%23%20for%20decoder%0A%20%20%20%20print(bits%5B0%3A10%5D)%0A%0A%0A%0A%20%20%20%20if%20False%3A%0A%20%20%20%20%20%20%20%20%23%20FOR%20TESTING%20USE%20KNOWN%20TO%20WORK%20BYTES%0A%20%20%20%20%20%20%20%20working_bits%20%3D%20np.fromfile('%2Fhome%2Fmarc%2FDownloads%2Fworking_bytes.uint8'%2C%20dtype%3Dnp.uint8)%0A%0A%20%20%20%20%20%20%20%20plt.figure(0)%0A%20%20%20%20%20%20%20%20plt.plot(bits)%0A%20%20%20%20%20%20%20%20plt.figure(1)%0A%20%20%20%20%20%20%20%20plt.plot(working_bits)%0A%20%20%20%20%20%20%20%20plt.show()%0A%20%20%20%20%20%20%20%20exit()%0A%20%20%20%20return%20(bits%2C)%0A%0A%0A%40app.cell%0Adef%20_(bits%2C%20np)%3A%0A%20%20%20%20%23%23%23%23%23%23%23%23%23%23%23%23%23%23%23%0A%20%20%20%20%23%20RDS%20Decoder%20%23%0A%20%20%20%20%23%23%23%23%23%23%23%23%23%23%23%23%23%23%23%0A%0A%20%20%20%20%23%20Writing%20notes%0A%20%20%20%20%23%20Remember%20to%20give%20credit%20to%20the%20GR%20block%0A%20%20%20%20%23%201187.5%20bps%20%2F%20104%20bits%20%3D%2011.4%20groups%2Fsec%2C%20or%2045.7%20blocks%2Fsec%20*%2F%0A%20%20%20%20%23%20Point%20out%20how%20the%20bad%20blocks%20towards%20the%20beginning%20are%20probably%20due%20to%20the%20time%20it%20takes%20to%20sync%0A%0A%20%20%20%20%23%20Constants%0A%20%20%20%20syndrome%20%3D%20%5B383%2C%2014%2C%20303%2C%20663%2C%20748%5D%0A%20%20%20%20offset_pos%20%3D%20%5B0%2C%201%2C%202%2C%203%2C%202%5D%0A%20%20%20%20offset_word%20%3D%20%5B252%2C%20408%2C%20360%2C%20436%2C%20848%5D%0A%0A%20%20%20%20%23%20see%20Annex%20B%2C%20page%2064%20of%20the%20standard%0A%20%20%20%20def%20calc_syndrome(x%2C%20mlen)%3A%0A%20%20%20%20%20%20%20%20reg%20%3D%200%0A%20%20%20%20%20%20%20%20plen%20%3D%2010%0A%20%20%20%20%20%20%20%20for%20ii%20in%20range(mlen%2C%200%2C%20-1)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20reg%20%3D%20(reg%20%3C%3C%201)%20%7C%20((x%20%3E%3E%20(ii-1))%20%26%200x01)%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20(reg%20%26%20(1%20%3C%3C%20plen))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20reg%20%3D%20reg%20%5E%200x5B9%0A%20%20%20%20%20%20%20%20for%20ii%20in%20range(plen%2C%200%2C%20-1)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20reg%20%3D%20reg%20%3C%3C%201%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20(reg%20%26%20(1%20%3C%3C%20plen))%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20reg%20%3D%20reg%20%5E%200x5B9%0A%20%20%20%20%20%20%20%20return%20reg%20%26%20((1%20%3C%3C%20plen)%20-%201)%20%23%20select%20the%20bottom%20plen%20bits%20of%20reg%0A%0A%20%20%20%20%23%20Initialize%20all%20the%20working%20vars%20we'll%20need%20during%20the%20loop%0A%20%20%20%20synced%20%3D%20False%0A%20%20%20%20presync%20%3D%20False%0A%0A%20%20%20%20wrong_blocks_counter%20%3D%200%0A%20%20%20%20blocks_counter%20%3D%200%0A%20%20%20%20group_good_blocks_counter%20%3D%200%0A%0A%20%20%20%20reg%20%3D%20np.uint32(0)%20%23%20was%20unsigned%20long%20in%20C%2B%2B%20(64%20bits)%20but%20numpy%20doesnt%20support%20bitwise%20ops%20of%20uint64%2C%20I%20don't%20think%20it%20gets%20that%20high%20anyway%0A%20%20%20%20lastseen_offset_counter%20%3D%200%0A%20%20%20%20lastseen_offset%20%3D%200%0A%0A%20%20%20%20%23%20the%20synchronization%20process%20is%20described%20in%20Annex%20C%2C%20page%2066%20of%20the%20standard%20*%2F%0A%20%20%20%20bytes_out%20%3D%20%5B%5D%0A%20%20%20%20for%20_i%20in%20range(len(bits))%3A%0A%20%20%20%20%20%20%20%20%23%20in%20C%2B%2B%20reg%20doesn't%20get%20init%20so%20it%20will%20be%20random%20at%20first%2C%20for%20ours%20its%200s%0A%20%20%20%20%20%20%20%20%23%20It%20was%20also%20an%20unsigned%20long%20but%20never%20seemed%20to%20get%20anywhere%20near%20the%20max%20value%0A%20%20%20%20%20%20%20%20%23%20bits%20are%20either%200%20or%201%0A%20%20%20%20%20%20%20%20reg%20%3D%20np.bitwise_or(np.left_shift(reg%2C%201)%2C%20bits%5B_i%5D)%20%23%20reg%20contains%20the%20last%2026%20rds%20bits.%20these%20are%20both%20bitwise%20ops%0A%20%20%20%20%20%20%20%20if%20not%20synced%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20reg_syndrome%20%3D%20calc_syndrome(reg%2C%2026)%0A%20%20%20%20%20%20%20%20%20%20%20%20for%20j%20in%20range(5)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20reg_syndrome%20%3D%3D%20syndrome%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20not%20presync%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lastseen_offset%20%3D%20j%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20lastseen_offset_counter%20%3D%20_i%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20presync%20%3D%20True%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20offset_pos%5Blastseen_offset%5D%20%3E%3D%20offset_pos%5Bj%5D%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_distance%20%3D%20offset_pos%5Bj%5D%20%2B%204%20-%20offset_pos%5Blastseen_offset%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_distance%20%3D%20offset_pos%5Bj%5D%20-%20offset_pos%5Blastseen_offset%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20(block_distance*26)%20!%3D%20(_i%20-%20lastseen_offset_counter)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20presync%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print('Sync%20State%20Detected')%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20wrong_blocks_counter%20%3D%200%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20blocks_counter%20%3D%200%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_bit_counter%20%3D%200%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_number%20%3D%20(j%20%2B%201)%20%25%204%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20group_assembly_started%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20synced%20%3D%20True%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20break%20%23%20syndrome%20found%2C%20no%20more%20cycles%0A%0A%20%20%20%20%20%20%20%20else%3A%20%23%20SYNCED%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20wait%20until%2026%20bits%20enter%20the%20buffer%20*%2F%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20block_bit_counter%20%3C%2025%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_bit_counter%20%2B%3D%201%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20good_block%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20dataword%20%3D%20(reg%20%3E%3E%2010)%20%26%200xffff%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_calculated_crc%20%3D%20calc_syndrome(dataword%2C%2016)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20checkword%20%3D%20reg%20%26%200x3ff%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20block_number%20%3D%3D%202%3A%20%23%20manage%20special%20case%20of%20C%20or%20C'%20offset%20word%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_received_crc%20%3D%20checkword%20%5E%20offset_word%5Bblock_number%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20(block_received_crc%20%3D%3D%20block_calculated_crc)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20good_block%20%3D%20True%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_received_crc%20%3D%20checkword%20%5E%20offset_word%5B4%5D%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20(block_received_crc%20%3D%3D%20block_calculated_crc)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20good_block%20%3D%20True%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20wrong_blocks_counter%20%2B%3D%201%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20good_block%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_received_crc%20%3D%20checkword%20%5E%20offset_word%5Bblock_number%5D%20%23%20bitwise%20xor%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20block_received_crc%20%3D%3D%20block_calculated_crc%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20good_block%20%3D%20True%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20wrong_blocks_counter%20%2B%3D%201%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20good_block%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20Done%20checking%20CRC%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20block_number%20%3D%3D%200%20and%20good_block%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20group_assembly_started%20%3D%20True%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20group_good_blocks_counter%20%3D%201%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_group%20%3D%20bytearray(8)%20%23%208%20bytes%20filled%20with%200s%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20group_assembly_started%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20not%20good_block%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20group_assembly_started%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20raw%20data%20bytes%2C%20as%20received%20from%20RDS.%208%20info%20bytes%2C%20followed%20by%204%20RDS%20offset%20chars%3A%20ABCD%2FABcD%2FEEEE%20(in%20US)%20which%20we%20leave%20out%20here%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20RDS%20information%20words%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23%20block_number%20is%20either%200%2C1%2C2%2C3%20so%20this%20is%20how%20we%20fill%20out%20the%208%20bytes%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_group%5Bblock_number*2%5D%20%3D%20(dataword%20%3E%3E%208)%20%26%20255%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20_group%5Bblock_number*2%2B1%5D%20%3D%20dataword%20%26%20255%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20group_good_blocks_counter%20%2B%3D%201%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23print('group_good_blocks_counter%3A'%2C%20group_good_blocks_counter)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20group_good_blocks_counter%20%3D%3D%205%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%23print(_group)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20bytes_out.append(_group)%20%23%20list%20of%20len-8%20lists%20of%20bytes%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_bit_counter%20%3D%200%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20block_number%20%3D%20(block_number%20%2B%201)%20%25%204%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20blocks_counter%20%2B%3D%201%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20blocks_counter%20%3D%3D%2050%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20if%20wrong_blocks_counter%20%3E%2035%3A%20%23%20This%20many%20wrong%20blocks%20must%20mean%20we%20lost%20sync%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(%22Lost%20Sync%20(Got%20%22%2C%20wrong_blocks_counter%2C%20%22%20bad%20blocks%20on%20%22%2C%20blocks_counter%2C%20%22%20total)%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20synced%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20presync%20%3D%20False%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20print(%22Still%20Sync-ed%20(Got%20%22%2C%20wrong_blocks_counter%2C%20%22%20bad%20blocks%20on%20%22%2C%20blocks_counter%2C%20%22%20total)%22)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20blocks_counter%20%3D%200%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20wrong_blocks_counter%20%3D%200%0A%20%20%20%20return%20(bytes_out%2C)%0A%0A%0A%40app.cell%0Adef%20_(bytes_out)%3A%0A%20%20%20%20%23%23%23%23%23%23%23%23%23%23%23%23%23%23%23%0A%20%20%20%20%23%20RDS%20Parser%20%20%23%0A%20%20%20%20%23%23%23%23%23%23%23%23%23%23%23%23%23%23%23%0A%0A%20%20%20%20%23%20Annex%20F%20of%20RBDS%20Standard%20Table%20F.1%20(North%20America)%20and%20Table%20F.2%20(Europe)%0A%20%20%20%20%23%20%20%20%20%20%20%20%20%20%20%20%20%20%20Europe%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20North%20America%0A%20%20%20%20pty_table%20%3D%20%5B%5B%22Undefined%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%22Undefined%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22News%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22News%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Current%20Affairs%22%2C%20%20%20%20%20%20%20%22Information%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Information%22%2C%20%20%20%20%20%20%20%20%20%20%20%22Sports%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Sport%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Talk%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Education%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%22Rock%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Drama%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Classic%20Rock%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Culture%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Adult%20Hits%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Science%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Soft%20Rock%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Varied%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Top%2040%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Pop%20Music%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%22Country%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Rock%20Music%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22Oldies%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Easy%20Listening%22%2C%20%20%20%20%20%20%20%20%22Soft%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Light%20Classical%22%2C%20%20%20%20%20%20%20%22Nostalgia%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Serious%20Classical%22%2C%20%20%20%20%20%22Jazz%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Other%20Music%22%2C%20%20%20%20%20%20%20%20%20%20%20%22Classical%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Weather%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Rhythm%20%26%20Blues%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Finance%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Soft%20Rhythm%20%26%20Blues%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Children%E2%80%99s%20Programmes%22%2C%20%22Language%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Social%20Affairs%22%2C%20%20%20%20%20%20%20%20%22Religious%20Music%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Religion%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Religious%20Talk%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Phone-In%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Personality%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Travel%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Public%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Leisure%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22College%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Jazz%20Music%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22Spanish%20Talk%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Country%20Music%22%2C%20%20%20%20%20%20%20%20%20%22Spanish%20Music%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22National%20Music%22%2C%20%20%20%20%20%20%20%20%22Hip%20Hop%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Oldies%20Music%22%2C%20%20%20%20%20%20%20%20%20%20%22Unassigned%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Folk%20Music%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22Unassigned%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Documentary%22%2C%20%20%20%20%20%20%20%20%20%20%20%22Weather%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Alarm%20Test%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%22Emergency%20Test%22%5D%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%5B%22Alarm%22%2C%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Emergency%22%5D%5D%0A%20%20%20%20pty_locale%20%3D%201%20%23%20set%20to%200%20for%20europe%20which%20will%20use%20first%20column%20instead%0A%0A%20%20%20%20%23%20page%2072%2C%20Annex%20D%2C%20table%20D.2%20in%20the%20standard%0A%20%20%20%20coverage_area_codes%20%3D%20%5B%22Local%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22International%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22National%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Supra-regional%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%201%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%202%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%203%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%204%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%205%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%206%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%207%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%208%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%209%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%2010%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%2011%22%2C%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%22Regional%2012%22%5D%0A%0A%20%20%20%20radiotext_AB_flag%20%3D%200%0A%20%20%20%20radiotext%20%3D%20%5B'%20'%5D*65%0A%20%20%20%20first_time%20%3D%20True%0A%20%20%20%20for%20group%20in%20bytes_out%3A%0A%20%20%20%20%20%20%20%20group_0%20%3D%20group%5B1%5D%20%7C%20(group%5B0%5D%20%3C%3C%208)%0A%20%20%20%20%20%20%20%20group_1%20%3D%20group%5B3%5D%20%7C%20(group%5B2%5D%20%3C%3C%208)%0A%20%20%20%20%20%20%20%20group_2%20%3D%20group%5B5%5D%20%7C%20(group%5B4%5D%20%3C%3C%208)%0A%20%20%20%20%20%20%20%20group_3%20%3D%20group%5B7%5D%20%7C%20(group%5B6%5D%20%3C%3C%208)%0A%20%20%0A%20%20%20%20%20%20%20%20group_type%20%3D%20(group_1%20%3E%3E%2012)%20%26%200xf%20%23%20here%20is%20what%20each%20one%20means%2C%20e.g.%20RT%20is%20radiotext%20which%20is%20the%20only%20one%20we%20decode%20here%3A%20%5B%22BASIC%22%2C%20%22PIN%2FSL%22%2C%20%22RT%22%2C%20%22AID%22%2C%20%22CT%22%2C%20%22TDC%22%2C%20%22IH%22%2C%20%22RP%22%2C%20%22TMC%22%2C%20%22EWS%22%2C%20%22___%22%2C%20%22___%22%2C%20%22___%22%2C%20%22___%22%2C%20%22EON%22%2C%20%22___%22%5D%0A%20%20%20%20%20%20%20%20AB%20%3D%20(group_1%20%3E%3E%2011%20)%20%26%200x1%20%23%20b%20if%201%2C%20a%20if%200%0A%0A%20%20%20%20%20%20%20%20%23print(%22group_type%3A%22%2C%20group_type)%20%23%20this%20is%20essentially%20message%20type%2C%20i%20only%20see%20type%200%20and%202%20in%20my%20recording%0A%20%20%20%20%20%20%20%20%23print(%22AB%3A%22%2C%20AB)%0A%0A%20%20%20%20%20%20%20%20program_identification%20%3D%20group_0%20%20%20%20%20%23%20%22PI%22%0A%0A%20%20%20%20%20%20%20%20program_type%20%3D%20(group_1%20%3E%3E%205)%20%26%200x1f%20%23%20%22PTY%22%0A%20%20%20%20%20%20%20%20pty%20%3D%20pty_table%5Bprogram_type%5D%5Bpty_locale%5D%0A%0A%20%20%20%20%20%20%20%20pi_area_coverage%20%3D%20(program_identification%20%3E%3E%208)%20%26%200xf%0A%20%20%20%20%20%20%20%20coverage_area%20%3D%20coverage_area_codes%5Bpi_area_coverage%5D%0A%0A%20%20%20%20%20%20%20%20pi_program_reference_number%20%3D%20program_identification%20%26%200xff%20%23%20just%20an%20int%0A%0A%20%20%20%20%20%20%20%20if%20first_time%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22PTY%3A%22%2C%20pty)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22program%3A%22%2C%20pi_program_reference_number)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(%22coverage_area%3A%22%2C%20coverage_area)%0A%20%20%20%20%20%20%20%20%20%20%20%20first_time%20%3D%20False%0A%0A%20%20%20%20%20%20%20%20if%20group_type%20%3D%3D%202%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%23%20when%20the%20A%2FB%20flag%20is%20toggled%2C%20flush%20your%20current%20radiotext%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20radiotext_AB_flag%20!%3D%20((group_1%20%3E%3E%204)%20%26%200x01)%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20radiotext%20%3D%20%5B'%20'%5D*65%0A%20%20%20%20%20%20%20%20%20%20%20%20radiotext_AB_flag%20%3D%20(group_1%20%3E%3E%204)%20%26%200x01%0A%20%20%20%20%20%20%20%20%20%20%20%20text_segment_address_code%20%3D%20group_1%20%26%200x0f%0A%20%20%20%20%20%20%20%20%20%20%20%20if%20AB%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20radiotext%5Btext_segment_address_code%20*%202%20%20%20%20%5D%20%3D%20chr((group_3%20%3E%3E%208)%20%26%200xff)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20radiotext%5Btext_segment_address_code%20*%202%20%2B%201%5D%20%3D%20chr(group_3%20%20%20%20%20%20%20%20%26%200xff)%0A%20%20%20%20%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20radiotext%5Btext_segment_address_code%20*4%20%20%20%20%20%5D%20%3D%20chr((group_2%20%3E%3E%208)%20%26%200xff)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20radiotext%5Btext_segment_address_code%20*%204%20%2B%201%5D%20%3D%20chr(group_2%20%20%20%20%20%20%20%20%26%200xff)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20radiotext%5Btext_segment_address_code%20*%204%20%2B%202%5D%20%3D%20chr((group_3%20%3E%3E%208)%20%26%200xff)%0A%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20%20radiotext%5Btext_segment_address_code%20*%204%20%2B%203%5D%20%3D%20chr(group_3%20%20%20%20%20%20%20%20%26%200xff)%0A%20%20%20%20%20%20%20%20%20%20%20%20print(''.join(radiotext))%0A%20%20%20%20%20%20%20%20else%3A%0A%20%20%20%20%20%20%20%20%20%20%20%20pass%0A%20%20%20%20%20%20%20%20%20%20%20%20%23print(%22unsupported%20group_type%3A%22%2C%20group_type)%0A%0A%0A%0A%0A%0A%0A%0A%20%20%20%20%23%20Once%20its%20all%20working%2C%20try%20decoding%20the%20pos%20and%20neg%20portions%20of%20RDS%20and%20compare%20BER%0A%20%20%20%20return%0A%0A%0Aif%20__name__%20%3D%3D%20%22__main__%22%3A%0A%20%20%20%20app.run()%0A