日出日落时间计算计算

import ephem
import datetime
'''局地经纬度时间计算日出日落的时间，以UTC的形式输出'''
def calculate_sunrise_sunset(latitude, longitude, date):
    # 创建观测者对象
    observer = ephem.Observer()
    observer.lat = str(latitude)  # 纬度 40.7128
    observer.lon = str(longitude)  # 经度 -74.0060
    observer.elevation = 0  # 海拔，设为0即为海平面

    # 计算日出和日落时间
    sun = ephem.Sun()
    sunrise = observer.next_rising(sun, start=date)
    sunset = observer.next_setting(sun, start=date) #'2024/04/26'
    
    sunrise=datetime.datetime(1899,12,31)+datetime.timedelta(days=float(sunrise))
    sunset=datetime.datetime(1899,12,31)+datetime.timedelta(days=float(sunset))
    
    return sunrise, sunset