library(terra)
data(africa)
p <- vect(system.file("ex/sa.json", package = "phyloregion"))
K <- ncol(africa$omega)
CLRS <- hcl.colors(K)
plot_spatial_membership(africa$omega, pol = p, col=CLRS, type="blend")
Run the code above in your browser using DataLab