Newer
Older
# tools.py
import re
import math
from datetime import datetime
import julian
DEGREE = u"\u00B0" # u"\N{DEGREE SIGN}"
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
def dateTime2julian(month, day, year, hour, minute, second):
# initialize datetime object
dt = datetime
# format date and time to "May 09 2011 16:33:53"
dateTime = "%s/%s/%s %s:%s:%s" % (
day, month, year, hour, minute, second)
# dateTime conversion to "09/05/2011 16:33:53"
dateTime = "%s" % \
(dt.strptime(dateTime,
"%d/%b/%Y %H:%M:%S").strftime("%d/%m/%Y %H:%M:%S"))
# conversion to "20110509163353"
epic_date = "%s" % \
(dt.strptime(dateTime,
"%d/%m/%Y %H:%M:%S").strftime("%Y%m%d%H%M%S"))
# conversion to julian day
julian = float((dt.strptime(dateTime, "%d/%m/%Y %H:%M:%S").strftime("%j"))) \
+ ((float(hour) * 3600.) +
(float(minute) * 60.) + float(second)) / 86400.
# we use julian day with origine 0
julian -= 1
return julian
def julian2dt(jd):
# see: https://en.wikipedia.org/wiki/Julian_day
# Julian Date 12h Jan 1, 4713 BC
# Modified JD 0h Nov 17, 1858 JD − 2400000.5
# CNES JD 0h Jan 1, 1950 JD − 2433282.5
jd = jd + JULIAN
dt = julian.from_jd(jd, fmt='mjd')
return dt
def dt2julian(dt):
jd = julian.to_jd(dt, fmt='mjd')
jd = jd - JULIAN
return jd
def dateTime2epic(month, day, year, hour, minute, second):
# initialize datetime object
dt = datetime
# format date and time to "May 09 2011 16:33:53"
dateTime = "%s/%s/%s %s:%s:%s" % (
day, month, year, hour, minute, second)
# dateTime conversion to "09/05/2011 16:33:53"
dateTime = "%s" % \
(dt.strptime(dateTime,
"%d/%b/%Y %H:%M:%S").strftime("%d/%m/%Y %H:%M:%S"))
# conversion to "20110509163353"
epic_date = "%s" % \
(dt.strptime(dateTime,
"%d/%m/%Y %H:%M:%S").strftime("%Y%m%d%H%M%S"))
return epic_date
# Dec2dmc convert decimal position to degree, mim with centieme string,
# hemi = 0 for latitude, 1 for longitude
def Dec2dmc(position, hemi):
if re.match('[EW]', hemi):
neg = 'W'
pos = 'E'
else:
neg = 'S'
pos = 'N'
if position < 0:
geo = neg
else:
geo = pos
# get integer and decimal part
dec, intg = math.modf(position)
dec = abs(dec)
intg = abs(intg)
if re.match('[EW]', hemi):
str = "{:0>3.0f}{:s}{:0>7.4f} {}".format(intg, DEGREE, (dec / 100 ) * 6000, geo)
str = "{:0>2.0f}{:s}{:0>7.4f} {}".format(intg, DEGREE, (dec / 100 ) * 6000,geo)
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
return str
# Dec2dms convert decimal position to degree, mim with second string,
# hemi = 0 for latitude, 1 for longitude
def Dec2dms(position, hemi):
if re.match('[EW]', hemi):
neg = 'W'
pos = 'E'
else:
neg = 'S'
pos = 'N'
if position < 0:
geo = neg
else:
geo = pos
# get integer and decimal part
dec, intg = math.modf(position)
# get integer and decimal part of min.sec
sec, min = math.modf((dec / 100) * 6000)
if re.match('[EW]', hemi):
str = "{:0>3.0f}{:s}{:2.4f} {}".format(intg, DEGREE, min + sec/100*60, geo)
str = "{:0>2.0f}{:s}{:2.4f} {}".format(intg, DEGREE, min + sec/100*60, geo)