Sonnenaktivität 2015/16 via SID

Hi Ralph,

<blockquote id="quote"><font size="1" face="Verdana, Arial, Helvetica" id="quote">Zitat:<hr height="1" noshade id="quote">Der Frequenzbereich des Sender ist kHz Bereich und damit ausserhalb der Spezifikation dieses Radioastronomieforums aber ich hoffe das geht auch noch durch :).<hr height="1" noshade id="quote"></blockquote id="quote"></font id="quote">

ich seh das recht locker. Die anderen Forenteilnehmer hoffentlich auch. Du beobachtest die Sonne mit Amateurfunkmitteln, dann paßt das schon. Ich freu mich über jede lesbare Info auf deutsch. Denn fast alles ist in englisch und damit tue ich mir schwer...

Irgendwann möchte ich die Sonne auch so beobachten können. Daher finde ich es gut und richtig, daß wir heir jetzt verschiedene Bereiche (Meteorscattering, Sonnenbeobachtung, ...) aufziehen. Könntest Du ein bissy auf die eingesetzte Hardware eingehen?

Z.B. ein Bild Deiner ?selbstgebauten? Antenne posten? Ich kann mir das mit den Wicklungen grad net so recht vorstellen...

Dann interessiert mich auch, was so ein "Super SID Monitor" eigentlich ist, wo man den herbekommt (oder selber bauen kann), was der kostet und welchen konkreten USB A/D Wandler Du da verwendest...

Ich bin kein Amateurfunker und verstehe daher viele für Funker triviale Aussagen nicht.

Grüße
Silvia

Hallo Ralph,

danke für Deinen Beitrag.

Ich teile die Meinung von Silvia; die Erforschung der Sonne ist ein wichtiger Teil der Astronomie und die dazu verwendete Wellenlänge unterliegt keinen Einschränkungen.

Weitere technische Infos sind bestimmt für mehrere Forumsteilnehmer von Interesse.

Gruß
Tim

Hi Ralph,

<blockquote id="quote"><font size="1" face="Verdana, Arial, Helvetica" id="quote">Zitat:<hr height="1" noshade id="quote">So hier ein Bild meines technischen Aufbaus: Antenne <hr height="1" noshade id="quote"></blockquote id="quote"></font id="quote">

jetzt hats bei mir geschnackelt. Merci für das Bild. [:)]

<blockquote id="quote"><font size="1" face="Verdana, Arial, Helvetica" id="quote">Zitat:<hr height="1" noshade id="quote">(mit langem schwarzen Kabel)<hr height="1" noshade id="quote"></blockquote id="quote"></font id="quote">

ist das das Kabel, mit dem Du die Wicklungen um den Rahmen machst oder mußt Du mit diesem langen Kabel eine größere Entfernung zwischen Antenne und SID, A/D Wandler und Raspberry Pi übewinden?

<blockquote id="quote"><font size="1" face="Verdana, Arial, Helvetica" id="quote">Zitat:<hr height="1" noshade id="quote">Den SuperSID Monitor und den Draht für die Antenne habe ich von http://solar-center.stanford.edu/SID/sidmonitor für $ 125,-. Für astronomische Verhältnisse alles echt günstig ;)<hr height="1" noshade id="quote"></blockquote id="quote"></font id="quote">

[:)]

<blockquote id="quote"><font size="1" face="Verdana, Arial, Helvetica" id="quote">Zitat:<hr height="1" noshade id="quote">Ich werde noch meine Anleitung zum Aufsetzen des Raspberry Pi aus dem Englischen ins Deutsche übersetzen und einen Link dazu posten.<hr height="1" noshade id="quote"></blockquote id="quote"></font id="quote">

[:)]

Grüße
Silvia

Hi Ralph,

