def convert_to_LHC(imt):
    """
    Converts from GMRotI50 to Larger of two horizontal components using
    global equation of:
    Boore, D and Kishida, T (2016). Relations between some horizontal-
    component ground-motion intensity measures used in practice.
    Bulletin of the Seismological Society of America, 107(1), 334-343.
    doi:10.1785/0120160250
    No standard deviation modification required.
    """
    # get period t
    if isinstance(imt, SA):
        t = imt.period
    else:
        t = 0.01

    T1 = 0.08
    T2 = 0.56
    T3 = 4.40
    T4 = 8.70
    R1 = 1.106
    R2 = 1.158
    R3 = 1.178
    R4 = 1.241
    R5 = 1.241

    Ratio = max(R1,
                max(min(R1+(R2-R1)/np.log(T2/T1)*np.log(t/T1),
                        R2+(R3-R2)/np.log(T3/T2)*np.log(t/T2)),
                    min(R3+(R4-R3)/np.log(T4/T3)*np.log(t/T3), R5)))
    SF = np.log(Ratio)

    return SF
Converts from GMRotI50 to Larger of two horizontal components using
    global equation of:
    Boore, D and Kishida, T (2016). Relations between some horizontal-
    component ground-motion intensity measures used in practice.
    Bulletin of the Seismological Society of America, 107(1), 334-343.
    doi:10.1785/0120160250
    No standard deviation modification required.