SkyMapMod 0.1.39__py3-none-any.whl → 0.1.41__py3-none-any.whl

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/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.dist-info → skymapmod-0.1.41.dist-info}/METADATA

@@ -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.dist-info → skymapmod-0.1.41.dist-info}/RECORD

@@ -1,5 +1,5 @@
 SkyMapMod/__init__.py,sha256=4NTV_OA0-1cPJVcsHxxPWWu5TjnzJxGvVD4Eohst0Bc,331
-SkyMapMod/airglow.py,sha256=tE2Te7WsGbpJtHQLCdzs3-6MoayVIbAqhXb8244BFbo,4473
+SkyMapMod/airglow.py,sha256=Z51b_N-BWUjuPVY2o_pXrFGvPdqSa8LLELiWjeMVyoo,4494
 SkyMapMod/airglow_spectrum.py,sha256=4pY4QSGjIDvjs_jvxrowsZ9WLtIUwR0DesoCK6moWjs,127314
 SkyMapMod/albedo_of_planets.py,sha256=EmK8jov_cprHFc71irm1hzeL1mkQs7MmNkKRWBZSrDM,66432
 SkyMapMod/band_V_data.py,sha256=fOj9H2IdZuqErAxnoLq4A02UXfYMgLNBiLjwrzjfWTs,15223
@@ -8,14 +8,14 @@ SkyMapMod/planets.py,sha256=KdStIERhzu0ThplaclKnHSFiJ02WCO_XnSWJBEZGGQA,28483
 SkyMapMod/solar_radio_flux.py,sha256=I02O6Mltv3lfqxif3BHUNAzeElQ5AZax_2L09A3cXb8,842
 SkyMapMod/solar_spectrum.py,sha256=LjBPSTI6Yq86pwF7haUja7RjRYH4C8HZL8SjyI7zp6g,24433
 SkyMapMod/star_catalogues.py,sha256=-JGHE3m6NCg4Ea57ZR6RuxfftkPJyjJHOaU5Wf7qvRc,4102
-SkyMapMod/total_background.py,sha256=Qi-5RnUZEwABdUo--P6Z5JKJSUhaJ8hras0pa8A3WHg,35291
+SkyMapMod/total_background.py,sha256=B-UERqAljsND31xX8LqvBYnhciSlXQlMXSC7dCE8LNc,35709
 SkyMapMod/transparency.py,sha256=HiSbbVX80neh3n284f64sWnrbxpRzPGKvFyid9sVo7I,376
 SkyMapMod/zodiac.py,sha256=mT08Cjh2qcL3FVEM0B4G52OeuzWPB9w9mpp89sWRVJo,40729
 SkyMapMod/data/__init__.py,sha256=_292qgM72hPa2vxDNSaI9QvPunvi-aSXUsZki99q63o,83
 SkyMapMod/data/load_data.py,sha256=p-Nq5St713HTTDunSWu9ANrF9V0Slaak17IhD0n_bSY,931
 SkyMapMod/data/star_brightness.npy,sha256=aTsrvGU-vl9pqcOLMz3kqsSA15piynS1aP19aqjT6fQ,51840128
 SkyMapMod/data/star_temperatures.npy,sha256=Yjp252bzqmKQz00K0ADTkstjzGFulRMgyW1b8kaDOPE,51840128
-skymapmod-0.1.39.dist-info/METADATA,sha256=TmUo6BIKnBR8z2v66IRqdl820UV_YwxsDBW81B0mXME,469
-skymapmod-0.1.39.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
-skymapmod-0.1.39.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
-skymapmod-0.1.39.dist-info/RECORD,,
+skymapmod-0.1.41.dist-info/METADATA,sha256=oLK1fH0YNW8I7qyUwDmAXu52XzAv9E1Ek2IPTP3_3GQ,469
+skymapmod-0.1.41.dist-info/WHEEL,sha256=qtCwoSJWgHk21S1Kb4ihdzI2rlJ1ZKaIurTj_ngOhyQ,87
+skymapmod-0.1.41.dist-info/licenses/LICENSE,sha256=eShhW1a77wkzEW9UASG0IaNNaJ3vnxmpcxp8tXh5WCA,1045
+skymapmod-0.1.41.dist-info/RECORD,,