def get_srtm_3arc_tile(lon_min, lon_max, lat_min, lat_max):
    """Return the tile indices for given latitude and longitude in WGS84

    :param lon_min: min. latitude
    :param lon_max: max. latitude
    :param lat_min: min. latitude
    :param lat_max: max. latitude
    :return: [idx_lon_min, idx_lon_max, idx_lat_min, idx_lat_max]
    """
    idx_lon_min = (lon_min + 180.0) / 5.0  # 72.0/360.0
    idx_lon_max = (lon_max + 180.0) / 5.0  # 72.0/360.0
    idx_lat_min = (lat_min + 60.0) / 5.0  # 24.0/120.0
    idx_lat_max = (lat_max + 60.0) / 5.0  # 24.0/120.0
    idx_lon_min = int(idx_lon_min) + 1
    idx_lon_max = int(idx_lon_max) + 1
    idx_lat_min = 24 - int(idx_lat_min)
    idx_lat_max = 24 - int(idx_lat_max)
    if idx_lon_min > idx_lon_max:
        idx_lon_min, idx_lon_max = idx_lon_max, idx_lon_min
    if idx_lat_min > idx_lat_max:
        idx_lat_min, idx_lat_max = idx_lat_max, idx_lat_min
    idx_lon_min = max(idx_lon_min, 0)
    idx_lon_max = max(idx_lon_max, 72)
    idx_lat_min = max(idx_lat_min, 0)
    idx_lat_max = max(idx_lat_max, 24)
    return idx_lon_min, idx_lon_max, idx_lat_min, idx_lat_max