similar to: calculating distance inland from coastline

Does anyone have any experience plotting complex polygons on a scatterplot3d display? I would like to present some spatially binned data in a 3-d type plot with background polygons (i.e.. a filled coastline map). I can get the coastline onto the plot window as a line type but cannot figure out how to fill the shape. I'm certain it may be possible but might require a bit of digging

Hi R-Help, I would like to make filled contour maps of ocean data overlaid by costlines from the map package. I can draw the 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

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

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

I've been able to do most things I'd like to using the map package, but I can't figure out how to fill a region without drawing a black outline. I can draw the outline in a specified colour, but as soon as I set fill to TRUE, I get a black outline in addition. My map has a lot of activity on the coastline so I don't want a black line obscuring what I do there. A pale colour to

Strange things are going on in R, if you reshape a matrix in R: > g=gretldata[1:2,] > g Empfang Versand Transit Inland Ausland SumS 1 787844.0 1307176.6 223395.4 1474726 16199.1 3809341 2 421473.1 306445.4 448801.2 1779402 14445.6 2970567 > dim(g) [1] 2 6 > as.vector(g) Empfang Versand Transit Inland Ausland SumS 1 787844.0 1307176.6 223395.4 1474726 16199.1

Hi Sherri, There are examples of topographic maps which you have been pointed to, however, I suspect that you want to know where you can obtain topographic data from rather than a canned example. There are quite a few intricacies to the process so I will go through them for you. (1) Topography files can be found in the geomapdata library. You will probably want to use the maps package too (if

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

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:

Hi I'm using R Version 1.4.0 Patched (2002-01-16) in a SuSE 7.1 box. I've a plot with the european coastline. After I want to plot a single point in the map using points, but it plots this point in the wrong place. Why is this happening ? Is there a known problem of the points function ? Regards EJ > plot(eucoast,type="l") > points(c(-15.0,57.5)) ps: shall I send the

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

Very helpful. One (possible) correction inline below... On Thu, 30 Mar 2006, Ray Brownrigg wrote: > > From: nhy303 at abdn.ac.uk > > > > 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) > > > >

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

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",

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()

Hello, I would like to retrieve data stored in MySQL database, so I installed RMySQL package. I can successfully connect with the my database using the following code > dvr<-dbDriver("MySQL") > con2<-dbConnect(dvr,group="exbardiv") > mysqlDescribeConnection(con2) <MySQLConnection:(972,0)> User: mmorag Host: localhost Dbname: exbardiv

Hi! Is there a simple way for an * newbie to force * to use different sip-trunks for different calls. I have 2 siptrunks, one for inland calls and one for international calls. All in country numbers starts with 0 and all international starts with 00. This I have configured in the outbound routing. But * always use the incountry trunk because the 0. dialpattern is also true for international

Hi- I am trying to find a way to diffuse GIS data on a European map. I have a dataset 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)

Dear all, I am modeling the incidence of recreational anglers along a stretch of coastline, and with a vary large proportion of zeros (>80%) have chosen to use a zero inflated negative binomial (ZINB) distribution. I am using the same variables for both parts of the model, can anyone help me with R code to compute overall marginal effects of each variable? My model is specified as follows:

Hello, I am trying to draw geographical maps with the maps package. However if I try to access the data following error occurs: > require(maps) > map() Fehler in zip.file.extract(file, "Rdata.zip") : 'destination' existiert nicht # Error in zip.file.extract(file, "Rdata.zip") : 'destination' does not exist > traceback() 4: zip.file.extract(file,