Le Swiss Ephemeris portate in Python

Le Swiss Ephemeris, routine di calcolo astronomico di alta precisione, verranno utilizzate in questo progetto solo per alcune funzioni essenziali, non utilizzando i metodi di interpolazione delle effemeridi, e limitatamente al sistema di domificazione più diffuso (Placidus) di cui è nota l'impossibilità di utilizzo per latitudini polari a causa della degenerazione del metodo stesso. Una versione per Python in ambiente linux è gia disponibile sotto il nome di pyswisseph.
Le librerie pyswisseph sono molto facili da installare usando alcuni semplici comandi.
Li riassumo brevemente nel seguito, per utenti Ubuntu come il sottoscritto; credo sia altrettanto facile importarle in altre distribuzioni. Per Windows si possono ottenere dal curatore delle librerie stesse "for a small fee", oppure possono essere compilate dai sorgenti, più avanti faremo qualche tentativo.
In breve, per installare le pyswisseph bisogna:
caricare le python-developer  e pip con il comando :

sudo apt-get install python-dev python-pip

pip è un package realizzzato per installare i moduli presenti nel repository principale dei python developers: Pypi

quindi bisogna installare la libreria pyswisseph:

sudo pip install pyswisseph

a questo punto si possono usare le librerie nel vostro codice python semplicemente richiamandole con l'istruzione:

import swisseph 

oppure, in alternativa:

from swisseph import *

la seconda forma è, di regola, sconsigliata perchè rischia di contaminare il namespace principale con variabili e nomi di funzione che si possono sovrapporre ad omonimi già esistenti. Personalmente preferisco usare sempre il primo metodo, anche se occorre digitare ogni volta il nome del pacchetto principale, seguito da un punto, prima delle funzioni chiamate.

Di seguito dimostro, con poche righe di codice e molti commenti, come il cuore del nostro progetto si riduca a poche righe, naturalmente senza una gestione automatica della localizzazione, della grafica e della simbologia planetaria. Tuttavia si vede facilmente come una buona libreria risolva gran parte dei problemi legati alla stesura di codice. Il codice, così com'è, può essere copiato in un file astro.py ed eseguito direttamente tramite l'interprete python.

# -*- coding: utf-8 -*-
import swisseph
# input dati essenziali, ora espressa in UTC, tempo coordinato universale
# (sostanzialmente l'ora di Greenwich), le longitudini orientali sono positive,

# quelle ad Ovest di Greenwich sono negative
anno = 1990
mese = 7
giorno = 23
ora = 19
minuto = 10
latitudine = 45.00
longitudine = 11.00

# chiamata alla funzione swisseph._julday, che restituisce il giorno giuliano standard
# riferito alla data inputata (l'underscore prima del nome è intenzionale,

# identifica una funzione differente da swisseph.julday)
giorno_giuliano=swisseph._julday(anno, mese, giorno, ora, minuto)

# chiamata alla funzione houses che restituisce una tupla di tuple
# la longitudine di ogni casa si trova facilmente con gli indici:
# case[0][0], case [0][1] e così via fino alla dodicesima casa[0][11]
case=swisseph.houses(giorno_giuliano, latitudine, longitudine)

# chiamata alla funzione calc, che accetta come input il giorno giuliano
# precedentemente trovato e un indice che individua il pianeta:
# 0=sole, 1=luna, 2=mercurio, 3=venere, 4=marte, 5=giove, 6=saturno
# 7=urano, 8=nettuno, 9=plutone, 10=nodo lunare medio, 11= nodo lunare vero
# e una serie di altri corpi celesti che non consideriamo.
# Per facilità e rapidità di scrittura utilizzo una list comprehension per
# indicizzare i pianeti richiesti. 

long_planet = [swisseph.calc(giorno_giuliano, x) for x in xrange(0,12)]

print "il giorno giuliano è:", giorno_giuliano
for i in range(0,12):
    print swisseph.get_planet_name(i), long_planet[i][0]
for i in range(0,12):
    print "cuspide ",i+1,"° casa:",case[0][i]

L'output è il seguente:

il giorno giuliano è: 2448096.29861
Sun 120.668415422
Moon 142.547706141
Mercury 141.295122231
Venus 94.3882989781
Mars 37.4190095096
Jupiter 114.409427893
Saturn 291.353521287
Uranus 276.657709825
Neptune 282.705341008
Pluto 224.96999813
mean Node 307.670417261
true Node 307.261137012
cuspide  1 ° casa: 308.094747498
cuspide  2 ° casa: 359.633528903
cuspide  3 ° casa: 36.8373098779
cuspide  4 ° casa: 61.8615329935
cuspide  5 ° casa: 82.0092465108
cuspide  6 ° casa: 102.058362387
cuspide  7 ° casa: 128.094747498
cuspide  8 ° casa: 179.633528903
cuspide  9 ° casa: 216.837309878
cuspide  10 ° casa: 241.861532994
cuspide  11 ° casa: 262.009246511
cuspide  12 ° casa: 282.058362387

nel prossimo post mostrerò come incasellare i risultati, resi più comprensibili, in una piccola interfaccia grafica