similar to: Do you know if you can map a large minmum spanning tree in R?

Thanks Mike. My data has longitude and latitude coords and I used distAB {clim.pact} then mst {ape} to calculate my minimum spanning tree. The nodes are telecoms sites from all over Australia. My goal is to determine the minimum cost of linking them via cabling, and I'm starting by calculating the distance "as the crow flies", but will probably eventually need to calculate the

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

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

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

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

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

Are there any functions available to calculate sunrise and sunset times for given latitude,longitude and dates? If not, I'll appreciatte any pointers to C code I could use/port. Thanks, Angel

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:

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

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

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,

Does anyone have suggestions about how to produce a map of British Colonial America 1775? At that time, Great Britain had 26 colonies in the Americas, including Bermuda, several Caribbean islands, "Quebec" (extending then almost to New Orleans), Nova Scotia, Newfoundland, Belize, and the 13 that declared independence in 1776. I've reviewed the "map"

Hi list, Well, this time I’ve a doubt with mapping generation. I was already able to read and plot shapefiles, plot point on this map. All this without any problems, but now I want to do something that I think, only Golden Software Surfer is capable of. I would like to plot a kriging result on the existing map (map script below). Well I looked for on the net, but I was not able to be sure of

Hello everybody, I'm an absolute newbe in R, so please be gentle to me. I am looking for a german map in the package maps by Becker and Wilks. After that I want to draw the German cities from world.cities of this package in the map and want to draw further a line from Berlin to Munich. Has anybody an idea? Thank you so much and sorry for the question! Best regards from Berlin, Maja!

Good morning, I am currently displaying a time series with the time on the Y axis, and the values on the X axis. The problem is that I wanted to reverse the Y axis (that is, to have the latest date the closest to the X axis). So I looked for a way to do this on this mailing list and I found a response: inverting the ylim parameter of the plot. I did it, the data are displayed correctly, but now

Hi list. I would like to know how do I insert a North arrow, stored as a png file in my computer, in a map? I found lots of post asking similar things, one of them mentioned the pixmap package. The map was done using map() and shapefiles (the code is below). I’m using the pixmap () and addlogo() functions. Well I can import the png with pixmap() function (I guess, once there’s no error message),

Hi R-user, I am new with R and I have a problem with the map and mapproj fonctions : with the following code : >test_proj=mapproject(LONG_d01_vec_t1, LAT_d01_vec_t1, projection="lambert", parameters=c(30,60)) >longitude_vec=test_proj$x >latitude_vec=test_proj$y >longitude_mat=matrix(longitude_vec,ncol=e_sn_d01) >latitude_mat=matrix(latitude_vec,ncol=e_sn_d01)

Dear everyone, as a newbie I fear I have made an easy mistake and am just too blind to see it. Can anyone help me seeing the (possibly) obvious? I have a command that creates a map of South America. In a separate txt-file I stored data about latitude and longitude of languages, and in a third column values that differ from language to language. The command should read that file and assign the

Does anyone know how to do hermite or Bezier splines in R? I can find Matlab routines, but really need to implement them in R. Failing that I'd be interested in other conformal splines. I need to smooth and interpolate animal tracking data. Antarctic Wildlife Research Unit School of Zoology University of Tasmania PO Box 252-05 Hobart TAS, 7001 Australia Phone: (0)3 6226 2645 Fax: