def find_points(iso3, grid_sample, tile_lookup, filename):
    '\n\n    '
    filename_2 = (filename + '.shp')
    directory = os.path.join(DATA_INTERMEDIATE, iso3, 'sampling_points')
    path_output = os.path.join(directory, filename_2)
    if os.path.exists(path_output):
        return gpd.read_file(path_output, crs='epsg:4326')
    output = []
    for (idx, grid_tile) in grid_sample.iterrows():
        (min_x, min_y, max_x, max_y) = grid_tile['geometry'].bounds
        geom = Point(random.uniform(min_x, max_x), random.uniform(min_y, max_y))
        output.append({'type': 'Feature', 'geometry': geom, 'properties': {}})
    output = gpd.GeoDataFrame.from_features(output, crs='epsg:4326')
    return output
def find_points(iso3, grid_sample, tile_lookup, filename):
    '\n\n    '
    filename_2 = (filename + '.shp')
    directory = os.path.join(DATA_INTERMEDIATE, iso3, 'sampling_points')
    path_output = os.path.join(directory, filename_2)
    if os.path.exists(path_output):
        return gpd.read_file(path_output, crs='epsg:4326')
    output = []
    for (idx, grid_tile) in grid_sample.iterrows():
        (min_x, min_y, max_x, max_y) = grid_tile['geometry'].bounds
        geom = Point(random.uniform(min_x, max_x), random.uniform(min_y, max_y))
        output.append({'type': 'Feature', 'geometry': geom, 'properties': {}})
    output = gpd.GeoDataFrame.from_features(output, crs='epsg:4326')
    return output