Switch the logic of density interpolation to match that of current.
The density plot is optimised to show sites faithfully. However this produces a less aesthetically pleasing result for large systems (especially with non-square lattices, and systems where site positions coincide.
I suggest to change the logic of the density interpolation to closer match that of the current plotter (potentially even using the same source kernel so that divergence of the plotted current is the plotted source density).