astropyで初日の出の時間を計算する🌅 ref: https://qiita.com/nanki/items/97191fee71ec5a3255f5

file0.txt

import astropy.coordinates
import astropy.units as u
from astropy.coordinates import EarthLocation, AltAz
from astropy.time import Time
from datetime import datetime
from scipy import optimize

tt = EarthLocation(lat=35.6585 * u.deg, lon=139.7455 * u.deg, height=250 * u.m)

def alt(timestamp):
    t = Time(timestamp, format='unix')
    sun = astropy.coordinates.get_sun(t)
    return sun.transform_to(AltAz(obstime=t, location=tt, pressure=1013 * u.hPa)).alt / u.deg + 16 / 60

hinode = optimize.bisect(alt, datetime(2018, 1, 1).timestamp(), datetime(2018, 1, 1, 9).timestamp())
print(datetime.fromtimestamp(hinode))

file1.txt

2018-01-01 06:52:44.560696