SkyMapMod 0.0.13__tar.gz → 0.0.14__tar.gz

{skymapmod-0.0.13 → skymapmod-0.0.14}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: SkyMapMod
-Version: 0.0.13
+Version: 0.0.14
 Project-URL: Homepage, https://istina.msu.ru/profile/aminevann/
 Author-email: Anna Amineva <amineva.2219@physics.msu.ru>
 License-File: LICENSE

{skymapmod-0.0.13 → skymapmod-0.0.14}/pyproject.toml

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

{skymapmod-0.0.13 → skymapmod-0.0.14}/src/SkyMapMod/__init__.py

@@ -1,4 +1,5 @@
 from .zodiac import *
 from .airglow import *
 from .airglow_spectrum import *
-from .modtran_default_kp_transparency import *
+from .modtran_default_kp_transparency import *
+from transparency import *

skymapmod-0.0.14/src/SkyMapMod/airglow.py

@@ -0,0 +1,19 @@
+#для работы с собственным свечением атмосферы. возможно, сюда можно досыпать поправку на радиопоток
+import numpy as np
+import math
+import matplotlib.pyplot as plt
+
+from .airglow_spectrum import wavelenght_kp, intensity_kp
+from .modtran_default_kp_transparency import wavelenght_modtran_kp, trancparency_modtran_kp
+
+def airglow(wavelenght_airglow = wavelenght_kp, intensity_airglow = intensity_kp, wavelenght_atmosphere = wavelenght_modtran_kp, transparency_atmosphere = trancparency_modtran_kp):
+    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]:
+                wavelenght.append(wavelenght_airglow[i])
+                intensity.append(intensity_airglow[i] / transparency_atmosphere[j])
+    return(np.array(wavelenght), np.array(intensity))
+#функция написана для спектра взятого от китт пик, прозрачность атмосферы взята из модтрана, параметры -- см. доклад, там прямо скрин. Прозрачность взята для 45 градусов, пересчитана в 30 градусов
+

skymapmod-0.0.14/src/SkyMapMod/transparency.py

@@ -0,0 +1,7 @@
+#для работы с прозрачностями атмосферы -- пересчет и тд
+import numpy as np
+from .zodiac import cos
+
+#на вход подаем спектр и два угла, на выходе получаем другой спектр
+def other_angle_transparency(initial_tr, initial_angle, angle):
+    return initial_tr ** (cos(initial_angle) / cos(angle))

skymapmod-0.0.14/src/SkyMapMod/transparency.py.ipynb

@@ -0,0 +1,6 @@
+{
+ "cells": [],
+ "metadata": {},
+ "nbformat": 4,
+ "nbformat_minor": 5
+}

{skymapmod-0.0.13 → skymapmod-0.0.14}/src/SkyMapMod/zodiac.py

@@ -2,6 +2,10 @@ import numpy as np
 import math
 import matplotlib.pyplot as plt
 
+#коэффициент пересчета из единиц S10 в количество фотонов от Солнца
+#вычислен на основе спектра АЧТ от звезды класса A0, солнечного спектра и факта, что 1 ед. S10 соответствует 100 фотонов от звезды класса A0
+k = 100.63910099562007
+
 #переопределяю тригонометрические функции, чтобы работать с градусами
 def cos(x):
     return math.cos(x * math.pi / 180)
@@ -134,3 +138,5 @@ def zodiacal_light(lmbd, beta, lmbd_sun):
     fourth = (8800 * exp((1 - sqrt(1 + ((abs_gm - 90)/3)**2))/10) - 1200) * exp(-eps/10)
     zod = first  + (second_1 + second_2) * second_mult + third + fourth + S(lmbd, beta, lmbd_sun, Omega) + D(lmbd, beta, lmbd_sun) + E(lmbd, beta) + F(lmbd, beta, lmbd_sun) + G(lmbd, beta, lmbd_sun)
     return zod/0.415 #деление на 0.415 осуществляет переход из ед. ADU в ед. S10
+
+

skymapmod-0.0.13/src/SkyMapMod/airglow.py

@@ -1,16 +0,0 @@
-import numpy as np
-import math
-import matplotlib.pyplot as plt
-
-from .airglow_spectrum import wavelenght_kp, intensity_kp
-from .modtran_default_kp_transparency import wavelenght_modtran_kp, trancparency_modtran_kp
-
-def airglow_custom(wavelenght_airglow = wavelenght_kp, intensity_airglow = intensity_kp, wavelenght_atmosphere = wavelenght_modtran_kp, transparency_atmosphere = trancparency_modtran_kp):
-    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]:
-                wavelenght.append(wavelenght_airglow[i])
-                intensity.append(intensity_airglow[i] / transparency_atmosphere[j])
-    return(np.array(wavelenght), np.array(intensity))