[R] Using map with filled.contour doesn't display map

Jeff Charlton jcharlto16 at gmail.com
Fri Nov 4 19:26:00 CET 2016


I'm trying to overlay a map on top of data showing temperatures across the
world. The code I'm using is:

      filledContour(tempdata, plot.axes={axis(1); axis(2); map("world2",

where tempdata is a raster file made from a netcdf file downloaded from the
NOAA website. The filledContour is a wrapper that translates the raster
image into something that can be used by filled.contour.

If I run the code:


I get the image I want, but when I add the plot.axes parameter, I get a
"plot.new has not been called yet" error in axis(1). If I run plot.new()
and the original filledContour call, I get the same result as if I just did
filledContour(tempdata) (i.e. no map is overlayed).

The kicker in all this is that this code worked a couple of years ago. I
ran into this error when I recently dusted of my code and installed the
latest version of R and its various packages. I have a feeling that
something has changed in the maps package but I can't find any reference to
it online.

Can anyone help?

Thank you,


Just in case, here's my session info:

R version 3.3.2 (2016-10-31) Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows 7 x64 (build 7601) Service Pack 1

locale: [1] LC_COLLATE=English_United States.1252 LC_CTYPE=English_United
States.1252 LC_MONETARY=English_United States.1252 [4] LC_NUMERIC=C
LC_TIME=English_United States.1252

attached base packages: [1] stats graphics grDevices utils datasets methods

other attached packages: [1] ncdf4_1.15 maps_3.1.1 XML_3.98-1.4
RCurl_1.95-4.8 bitops_1.0-6 raster_2.5-8 sp_1.2-3

loaded via a namespace (and not attached): [1] parallel_3.3.2 tools_3.3.2
Rcpp_0.12.7 grid_3.3.2 lattice_0.20-34

	[[alternative HTML version deleted]]

More information about the R-help mailing list