...
 
Commits (2)
......@@ -126,7 +126,7 @@ process in couple of minutes do have a strong signal. This indicates that the
software is looping faster though its loops.
See :download:`Timestamps wu</files/ON0EMEProcessTimeWU.txt>`
=======
In the below figure you can see the morse signal detected in wu 102 on -4,7 kHz.
In the below figure you can see the Morse signal detected in wu 102 on -4,7 kHz.
::
......@@ -143,7 +143,6 @@ After splitting the 2.5MHz raw data into 256 workunits, these workunits are
processed by seti_boinc or the setiathome_v8 program provided by Seti at Home.
Also the seti grafics is used to visualize the processed workunits.
Place the wu in any preferred directory and run seti_boinc in the same directory.
As soon this is running the seti grafics can be started.
......@@ -176,12 +175,15 @@ Recorded Broadband signal using new DT backend
Next to the main signal this figure depicts cross product artifacts, which
complicates the investigation on frequency mirroring. This signal is splited
using splitter_pfb into 256 workunits. All 256 workunits were processed by
SaH software. From each result.sah file the spikes were extracted. In below
figure you can see the same artifats as in previous figure.
using splitter_pfb into 256 workunits. All 256 workunits were processed by SaH
software. From each result.sah file the spikes were extracted. The
wus_data_extraction.py has been developed to run seti at home software and to
extract the signals from the result.sah file into a csv file which is plotted in
gnuplot, see data folder of this repository. In below figure you can see the
same artifacts as in previous figure.
.. figure:: /graphics/spikes.png
Extracted spikes from SaH result file.
Extracted spikes from SaH result file. Note: frequency dimension is incorrect
Workunit having more than 8 spikes are selected and plotted with gnuplot to
check the frequency as function of time. The overall idea is that the frequency
......@@ -208,11 +210,7 @@ Workunit 128
.. note::
Frequency axis is not well define. It is number divided by 10 in kHz
Frequency axis is not well defined. It is number/10 in kHz (delta)
Conlusion
---------
......
#!/home/marc/venvs/math/bin/python3.6
import wget
import subprocess
import xmltodict
import json
from shutil import copyfile
from xml.parsers.expat import ExpatError
url_wu = 'https://charon.camras.nl/public/setiatcamras/Sweep2/CAMRAS_DT_1573991026.24490.0.3.0.'
seti_at_home = 'setiathome_8.00_x86_64-pc-linux-gnu'
seti_graphics = 'setigraphics_8.00_x86_64-pc-linux-gnu'
for wu in range(1):
wu_file_name = wget.download(url_wu + str(wu+31))
copyfile(wu_file_name, 'work_unit.sah')
wu_process = subprocess.run([seti_at_home],check=True)
try:
with open('result.sah') as fd:
result_dict = dict(xmltodict.parse(fd.read()))
except ExpatError as e:
print(e)
print("Need to fix result.sah")