В настоящее время я создаю спектры с помощью pysynphot и знаком с командой sp (900), которая возвращает значение потока фотонов при 900 нм.
Я не писал никаких математических функций в python раньше, и мне любопытно, включил ли кто-нибудь интегральную функцию для команды потока? Я ищу, чтобы получить значение для интеграла между 900-1100 нм. Мой текущий код для скрипта:
import pysynphot as S
import matplotlib.pyplot as plt
from synphot import SourceSpectrum
from synphot.models import BlackBodyNorm1D
sp3 = S.Icat('ck04models', 13000, 0.0, 4.5)
plt.semilogy(sp3.wave, sp3.flux, 'r', label='B8')
plt.xlim(800, 2000)
plt.ylim(10e3, 10e11)
plt.xlabel(sp3.waveunits)
plt.ylabel(sp3.fluxunits)
plt.legend(loc='upper right')
sp(900)