pyephem 항성시 내가 처음 <code>ephem</code>을 사용하고 문제가 내가 시간 각도에서 태양 시간을 결정하기 위해 몇 가지 스크립트를 작성했습니다 oberver.sidereal_time의 출력()</p> <p>을 이해하는 데 문제가 예기치 않은 결과를

Question

을 제공합니다. 첫 번째 것은 에펜 (ephem)을 사용하여 적경을 계산하고 Meeus의 Astronomical Algorithms에서 그리니치 평균 항성석을 얻는 수식을 사용합니다.이 시간은 경도와 함께 지역 평균 항성으로 변환 될 수 있습니다.

import sys 
from datetime import datetime, time, timedelta 
import ephem 

def hour_angle(dt, longit, latit, elev): 
    obs = ephem.Observer() 
    obs.date = dt.strftime('%Y/%m/%d %H:%M:%S') 
    obs.lon = longit 
    obs.lat = latit 
    obs.elevation = elev 
    sun = ephem.Sun() 
    sun.compute(obs) 
    # get right ascention 
    ra = ephem.degrees(sun.g_ra) 

    # get sidereal time at greenwich (AA ch12) 
    jd = ephem.julian_date(dt) 
    t = (jd - 2451545.0)/36525 
    theta = 280.46061837 + 360.98564736629 * (jd - 2451545) \ 
      + .000387933 * t**2 - t**3/38710000 

    # hour angle (AA ch13) 
    ha = (theta + longit - ra * 180/ephem.pi) % 360 
    return ha 

def main(): 
    if len(sys.argv) != 6: 
     print 'Usage: hour_angle.py [YYYY/MM/DD] [HH:MM:SS] [longitude] [latitude] [elev]' 
     sys.exit() 
    else: 
     dt = datetime.strptime(sys.argv[1] + ' ' + sys.argv[2], '%Y/%m/%d %H:%M:%S') 
     longit = float(sys.argv[3]) 
     latit = float(sys.argv[4]) 
     elev = float(sys.argv[5]) 

    # get hour angle 
    ha = hour_angle(dt, longit, latit, elev) 

    # convert hour angle to timedelta from noon 
    days = ha/360 
    if days > 0.5: 
     days -= 0.5 
    td = timedelta(days=days) 

    # make solar time 
    solar_time = datetime.combine(dt.date(), time(12)) + td 
    print solar_time 

if __name__ == '__main__': 
    main() 

이

내가 일부 데이터를 연결하면 내가 기대하는 출력을 제공 : 내가 쓴

> python hour_angle_ephem.py 2012/11/16 20:34:56 -122.2697 37.8044 3.0 
2012-11-16 12:40:54.697115 

두 번째 스크립트는 같은 방법으로 적경을 계산하지만, 로컬 명백한을 얻을 수 ephem의 sidereal_time()를 사용 항시적인 시간.

import sys 
from datetime import datetime, time, timedelta 
import math 
import ephem 

def solartime(observer, sun=ephem.Sun()): 
    sun.compute(observer) 
    # sidereal time == ra (right ascension) is the highest point (noon) 
    t = observer.sidereal_time() - sun.ra 
    return ephem.hours(t + ephem.hours('12:00')).norm # .norm for 0..24 

def main(): 
    if len(sys.argv) != 6: 
     print 'Usage: hour_angle.py [YYYY/MM/DD] [HH:MM:SS] [longitude] [latitude] [elev]' 
     sys.exit() 
    else: 
     dt = datetime.strptime(sys.argv[1] + ' ' + sys.argv[2], '%Y/%m/%d %H:%M:%S') 
     longit = float(sys.argv[3]) 
     latit = float(sys.argv[4]) 
     elev = float(sys.argv[5]) 

    obs = ephem.Observer() 
    obs.date = dt.strftime('%Y/%m/%d %H:%M:%S') 
    obs.lon = longit 
    obs.lat = latit 
    obs.elevation = elev 
    solar_time = solartime(obs) 
    print solar_time 

if __name__ == '__main__': 
    main() 

예상치 못한 결과가 나옵니다.

python hour_angle_ephem2.py 2012/11/16 20:34:56 -122.2697 37.8044 3.0 
9:47:50.83 

AFAIK 두 스크립트 사이의 유일한 차이는 로컬 평균 항성시의 첫번째 염기 시간 각도 고려 지구의 장동 소요 로컬 명백 항성시 두 번째 염기 시간 각도 반면 , 나는 그것이 매우 작은 요인이어야한다고 생각한다. 대신 약 3 시간의 차이가 있습니다. 아무도 내게 무슨 일이 일어 났는지 설명 할 수 있습니까?

Answer 1

PyEphem에 각도를 예상하는 원시 부동 소수점 숫자를 제공하면 먼저 부동 소수점 각도를 라디안으로 처리하므로 일관성을 유지하기 위해 각도를 라디안으로 먼저 변환 한 것으로 신뢰합니다. 하지만 두 번째 스크립트에서는 경도와 위도가도 단위로 표시되고 PyEphem에 라디안 단위로 제공됩니다. 당신처럼 Observer 모양의 .lon 및 .lat 속성이 무엇인지 확인하기 위해 print 문 또는 두 가지를 추가하면 당신은 결과를 볼 수 있습니다

print observer.lon #--> -7005:32:16.0 
print observer.lat #--> 2166:01:57.2 

내가 무엇을 대신 단순히 원시 경도를 제공한다하고 싶은 것을 생각하고 위도 문자열을 PyEphem에 추가하여 두 번째 스크립트에서 float() 번의 호출을 argv[3] 및 argv[4] 주변에서 제거하여 사람이 읽을 수있는 각도로 해석 할 수 있습니다.

$ python tmp11.py 2012/11/16 20:34:56 -122.2697 37.8044 3.0 
12:40:55.59