Forum der Vereinigung der Sternfreunde

Liebe Sternfreunde,

in dieser Woche habe ich das VdS-Journal Nr. 72 erhalten. Für die Fachgruppe Astrophysik und Algorithmen ist darin der letzte Aufsatz zur kleinen Serie zu Kartenprojektionen enthalten, azimutale Karten. Der Aufsatz wist wie fast immer von einem kleinen Python-Programm begleitet. Das Programm selbst ist klein, aber die Koordinatenliste ist lang. Deshlabverweise ich auf die Internetseite der Fachgruppe und drucke das hier nicht ab:

http://fg-astrophysik.vdsastro.de/prg72.html

Ich hatte übrigens zur Mitgliederversammlung einen kleinen Vortrag über diese Projektionen gehalten, auch dazu gibt es einen kurzen Kommentar.

~

Ab der nächsten Ausgabe widme ich mich in einer kleinen Serie dem Gravitationsgesetz, genauer dem sog. Mehrkörperproblem. Damit lassen sich viele interessante Fragestellungen der Astronomie simulieren, mit den heutigen Computern kann das jeder von zu Hause aus.

~

Zur Würzburger Frühjahrstagung gibt es einen Beitrag der Fachgruppe: Schwarze Löcher: Wo Licht ist, ist auch starker Schatten. Es geht um die Nachrechnung des im vorigen Jahr erstmals abgebildeten sog. Schatten eines Schwarzen Lochs. . Ihr wißt doch: Unsere Computerprogramme dienen nicht zuerst den Ergebnissen, sondern den Erkenntnissen.

_________________
Uwe Pilz, Fachgruppen Kometen und Astrophysik/Algorithmen.
Oft benutzte Instrumente: Fujinon 16x70 FMT-SX-2, TMB Apo 105/650, Skywatcher Evostar 120/900 ED, Ninja Dobson 320/1440, Smartscope Dwarf-3

Ein Sternfreund frug mich nach einer Kartenprojektion, welche oft für alte Weltkarten benutzt wird.

Dies ist eine Abart der stereographischen Projektion, welche ich in Heft 72 mit erläutert habe. Der Zentralpunkt im Aufsatz waren der Nord- oder der Südpol. Für die Weltkarten wird aber der Äquator benutzt. Ich habe die Formeln so umgestellt, dass man eine äquatorständige Projektion berechnen kann:

# Hauptprogramm
z=len(k)
tracer(0,0) # Bildschirmaktualisierung aus
# Zentralmeridian einstellen (in Grad)
L0=-110
L0=L0*pi/180
for i in range(z):
    l= k[-i][0]-L0
    f= k[-i][1]
    # lambda wieder normalisieren
    if (l>pi): l=l-2*pi
    if (l<-pi): l=l+2*pi
    # transversale stereographische Projektion
    delta=acos(cos(f)*cos(l))
    alpha=acos(sin(f)/sin(delta))
    r=tan(delta/2)
    y=r*cos(alpha)
    x=r*sin(alpha)
    if (l<0): x=-x
    if abs(delta)<pi/2:
        plot(400*x,400*y)
# noch einen Kreis drumrum
for j in range(1000):
    plot(400*sin(j*2*pi/1000),400*cos(j*2*pi/1000))
update() # Bildschirm anzeigen  
done()

Das gesamte Programm mit der Koordinatenliste habe ich in unsere Fachgruppenseiten gelegt. Ich habe mit diesem Formelsatz ein Kartenpaar mit den Zentralmeridianen von -110° und 70° berechnet:

_________________
Uwe Pilz, Fachgruppen Kometen und Astrophysik/Algorithmen.
Oft benutzte Instrumente: Fujinon 16x70 FMT-SX-2, TMB Apo 105/650, Skywatcher Evostar 120/900 ED, Ninja Dobson 320/1440, Smartscope Dwarf-3