本文整理汇总了Python中sunpy.time.julian_centuries函数的典型用法代码示例。如果您正苦于以下问题:Python julian_centuries函数的具体用法?Python julian_centuries怎么用?Python julian_centuries使用的例子?那么恭喜您, 这里精选的函数代码示例或许可以为您提供帮助。
在下文中一共展示了julian_centuries函数的20个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于我们的系统推荐出更棒的Python代码示例。
示例1: apparent_longitude
def apparent_longitude(t=None):
"""Returns the apparent longitude of the Sun."""
T = julian_centuries(t)
omega = 259.18 - 1934.142 * T
true_long = true_longitude(t)
result = true_long - 0.00569 - 0.00479 * np.sin(np.radians(omega))
return result
开发者ID:khughitt,项目名称:sunpy,代码行数:7,代码来源:sun.py
示例2: mean_anomaly
def mean_anomaly(t='now'):
"""Returns the mean anomaly (the angle through which the Sun has moved
assuming a circular orbit) as a function of time."""
T = julian_centuries(t)
result = 358.475830 + 35999.049750 * T - 0.0001500 * T ** 2 - 0.00000330 * T ** 3
result = result * u.deg
return Longitude(result)
开发者ID:Alex-Ian-Hamilton,项目名称:sunpy,代码行数:7,代码来源:sun.py
示例3: apparent_longitude
def apparent_longitude(t='now'):
"""Returns the apparent longitude of the Sun."""
T = julian_centuries(t)
omega = (259.18 - 1934.142 * T) * u.deg
true_long = true_longitude(t)
result = true_long - (0.00569 - 0.00479 * np.sin(omega)) * u.deg
return Longitude(result)
开发者ID:Alex-Ian-Hamilton,项目名称:sunpy,代码行数:7,代码来源:sun.py
示例4: true_obliquity_of_ecliptic
def true_obliquity_of_ecliptic(t='now'):
"""
Returns the true obliquity of the ecliptic.
"""
T = julian_centuries(t)
result = 23.452294 - 0.0130125 * T - 0.00000164 * T**2 + 0.000000503 * T**3
return Angle(result, u.deg)
开发者ID:DanRyanIrish,项目名称:sunpy,代码行数:8,代码来源:sun.py
示例5: equation_of_center
def equation_of_center(t=None):
"""Returns the Sun's equation of center (in degrees)"""
T = julian_centuries(t)
mna = mean_anomaly(t)
result = ((1.9194600 - 0.0047890 * T - 0.0000140 * T
** 2) * np.sin(np.radians(mna) + (0.0200940 - 0.0001000 * T) *
np.sin(np.radians(2 * mna)) + 0.0002930 * np.sin(np.radians(3 * mna))))
return result
开发者ID:khughitt,项目名称:sunpy,代码行数:8,代码来源:sun.py
示例6: eccentricity_SunEarth_orbit
def eccentricity_SunEarth_orbit(t='now'):
"""
Returns the eccentricity of the Sun Earth Orbit.
"""
T = julian_centuries(t)
result = 0.016751040 - 0.00004180 * T - 0.0000001260 * T**2
return result
开发者ID:DanRyanIrish,项目名称:sunpy,代码行数:8,代码来源:sun.py
示例7: equation_of_center
def equation_of_center(t='now'):
"""Returns the Sun's equation of center (in degrees)"""
T = julian_centuries(t)
mna = mean_anomaly(t)
result = ((1.9194600 - 0.0047890 * T - 0.0000140 * T**2) * np.sin(mna) +
(0.0200940 - 0.0001000 * T) * np.sin(2 * mna) + 0.0002930 * np.sin(3 * mna))
result = result * u.deg
return Angle(result)
开发者ID:Hypnus1803,项目名称:sunpy,代码行数:8,代码来源:sun.py
示例8: mean_ecliptic_longitude
def mean_ecliptic_longitude(t='now'):
"""
Returns the mean ecliptic longitude.
"""
T = julian_centuries(t)
result = 279.696680 + 36000.76892 * T + 0.0003025 * T**2
result = result * u.deg
return Longitude(result)
开发者ID:DanRyanIrish,项目名称:sunpy,代码行数:9,代码来源:sun.py
示例9: geometric_mean_longitude
def geometric_mean_longitude(t='now'):
"""
Returns the geometric mean longitude (in degrees).
"""
T = julian_centuries(t)
result = 279.696680 + 36000.76892 * T + 0.0003025 * T**2
result = result * u.deg
return Longitude(result)
开发者ID:DanRyanIrish,项目名称:sunpy,代码行数:9,代码来源:sun.py
示例10: eccentricity_sun_earth_orbit
def eccentricity_sun_earth_orbit(t='now'):
"""
Returns the eccentricity of the Sun Earth Orbit.
Parameters
----------
t : {parse_time_types}
A time (usually the start time) specified as a palongitude_Sun_perigee, true_latitude and apparent_latitude need to be fixedrse_time-compatible
time string, number, or a datetime object.
"""
T = julian_centuries(t)
result = 0.016751040 - 0.00004180 * T - 0.0000001260 * T**2
return result
开发者ID:Cadair,项目名称:sunpy,代码行数:13,代码来源:sun.py
示例11: true_obliquity_of_ecliptic
def true_obliquity_of_ecliptic(t='now'):
"""
Returns the true obliquity of the ecliptic.
Parameters
----------
t : {parse_time_types}
A time (usually the start time) specified as a parse_time-compatible
time string, number, or a datetime object.
"""
T = julian_centuries(t)
result = 23.452294 - 0.0130125 * T - 0.00000164 * T**2 + 0.000000503 * T**3
return Angle(result, u.deg)
开发者ID:Cadair,项目名称:sunpy,代码行数:13,代码来源:sun.py
示例12: mean_ecliptic_longitude
def mean_ecliptic_longitude(t='now'):
"""
Returns the mean ecliptic longitude.
Parameters
----------
t : {parse_time_types}
A time (usually the start time) specified as a parse_time-compatible
time string, number, or a datetime object.
"""
T = julian_centuries(t)
result = 279.696680 + 36000.76892 * T + 0.0003025 * T**2
result = result * u.deg
return Longitude(result)
开发者ID:Cadair,项目名称:sunpy,代码行数:14,代码来源:sun.py
示例13: apparent_longitude
def apparent_longitude(t='now'):
"""
Returns the apparent longitude of the Sun.
Parameters
----------
t : {parse_time_types}
A time (usually the start time) specified as a parse_time-compatible
time string, number, or a datetime object.
"""
T = julian_centuries(t)
omega = (259.18 - 1934.142 * T) * u.deg
true_long = true_longitude(t)
result = true_long - (0.00569 - 0.00479 * np.sin(omega)) * u.deg
return Longitude(result)
开发者ID:Cadair,项目名称:sunpy,代码行数:15,代码来源:sun.py
示例14: solar_north
def solar_north(t=None):
"""Returns the position of the Solar north pole in degrees."""
T = julian_centuries(t)
ob1 = true_obliquity_of_ecliptic(t)
# in degrees
i = 7.25
k = 74.3646 + 1.395833 * T
lamda = true_longitude(t) - 0.00569
omega = apparent_longitude(t)
lamda2 = lamda - 0.00479 * np.sin(np.radians(omega))
diff = np.radians(lamda - k)
x = np.degrees(np.arctan(-np.cos(np.radians(lamda2)*np.tan(np.radians(ob1)))))
y = np.degrees(np.arctan(-np.cos(diff) * np.tan(np.radians(i))))
result = x + y
return result
开发者ID:khughitt,项目名称:sunpy,代码行数:15,代码来源:sun.py
示例15: mean_anomaly
def mean_anomaly(t='now'):
"""
Returns the mean anomaly (the angle through which the Sun has moved
assuming a circular orbit) as a function of time.
Parameters
----------
t : {parse_time_types}
A time (usually the start time) specified as a parse_time-compatible
time string, number, or a datetime object.
"""
T = julian_centuries(t)
result = 358.475830 + 35999.049750 * T - 0.0001500 * T**2 - 0.00000330 * T**3
result = result * u.deg
return Longitude(result)
开发者ID:Cadair,项目名称:sunpy,代码行数:15,代码来源:sun.py
示例16: solar_north
def solar_north(t='now'):
"""Returns the position of the Solar north pole in degrees."""
T = julian_centuries(t)
ob1 = true_obliquity_of_ecliptic(t)
# in degrees
i = 7.25 * u.deg
k = (74.3646 + 1.395833 * T) * u.deg
lamda = true_longitude(t) - (0.00569 * u.deg)
omega = apparent_longitude(t)
lamda2 = lamda - (0.00479 * np.sin(omega)) * u.deg
diff = lamda - k
x = np.arctan(-np.cos((lamda2) * np.tan(ob1)))
y = np.arctan(-np.cos(diff) * np.tan(i))
result = x + y
return Angle(result.to(u.deg))
开发者ID:Alex-Ian-Hamilton,项目名称:sunpy,代码行数:15,代码来源:sun.py
示例17: equation_of_center
def equation_of_center(t='now'):
"""
Returns the Sun's equation of center (in degrees).
Parameters
----------
t : {parse_time_types}
A time (usually the start time) specified as a parse_time-compatible
time string, number, or a datetime object.
"""
T = julian_centuries(t)
mna = mean_anomaly(t)
result = ((1.9194600 - 0.0047890 * T - 0.0000140 * T**2) * np.sin(mna) +
(0.0200940 - 0.0001000 * T) * np.sin(2 * mna) + 0.0002930 * np.sin(3 * mna))
result = result * u.deg
return Angle(result)
开发者ID:Cadair,项目名称:sunpy,代码行数:16,代码来源:sun.py
示例18: heliographic_solar_center
def heliographic_solar_center(t='now'):
"""Returns the position of the solar center in heliographic coordinates."""
jd = julian_day(t)
T = julian_centuries(t)
# Heliographic coordinates in degrees
theta = ((jd - 2398220)*360/25.38) * u.deg
i = 7.25 * u.deg
k = (74.3646 + 1.395833 * T) * u.deg
lamda = true_longitude(t) - 0.00569 * u.deg
diff = lamda - k
# Latitude at center of disk (deg):
he_lat = np.degrees(np.arcsin(np.sin(diff)*np.sin(i)))
# Longitude at center of disk (deg):
y = -np.sin(diff)*np.cos(i)
x = -np.cos(diff)
rpol = (np.arctan2(y, x))
he_lon = rpol - theta
return [Longitude(he_lon), Latitude(he_lat)]
开发者ID:Alex-Ian-Hamilton,项目名称:sunpy,代码行数:18,代码来源:sun.py
示例19: heliographic_solar_center
def heliographic_solar_center(t=None):
"""Returns the position of the solar center in heliographic coordinates."""
jd = julian_day(t)
T = julian_centuries(t)
# Heliographic coordinates in degrees
theta = (jd - 2398220)*360/25.38
i = 7.25
k = 74.3646 + 1.395833 * T
lamda = true_longitude(t) - 0.00569
#omega = apparent_longitude(t)
#lamda2 = lamda - 0.00479 * math.sin(np.radians(omega))
diff = np.radians(lamda - k)
# Latitude at center of disk (deg):
he_lat = np.degrees(np.arcsin(np.sin(diff)*np.sin(np.radians(i))))
# Longitude at center of disk (deg):
y = -np.sin(diff)*np.cos(np.radians(i))
x = -np.cos(diff)
rpol = cmath.polar(complex(x,y))
he_lon = np.degrees(rpol[1]) - theta
he_lon = he_lon % 360
if he_lon < 0:
he_lon = he_lon + 360.0
return [he_lon, he_lat]
开发者ID:khughitt,项目名称:sunpy,代码行数:24,代码来源:sun.py
示例20: geometric_mean_longitude
def geometric_mean_longitude(t=None):
"""Returns the geometric mean longitude (in degrees)"""
T = julian_centuries(t)
result = 279.696680 + 36000.76892 * T + 0.0003025 * T ** 2
result = result % 360.0
return result
开发者ID:khughitt,项目名称:sunpy,代码行数:6,代码来源:sun.py
注:本文中的sunpy.time.julian_centuries函数示例由纯净天空整理自Github/MSDocs等源码及文档管理平台,相关代码片段筛选自各路编程大神贡献的开源项目,源码版权归原作者所有,传播和使用请参考对应项目的License;未经允许,请勿转载。 |
请发表评论