<blockquote id="quote"><font size="1" face="Verdana, Arial, Helvetica" id="quote">Zitat:<hr height="1" noshade id="quote">um den Holzrahmen ist ein dünner, gelb isolierter Draht gewickelt, der als Antenne dient - wenn du genau hinschaust siehst du wie ich ihn vom Rahmen wegführe. Das langen schwarze Koaxialkabel verbindet die Antenne (= dünner gelben Draht) mit dem SID Monitor. Die weitere Signalverarbeitung geht dann vom SID Monitor via A/D Wandler zum Raspberry Pi.<hr height="1" noshade id="quote"></blockquote id="quote"></font id="quote">

ok, jetzt hab ich den Aufbau komplett verstanden. Danke.

Grüße
Silvia

Die wichtigste Information zur Rahmenantenne fuer den Nichtfachmann fehlt leider: Drahtart und Windungszahl. Interessant waeren auch die Frequenzen der Sender GQD und ICV.
Darum sei hier nachgetragen: Der Drahtdurchmesser ist egal,
die Windungszahl abhaengig vom Eingangswiderstand der Elektronik.

Zum Thema VLF gibt es im Netz haufenweise Informationen, auch zum
Antennen- und Verstaerkerbau. $125 sind astronomisch billig, nicht aber funktechnisch, abgesehen vom umstaendlichen Import aus den USA. Alles was man braucht, kriegt man preiswert und uebernacht bei deutschen Elektronikhaendlern.

Die allerbeste Webseite zum Thema: http://www.vlf.it

Bedauerlicherweise ist der letzte Post schon 3 Jahre alt. Sind die Radioastronomen ausgestorben?
Darum diese Zeilen zur Ergaenzung und Belebung. Ich finde es aeusserst interessant, das Thema Sonne von einer ganz anderen Seite anzugehen und werde dazu noch einiges sagen. Gruss, eric2

Hallo Eric,

bist du auf dem Gebiet noch aktiv?

Hast du schon ein paar SIDs (Sudden Ionoshperic Disturbance) aufgezeichnet?

Für mich ist das Thema ziemlich spannend, aber neu. Den Sonnenaufgang auf Sizilien und in Gescher konnte ich gestern schon loggen.

Viele Grüße,

Wilhelm

Hallo Wilhelm,

eine 120-m-Antenne habe ich mir nun gebaut und frage mich, was genau und wie Du in dem Diagramm aufgezeichnet hast. Kann ich das mit HDSDR oder was brauche ich, evtl. auch an zusätzlich Hardware.

Mit den vorhandenen Anzeigen in HDSDR konnte ich keine zeitliche Funktion der Signalstärke (der Trägerwelle eines KW-Senders) erstellen.

Freue mich auf Tipps.

Hallo Erik,

das Projekt steht noch ziemlich am Anfang.

Zum Aufzeichnen der Langwellensender und zum Aufzeichnen von Meteorströmen verwende ich einen Nvidia Jetson Xavier NX. Eine preiswerte Variante, der Nvidia Jetson Nano oder ein Raspberry Pi sollten aber auch ausreichen.

Dieser Kleincomputer ist über einen HDMI zu USB3 Adapter mit dem HDMI Ausgang meines Notebooks verbunden. Der Jetson sieht dann den PC Bildschirm als USB3 Kamera an. Ein Python3 Programm liest die umrandeten Bereiche des Wasserfalls (s. Screenshot) aus und summiert die Werte auf.

Meine Python3 Programme zum loggen und plotten gebe ich gerne weiter. Die momentane Version des SID Aufnahmeprogramm ist angehängt. Es ist noch sehr roh…

Viele Grüße,

Wilhelm

# (C) WiSi-Testpilot 2021
import datetime
import numpy as np
import cv2
import time
import math as m

RED = (0, 0, 255)
YELLOW = (0, 255, 255)
light_BLUE = (200, 100, 100)
BLUE = (255, 0, 0)
GREEN = (0, 255, 0,)
MAGENTA = (255, 255, 0)
WHITE = (255, 255, 255)

bg = np.zeros((740, 1500, 3), np.uint8)             # Result Screen
bg = cv2.line(bg,(1,1),(1,700),(WHITE),2)
bg = cv2.line(bg,(1,699),(1500,699), WHITE,2)       # Skala

