#Figure 3 - having decided which to plot
#Fig 3 – Trees
shades = shading(breaks=c(0.9470, 0.9472, 0.9475, 0.9477),cols=brewer.pal(5,'OrRd')) # for T
if (.Platform$GUI == "AQUA") {
quartz(w=8,h=7) } else {
x11(w=8,h=7) }
par(mar=c(1,1,2,1))
level.plot(gwr.Bool.res.T,shades)
sub.masker = poly.outer(gwr.Bool.res.T, lib, extend = 0)
add.masking(sub.masker)
plot(lib, add = T)
level.plot(gwr.Bool.res.T,shades)与plot(lib, add = T)无重叠部分,本来应该有的,问题就在这