similar to: SpatialLines

Hi Ashraf, It is not obvious to me what your structures are but one problem in your function is the assignment tt1 <- SpatialLines(list(tt[[i]])). This will set tt1 to just have one item. Consider the following test.func <- function(x) { tt1 <- list() for ( i in ... ) { ... tt1[[i]] <- SpatialLines(tt[[i]]) } return(tt1) } HTH, Eric On Mon, Sep 25,

Hi Eric, Thanks for the help.But this will not solve the problem as it will generate a list and what I need is an object of class sp using SpatialLine function from sp package.So, I need to convert each matrix to coordinates and then to a line and then to a spatial line as figured in the code. My data structure is a list of 141 matrices.Each matrix represents coordinates of the river lines

Hi all,I'm trying to build a SpatialLines object from a list that contains 124 river segments. Each segment in the list contains the x,y coordinates. I'm using the following code to create the SpatialLines object, but it just retrieves one segment. Any suggestions? test.func = function(x){ ??? for (i in 1:length(x)) {??????? tt[[i]] <- x[i]; tt[[i]]? = Line(tt[[i]]); tt[[i]]? =

Hi Ashraf, In that case I think you may need to structure the code to first build the list and only at the end supply that to the SpatialLines function, something like test.func <- function(x) { tt <- list() for ( i in ... ) { ... tt[[i]] <- (whatever) } return(SpatialLines(tt)) } Eric On Tue, Sep 26, 2017 at 12:36 PM, Ashraf Afana <asafaneh at

Hi! Is there a way to manually costumize the color legend in an spplot() - especially where to draw ticks and labels for the ticks? The reason I'm asking: Usually spplot() automatically divides the data into fitting slices and makes a color legend (also automatically). I want to assign the slices myself and have a fixed scale instead of an automatic/dynamic scale. I think what I want

Hi! How does one change the background color of the map-panel in spplot()? Example: library(sp) data(meuse.grid) gridded(meuse.grid) = ~x+y spplot(meuse.grid, "part.a") How would I get another background-color for the map-panel (but not for the whole plot) here? Thank you! Marcel

Hi, I am very new to R and I am trying to write KML files. I am running R-2.13.1 and I have the following packages: - AKIMA - RGDAL - RSAGA - SP I was trying to run the simple example that I found at the following web page: http://spatial-analyst.net/wiki/index.php?title=Export_maps_to_GE The lines that I can run are: data(meuse.grid)

SUMMARY: Specific problem: I'm regridding biomass-burning emissions from a global/unprojected inventory to a regional projection (LCC over North America). I need to have boundaries for Canada, Mexico, and US (including US states), but also Caribbean and Atlantic nations (notably the Bahamas). I would also like to add Canadian provinces and Mexican states. How to put these together? General

Dear forum,   I am working in R 2.9.1 and I am trying to sample locations from a network file. Reading in nor plotting is a problem, however when I am trying to sample from the file I get the following message:   nwlim<-readShapeLines("C:/Limburg_nwshape", proj4string=CRS("+init=epsg:31300")) plot(nwlim)   randacc<-spsample(nwlim,n=1000,"random")

hi all, I have a plot to which i would like to add text labels. And i cant find a way...here is the code : enaD2<-idw(D2~1, loca=dva, newdata=grd) pts = list("sp.points", dva, pch = 20, cex=1.5, col = "darkred spplot(enaD2, "var1.pred",sp.layout=pts, main = "globina 60 cm", sub="D2",

Hi ALL, With the pakage gstat How can one create a grid of points as in the meuse.grid example, from measuring sites as in the Meuse example Thanks in advance! Alain PONSERO Conservateur Réserve Naturelle de la Baie de Saint-Brieuc site de l'Etoile 22120 HILLION tel/fax : 02.96.32.31.40 <mailto:aponsero@cabri22.com> aponsero@cabri22.com

Hi, please help!! I've been trying to fit variogram models using the REML method in the gstat package. Every time the Windows GUI crashes. For example library(gstat) data(meuse) x <- variogram(zinc ~ 1, ~x + y, meuse) v <- vgm(140000, "Sph", 800, nug = 10000) plot(x, model = fit.variogram(x, model = v, fit.method=5)) Other fit methods are non problematic (eg. fit.method=7

One more thing I forgot to mention. The chapter calls for editing /etc/resolv.conf, but in my case it won't stay edited--it keeps getting set back to an original form (for searching my ISP) on every reboot. Thanks again. Eric Hines -------- Original Message -------- Subject: Samba-3 By Ex Chapt 3 Date: Tue, 05 Jul 2005 17:47:09 -0500 From: Eric Hines <eehines@comcast.net> To:

Dumb question time; I can't find the answer in the Samba-3 or HOWTO docs: how do you add users to the passdb.tdb (tdbsam's db)? Is pdbedit the only way? I ask because the Chapt 3 Samba-3 example has passdb = tdbsam in the samba config file, but the instructions for adding users are to use useradd and smbpasswd, which leave passdb.tdb empty (except for root--I have no idea how that

Hi ALL, Can you please help me on how to determine the prediction error for ordinary kriging?Below are all the commands i used to generate the OK plot: rsa2 <- readShapeSpatial("residentialsa", CRS("+proj=tmerc +lat_0=39.66666666666666 +lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs")) x2 <- readShapeSpatial("ptna2",

Folks, This is a strange problem, and I'm not sure this is the right place for it, so feel free point me to a better location. I'm running FC3 and Samba 3.0.14a. The problem is this: when I boot up my Linux box, and it gets to the point of turning on my NIC (eth0), where the PC used to turn it right on and continue, now there's an inordinately long delay in getting eth0 to

Hi! Is there a way to manually costumize the color legend in an spplot() - especially where to draw ticks and labels for the ticks? The reason I'm asking: Usually spplot() automatically divides the data into fitting slices and makes a color legend (also automatically). I want to assign the slices myself and have a fixed scale instead of an automatic/dynamic scale. I think what I want gets

Soon on CRAN a new version of package gstat will emerge, which has a few minor changes and possible incompatibilities w.r.t. the previous version(s). The new gstat (0.9-23) now: + depends on sp, and uses internally with Spatial* classes from sp if data are provided in the old-fashoned way (as data.frame) + has a vignette to get you started with the classes in sp + defines krige as a generic;

Soon on CRAN a new version of package gstat will emerge, which has a few minor changes and possible incompatibilities w.r.t. the previous version(s). The new gstat (0.9-23) now: + depends on sp, and uses internally with Spatial* classes from sp if data are provided in the old-fashoned way (as data.frame) + has a vignette to get you started with the classes in sp + defines krige as a generic;

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.