Pulsar Voices

<p><strong>The data is sourced from CSIRO Parkes ATNF.</strong></p><p><strong>eg </strong>http://www.atnf.csiro.au/research/pulsar/psrcat/</p><p><strong>Feel the pulse of the universe</strong></p><p>We're taking signal data from astronomical "pulsar" sources and creating a way to listen to their signals audibly.</p><p>Pulsar data is available from <a href="http://www.atnf.csiro.au/outreach/education/everyone/pulsars/index.html">ATNF at CSIRO.au</a>. Our team at <a href="http://sciencehackau.github.io/melbourne/">#SciHackMelb</a> has been working on a #datavis to give researchers and others a novel way to explore the Pulsar corpus, especially through the sound of the frequencies at which the Pulsars emit pulses.</p><p>Link to project page at #SciHackMelb - http://www.the-hackfest.com/events/melbourne-science-hackfest/projects/pulsar-voices/</p><p>The files attached here include: source data, project presentation, data as used in website final_pulsar.sql, and other methodology documentation. Importantly, see the Github link which contains data manipulation code, html code to present the data, and render audibly, iPython Notebook to process single pulsar data into an audible waveform file. Together all these resources are the Pulsar Voices activity and resulting data.</p><p><strong>Source Data;</strong></p><p>* RA - east/west coordinates (0 - 24 hrs, roughly equates to longitude) [theta; transforms RA to 0 - 360*]<br>* Dec - north/south coordinates (-90, +90 roughly equates to latitude i.e. 90 is above north pole, and -90 south pole)<br>* P0 - the time in seconds that a pulsar repeats its signal<br>* f - 1/P0 which ranges from 700 cycles per sec, to some which pulses which occur every few seconds<br>* kps - distance from Earth in kilo-parsecs. 1 kps = 3,000 light years. The furthest data is 30 kps. The galactic centre is about 25,000 light years away i.e. about 8kps.</p><p>psrcatShort.csv = 2,295 Pulsars all known pulsars with above fields; RA, Dec, Theta</p><p>psrcatMedium.csv - add P0 and kps, only 1428 lines - i.e. not available for all 2,295 datapoint</p><p>psrcatSparse.csv - add P0 and kps, banks if n/a, 2,295 lines</p><p>short.txt - important pulsars with high levels of observation (** even more closely examined)</p><p>pulsar.R - code contributed by Ben Raymond to visualise Pulsar frequency, period in histogram</p><p><b>Added to the raw data:</b></p><p>- Coordinates to map RA, Dec to screen width(y)/height(x)</p><p>y = RA[Theta]*width/360; x = (Dec + 90)*height/180</p><p>- audible frequency converted from Pulsar frequency (1/P0)</p><p>Formula for 1/P0(x) -> Hz(y) => y = 10 ^ (0.5 log(x) + 2.8)</p><p>Explanation in text file; Convert1/P0toHz.txt</p><p>Tone generator from: http://www.softsynth.com/webaudio/tone.php</p><p>- detailed waveform file audible converted from Pulsar signal data, and waveform image (and python notebook to generate; available):</p><p>The project source is hosted on github at:<br></p><p><a href="https://github.com/gazzar/pulsarvoices">https://github.com/gazzar/pulsarvoices</a></p><div>An IPython/Jupyter notebook contains code and a rough description of the method used to process a psrfits .sf file<br>downloaded via the CSIRO Data Access Portal at <a href="http://doi.org/10.4225/08/55940087706E1" target="_blank">http://doi.org/10.4225/08/55940087706E1</a><br>The notebook contains experimental code to read one of these .sf files and access the contained spectrogram data, processing it to generate an audible signal.<br>It also reads the .txt files containing columnar pulse phase data (which is also contained in the .sf files) and processes these by frequency modulating the signal with an audible carrier.<br>This is the method used to generate the .wav and .png files used in the web interface.</div><ul><li>https://github.com/gazzar/pulsarvoices/blob/master/ipynb/hackfest1.ipynb </li></ul><div>A standalone python script that does the .txt to .png and .wav signal processing was used to process 15 more pulsar data examples. These can be reproduced by running the script.</div><ul><li>https://github.com/gazzar/pulsarvoices/blob/master/data/pulsarvoices.py</li></ul><p>Processed file at: https://github.com/gazzar/pulsarvoices/tree/master/web</p><p>https://github.com/gazzar/pulsarvoices/blob/master/web/J0437-4715.png</p><p><a href="https://github.com/gazzar/pulsarvoices/blob/master/web/J0437-4715.wav">J0437-4715.wav</a> | <a href="https://github.com/gazzar/pulsarvoices/blob/master/web/J0437-4715.png">J0437-4715.png</a></p><p>)</p><p>#Datavis online at: http://checkonline.com.au/tooltip.php. Code at Github linked above. See especially:</p><p>https://github.com/gazzar/pulsarvoices/blob/master/web/index.php</p><p>particularly, lines 314 - 328 (or search: "SELECT * FROM final_pulsar";) which loads pulsar data from DB and push to screen with Hz on mouseover.</p><p><b>Pulsar Voices webpage Functions</b>:<br>1.There is sound when you run the mouse across the Pulsars. We plot all known pulsars (N=2,295), and play a tone for pulsars we had data on frequency i.e. about 75%.<br>2. In the bottom left corner a more detailed Pulsar sound, and wave image pops up when you click the star icon. Two of the team worked exclusively on turning a single pulsars waveform into an audible wav file. They created 16 of these files, and a workflow, but the team only had time to load one waveform. With more time, it would be great to load these files.<br>3. If you leave the mouse over a Pulsar, a little data description pops up, with location (RA, Dec), distance (kilo parsecs; 1 = 3,000 light years), and frequency of rotation (and Hz converted to human hearing).<br>4.If you click on a Pulsar, other pulsars with similar frequency are highlighted in white. With more time I was interested to see if there are harmonics between pulsars. i.e. related frequencies.<br></p><p><b>The Team</b></p><p>Michael Walker is:  <a href="http://orcid.org/0000-0003-3086-6094" target="_blank">orcid.org/0000-0003-3086-6094</a> ; Biosciences PhD student, Unimelb, Melbourne.</p><p>Richard Ferrers is: orcid.org/<a href="http://orcid.org/0000-0002-2923-9889">0000-0002-2923-9889</a> ; ANDS Research Data Analyst, Innovation/Value Researcher, Melbourne.</p><p>Sarath Tomy is:  <a href="http://orcid.org/0000-0003-4301-0690" target="_blank">http://orcid.org/0000-0003-4301-0690</a> ; La Trobe PhD Comp Sci, Melbourne.</p><p>Gary Ruben is: <a href="http://orcid.org/0000-0002-6591-1820" target="_blank">http://orcid.org/0000-0002-6591-1820</a> ; CSIRO Postdoc at Australian Synchrotron, Melbourne.</p><p>Christopher Russell is: Data Manager, CSIRO, Sydney.</p><p>https://wiki.csiro.au/display/ASC/Chris+Russell</p><p>Anderson Murray is:  <a href="http://orcid.org/0000-0001-6986-9140" target="_blank">orcid.org/0000-0001-6986-9140</a>; Physics Honours, Monash, Melbourne.</p><p>Contact: richard.ferrers@ands.org.au for more information.</p><p><br></p><p><strong>What is still left to do?</strong></p><p>* load data, description, images fileset to figshare :: DOI ; DONE except DOI<br>* add overview images as option eg frequency bi-modal histogram<br>* colour code pulsars by distance; DONE<br>* add pulsar detail sound to Top three Observants; 16 pulsars processed but not loaded<br>* add tones to pulsars to indicate f; DONE<br>* add tooltips to show location, distance, frequency, name; DONE<br>* add title and description; DONE<br>* project data onto a planetarium dome with interaction to play pulsar frequencies.DONE see youtube video at https://youtu.be/F119gqOKJ1U<br>* zoom into parts of sky to get separation between close data points - see youtube; function in Google Earth #datavis of dataset. Link at youtube.</p><p>* set upper and lower tone boundaries, so tones aren't annoying</p><p>* colour code pulsars by frequency bins e.g. >100 Hz, 10 - 100, 1 - 10, <1 Hz</p>