# Color plots on surface of sphere

I would like to use Sage to draw colormap plots of spherical functions on the surface of spheres. I want to produce something like this example.

Are such plots currently possible in Sage, or using optional toolkits? I would like to have a color plots of arbitrary functions f(theta,phi).

Moreover, is it possible to have contour plots of such functions as well?

If it is possible using scipy and matplotlib, then you can use those directly in sage. I don't think it is possible by using some Sage function directly.