import math


def _find_utm_crs(lat, lon):
    """Find the UTM CRS based on lat/lon coordinates.

    Parameters
    ----------
    lat : float
        Decimal latitude.
    lon : float
        Decimal longitude.

    Returns
    -------
    crs : dict
        Corresponding UTM CRS.
    """
    utm_zone = (math.floor((lon + 180) // 6) % 60) + 1
    if lat >= 0:
        pole = 600
    else:
        pole = 700
    epsg = 32000 + pole + utm_zone
    return {'init': f'epsg:{epsg}'}