search for: worldhires

...he filled contours and the coastlines om the same plot, but the filled contour also covers part of the land. To get rid of that I tried to draw a filled coastline map on top of the filled contour, but the filled map only draws the closed contours - so most of the land is missing. Example: map("worldHires",xlim=c(0,15),ylim=c(50,60)) #Draw relevant region (North Sea and Denmark waters) map("worldHires",xlim=c(0,15),ylim=c(50,60),fill=TRUE) # This only draws Denmark and Holland (Sweden, uk, Germany etc. disappears because they are not closed polygons). Any hint on how to fix this prob...

Hi R-help I am trying to mask the ocean from an image plot I have made. Here is some example code: library(mapdata) image(x=110:155, y =-40:-10, z = outer(1:45, 1:30, "+"), xlab = "lon", ylab = "lat") outline <- map("worldHires", plot=FALSE) # returns a list of x/y coords xrange <- range(outline$x, na.rm=TRUE) # get bounding box yrange <- range(outline$y, na.rm=TRUE) xbox <- xrange + c(-2, 2) ybox <- yrange + c(-2, 2) # create the grid path in the current device polypath(c(outline$x, NA, c(xbox, rev...

...body know of any existing functions that will calculate distance inland from a coastline? It's possible to test if a lon,lat location is land or sea using map.where(), but I need to add a buffer to this of say 2km, to allow for points that are just on the coast, and below the resolution of the worldHires database. I'm working with a marine mammal satellite telemetry dataset and wish to filter out spurious locations on land. On another issue, does anybody know of any free vector map datasets that are more up to date and have a higher resolution than the worldHires database from the 'maps...

