f_plot_sp: Plot spatial data with overlap shpfile
In ln1267/dWaSSI: dWaSSI-C model

Plot spatial data with overlap shpfile

1
2
3

f_plot_sp(da, filename, colstyle = "RdYlGn", pretty = T,
  margin = list(), shpname = NA, varnames = NA, cuts = NA,
  ranges = NA, width = 7, height = 7, plot = T)

`da`	The input nc file.
`filename`	The input polygon.
`colstyle`	"RdYlGn",
`pretty`	Logic
`margin`	list or logic. This is the Statistic information for the data showed on the margin
`shpname`	It can be any sting if a shp will be ploted
`varnames`	This is the var names for the plot variable
`cuts`	This is how many segments of the plot
`ranges`	This is the min and max of the data
`width`	This is the width of the pdf in "inch"
`height`	This is the height of the pdf in "inch"

# the shpfile has to be loaded before the plot
library(maptools)
shp<-readShapeLines("data/shp/AU_STATES.shp")
f_plot_sp(da = WTD,
       filename = "/Dataset/www/images/WUE_AU/tt1.pdf",
       width=7,height=7,
       cuts = 10,
       ranges = c(3,100),
       shpname = "shp")