similar to: mapping bathymetries and species abundances

System: Linux Ubuntu 7.10 Gibbon kernel 2.6.22-14-generic Emacs ver. 22.1.1 ESS ver. 5.3.0 R version 2.6.0 (2007-10-03) ------------------------------------------------ Colleagues This is a follow-up note to my earlier post under this header, giving my solution. First, with some trepidation I upgraded my ubuntu distro, and this time did a full reinstall, so I did not break my system, as I have

Dear List, I am trying to?plot?bathymetry contours around the Hawaiian Islands using the package rgdal and PBSmapping.? I have run into a memory limit when trying to combine two fairly small objects using cbind().? I have increased the memory to 4GB, but am being told I can't allocate a vector of size 240 Kb.? I am running R 2.11.1 on a Dell Optiplex 760 with Windows XP.? I have pasted

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, i have a longitude vector (x) a latitude vector (y) and a matrix of bathymetry (z) with the dimensions (x,y). I have already succeeded in plotting it with the image.plot (package 'field') and the contour functions. But now, I want to make a grid in order to extract easily the bathymetry corresponding to a couple of longitude, latitude coordinates. Do you know a function or a

Dear R Users, Using the following R-script I created the first image > require(akima) > require(spatial) > dep <- interp(long, lat, depth, xo=seq(1,990,10), yo=seq(1,990,10), + extrap=FALSE, ncp=0,duplicate = "mean", dupfun = NULL) http://r.789695.n4.nabble.com/file/n3431391/Rpics.bmp Where "long" are x-coordinates between 1 and 1000, "lat" are

I am trying to create a NetCDF file with bathymetry (a matrix) and then to add a grid type (an integer) as variables. This is the relevant part of the code: library(ncdf) wfile="my_file.nc" msv=-10 grdtp=2 # Dimensions d1=dim.def.ncdf("lon","degrees",as.double(lon)) d2=dim.def.ncdf("lat","degrees",as.double(lat)) # Variables

Hi I am trying to estimate bottom temperatures over a particular depth range, based on one dataset containing synthetic temperature profiles (for set depths) and another that contains information on bathymetry. I need to do this for multiple regions and thus would ideally like an automated solution that is quicker and more accurate than the manual option I currently have! I have not had much luck

Using R 2.5.1 on Windows XP Professional, and PBSmapping package version 2.51, I have encountered some behavior which puzzles me. I am including the package's listed maintainer on this email but also seek the thoughts of the R-help community. I have a set of EventData, which I want to plot as points, and to color the points according to some criterion. It turns out that some of my

Hi everybody. i have 2 datasets : one with a long grid a lat grid and a bathymetry grid the second one only has the long and lat coordinates. So i want to know the bathymetry associated to the second dataset. I thought it was a good idea to do an interpolation between the two datasets. I find a script doing that on matlab but i want to do it on R. Do you how to do it? for ii =

Hi folks, I've been using 'PBSmapping' to make a map of Europe with some labels. I've been using the 'worldLL' PolyData, as my computer is too slow to make my own from the GSHHS files. The only problem is this PolyData does not seem to include Ireland. I have no idea why this should be so, other European islands such as Sardinia etc. are included.

Hi All- I would like to create a 3d topographic map using lat/lon and z(height). I have been scouring the R help pages and have not located the package I am looking for. Does anyone have a suggestion of package that will work for this? thanks- sherri

Dear List, I am trying to use the PBSmapping package to import and plot the shoreline of Hawaii. I am having problems importing and plotting the specific regions that I would like to plot. Specifically, I can't get it to import the range of x variables that I would like. I think my problem is with a variable called xoff, but I have no idea what it is, what it does, or how I am supposed to

Hello, List, I am a new user of the R project, and I need some help in plotting a legend. I am using the PBSmapping library to plot map of Ohio and heat color it with the count of employees in each county. As a guide, I am using "Data Mashups in R." I am able to plot the map with the colors; however, I would like to put a legend with a single box full of these colors, but only show

Hi, In another thread ("PBSmapping and shapefiles") I asked for an easy way to read "shapefiles" and transform them in data that PBSmapping could use. One person is exploring some ways of doing this, but it is possible I'll have to do this "manually". With package "maptools" I am able to extract the information I need from a shapefile but it is

Hi R helpers, I tried to convert a list of LatLong coordinates (in DD format) into UTM zone 11 NAD 27. I first tried this from PBSmapping: library(PBSmapping) LatLong<-cbind(c(56.85359, 56.85478),c(-118.4109, -118.4035)) colnames(LatLong)<-c("X","Y") attr(LatLong, "projection") <- "UTM" attr(LatLong, "zone") <- 11

Hola, Estoy tratando de sobreponer un mapa sobre otro, uno que sirva de ?background? y otro con datos. Mi primer intento lo hice usando el paquete PBSmapping. Aquí lo que hice: library(PBSmapping) estados <- importShapefile(?Entidades_2013.shp") #Aquí cargamos las entidades federativas (32 entidades) zonas <- importShapefile(?zm.shp?) #Aquí cargamos las zonas metropolitanas del país

I have a spatialpixelsdataframe containing bathymetry data. Is it possible to save this as something similar to the meuse.grid so it can be loaded without having to run the spatial interpolation each time the script is called? I don't see a straightforward method in the sp documentation.

Dear all, I want make correction depth of a bathymetric data set. To do so, I have the depth data set sample every second (a depth at each second) in one hand, and in the other hand, I have a tide variation level data set sample every 250 ms. The time register in each data sets (tide and bathymetric) is express in seconds followinf this format : hh:mm:ss.ss I would like to rectify the depth

Dear UseRs, I would like to know which function is the most efficient in finding convex hull of points in 3(or 2)-dimensional case? Functions for finding convex hull is the following: convex.hull (tripack), chull (grDevices), in.chull (sgeostat), convhulln (geometry), convexhull.xy (spatstat), calcConvexHull (PBSmapping). I also would like to know if there is a function that can be used

I have a dataset giving traffic between pairs of ports, and the lat/lon of each port, roughly as follows: sample data as follows > df2VisitTrips[1:4,] Origin Destination NumberOfTrips OriginLatitude OriginLongitude DestinationLatitude DestinationLongitude 1 P1 P16 1 39.45965 -80.15633 40.76111 -79.54583 2 P1 P3