Extraire le mois d'une date au format juliandate

Version imprimable

01/12/2017, 09h02
habasque

Extraire le mois d'une date au format juliandate
Bonjour,

Débutant en python, je ne parviens pas à extraire le mois (01, 02...,12) d'une date au format juliandate.
Je récupère la date d'un fichier texte puis je la convertis:
Code:

1 2 #convert dates given into true Julian days, add 2451544.5 date_first_juld = date_first + 2451544.5
Si je fais un print du champ 'date_first_juld' voici le contenu:

Citation:

[ 2450659.34 2450670.311 2450681.351 ..., 2455192.144 2455173.01]

A partir de là j'ai fait des tentatives avec 'datetime.datetime.strptime', etc. mais sans succès.

Si quelqu'un a une idée, je suis preneur.

Merci d'avance,
Cdt
01/12/2017, 10h53
Sve@r

Bonjour

Citation:

Envoyé par habasque

Si je fais un print du champ 'date_first_juld'...

Surtout ne nous donne pas le contenu de "date_first", ça pimente un peu le challenge !!!

PS: le type aussi serait aussi une bonne idée mais je suis sûr que tu y aurais pensé tout seul...
01/12/2017, 12h16
habasque

Voici le contenu de date_first:

Citation:

[ -885.16 -874.189 -863.149 ..., 3647.644 3628.51 3638.193]

Le type : <type 'numpy.ndarray'>
01/12/2017, 14h10
Sve@r

Citation:

Envoyé par habasque

Le type : <type 'numpy.ndarray'>

Mouais, ça m'aide encore moins (mais bon, là ce n'est pas ta faute...)

Citation:

Envoyé par habasque

Je récupère la date d'un fichier texte puis je la convertis

Tu voudrais bien nous montrer sous quel format est écrit cette date dans le fichier texte... parce que partir d'une date (que j'imagine assez naïvement jj/mm/aaaa) et finir avec une suite de nombres "aléatoires?" dans un "ndarray" de numpy...:weird:

Accessoirement tu devrais aller voir cette discussion...

date_first, c'est le contenu du fichier que tu lis ?
Si tu le capture comme ça sous la forme d'un array, pour moi il contient plusieurs dates. Chaque nombre en représentant une.

Si je traite avec un seul nombre, il n'y a pas trop de problème :
Code:

1 2 3 4 5 6 7 from datetime import datetime date_first = -885.16 date_first_juld = date_first + 2451544.5 print( datetime.fromtimestamp(date_first_juld) ) print( datetime.fromtimestamp(date_first_juld).month ) ### affiche juste le mois (un entier) print( datetime.fromtimestamp(date_first_juld).strftime('%m') ) ### ou comme ca si tu le prefère dans une chaine de caractère
Et du coup si date_first est un array, et bien tu itères dessus, et tu fais l'opération sur chacun de ses éléments.

01/12/2017, 15h17
habasque

Merci on approche de la solution ! En revanche, peu importe le contenu du champ 'date_first', la date de sortie est:
1970-01-29
Seule l'heure change...

Pour Sve@r:

Les dates sont au format indiqué plut haut:
[ -885.16 -874.189 -863.149 ..., 3647.644 3628.51 3638.193]
La date de référence est January 1 2000 at 12h UTC soit 2451545 en nombre jour julien.
C'est pour cela que l'on a des valeurs négatives.
01/12/2017, 16h01
wiztricks

Salut,

Citation:

Envoyé par habasque

Merci on approche de la solution ! En revanche, peu importe le contenu du champ 'date_first', la date de sortie est:

.fromtimestamp se contente de convertir un nombre flottant en date... Mais le zéro est au 1/1/1970.
Vous ce que vous cherchez à convertir c'est la représentation d'un jour Julien qui est un calendrier astronomique qui commence en 4713 avant Jesus Christ. Pour çà il faut trouver un algo. spécifique (Wikipedia aide bien pour çà) et le coder en Python.

- W
01/12/2017, 18h38
Sve@r

Citation:

Envoyé par habasque

La date de référence est January 1 2000 at 12h UTC soit 2451545 en nombre jour julien.

Déjà si on divise 2451545 par 365.2425 (division par 4 mais pas par 100 ou par 400) ça donne 6712.10. On les retranche à 2000 et on tombe sur -4712.10 avant JC ce qui est très très proche de la date donnée par wiztricks. Donc on peut considérer que 2451545 est le nombre de jours depuis le début du calendrier.
Maintenant si on dit que 2451545 c'est le 1/1/2000 (en nombre de jours), et que ce "-885.16" (premier nombre de ton tableau) est un nombre de jours négatifs à partir du 1/1/2000, alors on peut essayer jour=datetime(2000, 1, 1) + timedelta(-885.16). Ca donnera le 29/7/1997 20h09...

Bonjour,

Merci Sve@r !
Finalement voici comment j'ai procédé pour extraire le mois de mon champ "date_first" qui est au format juliandate:

Citation:

[ -885.16 -874.189 -863.149 ..., 3647.644 3628.51 3638.193]

Ci-dessous le code utilisé pour récupérer le mois de chaque date:
Code:

1 2 3 4 5 6 #calcul du mois pour chaque date month_vel = np.empty((np.size(date_first))) date_number = np.linspace(1,np.size(date_first),np.size(date_first),dtype='int32') for i in date_number: jour=datetime(2000, 1, 1) + timedelta(date_first[list(date_number).index(i)]) month_vel[list(date_number).index(i)]= jour.month