map.stippling does not work with station data

[rmanzanas]

The map.stippling function gives an error when applied on station data (it works well for gridded data). This would be an example:

pval = climatology(VALUE_Iberia_tas, clim.fun = list(FUN = "trend.1D", dates = getRefDates(VALUE_Iberia_tas), method = "kendall", return.pvalue = T))

sig = map.stippling(clim = pval, threshold = 0.05, condition = "LT", pch = 19, cex = .1, col = "black") # points exhibiting significant trends

[jbedia]

Hi, the stippling is not in principle intended to work with point data. You can use the index of significant/non-significant stations to highlight them at your convenience. Is it the best way to overlay one point layer on top of another? It can be easily implemented if this is what you want

[rmanzanas]

Hi, @jbedia. Perhaps implementing the function in a way that allows working indistinctly with gridded and point-wise data would be ideal in terms of consistency (just a suggestion, this is nothing critical...).