def create_square(lat, lon, radius=0.0001):
    """
    Create the a geojson square polygon

    Args:
        lat: the center latitude for the polygon
        lon: the center longitude for the polygon
        radius (int): half of the length of the edge of the square

    Returns:
        list: a list of lat/lon points defining a square polygon
    """
    return [
        [round(lon + radius, 7), round(lat + radius, 7)],
        [round(lon + radius, 7), round(lat - radius, 7)],
        [round(lon - radius, 7), round(lat - radius, 7)],
        [round(lon - radius, 7), round(lat + radius, 7)]
    ]