Бесплатный авторский астрософт

[Freddykrug]

После почти десятилетия заброшенности разработчик Winstars вернулся к ней и -судя по обновлениям версии - стал активно разрабатывать. https://winstars.net/ru/
Программа бесплатна для Линукса и платна для Винды и Андроида... В интернете можно поискать первую версию, которая распаковывалась в любую папку, которую можно было таскать (=портейбл).

[newerow1989]

Здравствуйте!

Ранее я публиковал свою программу "Атлас (карта) звездного неба" на этом форуме Бесплатный авторский астрософт.
Я ее обновил. Добавил следующие опции:
- моменты вычисления конфигурации планет, перигелия/афелия,
- промежуточные значения координат небесных объектов,
- угловое удаление планет от Солнца,
- вычисление звездных величин,
- атмосферная рефракция (в горизонтальной системе координат).

Прямая ссылка: http://newerow1989.ru/my_program/atlaszwezd/AtlasZwezd.rar

Жду ваши отзывы, пожелания. Если есть недоработки моей программы, тоже пишите!

[Toth]

Для расчета транзитов/прохождений МКС и вообще спутников по Солнцу, Луне, планетам, и неск. дипскаям.
https://disk.yandex.ru/d/cTvxrjpuxv6eSQ

[Toth]

Тут пару раз такое просили.
https://disk.yandex.ru/d/43K20UuFK5ee6Q

[GraY25]

Тут пару раз такое просили.
https://disk.yandex.ru/d/43K20UuFK5ee6Q

Круто!
А просто в виде формулы/функции есть возможность выложить?
Т.к. мне не для прямого использования а для интеграции в другой некоммерческий софт (сиинг-бот)

[Toth]

А просто в виде формулы/функции есть возможность выложить?

Там у меня много формул , я давно наделал всяких модулей. И там избыточная точность в некоторых местах, поэтому слишком много.
Проще вот тут - Помогите рассчитать азимут и высоту (угол) Солнца над горизонтом? Pluto выложил исходник на Си , там вроде все можно понять.

[GraY25]

Там у меня много формул , я давно наделал всяких модулей. И там избыточная точность в некоторых местах, поэтому слишком много.
Проще вот тут - Помогите рассчитать азимут и высоту (угол) Солнца над горизонтом? Pluto выложил исходник на Си , там вроде все можно понять.

То что надо!
Если кому надобно переделал под Python:

(кликните для показа/скрытия)

import math

def calc_sun_az_alt(year, month, day, ut, lon, lat):
    """
    Calculate sun's azimuth and altitude by UT date|time and coords lon lat
    :param: year, month, day, ut, lon, lat
    :return: azimut, altitude
    """
    def rev(q):
        return q - math.floor(q / 360.0) * 360.0
    rad_deg = 180.0 / math.pi
    deg_rad = math.pi / 180.0
    # Расчет количества целых суток прошедших от 01 января 0.0ч 2000г
    data = 367 * year - (7 * (year + (month + 9) / 12)) / 4 + (275 * month) / 9 + day - 730530
    data1 = data * 1.0 + ut / 24.0  # расчет даты с учетом текущего времени UT
    # Элементы орбиты Солнца (Земли)
    w = 282.9404 + 0.0000470935 * data1  # долгота перигелия
    e = 0.016709 - 0.000000001151 * data1  # эксцентриситет орбиты
    m = 356.047 + 0.9856002585 * data1  # средняя аномалия
    m = rev(m)
    oblecl = 23.4393 - 0.0000003563 * data1  # наклон плоскости эклиптики
    l = w + m
    l = rev(l)
    ee = m + (180.0 / math.pi) * e * math.sin(m * deg_rad) * (1 + e * math.cos(m * deg_rad))  # эксцентрическая аномалия
    x = math.cos(ee * deg_rad) - e
    y = math.sin(ee * deg_rad) * math.sqrt(1 - e * e)
    r = math.sqrt(x * x + y * y)
    v = math.atan2(y, x) * rad_deg
    longt = v + w  # эклиптическая долгота Солнца
    longt = rev(longt)
    x = r * math.cos(longt * deg_rad)
    y = r * math.sin(longt * deg_rad)
    z = 0.0
    xequat = x
    yequat = y * math.cos(deg_rad * oblecl) - z * math.sin(deg_rad * oblecl)
    zequat = y * math.sin(deg_rad * oblecl) + z * math.cos(deg_rad * oblecl)
    # экваториальные координаты Солнца
    ra = math.atan2(yequat, xequat) * rad_deg  # прямое восхождение
    dec = math.asin(zequat / r) * rad_deg  # склонение
    # расчет азимутальных координат
    gmst0 = l / 15.0 + 12.0  # звездное время для гринвичского меридиана
    sidtime = gmst0 + ut + lon / 15.0  # Звездное время
    if sidtime < 0:
        sidtime = sidtime + 24.0
    if sidtime > 24:
        sidtime = sidtime - 24.0
    ha = sidtime * 15.0 - ra  # Часовой угол
    xx = math.cos(ha * deg_rad) * math.cos(dec * deg_rad)
    yy = math.sin(ha * deg_rad) * math.cos(dec * deg_rad)
    zz = math.sin(dec * deg_rad)
    xhor = xx * math.sin(lat * deg_rad) - zz * math.cos(lat * deg_rad)
    yhor = yy
    zhor = xx * math.cos(lat * deg_rad) + zz * math.sin(lat * deg_rad)
    azimut = math.atan2(yhor, xhor) * rad_deg + 180  # азимут Солнца [град]
    altitude = math.asin(zhor) * rad_deg  # высота Солнца над горизонтом [град]
    return azimut, altitude

[Toth]

переделал под Python:

А я в Excel -
Там кстати коммент один неправильный - " Расчет количества целых суток прошедших от 01января 0.0ч 2000г"
На самом деле считается от 31.12.1999 00:00 . Но считается все правильно.
Только лишние преобразования градусы <-> радианы для некоторых промежуточных значений.