for j in range (1,7):
    bg = cv2.line(bg,(0, j * 100),(15, j * 100), WHITE,2)   # Y- Ticks

cap = cv2.VideoCapture(0)                           # /dev/video0, Spektrum Lab live Screen Input
cap.set(3,1280)                                     # HD  1280  720
cap.set(4,720)

font = cv2.FONT_HERSHEY_SIMPLEX
counter = 0    # Debug Info
altsek = 0

warda = 0
neunhundert = 110 # 100 Vorschub
#neunhundert = 5

yo = 100      #Positionen der Rahmen
yu = 300

xl = 900
xr = 1150

bxl = 200
bxr = 450

cxl = 520
cxr = 770  #770

teiler =  (xr - xl)   * (yu - yo)
bteiler = (bxr - bxl) * (yu - yo)
cteiler = (cxr - cxl) * (yu - yo)

print()
print('   Start: ',time.strftime("%T"),' UTC') 
print()
while cap.isOpened():                                              # diese Schleife wird 25 x in der Sekunde durchlaufen
    ret, streamout = cap.read()                                    # Einlesen des Frame 
    save = cv2.cvtColor(streamout, cv2.COLOR_BGR2GRAY)

    streamout = cv2.rectangle(streamout, (xl, yo), (xr, yu), GREEN, 2)         # plottet grünen Rahmen    
    streamout = cv2.rectangle(streamout, (bxl, yo), (bxr, yu), YELLOW, 2)      # plottet gelben Rahmen
    streamout = cv2.rectangle(streamout, (cxl, yo), (cxr, yu), RED, 2)         # plottet roten Rahmen
   
    # drei Summierer
    signal = 0     
    for i in range (xl, xr):
        for j in range (yo, yu):
            bw = save[j,i]
            signal = signal + bw
    fl = 700 - int (signal / teiler * 6)    
    # print('Ch1 ',time.strftime("%T"),' ',str(signal).zfill(8))

    bsignal = 0
    for i in range (bxl, bxr):
        for j in range (yo, yu):
            bw = save[j,i]
            bsignal = bsignal + bw
    bfl = 700 - int (bsignal / bteiler * 6)  
    # print('Ch2 ',time.strftime("%T"),' ',str(bsignal).zfill(8))   

    csignal = 0
    for i in range (cxl, cxr):
        for j in range (yo, yu):
            bw = save[j,i]
            csignal = csignal + bw
    cfl = 700 - int (csignal / cteiler * 6)  
    print(datetime.datetime.now(),' ',str(signal).zfill(9),' ',str(bsignal).zfill(9),' ',str(csignal).zfill(9))
            
    y_t_color = GREEN
    y_t_fuellung = -1            
    dicke = 1
    y_t = (int(counter/neunhundert), fl)   
    bg = cv2.circle(bg, y_t, dicke , y_t_color, y_t_fuellung)

    y_t_color = YELLOW
    y_t = (int(counter/neunhundert), bfl)   
    bg = cv2.circle(bg, y_t, dicke , y_t_color, y_t_fuellung) 

    y_t_color = RED
    y_t = (int(counter/neunhundert), cfl)   
    bg = cv2.circle(bg, y_t, dicke , y_t_color, y_t_fuellung) 

    Hora = int (time.strftime("%H"))
    Minute = int (time.strftime("%M"))
    Sekunde = int (time.strftime("%S"))
    fontScale = 0.8
    fthickness = 1
    Timestamp = str(Hora)
    Timestr = (Timestamp + 'UTC')

    org = (int(counter/neunhundert), 735)
    if (Minute == 0) and (Sekunde == 0) and (warda == 0):
        bg = cv2.putText(bg, Timestr, org, font, fontScale, WHITE, fthickness) 
        bg = cv2.line(bg,(int(counter/neunhundert),700),(int(counter/neunhundert), 715), WHITE, 2)  # X-Ticks
        print('1 TimeTi ',time.strftime("%T")) 
        warda = 50

    if ((Minute == 15) or (Minute == 30) or (Minute == 45)) and (Sekunde == 0) and (warda == 0):
        bg = cv2.line(bg,(int(counter/neunhundert),700),(int(counter/neunhundert), 710), WHITE, 2)  # X-Ticks
        print('1 TimeTi ',time.strftime("%T")) 
        warda = 50

    if (warda > 0):
        warda -= 1
       
    cv2.imshow('Output', bg)
    cv2.imshow("Auslesebildschirm", streamout)

    counter += 1   # zum messen der Framerate und Timing
  #  print("  Sekunde: ", time.strftime("%S"),"Count:",counter)
  #  print('   Start: ',time.strftime("%T"),' UTC', '  Signal: ', int (signal / teiler * 5)) 
    keyCode = cv2.waitKey(30) & 0xFF
    if keyCode == 27:
        cap.release()
        cv2.destroyAllWindows()
        break

