From d3fdfa587060d1749cfe9d75b88a5b3146acafc2 Mon Sep 17 00:00:00 2001 From: Thomas Else <56722586+tomelse@users.noreply.github.com> Date: Mon, 20 Jan 2025 14:31:22 +0000 Subject: [PATCH] Correct the formula for Rayleigh scattering --- simpa/utils/libraries/spectrum_library.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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)