Very similar to this question
except instead of integers we have multidimensional feature vectors with an arbitrary distance function for measuring their similarity.
The person who posted that question claimed that they had found a polynomial time solution to their problem that would appear to also work for my problem, but I couldn't understand their answer and it wasn't clear if it was actually correct.
My problem is also equivalent to this question with the weights of the graph equal to the similarity function's outputs.
max-weight k-clique in a complete k-partite graph
But it's unclear if the k-partiteness makes this problem any easier than the general max-weight clique problem.
edit:
From David Eisenstat's answer, it looks the original question isn't solvable in polynomial time. But does anything change if the objective is minimize the maximum pairwise distance between the selected points?