search for: bathymetry

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

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 fun...

...tempted this using the following script, which generated the second image. >dep <- interp(long, lat, depth, xo=seq(1,990,10), yo=seq(1,990,10), + linear=FALSE,extrap=TRUE, ncp=NULL,duplicate = "mean", dupfun = NULL) However, as seen in the image, this results in a messy looking bathymetry, with erros. Any suggestions on what went wrong or other packages or functions I could try to would be greatly appreciated? I played around with the krig.image function in the "fields" package, but didnt have much luck. Sincerly Jens -- View this message in context: http://r.789695...

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 lon -52W -53W. Besides, it seems that the shallower isobath is 500 m. Or maybe I am wrong here? 2) import a shapefile with the bathymetry,...

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 i...

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 and was wondering if someone could help? I understand that this may be a bit of an ask, but any help/ suggest...

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

...G file. This is good where the journal wants your graphics files separately and may not be able to deal with SVG files, or very large files. I have one caveat. It seems to me that RSvgDevice had problems writing to file when the data were unreasonably voluminous. I produced another graphic using bathymetry and coastlines recommended in the PBSmapping package. These files are high resolution and large (85 and 64 MB respectively). The plot took too long and the resulting SVG could not be loaded by Inkscape, nor converted to EPS using Image Magick (suggesting failure somewhere?). I'm not criticizin...

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

...et 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 (bathymetry) by substrating the tide at each depth time. I tried to construct time series using pastecs package but didn't succed Example of the tide data set: Time Correction 11:55:07.58 -2.64 11:55:07.68 -2.64 11:55:07.88 -2.64 11:55:08.28 -2.64 11:55:08.38 -2.64 Example of the bathym...

Dear members, How can I read and plot an ArcGIS raster file into R ? The file has extension .aux and contains floating point bathymetry data. The purpose is to create a spatial model in R that uses ArcGIS map data. I have managed to read and plot various shape files into my R project, but I am stuck with this now. I am new to this list and also to R, so any help would be much appreciated. Many thanks and best wishes, Juliane...

...unction myquadrature inside the script so that it returns the quadrature of descrete data using the trapezoidal rule. Modify the call to the function at the bottom of the script so that is uses your modifies routine to compute the total volume flux through a transect using descrete observations of bathymetry and vertically-averaged velocity. The professor supplied us with some sample code and i'm pretty sure it outlines the whole problem. My issue is, I have no idea how to write the sum code for the trapezoid rule to estimate the area under the curve: Here is the code provided: myquadrature = fu...

..., if you want to multiple two series of numbers together % for example c = a*b where c(1) = a(1)*b(1), c(2) = a(2)*b(2), etc. % you must tell Matlab you want "element by element" multiplication % e.g.: c = a.*b % note the "." % h = 10.*(cos(((2*pi)/2000)*(x-mean(x)))+1); %bathymetry u = 1.*(cos(((2*pi)/2000)*(x-mean(x)))+1); %vertically-averaged cross-transect velocity plot(x,-h) % set begin and end points for the integration a = x(1); b = x(end); % call your quadrature function. Hint, the answer should be 30000. f=u.*h; val = myquadrature(f,a,b); fprintf('the solution...

Dear all; Thank you for your reply. David has explained an interesting method. David I have DEM file of the region and I have extracted the xyz data from that. Also I can extract bathymetry data as xyz file. I have calculated the storage (volume) of reservoir at the current elevation. But the method I have used to calculate the volume is different from your method. I have crop DEM by the reservoir boundary and then calculate the volume. I would be more than happy if you please please...

...he 'method' parameter doesn't work very well (or at all?) when xlim and ylim are also set. The following code should recreate and demonstrate the problem: --------------------------------------------------------------------------------------- library(marelac) data(Bathyometry) B <- Bathymetry xrange <- c(-150, -120) yrange <- c(40,60) par(mfrow=c(1,2)) ### using xlim and ylim #### contour(B, xlim=xrange, ylim=yrange, method='edge', main='x/y limited') #### using subsetting #### x <- !is.na(cut(B$x,breaks=xrange)) y <- !is.na(cut(B$y,breaks=yrange)) contou...

Water engineer here. The standard approach is to 1) get the storage vs. elevation data from the designers of the reservoir or, barring that, 2) get the bathymetry data from USBR or state DWR, or, if available, get the DEM data from USGS if the survey was done before the reservoir was built or 3) get a boat+sonar with GPS? +lots of time and survey the bottom elevation yourself. Put the xyz data into ArcGIS and have it create the bottom surface, then, with...

Hi, As already mentioned, this is a great topic for R-sig-geo, where you'd probably get specialist answers like the lake morphology package https://cran.r-project.org/web/packages/lakemorpho/index.html which is explicitly designed for this kind of question. Sarah On Tue, Dec 5, 2023 at 11:13?PM javad bayat <j.bayat194 at gmail.com> wrote: > > Dear all; > I am trying to

System: Linux kernel 2.6.15 Ubuntu dapper R version 2.5.1 ESS 5.2.11 on Emacs 21.4.1 ------------------------- Colleagues Having read the posts on producing perfect graphs in R, or using inkscape to edit R graphics output, I have a related question. Lately I am publishing in a journal that is very picky about their graphics formats. It is more efficient for me to get the final picky details

Hello everyone (and Hadley in particular), I often need to plot data from multiple datasets on the same graph. A common example is when mapping some values: I want to plot the underlying map and then add the points. I currently do it with base graphics, by recording the maximum region in which my map+point will fit, plotting both with these xlim and ylim parameters, adding par (new=T)