What is a straightforward way to produce a heat map with the heat map strictly constrained within but filling shape file borders? If someone can illustrate an example in R I would be most grateful. Best wishes Graham