I have just learned how to play with map, but something weird (or not) is happening. Suppose I want to draw a map of two countries (that have disconnected components), like Argentina and Brazil. If I command: library(maps) library(mapdata) map("worldHires", c("Argentina", "Brazil")) It works fine. However, if I want to _colour_ the interior: map("worldHires", c("Argentina", "Brazil"), c("cyan", "green"), fill=T) Then the colors will be assigned to the islands (Marajo in Br...

R users, In S, there was a function called usa() that would draw the map of the United States, plus it had other options for graphics. I have looked but I can't find the equivalent in R. Is there one? Thanks, Jason

Hello, The other day Justin Peter presented a mini program to plot a topographic map with an overlay of the worldHires. I seemed interesting so I checked the ETOPO5 site and find that there is a new file ETOPO1 with a 1 minute grid. I downloaded it and tried a similar procedure. Now the ETOPO1.gz is 1 Gb and the uncompressed file is 5 Gb. They do not fit into my laptop. I tried the following program that works, but...

Dear all, As a newbie in R I would like to do the following (simple?) thing: to plot a filled.contour plot over a map showing country boundaries (e.g. for Europe) What i do is: map('worldHires',xlim=c(-10,40),ylim=c(35,70),boundary = TRUE,border=0.1) map.axes() filled.contour(mslp, zlim=c(1000,1020),color.palette = colorRampPalette(c("blue", "white", "red")),main="Avegared MLSP (hPa) ERA40 JJA [1996-2002]", xlab="Longitude",ylab=&quot...

Dear all, I would like to employ a European map in a presentation. My idea was to use: library(mapdata) map("worldHires", c("Austria", "Switzerland", "Germany")) where I include all countries from my analysis as a vector of character strings like in the example above. Unfortunately, I was unable to specify the Baltic States (Latvia, Lithuania, Estonia) or the Czech Republic. Init...

dear all: Im trying to obtain maps on R, under mapdata library, but I cannot define the labels. Here is an example: library(mapdata) map("worldHires", c("portugal","spain"),ylim=c(34,46),xlim=c(-14,3.5)) axis(1,at=seq(-12,3,3)) axis(2) when I try the parameter xlab or ylab in axis, i get: parameter "ylab" couldn't be set in high-level plot() function  any suggestion?  thanks in advance for your always u...

...re out how to fill the shape. I'm certain it may be possible but might require a bit of digging into the scatterplot3d code and rotational method. Any and all help or suggestions appreciated. EX. library(scatterplot3d) # get a coastline library(maps) library(mapdata) coast<-map('worldHires', c('usa','canada'),xlim=c(-72,-67), ylim=c(40,45),fill=TRUE, col='gray') # note: substitute a simple polygon if you don't have maps and mapdata # coast<-list() # coast$x<-c(-70.5,-70.5,-70,-70) # coast$y<-c(41.2,43,43,41.2) lat<-c(41,41.5,42,42.5,43)...

...consisting of particular locations scattered across Europe, along with magnitude and value information. I can plot these as discrete points with something like the following: "geocode" is a dataframe with four columns: LAT; LONG; MAGNITUDE;VALUE. library(maps) library(mapdata) map("worldHires", regions=c("Germany", "Belgium", "Netherlands")) points(geocode$LONG, geocode$LAT, cex=geocode$MAGNITUDE / 2500, col=rainbow(length(geocode$VALUE), start=0, end=.4)[rank(geocode$VALUE)]) This gives me a map of Europe with my datapoints highlighted in two ways: m...

...le.extract(file, "Rdata.zip") : 'destination' does not exist > traceback() 4: zip.file.extract(file, "Rdata.zip") 3: data(list = dbname) 2: maptype(database) 1: map() To check if I have principle problems with my system I tried: > require(mapdata) > map('worldHires', col=1:10) # Works I am running R 2.6.0 and a current version of maps. But I had a similar problem with R 2.5.1 Is this a known issue? thanks Lukas

Would anyone with experience with the map functions know how to divide Czechoslovakia into the Czech Republic and Slovakia. They have been two separate countries for some time now. I'm thinking about the worldhires map database in particular. -- Dean Sonneborn, MS Programmer Analyst Department of Public Health Sciences University of California, Davis (530) 754-9516

Dear All, This is probably a very basic question but: I have plotted a map of the Barents Sea and surrounding coastline using: map('worldHires',ylim=c(50,85),xlim=c(5,65),fill=T,resolution=0) map.axes() map.scale(x=30,metric=T) Next, I imported a shapefile with depth contours for the sea: contours<-read.shape("D://My Documents/BarentsSea.shp",dbf.data=T) (This is in mercator projection). Despite extensive searches of...

I've generated a levelplot showing the density distribution of a species derived from survey transects, with lon, lat co-ordinates. I'd like to overlay this on a map of the study region specified by: map('worldHires', xlim = range(mlon), ylim = range(mlat)), where mlon, mlat specifies the study region. I've tried plotting the map first, then using the options: 'add=TRUE, axes=FALSE' in the levelplot, but this doesn't work - a new plot is generated for the level plot despite the 'add...

...quot;#7800FF","#B400FF","#FF0096") zprcp.quants <- rev(quantile(zprcp,na.rm=T,probs=c(1,0.98,0.9,0.8,0.7,0.6,0.5,0.4,0.3,0.2,0.1))) zprcp.breaks <-c(0,10*ceiling(zprcp.quants/10)) filled.contour(xlon,ylat,zprcp,ylim=c(20,50),xlim=c(10,60), asp=1.0, plot.axes=map('worldHires',xlim=c(10,60),ylim=c(20,50), border=0.9, add =TRUE),levels=zprcp.breaks,col=zprcp.colour,key.axes = axis(4,zprcp.breaks)) I would like the colour bars to be even (and the labels to represent the actual quantile values). I tried to modify the key.axes=axis(..) to force an evenly spaced colou...

...n’t look good. Any ideas! Bellow is my plot code without the above mentioned attempts. I’d like to add this effect at the second text() command. png('BaseMap.png',width=5.2,height=3.6,pointsize=5,units='in',bg=NA,restoreC onsole=T,res=1200) par(mar=c(5,0,0,0)) map('worldHires','brazil',ylim=c(-5.15,-4.55),xlim=c(-37,-36.1),type='n') rect(-37.1,-5.25,-36,-4.9,density=NULL,angle=45,col='#dbf4ff',border=F) plot(Batimetria,ylab='',xlab='',border=F,col=Cor$Cor,add=T) plot(Municipios,ylab='',xlab='',col=rgb(245,23...

...=c(-47,-7)) There are plenty of other options, if you have your own map data (or want to use another source). Feel free to provide more detail about your current plotting methods and coordinate system. Also, the package "mapdata" contains a high resolution continental dataset -"worldHires" Hope that helps, Mike. ############################################### Michael Sumner - PhD. candidate Maths and Physics (ACE CRC & IASOS) and Zoology (AWRU) University of Tasmania Private Bag 77, Hobart, Tas 7001, Australia Phone: 6226 1752

...('Pontos.xls',sheet=1,rowNames=T) Prof<-read.xls('Pontos.xls',sheet=2,rowNames=T) Batimetria<-readShapeSpatial(‘Batimetria.shp') Municipios<-readShapeSpatial('Municipios.shp') Recifes<-readShapeSpatial('Recifes.shp') par(mar=c(4,0,0,0)) map('worldHires','brazil',ylim=c(,),xlim=c(,),type='n') plot(Batimetria,ylab='',xlab='',border='lightgray',add=T) plot(Recifes,ylab='',xlab='',lty='dashed',add=T) plot(Municipios,ylab='',xlab='',col=rgb(229,229,229,max=255),add=...

...’m trying to substitute the points() functions with one bubble() function. for (i in 1:length(colnames(Spp10ml))) { png(paste(colnames(Spp10ml),'.png',sep='')[i],width=5.2,height=3.6,pointsize =5,units='in',bg=NA,restoreConsole=T,res=1200) par(mar=c(5,0,0,0)) map('worldHires','brazil',ylim=c(-5.15,-4.55),xlim=c(-37,-36.1),type='n') rect(-37.1,-5.25,-36,-4.9,density=NULL,angle=45,col='#dbf4ff',border=F) plot(Batimetria,ylab='',xlab='',border=F,col=Cor$Cor,add=T) plot(Municipios,ylab='',xlab='',col=rgb(245,23...