Source code for mosqito.sq_metrics.loudness.utils.sone2phone

# -*- coding: utf-8 -*-
from numpy import log10




def sone2phone(sones):
    """Conversion between sones and phones, based on the work done by Hugo Fastl and Eberhard Zwicker in the article
    "Psychoacoustics: Facts and Models".

    Parameters
    ----------
    sones: float
        One of the units in which loudness is measured. 1 sone equals to 40 phones. The conversion from sones to phones
        is not linear.

    Returns
    -------
    phons: float
        Logarithmic unit of loudness (phon = Ln). The dB SPL values of the signal at 1 kHz have the same value in
        phons. For example, 50 db SPL @ 1 kHz corresponds to 50 phones. The rest of the frequencies does not have
        the same co-relation between each other.

    """

    if sones >= 1:
        phons = 40 + ((10 * log10(sones)) / log10(2))

    else:
        phons = 40 * ((sones + 0.0005) ** 0.35)
    return phons