diff --git a/simpa/utils/libraries/spectrum_library.py b/simpa/utils/libraries/spectrum_library.py index 6a1af7ee..be1cdddd 100644 --- a/simpa/utils/libraries/spectrum_library.py +++ b/simpa/utils/libraries/spectrum_library.py @@ -247,7 +247,7 @@ def scattering_from_rayleigh_and_mie_theory(name: str, mus_at_500_nm: float = 1. :return: A Spectrum instance based on Rayleigh and Mie scattering theory. """ wavelengths = np.arange(450, 1001, 1) - scattering = (mus_at_500_nm * (fraction_rayleigh_scattering * (wavelengths / 500) ** 1e-4 + + scattering = (mus_at_500_nm * (fraction_rayleigh_scattering * (wavelengths / 500) ** (-4) + (1 - fraction_rayleigh_scattering) * (wavelengths / 500) ** -mie_power_law_coefficient)) return Spectrum(name, wavelengths, scattering)