Assuming that your ground truth is in a shapefile, the basic idea is that each labelled pixel is matched with a corresponding embedding vector; your labels plus those vectors then beocome a standard supervised-learning problem. tessera-eval handles the details, fetching and mosaicking the Tessera tiles, sampling pixels from your polygons, and scoring models. It supports both classification (habitat class, land cover) and regression (canopy height, biomass, percent cover).

It supports:

• k-fold cross-validation — a robust accuracy estimate (macro-F1, or R²/RMSE/MAE for regression) with the fold-to-fold spread.
• Learning curves — accuracy versus training-set size, so you can see whether collecting more labels would still help.
• Confusion matrices, and predicted-vs-actual scatters for regression.
• Spatial train/test splits — hold out a whole region instead of random pixels, so nearby look-alike pixels can’t inflate the score.
• Spatial-context models — MLPs over a 3×3 or 5×5 neighbourhood, alongside per-pixel k-NN, random forest, XGBoost and MLP.

Here is a minimal run. austria_crops.geojson is 349 field parcels near Vienna labelled by crop type:

import geopandas as gpd
from geotessera import GeoTessera
from tessera_eval import load_embeddings_for_shapefile, run_kfold_cv

gdf = gpd.read_file("austria_crops.geojson").to_crs(4326)   # polygons with a "crop" column

# One embedding per 10 m pixel inside your polygons (downloads tiles on first run)
vectors, labels, class_names, stats = load_embeddings_for_shapefile(
    gdf, field="crop", year=2024, gt_instance=GeoTessera()
)
print(f"{len(labels):,} labelled pixels across {len(class_names)} classes")

# 5-fold cross-validation, random forest on the raw embeddings
for event in run_kfold_cv(vectors, labels, ["rf"], k=5, task="classification"):
    if event["type"] == "aggregate":
        m = event["models"]["rf"]
        print(f"macro-F1: {m['mean_f1']:.3f} ± {m['std_f1']:.3f}")

That prints macro-F1: 0.796 ± 0.004, a reasonable result for 10-way crop classification from embeddings alone. Replace [“rf”] with [“nn”, “rf”, “xgboost”, “mlp”] to compare all four per-pixel models at once. For a continuous target (canopy height, biomass, percent cover), pass task=“regression” with the regressor names — [“nn_reg”, “rf_reg”, “xgboost_reg”, “mlp_reg”] to get R²/RMSE/MAE instead of F1. One caveat: random k-fold on satellite data is optimistic because neighbouring pixels are correlated. For a realistic comparison, use a spatial split.

Install with pip install "tessera-eval[geotessera,xgboost]". A follow-up post will cover the command-line interface, which runs the same evaluations without writing any Python.

To find out more, here is the link to the git repo.