similar to: smoothing bathymetry

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

I'm new at this, so please bear with me. I am trying to get Krig (actually tps) results into a data frame at a set of points. Frankly, all I really need is to query the Krig on a specific XY pair and return a result. I can carry the rest. I am using tps from the fields package. Here's the code I use to create the Krig object: =========== Snip! =============================== dCurr

Hello, I have a long list of x-, y- and z-data and try to generate a heatmap. Obviously there are several data with identical x- and y-values. I get the following error message: Error in interp.old(x, y, z, xo = xo, yo = yo, ncp = 0, extrap = extrap, : duplicate data points: need to set 'duplicate = ..' Unfortunately there seems no help screen on "duplicate". I'd prefer

Dear all, I am new to R and even newer to mapping with R. I have installed PBSMapping, maptools and maps, and I am struggling to produce a map of fish abundance in the Grand Banks. I have tried several approaches to plot the bathymetry in my study area: 1) download the data from http://topex.ucsd.edu/cgi-bin/get_data.cgi (42N46N313W307W) To start with, the map I get is blank between

DeaR list, I'm running an external program that computes some electromagnetic response of a scattering body. The numerical scheme is based on a discretization with a characteristic mesh size "y". The smaller y is, the better the result (but obviously the computation will take longer). A convergence study showed the error between the computed values and the exact solution

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

Dear all, I would like to estimate rainfall with the krige function. First, I produce a polygon of my region of interest (where poly_test.txt is a x,y suite of points defining a polygon, obtained with the grass v.out.ascii command) : user>border_limite<-read.table("/home/lionel/rwork/poly_test.txt",header=FALSE) user>polygone<-list(x=border_limite[,1],y=border_limite[,2])

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

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 =

Dear R users, I want to do a fitted.contour plot of selected columns of a dataframe M with M$AM and M$Irradiance as x and y axes respectively. The level of the contour shall be determined by M$PR. Some words on my data first. Dataframe M looks like: head(M$Irradiance) [1] 293 350 412 419 477 509 head(M$AM) [1] 2.407 2.161 1.964 1.805 1.673 1.563 head(M$PR) [1] 70.102 72.600 75.097 80.167

This script below pulls yahoo data via a function in quantmod, then massages the data around to forumalate a 3D graph with RGL library, attached is a ggplot to show the data i'm trying to create a surface with in separate line geoms . the issue is that the 3D graph looks very ugly and cut up because of the limited quantities of points on the front month expirations.. can anyone tell me whats

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

Hi, I want to make a surface3D plot of a landscape. I have cordinates (x, y, z) recorded with a GPS. The datapoints are not evenly distributed within the rectangular area. To do a fast 3D plot I used following. > library(grid) > library(lattice) > v <- read.table("clipboard") > names(v) <- c("x", "y", "z") > wireframe(z ~ x * y, data =

Full_Name: Jonathan Lees Version: 2.0.1 OS: linux-gnu Submission from: (NULL) (152.2.75.65) there is a problem with calculating the convex hull in 2-D interpolation using the codes interp fromt eh akima package: x =c(0.6505304, -1.1821562, -0.2600792, 0.7913716) y = c(1.0424226, 0.1754048, -1.4523334, 0.2349112) z = c(0.000, 3.042, 0.370, 0.122) EX = seq(from=min(x), to=max(x),

Dear R-helpers, I'm having trouble compiling R packages from source on Win2K. I installed R 1.7.1beta [1] on my win2k machine [2], downloaded the fields package as source [3] and tried but failed to install the package [4]. I get the same problem with other packages, so it doesn't appear to be limited to fields. Installation of precompiled packages seems to work fine. I've been

In the code to follow, I'm trying to label points with their corresponding values but have been unsuccessful after many variations to the following code. The code below creates the plot I want, I simply cannot get the black points ("+") to display the actual value. I'm guessing the problem is somewhere in the second to last line of code (starts with "pts<-..."). I

Dear List, I'm using interp() to prepare 3d data for plotting with the contour() function. If have x,y and z data. All are arrays. X and Y are sampled in an orderly fashion on a grid (a circular sub-area of a grid - see plot). I'm trying to use interp() to get x and y arrays and a z matrix that can be fed to contour(). This is the command:

I am returning to using R after a little time away. Saw some great new functions in the fields package and decided to upgrade. Everything seems to go fine until I tried to actually load the package, when I get an error: >library(fields) Error in parse(file, n, text, prompt) : syntax error on line 227 And (of course) I am unable to access any of the functions in the fields package. Below is

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 list, I get the following error message when installing the package GRASS. missing link(s): interp.new missing link(s): vgram Krig predict.Krig predict.se.Krig See also below for the full installation message. What does this mean exactly? Regards, Ulrich ------------------------------------------------------------------ > R CMD INSTALL