SkyMapMod 0.1.21__tar.gz → 0.1.22__tar.gz

{skymapmod-0.1.21 → skymapmod-0.1.22}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: SkyMapMod
-Version: 0.1.21
+Version: 0.1.22
 Summary: for modeling
 Project-URL: Homepage, https://istina.msu.ru/profile/aminevann/
 Author-email: Anna Amineva <amineva.2219@physics.msu.ru>

{skymapmod-0.1.21 → skymapmod-0.1.22}/SkyMapMod/__init__.py

@@ -8,4 +8,5 @@ from .band_V_data import *
 from .solar_radio_flux import *
 from .star_catalogues import *
 from .albedo_of_planets import *
-from .planets import *
+from .planets import *
+from .total_background import *

{skymapmod-0.1.21 → skymapmod-0.1.22}/SkyMapMod/planets.py

@@ -18,7 +18,7 @@ def coordinates_of_planet(body, date, time): #'mars', 'venus', 'saturn', 'jupyte
     return(l, b)
 
     
-def magnitude_of_planet(body, date, time): #'mars', 'venus', 'saturn', 'jupyter'; 2023-11-03; 12:00:00';
+def magnitude_of_planet(body, date, time): #'mars', 'venus', 'saturn', 'jupiter'; 2023-11-03; 12:00:00';
     if body == 'venus':
         planet = ephem.Venus()
     elif body == 'mars':

{skymapmod-0.1.21 → skymapmod-0.1.22}/SkyMapMod/star_catalogues.py

@@ -7,7 +7,12 @@ from.zodiac import black_body_phot
 #на вход подаем галактические координаты (l -- долгота, b -- широта) и массив длин волн, по которым строим спектр
 #РАЗОБРАТЬСЯ С КОЭФФИЦИЕНТАМИ А -- НАДО ЛИ НОРМИРОВАТЬ НА C (КАЖЕТСЯ НАДО) А ЕЩЕ ВРОДЕ Я ЕГО ОТНОРМИРОВАЛА И ПОДГРУЗИЛА -- ПРОВЕРИТЬ
 def star_spectrum(l, b, lmbd):
-    j = int((1800 - round(l*10)) + 3600) #пересчет из галактических к позиции в массиве
+    #пересчет из галактических к позиции в массиве, ЕСЛИ ПРОБЛЕМЫ ВОЗНИКЛИ, ПЕРЕПРОВЕРИТЬ
+    if round(l*10) >= 1800:
+        j = int((1800 - round(l*10)) + 3600)
+    elif round(l*10) < 1800:
+        j =  int(1800 - round((l+10)*10))
+     
     i = int(900 + round(b*10))
     star_brightness = load_star_brightness()
     star_temperatures = load_star_temperatures()    

skymapmod-0.1.22/SkyMapMod/total_background.py

