与えられた日時のある場所で太陽の位置を計算するにはどうすればいいですか？

Question

UTC時刻（時、分、秒）、経度（deg E）、緯度（deg N）、および日付があります。 Python 2.7で太陽天頂角を計算するコードを誰かに教えてもらえますか？

Answer 1

これは面白い問題ですが、私は良い答えがあると思います。

素晴らしいastropy packageをご覧ください。私はあなたがcoordinates moduleを使用する必要があると信じています。これらの線に沿って

何か：ここ

import astropy.coordinates as coord 
from astropy.time import Time 
import astropy.units as u 


loc = coord.EarthLocation(lon=0.1 * u.deg, 
          lat=51.5 * u.deg) 
now = Time.now() 

altaz = coord.AltAz(location=loc, obstime=now) 
sun = coord.get_sun(now) 

print(sun.transform_to(altaz).alt) 

は、我々が現在の時刻に0.1度経度と緯度51.5場所の地平線上に太陽の角度を取得しています。

FYI、.zen天頂角を与えるでしょう。 alecxeの答え@

Answer 2

は素晴らしいですが、私は元の質問が

from astropy.coordinates import get_sun, AltAz, EarthLocation 
from astropy.time import Time 

sun_time = Time('2017-12-6 17:00') #UTC time 
loc = EarthLocation.of_address('Baltimore, MD') # anything the google geocoding API resolves 
altaz = AltAz(obstime=sun_time, location=loc) 

zen_ang = get_sun(sun_time).transform_to(altaz).zen 

zen_ang（特定の時間の天頂角）を求めているものに少し近いです若干の修正を追加しようと思いましたAngleオブジェクトです - the docsの詳細については、numpyスカラーと "degrees"の関連単位のように機能します。