SkyMapMod 0.1.39__tar.gz → 0.1.41__tar.gz

{skymapmod-0.1.39 → skymapmod-0.1.41}/PKG-INFO

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

{skymapmod-0.1.39 → skymapmod-0.1.41}/SkyMapMod/airglow.py

@@ -23,7 +23,7 @@ def radioflux(date, time): #date -- в формате строки 'дд.мм.г
     return float(fluxobsflux[np.argmin(box)])
 
 
-def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = intensity_kp, wavelenght_atmosphere = wavelenght_modtran_kp, transparency_atmosphere = trancparency_modtran_kp):
+def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = intensity_kp, wavelenght_atmosphere_kp = wavelenght_modtran_kp, transparency_atmosphere_kp = trancparency_modtran_kp):
     """
     Возвращает длины волн и спектр собственного свечения за атмосферой.
 
@@ -35,9 +35,9 @@ def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = int
             По умолчанию — wavelenght_kp.
         intensity_airglow (np.ndarray, optional): Интенсивность собственного свечения, Рэлей/Ангстрем.
             По умолчанию — intensity_kp.
-        wavelenght_atmosphere (np.ndarray, optional): Длины волн атмосферных данных, нм.
+        wavelenght_atmosphere_kp (np.ndarray, optional): Длины волн атмосферных данных, нм.
             По умолчанию — wavelenght_modtran_kp.
-        transparency_atmosphere (np.ndarray, optional): Прозрачность атмосферы (от 0 до 1).
+        transparency_atmosphere_kp (np.ndarray, optional): Прозрачность атмосферы (от 0 до 1).
             По умолчанию — transparency_modtran_kp.
 
     Возвращает:
@@ -53,10 +53,10 @@ def airglow_spectrum(wavelenght_airglow = wavelenght_kp, intensity_airglow = int
     wavelenght = []
     intensity = []
     for i in range(wavelenght_airglow.shape[0]):
-        for j in range(wavelenght_atmosphere.shape[0]):
-            if wavelenght_airglow[i] == wavelenght_atmosphere[j]:
+        for j in range(wavelenght_atmosphere_kp.shape[0]):
+            if wavelenght_airglow[i] == wavelenght_atmosphere_kp[j]:
                 wavelenght.append(wavelenght_airglow[i])
-                intensity.append(intensity_airglow[i] / transparency_atmosphere[j])
+                intensity.append(intensity_airglow[i] / transparency_atmosphere_kp[j])
     intensity = np.array(intensity)
     intensity = intensity * 10**11 / (4 * np.pi) * np.pi**2 / (1800 * 1800)
     wavelenght = np.array(wavelenght)

{skymapmod-0.1.39 → skymapmod-0.1.41}/SkyMapMod/total_background.py

@@ -7,6 +7,7 @@ from astropy import 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 *
@@ -46,10 +47,7 @@ def galactic_to_equatorial(l, b):
         - `astropy.coordinates.Galactic`: Для представления галактических координат.
         - `astropy.coordinates.GeocentricMeanEcliptic`: Для представления эклиптических координат.
     """
-#     import astropy.units as u
-
-
-
+    import astropy.units as u
 # Входные данные: галактические координаты и эпоха
     galactic_l = l * u.deg  # Галактическая долгота
     galactic_b = b * u.deg # Галактическая широта
@@ -151,7 +149,7 @@ def sum_of_spectrums(array_1, meaning_1, array_2, meaning_2):
 #собираем за атмосферой все 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):
+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_kp = wavelenght_modtran_kp, transparency_atmosphere_kp = 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):
     """
     Вычисляет суммарный спектр фонового излучения для заданных галактических координат, даты и времени.
 
@@ -178,9 +176,9 @@ def total_background(l, b, date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_s
                                                    По умолчанию используется `wavelenght_kp`.
         intensity_airglow (ndarray, optional): Массив интенсивностей собственного свечения атмосферы.
                                                   По умолчанию используется `intensity_kp`.
-        wavelenght_atmosphere (ndarray, optional): Массив длин волн для модели атмосферы (в нм).
+        wavelenght_atmosphere_kp (ndarray, optional): Массив длин волн для модели атмосферы (в нм).
                                                       По умолчанию используется `wavelenght_modtran_kp`.
-        transparency_atmosphere (ndarray, optional): Массив значений прозрачности атмосферы.
+        transparency_atmosphere_kp (ndarray, optional): Массив значений прозрачности атмосферы.
                                                         По умолчанию используется `trancparency_modtran_kp`.
         venus_albedo_wl (ndarray, optional): Массив длин волн спектра альбедо Венеры (в нм).
                                                 По умолчанию используется `venus_alb_wl`.
@@ -230,7 +228,7 @@ def total_background(l, b, date, time, Sun_sp_wl = wavelenght_newguey2003, Sun_s
 #     return(zodiac_wl, zodiac_spec)
 
     #собственное свечение достаем
-    airglow_wl, airglow_spec = airglow_spectrum(wavelenght_airglow, intensity_airglow, wavelenght_atmosphere, transparency_atmosphere)
+    airglow_wl, airglow_spec = airglow_spectrum(wavelenght_airglow, intensity_airglow, wavelenght_atmosphere_kp, transparency_atmosphere_kp)
     print('ФОТОНОВ СОБСТВЕННОГО:', integral(airglow_wl, airglow_spec))
 
 #     return(airglow_wl, airglow_spec)
@@ -360,7 +358,7 @@ def transparency(h, zenith_angle): #на выходе дает расчитан
     
     return(wl, tr)
 
-def total_background_on_Earth(lat_obs, lon_obs, h, zenith_angle, azimuth_angle, date, time, wl_tr = [None], sp_tr = [None], 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):
+def total_background_on_Earth(lat_obs, lon_obs, h, zenith_angle, azimuth_angle, date, time, wl_tr = [None], sp_tr = [None], 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_kp = wavelenght_modtran_kp, transparency_atmosphere_kp = 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):
     """
     Вычисляет суммарный спектр фонового излучения для заданных условий наблюдения на Земле
     для квадрата размером 0.1 x 0.1 градуса по галактическим координатам.
