similar to: Concatening two maps/shapefiles...

Dear R'helpers, I am reading lines in a .txt file Each line is stocked into a n elements object, as this: [958] " 422 287 339 31 203 602 547 1026 500 366 346 227" [959] " 410 67 11 220 110 451 562 598 732 163 163 220" [960] " 179 513 95 186 102 595 333 1289 804 210 294 459" [961] " 276 153 307 138 126 233 623 739 521 421 209 75" [962] " 64

Hi! I am trying to install in my R environment the rhdf5 package and library but it seems to have vanished from either the CRAN or BioConductors sites. Can you tell me where it would be possible to find it or any R library (or function) able to read hdf files? Sincerely, Nicolas Degallier UMR 7159 / IRD UR182 Laboratoire d'Oc?anographie et du Climat, Exp?rimentation et Approches

Dear R'helpers, I guess the solution is trivial but despite some hours of testing and looking at the doc, I was unable to fix the following problem. I want to do either: for (i in 1:lat) a[i,which(a[i,]==0] <- NA or: for (i in 1:lat) a[i,which(a[i,]==-9999|a[i,]==9999|a[i,]==9998|a[i,] ==-9998)] <- NA if the word "Time" is or is not present in variable

Dear R'Helpers and Colleagues, I have looked in the documentation, asked to some colleagues in the lab, but was unable to solve the following difficulty. I am running R on a iMac G5 with OS 10.4. The file below (73 rows x 144 col) shows the values of a climatic index on the globe with a grid of 2,5 ? x 2,5 ? (NA = no value): ? With image() and map() and running the following

My aim is to create a figure consisting of three maps: Alaska in the upper left corner, the 48 contiguous US states in the center right, and Hawaii in the lower left corner. In some ways the figure I'm trying to create is analogous to figure 1.5 in Paul Murrell's excellent "R Graphics", which combines a map of New Zealand with a map of the world. I tried adapting the code

My goal is to prepare a thematic map of the US, with states shaded according to their values for a variable of interest. I would like to include an inset for Alaska in the upper left and an inset for Hawaii in the lower left. If possible, I'd like to use Albers conic projection, or something similar. Thus far I have tried using the maps package with its state database (which omits

Hi, I'm trying to plot a map of the pacific ocean, centered on the dateline, using the maps package. library(maps) # Basic library to draw maps library(mapdata)# Library with specialized maps library(mapproj) map(database = "world", fill = TRUE, col = 1, plot = TRUE,add=F, xlim = c(120,300), ylim = c(-20,40)) I tried several different combinations of coordinates, but whatever

I just uploaded two packages to CRAN. shapefiles_0.1.tar.gz - functions to read and write ESRI shapefiles (including dbfs) emme2_0.1.tar.gz - functions to read binary data from an EMME/2 databank data (EMME/2 is a transportation modeling program) Please let me know if you find any bugs or have some suggestions. Thanks. Regards, Benjamin Stabler Transportation Planning Analysis Unit Oregon

cheers , I need help in shapefile manipulations. I have two shapefiles of my country Cameroon. The first contain 10 provinces and each province contains a certain number of administratives units. I dont have ESRI Arc view. I want to add a admistrative unit of the center province to the map with 10 province. Is it possible to to take a peace of shapefile and combining it to another ? Thanks !

Dear List, I have examined the sp, maptools, and shapefiles packages (Windows) but cannot seem to figure out a solution for writing shapefiles for 'curvy shapes'. The scenario is that I generate a series of polygons, circles, and ellipses on a plot and would like to convert them to shapefiles. Since the circles and ellipses are do not return a coordinate list during the plotting

Hi, I have overlayed 2 projected shapefiles using the "plot" function. When I plot them the numbers next to the x-axis ticks are printed in scientific format (e.g. -4e+05, -3e+05,...,), and the problem is these dont match up to the projected point locations I also wish to overlay from my dataset. For example, according to the x-axis ticks the prime meridian (0?) is in totally the

Now available on CRAN is a revised version of the shapefiles package for reading and writing shapefiles in R. New additions, courtesy of others, include the ability to convert a simple R data frame of points, polylines or polygons to a shp format list, which can then be written out to a shapefile with write.shp. There is also a function to convert the read.shp shp format list to a simple data

I am attempting to write what MATLAB does as : [S,A]= shaperead() which returns an N by 1 geographic array structure S containing geometric information, and a paralell N by 1 attribute structure array, A containing feature attribute information. How can this be done with R code I tried using shp<-readShapePoints() but this doesnt return the same structure that I am looking for. Please

Hi, I'm a SAS veteran and R newbie. Can anyone point me to where I can best find examples on reading in shapefiles with R, and then merging in measurement info for the different polygons? I'm struggling a bit finding a good toehold in the documentation, as I am not yet familiar with navigating it yet. Thanks, Matt [[alternative HTML version deleted]]

I'm hoping that someone could guide me in how to extract data from shapefiles. I want to extract data from a shapefile (classed as "SpatialPolygonsDataFrame") and more specifically, the data is contained within the slot called "coords" from the class "polygons" within this file and despite my best efforts (and much reading!), I'm at a lost. Many thanks, Kate

Hi, Is there a way, preferably with R, to read shapefiles and transform them in a format that I could then use with package PBSmapping? I have been able to read such files into R with maptools' read.shape and plot it with plot.Map, but I'd like to bring the data to PBSmapping and plot from there. I also looked at the package shapefile, but it does not seem to do what I want

Hi all, I'm looking for zip code shapefiles for Italy and France.Does anyone know where I might be able to find these? thank you! [[alternative HTML version deleted]]

Hi, I meet a strange trouble writing shapefiles with write.pointShape() from maptools. > write.pointShape(mylines[,5:9], file="Sxtplinshp", mycoordspt) Error in write.pointShape(mylines[, 5:9], file = "Sxtplinshp", mycoordspt) : shapefile names must conform to the 8.3 format It seems that the function request a special shapefile name spelling (conformed to 8.3

Dear List, I'm trying to import a multiple polygon shapefile in splancs. I've found the following command but it seems to work only for single polygons: Name_of_splancs_polygon<-getPolygonCoordsSlot (getPolygonsPolygonsSlot(getSpPpolygonsSlot(name_of_sp_polygon) [[1]]) [[1]]) I've tried to change the numbers but it seams to choose polygons instead of importing them

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