Re: [R-sig-Geo] maptools/kmlPolygon example: two SPDFs, different results

2009-08-27 Thread Robert J. Hijmans
Rick, I had the same problem. Perhaps I simply misunderstood the maptools function. I wrote the script below and that worked for me (it includes all polygons). It probably has some limitations. Robert polykml <- function(sw, filename, kmlname, namefield, kmldesc="Made with R" ) { out <- s

[R-sig-Geo] help on Initializing the matern at point zero

2009-08-27 Thread TANKISO THEJANE
R users Let me refine my request for help. I have a matern as per the code below. What i want to do is to initialize the matern at point/distance zero to be equal to 1.1 not 1.0 as is at present. I would appreciate your help very much.   library(geoR) curve(matern(x, phi= 1.0, kappa = 1.0),from

Re: [R-sig-Geo] Preparing KML files for GoogleMaps/Earth: Why no kmlPoint() in maptools package?

2009-08-27 Thread Robert J. Hijmans
Rick, I wondered the same thing a while ago and wrote this script; I worked for my purposes but never really finished it (e.g. SpatialPoints objects support; and I wanted to add folders). Perhaps you or someone else can finish it (or do a better one) and send it to the maptools maintainers? Robert

[R-sig-Geo] maptools/kmlPolygon example: two SPDFs, different results

2009-08-27 Thread rick reeves
Hello List: I'm using R 2.9 on Windows XP/64 for this example, with recent builds of maptools and rgdal packages. The shape file used in this example can be obtained from: http://nceas.ucsb.edu/scicomp/upload/UsCountiesGeogProj.zip While studying the example for kmlPolygon() in the maptools.

[R-sig-Geo] Lat/Lon Polygons

2009-08-27 Thread Matteo Mattiuzzi
Hello all! I've got a file with Lat/Lon organised in the following way: NAME | LATp1 | LONp1 | LATp2 | LONp2 | ... | LONp4 Each row correspond to an area coveres by the satellite scene (4 LAT/LON pairs). Two things are important for my now: 1. I need the line in the middle of the polygo

Re: [R-sig-Geo] Creating grid from xyz file

2009-08-27 Thread Paul Hiemstra
Mehdi Khan wrote: Sorry about that, this is what I did: wald<-read.table ("C:/Documents and Settings/MK/My Documents/USGS/R Data/wald/california.xyz") names(wald)[1]<- "LON" names(wald)[2]<- "LAT" names(wald)[3]<- "vs30" coordinates(wald)<- c("LON", "LAT") gridded(wald)<- TRUE suggested