In the similar way we implement zonal_stats for variable geometry, we should add support to extract_points.