This function is adapted to plotting band structure in the case
where the translation symmetry vectors are not orthogonal
(e.g. graphene). It can only handle a very limited subset of
systems (2D lattices in 2D realspace), and most of the time
use of 'plotter.spectrum' should be preferred.