在这里,“ OK ”当然意味着 AYOR(风险自负),但如果避免与现有属性名称发生明显冲突,则不会出现可预见的问题。
Skyfield物体——尤其是行星——通常具有有限数量的属性。我经常编写简短的脚本来提取我保存为文本并稍后使用的数字数据。这些本质上是“一次性”脚本,因为我很少使用它们超过一次或两次,并且从不共享它们。
当我编写更持久的代码时,我当然会创建自己的容器对象。
我的问题:它似乎对我来说效果很好,所以在这个特定的上下文中,除了属性名称冲突之外,还有什么可能出错的吗?
from skyfield.api import load
import numpy as np
import matplotlib.pyplot as plt
from mpl_toolkits.mplot3d import Axes3D
eph = load('de421.bsp')
earth = eph['earth']
sun = eph['sun']
ts = load.timescale()
t = ts.utc(2016, 1, np.linspace(0, 366, 1000))
# SLOPPY WAY: just add them directly
earth.pos = earth.at(t).position.km
sun.pos = sun.at(t).position.km
earth.r = np.sqrt(((earth.pos-sun.pos)**2).sum(axis=0))
earth.peri = earth.r.min()
earth.apo = earth.r.max()
print earth.peri, earth.apo, earth.pos.shape
# BETTER WAY: tedious but more cautious
uhoh = dict()
ep = earth.at(t).position.km
sp = sun.at(t).position.km
r = np.sqrt(((ep-sp)**2).sum(axis=0))
uhoh['pos'] = ep
uhoh['r'] = r
uhoh['peri'] = r.min()
uhoh['apo'] = r.max()
earth.uhoh = uhoh
print earth.uhoh['peri'], earth.uhoh['apo'], earth.uhoh['pos'].shape
返回:
147100175.99 152103762.948 (3, 1000)
147100175.99 152103762.948 (3, 1000)