@@ -375,11 +373,11 @@ def total_background_on_Earth(lat_obs, lon_obs, h, zenith_angle, azimuth_angle,
     для квадрата размером 0.1 x 0.1 градуса по галактическим координатам.
 
     Параметры:
-        lat_obs (float): Широта наблюдателя (в градусах).
-        lon_obs (float): Долгота наблюдателя (в градусах).
+        lat_obs (float): Широта наблюдателя (в градусах). Должна находиться в диапазоне [-90, 90]
+        lon_obs (float): Долгота наблюдателя (в градусах). Должна находиться в диапазоне [0, 360]
         h (float): Высота наблюдателя над уровнем моря (в метрах).
-        zenith_angle (float): Зенитный угол направления наблюдения (в градусах).
-        azimuth_angle (float): Азимутальный угол направления наблюдения (в градусах, от севера по часовой стрелке).
+        zenith_angle (float): Зенитный угол направления наблюдения (в градусах). Должен находиться в диапазоне [0, 90]
+        azimuth_angle (float): Азимутальный угол направления наблюдения (в градусах, от севера по часовой стрелке). Должен находиться в диапазоне [0, 360]
         date (str): Дата наблюдения в формате 'YYYY-MM-DD' (например, '2023-11-03').
         time (str): Время наблюдения в формате 'HH:MM:SS' (например, '12:00:00').
         wl_tr (ndarray, optional): Массив длин волн для прозрачности атмосферы (в нм). 
@@ -398,9 +396,9 @@ def total_background_on_Earth(lat_obs, lon_obs, h, zenith_angle, azimuth_angle,
                                                    По умолчанию используется `wavelenght_kp`.
         intensity_airglow (ndarray, optional): Массив интенсивностей собственного свечения атмосферы.
                                                   По умолчанию используется `intensity_kp`.
-        wavelenght_atmosphere (ndarray, optional): Массив длин волн для модели атмосферы (в нм).
+        wavelenght_atmosphere_kp (ndarray, optional): Массив длин волн для модели атмосферы (в нм).
                                                       По умолчанию используется `wavelenght_modtran_kp`.
-        transparency_atmosphere (ndarray, optional): Массив значений прозрачности атмосферы.
+        transparency_atmosphere_kp (ndarray, optional): Массив значений прозрачности атмосферы.
                                                         По умолчанию используется `trancparency_modtran_kp`.
         venus_albedo_wl (ndarray, optional): Массив длин волн спектра альбедо Венеры (в нм).
                                                 По умолчанию используется `venus_alb_wl`.
@@ -443,6 +441,7 @@ def total_background_on_Earth(lat_obs, lon_obs, h, zenith_angle, azimuth_angle,
     if wl_tr[0] == None and sp_tr[0] == None:
         wl_tr, sp_tr = transparency(h, zenith_angle)
     
+    import astropy.units as u
     lat_obs = lat_obs * u.deg  # Широта наблюдателя
     lon_obs = lon_obs * u.deg  # Долгота наблюдателя
     h = h * u.m  # Высота над уровнем моря
@@ -469,8 +468,9 @@ def total_background_on_Earth(lat_obs, lon_obs, h, zenith_angle, azimuth_angle,
 
     l = round(galactic_coords.l.value, 1)
     b = round(galactic_coords.b.value, 1)
+    print('Галактические координаты: ', l, b)
 
-    tot_gal_wl, tot_gal_sp = total_background(l, b, date, time, Sun_sp_wl, Sun_sp_fx, V_wl, V_tr, wavelenght_airglow, intensity_airglow, wavelenght_atmosphere, transparency_atmosphere, venus_albedo_wl, venus_albedo_rf, mars_albedo_wl, mars_albedo_rf, jupiter_albedo_wl, jupiter_albedo_rf, saturn_albedo_wl, saturn_albedo_rf)
+    tot_gal_wl, tot_gal_sp = total_background(l, b, date, time, Sun_sp_wl, Sun_sp_fx, V_wl, V_tr, wavelenght_airglow, intensity_airglow, wavelenght_atmosphere_kp, transparency_atmosphere_kp, venus_albedo_wl, venus_albedo_rf, mars_albedo_wl, mars_albedo_rf, jupiter_albedo_wl, jupiter_albedo_rf, saturn_albedo_wl, saturn_albedo_rf)
     
     result_wl, result_sp = convolution(tot_gal_wl, tot_gal_sp, wl_tr, sp_tr)
     

{skymapmod-0.1.39 → skymapmod-0.1.41}/pyproject.toml

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