Не полный ответ на ваш вопрос. Просто отправная точка.
tz = Entity["TimeZone", "Australia/Sydney"];
date = DateObject[{2019, 07, 14}, TimeZone -> tz];
Sunrise[Entity["City", {"Sydney", "NewSouthWales", "Australia"}], date, TimeZone -> tz]
Аналогично для Sunset
.
sunset = Sunset[
Entity["City", {"Sydney", "NewSouthWales", "Australia"}],
DateObject[{2019, 07, 14}, TimeZone -> tz], TimeZone -> tz]
Вы можете прибавить / вычесть 2 или 3 часа из результата, чтобы получить «астрономическое» значение.
Фаза луны
MoonPhase[date]
(* 0.8913 *)
MoonPhase[date, "Name"]["Name"]
(* waxing gibbous moon *)
Положение Луны на закате
MoonPosition[Entity["City", {"Sydney", "NewSouthWales", "Australia"}], sunset]
(* {98.67\[Degree],25.03\[Degree]} *)
Взгляните на это видео для обзора астрономических и космических данных об объектах в WL.