def select_trees(gdf, subplot):
    """
    Args:
        gdf: a geopandas dataframe
        subplot: a shapely box
    Returns:
        selected_trees: pandas dataframe of trees
    """
    selected_trees = gdf[gdf.intersects(subplot)]
    return selected_trees