Multiple plot utility

[Robert Benedict]

Here is a program to plot multiple raw files from different days and/or different nodes. It is an extensive mod of the WWV_plt2.py program by dkazdan & jgibbons. It can plot Doppler or Power for up to ten data files and adds an average curve. Source is at

https://github.com/rkayakr/PSWSmultiPlot

   73

     Bob  KD8CGH

[Robert Benedict]

Plots of one day with data from my Grape1 and that of David Waugh, KE8QEP, who lives nearby. My data is the first curve, his is the second.

  73

    Bob  KD8CGH

[Robert Benedict]

Version v1.1 updated to include plotting sunrise and sunset times for the location in the first rawdata file.

https://github.com/rkayakr/PSWSmultiPlot/releases/tag/v1.1

   Bob  KD8CGH

[William Smith]

I'm getting  a few errors running this on the RPi under ssh:

Unable to init server: Could not connect: Connection refused

Unable to init server: Could not connect: Connection refused

(multiPlot.py:25107): Gdk-CRITICAL **: 21:56:47.154: gdk_cursor_new_for_display: assertion 'GDK_IS_DISPLAY (display)' failed

(multiPlot.py:25107): Gdk-CRITICAL **: 21:56:47.163: gdk_cursor_new_for_display: assertion 'GDK_IS_DISPLAY (display)' failed

Though it does seem to work (with python3 and pip3 install suntime and some judicous editting for the Pi versus Windows),  so not 'real' errors.

Thanks for this!  Also, does my data look sane?

73, Willie N1JBJ

On Aug 5, 2021, at 3:46 PM, Robert Benedict <rka...@gmail.com> wrote:

Version v1.1 updated to include plotting sunrise and sunset times for the location in the first rawdata file.

https://github.com/rkayakr/PSWSmultiPlot/releases/tag/v1.1

   Bob  KD8CGH

<multi Power_N0000022_G1_EN91he_WWV10.png>

On Sunday, August 1, 2021 at 7:08:17 PM UTC-4 Robert Benedict wrote:

Plots of one day with data from my Grape1 and that of David Waugh, KE8QEP, who lives nearby. My data is the first curve, his is the second.

  73

    Bob  KD8CGH

On Sunday, August 1, 2021 at 5:46:02 PM UTC-4 Robert Benedict wrote:

Here is a program to plot multiple raw files from different days and/or different nodes. It is an extensive mod of the WWV_plt2.py program by dkazdan & jgibbons. It can plot Doppler or Power for up to ten data files and adds an average curve. Source is at

https://github.com/rkayakr/PSWSmultiPlot

   73

     Bob  KD8CGH

-- 
You received this message because you are subscribed to the Google Groups "hamsci-grape" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hamsci-grape...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/hamsci-grape/5b2565f3-84bc-4d64-ae69-dcb94eb43a5an%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
<multi Power_N0000022_G1_EN91he_WWV10.png>

[Robert Benedict]

Willie

   It looks like you are missing a lot of data. Are you running this on rawdata files? They should be about 4.9 MiB.

  Look at the output. You should see something like

      Ready to start processing records
      nf  0 len hours 84375 

  this is the length of the array of times for the first (zero index) file.

     Bob

[William Smith]

Hmm, I've really only got a couple of days of data so far:

-rw-r--r--  1 pi pi 845K Aug  3 23:59 analysis210803.csv

-rw-r--r--  1 pi pi 3.1M Aug  5 00:56 analysis210804.csv

-rw-r--r--  1 pi pi 4.8M Aug  5 23:59 analysis210805.csv

-rw-r--r--  1 pi pi 8.9K Aug  6 00:02 analysis210806.csv

The 3rd didn't start till amost 20:00, the 4th stopped recording just before 1500, and the 5th is my first full day.

Should I delete the partial days and just go with the full days, or is it worth leaving those in?

I got:

Ready to start processing records

nf  0 len hours 14158

Doppler min:  -0.92 ; Doppler max:  0.525

Vpk min:  0.000126 ; Vpk max:  0.10092

dB min:  -77.97 ; dB max:  -19.92

/usr/lib/python3/dist-packages/scipy/signal/_arraytools.py:45: FutureWarning: Using a non-tuple sequence for multidimensional indexing is deprecated; use `arr[tuple(seq)]` instead of `arr[seq]`. In the future this will be interpreted as an array index, `arr[np.array(seq)]`, which will result either in an error or a different result.

  b = a[a_slice]

2021-08-03

process file  1 Srawdata/analysis210804.csv

process file  2 Srawdata/analysis210805.csv

process file  3 Srawdata/analysis210806.csv

Final Plot for Decoded 10MHz WWV Beacon

Plot File: multi_Power_N0000000_G1_FN42kn_WWV10.png

At the very least it's showing the difference in signal 'power' between my 40M dipole and my YouLoop Magnetic Loop.  8*)

73, Willie N1JBJ

To view this discussion on the web visit https://groups.google.com/d/msgid/hamsci-grape/159a29f2-54c3-49d9-b5da-4ef37a040fb8n%40googlegroups.com.

[Robert Benedict]

Since you have partial data I think that the program is plotting the data you have correctly. Try collecting a few days of uninterrupted data and then take a look. It should look something like the plots posted in the first posts in the thread.

    73

      Bob

[Robert Benedict]

multiPlot version V1.2 cleans code up to make Pi and Windows use easier

    Bob  KD8CGH