@@ -0,0 +1,116 @@
+from astropy.coordinates import Galactic, ICRS, GeocentricMeanEcliptic, get_sun, get_body
+from astropy.time import Time
+import astropy.units as u
+
+from .zodiac import *
+from .airglow import *
+from .transparency import *
+from .airglow_spectrum import *
+from .modtran_default_kp_transparency import *
+from .solar_spectrum import *
+from .band_V_data import *
+from .solar_radio_flux import *
+from .star_catalogues import *
+from .albedo_of_planets import *
+from .planets import *
+
+def galactic_to_equatorial(l, b):
+# Входные данные: галактические координаты и эпоха
+    galactic_l = l * u.deg  # Галактическая долгота
+    galactic_b = b * u.deg # Галактическая широта
+    epoch = Time('J2000.0')  # Эпоха (стандартная для астрономических расчетов)
+
+    # Создаем объект галактических координат
+    galactic_coords = Galactic(l=galactic_l, b=galactic_b)
+
+    # Преобразуем галактические координаты в экваториальные (ICRS)
+    ecliptic_coords = galactic_coords.transform_to(GeocentricMeanEcliptic(equinox=epoch))
+
+
+    print("Техническая проверка")
+    print(f"Эклиптическая долгота (λ): {ecliptic_coords.lon.to(u.deg)}")
+    print(f"Эклиптическая широта (β): {ecliptic_coords.lat.to(u.deg)}")
+    return(ecliptic_coords.lon.to(u.deg).value, ecliptic_coords.lat.to(u.deg).value)
+
+def Sun_ecl_lon(date, time):
+    datetime = date + ' ' + time
+    observation_time = Time(datetime)
+
+    # Получаем координаты Солнца в геоцентрической системе координат
+    sun_position = get_sun(observation_time)
+
+    # Преобразуем координаты Солнца в эклиптическую систему (геоцентрическую)
+    ecliptic_coords = sun_position.transform_to(GeocentricMeanEcliptic())
+
+    # Извлекаем эклиптическую долготу (lon) и широту (lat)
+    ecliptic_longitude = ecliptic_coords.lon
+    ecliptic_latitude = ecliptic_coords.lat
+
+    # Выводим результат
+    print("Техническая проверка:")
+    print(f"Эклиптическая долгота Солнца: {ecliptic_longitude}")
+    print(f"Эклиптическая широта Солнца: {ecliptic_latitude}")
+    return(ecliptic_longitude.value)
+
+
+#собираем за атмосферой все 4 компоненты
+#принимает на вход галактические координаты
+#принимает на вход дату и время наблюдения по UTC
+def total_background(l, b, date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_sp_fx = flux_newguey2003, V_wl = wavelenght_band_V, V_tr = trancparency_band_V, wavelenght_airglow = wavelenght_kp, intensity_airglow = intensity_kp, wavelenght_atmosphere = wavelenght_modtran_kp, transparency_atmosphere = trancparency_modtran_kp, venus_albedo_wl = venus_alb_wl, venus_albedo_rf = venus_alb_rf, mars_albedo_wl = mars_alb_wl, mars_albedo_rf = mars_alb_rf, jupiter_albedo_wl = jupiter_alb_wl, jupiter_albedo_rf = jupiter_alb_rf, saturn_albedo_wl = saturn_alb_wl, saturn_albedo_rf = saturn_alb_rf):
+    #получаем эклиптическую долготу Солнца
+    lmbd_Sun = Sun_ecl_lon(date, time)
+    #переводим галактические координаты в эклиптические геоцентрические
+    lmbd, beta = galactic_to_equatorial(l, b)
+    #получаем спектр зодиакального света (нм, фот / (м^2 сек нм ср))
+    zodiac_wl, zodiac_spec = zodiacal_spectrum(lmbd, beta, lmbd_Sun, Sun_sp_wl, Sun_sp_fx, V_wl, V_tr)
+#     return(zodiac_wl, zodiac_spec)
+
+    #собственное свечение достаем
+    airglow_wl, airglow_spec = airglow_spectrum(wavelenght_airglow, intensity_airglow, wavelenght_atmosphere, transparency_atmosphere)
+
+#     return(airglow_wl, airglow_spec)
+    #зод. свет и собств. свечение -- окей. Нужно допилить звездные каталоги и планеты.
+    
+    #звездные каталоги
+    star_cat_wl = zodiac_wl
+    star_cat_spec = star_spectrum(l, b, star_cat_wl)
+#     return star_cat_wl, star_cat_spec
+
+    #планеты
+    #Венера
+    venus_l, venus_b = coordinates_of_planet('venus', date, time)
+    print('Техническая проверка координат Венеры:', round(venus_l, 1), round(venus_b, 1))
+    if round(venus_l, 1) == l and round(venus_b, 1) == b:
+        venus_wl, venus_sp = venus_spectrum(date, time, Sun_sp_wl, Sun_sp_fx, venus_albedo_wl, venus_albedo_rf, V_wl, V_tr)
+    else:
+        venus_wl = zodiac_wl
+        venus_sp = np.zeros(venus_wl.shape[0])
+    
+    #Марс
+    mars_l, mars_b = coordinates_of_planet('mars', date, time)
+    print('Техническая проверка координат Марса:', round(mars_l, 1), round(mars_b, 1))
+    if round(mars_l, 1) == l and round(mars_b, 1) == b:
+        mars_wl, mars_sp = mars_spectrum(date, time, Sun_sp_wl, Sun_sp_fx, mars_albedo_wl, mars_albedo_rf, V_wl, V_tr)
+    else:
+        mars_wl = zodiac_wl
+        mars_sp = np.zeros(mars_wl.shape[0])
+        
+    #Юпитер
+    jupiter_l, jupiter_b = coordinates_of_planet('jupiter', date, time)
+    print('Техническая проверка координат Юпитера:', round(jupiter_l, 1), round(jupiter_b, 1))
+    if round(jupiter_l, 1) == l and round(jupiter_b, 1) == b:
+        jupiter_wl, jupiter_sp = jupiter_spectrum(date, time, Sun_sp_wl, Sun_sp_fx, jupiter_albedo_wl, jupiter_albedo_rf, V_wl, V_tr)
+    else:
+        jupiter_wl = zodiac_wl
+        jupiter_sp = np.zeros(jupiter_wl.shape[0])
+        
+    #Сатурн
+    saturn_l, saturn_b = coordinates_of_planet('saturn', date, time)
+    print('Техническая проверка координат Сатурна:', round(saturn_l, 1), round(saturn_b, 1))
+    if round(saturn_l, 1) == l and round(saturn_b, 1) == b:
+        saturn_wl, saturn_sp = saturn_spectrum(date, time, Sun_sp_wl, Sun_sp_fx, saturn_albedo_wl, saturn_albedo_rf, V_wl, V_tr)
+    else:
+        saturn_wl = zodiac_wl
+        saturn_sp = np.zeros(saturn_wl.shape[0])
+        
+    #теперь бы объединить эти спектры

{skymapmod-0.1.21 → skymapmod-0.1.22}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "hatchling.build"
 
 [project]
 name = "SkyMapMod"
-version = "0.1.21"
+version = "0.1.22"
 authors = [
  { name="Anna Amineva", email="amineva.2219@physics.msu.ru" },
 ]