Python 星体太阳计算

我是python的新手，我想从数据帧的datetime列计算太阳的位置（方位角和天顶角）

我为这项任务找到了Astral模块，并使用各个时间戳进行了一些正确的计算。 但是当我试图处理整个datetime列时，我得到了一个错误

我还尝试将dataframe（datetime列）本地化为特定的时区，但错误保持不变

谢谢你的帮助

代码如下：

import datetime
from datetime import datetime
import pytz
pytz.all_timezones
from astral import Astral

d1=pd.to_datetime(d1)
d1=d1.dt.tz_localize("Europe/Berlin")

输出如下所示：

0   2015-02-04 09:10:21
1   2016-03-05 11:30:25
dtype: datetime64[ns]

方位角计算：

print(Astral().solar_azimuth(d1, 52.52, -13.40))

错误：

---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-95-2a14f69465f4> in <module>
----> 1 print(Astral().solar_azimuth(d1, 52.52, -13.40))

~\Anaconda3\lib\site-packages\astral.py in solar_azimuth(self, dateandtime, latitude, longitude)
   2276             latitude = -89.8
   2277 
-> 2278         if dateandtime.tzinfo is None:
   2279             zone = 0
   2280             utc_datetime = dateandtime

~\Anaconda3\lib\site-packages\pandas\core\generic.py in __getattr__(self, name)
   5137             if self._info_axis._can_hold_identifiers_and_holds_name(name):
   5138                 return self[name]
-> 5139             return object.__getattribute__(self, name)
   5140 
   5141     def __setattr__(self, name: str, value) -> None:

AttributeError: 'Series' object has no attribute 'tzinfo'

---------------------------------------------------------------------------
AttributeError回溯（最近一次呼叫上次）
在里面
---->1个打印（Astral（）.太阳方位角（d1，52.52，-13.40））
~\Anaconda3\lib\site packages\astral.py太阳方位角（self、dateandtime、纬度、经度）
2276纬度=-89.8
2277
->2278如果dateandtime.tzinfo为无：
2279区域=0
2280 utc_datetime=日期和时间
~\Anaconda3\lib\site packages\pandas\core\generic.py in\uuuuuu getattr\uuuuuu（self，name）
5137如果自身信息轴可以保存标识符，并且保存名称（名称）：
5138返回自我[姓名]
->5139返回对象。\uuuu getattribute\uuuuu（self，name）
5140
5141定义设置属性（self，名称：str，值）->无：
AttributeError:“Series”对象没有属性“tzinfo”

不确定您的代码示例是如何工作的，但基本上，您可以

将
太阳方位角
功能应用于
熊猫系列
，例如

import pandas as pd
from astral.location import Location, LocationInfo

# define your location
l = LocationInfo('name', 'region', 'timezone/name', 52.52, -13.40)

# assuming input datetimes are in the form of a pandas Series
d1 = pd.to_datetime(["2015-02-04 09:10:21", "2016-03-05 11:30:25"]).to_series().dt.tz_localize('Europe/Berlin')

azim = d1.apply(Location(l).solar_azimuth)

print(azim)
# 2015-02-04 09:10:21    112.190070
# 2016-03-05 11:30:25    137.505145
# dtype: float64

请注意，您也可以在没有熊猫的情况下执行此操作，例如

from datetime import datetime
from zoneinfo import ZoneInfo # Python 3.9

dt = ["2015-02-04 09:10:21", "2016-03-05 11:30:25"]
# set time zone
dt = [datetime.fromisoformat(t).replace(tzinfo=ZoneInfo('Europe/Berlin')) for t in dt]
# azimut calc
azim = [Location(l).solar_azimuth(t) for t in dt]

print(azim)
# [112.19006967766785, 137.50514457536954]

谢谢这就成功了。我忘了提到我以前使用过旧版本的Astral，但现在升级到了第2版。@mixi很高兴能帮上忙！一个小的补充（因为您说您是Python新手）：我建议使用它来处理时区
astral
但是在内部使用了
pytz
，Python3.9版本不推荐使用pytz。因此，如果您遇到问题，您可能必须返回到使用
pytz
。再次感谢。你帮了大忙。我也会调查zoneinfo。亲切的问候！