Compute x/y positions for phylopic silhouette overlays on a functional space plot.
compute_phylopic_positions <- function(scores, species, radius_frac = 0.25,