Путаница с использованием dec / ra для вычисления сублунного местоположения - PullRequest
Новая
       29

Путаница с использованием dec / ra для вычисления сублунного местоположения

0 голосов
/ 14 сентября 2018

Мне нужно рассчитать кажущийся азимут и углы места, а также сублунный широту / долготу для данной даты / времени. Я получаю, что углы аз / эль обычно согласуются с другими источниками (MoonCalc.org, Горизонты и т. Д.), Но нет хороших источников сравнения для сублунного широты / долготы. Что еще более важно, я сомневаюсь, что широта / долгота, которую я получаю, используя значения dec / ra, потому что ra почти не меняется в течение длительных периодов времени.

Вот основной звонок, который я делаю: 

roc.date='2018/1/1 01:00:00'
moon=ephem.Moon(roc)
print('rocMoonTest: %s UTC-4, lat/lon = %0.4f [+N], %0.4f [+E]' %
    (roc.date, math.degrees(roc.lat), math.degrees(roc.lon)))
print('Moon dec/ra     = %s [+N], %s [+W]' % (moon.dec, moon.ra ))
print('Moon a_dec/a_ra = %s [+N], %s [+W]' % (moon.a_dec, moon.a_ra ))
print('Moon g_dec/g_ra = %s [+N], %s [+W]' % (moon.g_dec, moon.g_ra ))
print('Moon az/el      = %0.4f, %0.4f' %
    (math.degrees(moon.az), math.degrees(moon.alt)))

А потом я повторяю это каждые 3 часа. Ниже вывод: 

rocMoonTest: 2018/1/1 01:00:00 UTC-4, lat/lon = 43.0000 [+N], -78.0000 [+E]
Moon dec/ra     = 18:53:07.1 [+N], 5:43:03.33 [+W]
Moon a_dec/a_ra = 19:22:21.3 [+N], 5:39:38.43 [+W]
Moon g_dec/g_ra = 19:22:44.7 [+N], 5:40:41.41 [+W]
Moon az/el      = 105.3953, 43.0670

rocMoonTest: 2018/1/1 04:00:00 UTC-4, lat/lon = 43.0000 [+N], -78.0000 [+E]
Moon dec/ra     = 19:07:55.4 [+N], 5:49:00.24 [+W]
Moon a_dec/a_ra = 19:32:24.2 [+N], 5:47:42.22 [+W]
Moon g_dec/g_ra = 19:32:35.1 [+N], 5:48:45.29 [+W]
Moon az/el      = 169.5907, 65.8406

rocMoonTest: 2018/1/1 07:00:00 UTC-4, lat/lon = 43.0000 [+N], -78.0000 [+E]
Moon dec/ra     = 19:13:15.7 [+N], 5:54:49.89 [+W]
Moon a_dec/a_ra = 19:41:07.2 [+N], 5:55:47.50 [+W]
Moon g_dec/g_ra = 19:41:05.5 [+N], 5:56:50.65 [+W]
Moon az/el      = 246.5737, 49.4664

Как и ожидалось, и согласно аз / углам, Луна поворачивается с Востока на Запад, когда Земля вращается и достигает пиковой высоты где-то в течение периода. Однако ни одно из различных значений dec / ra существенно не изменяется. За этот 6-часовой промежуток я бы ожидал увидеть примерно 6-часовое изменение в РА. Очевидно, что когда я использую любое из этих значений ra для вычисления долготы, я получаю неправильный ответ. Похоже, система отсчета для dev / ra не вращается вместе с землей. Тем не менее, документы указывают, что я должен ожидать этого.

Кто-нибудь хочет объяснить, где я ошибся в своем понимании различных переменных правильного вознесения и какой самый прямой способ вычислить сублунный широту / долготу? Обратите внимание, что я бы предпочел не использовать подход, который поворачивает кажущуюся позицию аз / эл в геодезическую широту / долготу.

Ответы [ 2 ]

0 голосов
/ 23 апреля 2019

относительно этой части вашего вопроса

Какой самый прямой способ вычислить сублунный широту / долготу?

Вот мой код для расчета подлунной точки. 

greenwich = ephem.Observer()
greenwich.lat = "0"
greenwich.lon = "0"
greenwich.date = datetime.utcnow()

#add Moon Sub Solar Point
moon = ephem.Moon(greenwich)
moon.compute(greenwich.date)
moon_lon = math.degrees(moon.ra - greenwich.sidereal_time() )
# map longitude value from -180 to +180 
if moon_lon < -180.0 :
  moon_lon = 360.0 + moon_lon 
elif moon_lon > 180.0 :
  moon_lon = moon_lon - 360.0

moon_lat = math.degrees(moon.dec)
print "moon Lon:",moon_lon, "Lat:",moon_lat

Надеюсь, это поможет. Я также использую тот же подход для расчета солнечной точки. Прекрасно работает для меня.

РЕДАКТИРОВАТЬ: Да ... широта Гринвича установлена ​​на ноль, потому что это не имеет никакого значения для этого расчета.

0 голосов
/ 16 сентября 2018

Измерение «Правильное Вознесение» проводится не против вращающейся поверхности Земли, а против неподвижных звезд неба - это своего рода долгота, но источником которой является точка на звездной карте, где два великих «экватора неба, экватора Земли и экватора Солнечной системы, эклиптики (которая не является истинно экватором Солнечной системы, потому что это плоскость орбиты Земли, а не средневзвешенное значение для всехпланетные орбиты), крест.

Поскольку точка их пересечения движется с течением времени, система Правильного Вознесения каждый год очень немного отличается и сильно различается на протяжении веков и тысячелетий. Таким образом, Правильное Вознесение и склонение(небесная широта) всегда нужно указывать относительно некоторой даты, определенного точного момента, например, B1950 или J2000.

Теперь существует фиксированная система координат для RA и dec, которая не перемещается, ICRS, которая являетсяориентируется как J2000, но определяется с помощью позиций of квазаров, которые (как мы предполагаем) не будут измеримо двигаться в течение жизни нашего вида.

...