Sample all existing points of a polygon in Python

Question

I'm trying to sample raster values to points based on the geometries of a shapefile using the below code but the raster values are only sampled at the coordinates of the polygon geometry. Is there a way to sample for all the points along the boundary and inside the polygon.

multi_values_points = pd.Series()
 with fiona.open(inputShape) as shp:
     for feature in shp:
        Name = feature['properties']['Name']
        coords = feature['geometry']['coordinates']
        with rio.open(inputRaster) as stack_src:
           value = [v for v in stack_src.sample([coords])]
        multi_values_points.loc[Name] = value

Answer 1

I tried to generate points inside my various polygon boundaries of the shapefile by using below code:

import ogr
import random
from osgeo import ogr
from shapely.geometry import Point
from shapely.prepared import prep
from shapely.wkb import loads

g = ogr.Open(inputShapeFile)# Open the shapefile
for layer in g.GetLayer(0):
    field = loads(layer.GetGeometryRef().ExportToWkb())
    fieldP = prep(field)
    lon_min, lat_min, lon_max, lat_max = field.bounds
    n_points = 0 # A counter
    field_lonlat = []
    while n_points < 10:
        r_lon = random.uniform(lon_min, lon_max)
        r_lat = random.uniform(lat_min, lat_max)
        pt = Point(r_lon, r_lat)
        if fieldP.contains(pt):
            field_lonlat.append([r_lon, r_lat])
            n_points += 1