Hallo Wilhelm,

so hatte ich mir das beinahe gedacht. Für mich wäre das erst einmal (noch) nicht der favorisierte Weg. Ich werde Deine interessante Lösung aber im Gedächtnis behalten. Vielen Dank für den Python-Code. Er wird zu gegebener Zeit sicherlich ein gute Starthilfe sein.

Guten Abend.

Gestern, am 7.5. 2021, gab es einen M3.9 Ausbruch, der auf GBZ (19.58 kHz) und GQD (22.1 kHz) sehr schön zu sehen ist.

https://sidstation.loudet.org/data-en.xhtml

https://www.swpc.noaa.gov/products/goes-x-ray-flux

Jetzt gerade, am 8.5. 18:40 UTC, findet wieder ein Ausbruch statt. Er ist C8.6 stark und auch bei beiden Frequenzen zu sehen.

Auch der C1.5 von 13 Uhr UTC ist noch gut zu erkennen.

Meine Anlage war leider nicht in Betrieb….

Viele Grüße,

Wilhelm

Hallo Erik,

meine Software habe ich nun auf Win10 und Postprocessing umgestellt. Im schwarzen Forum gibt es einen Thread dazu.

Auf den Bildern ist ein C4.8 Ausbruch von gestern, den 9.7.2021, um 17:25 UTC zu sehen, der auf 6 Langwellenfrequenzen eine Änderung bewirkte. In manchen Fällen gibt es einen Ausschlag nach unten, wie z.B. hier an der weißen Kurve zu sehen ist. Vielleicht bekommen wir gemeinsam heraus, was die Ursache dafür ist.

Ein schönes Wochenende wünsche ich,

viele Grüße,

Wilhelm

Hallo Ralph,

danke für deine Erklärung.

Ja, ich meinte diesen negativen Ausschlag in deiner oberen (DHO) Kurve etwas links von der Mitte. Es passiert wohl seltener, zumindest habe ich den Ausschlag nach unten noch nicht oft gesehen. In deinem Plot ist er ja sehr ausgeprägt.

Mein DHO38 Signalverlauf ist sehr flach im Vergleich zu deiner Kurve, s. die obere weiße Spur. Auch der C4.8 Ausbruch vom 9.7.2021 um 17:25 UTC ist nur schwach zu sehen. Hier im westlichen Münsterland empfange ich DHO offensichtlich direkt.

Viele Grüße,

Wilhelm

Hallo zusammen,

mich hat Caros Artikel "Suche nach langwelligen Gravitationswellensignalen" getriggert mal wieder zu schauen was mit SDR Empfängern sinnvoll möglich wäre.

Dies ist nun zwar nicht direkt für SDR, scheint mir ein gutes Projekt zu sein, da ich gerade aktiv mit der Sonnenbeobachtung in H-Alpha, und demnächst auch WL, beginne. Ich muss mir den Thread in den nächsten Tagen mal intensiv von oben nach unten durchlesen

Gibt es aus den vorhergehenden Experimenten weitere Ergebnisse, die ihr berichten könnt?

CS & VG

Stefan