Re: [R-sig-Geo] Time vs. Longitude (Hovmueller Diagram)

Hi Michael, Thanks for the reply. The data file can be rechieved from the following link: ftp://ftp.wiley.com/public/sci_tech_med/spatio_temporal_data SST011970_032003.dat SSTdata_descrip.txt

Re: [R-sig-Geo] Time vs. Longitude (Hovmueller Diagram)

At a guess, for a given latitude you could do this where "x" is a data.frame with columns Longitude, Latitude, Time1, ..., TimeN x1 <- x["Latitude" == -29, ] ## very simplistically ## drop Latitude, and order by Longitude (assuming Time columns are in the right order) x2 <- as.matrix(x1[order(x1

[R-sig-Geo] Time vs. Longitude (Hovmueller Diagram)

Dear Experts, I am trying to plot the Hovmueller Diagram (example below), which is time vs. longitude in R. http://www.esrl.noaa.gov/psd/cgi-bin/hovmoller/timeplot.pl?var=air&level=1000&mon1=1&mon2=1&dy1=1&dy2=31&yr1=2010&yr2=2011&datatype=operational&type=mean&fxdlon=no&postscript=no&lon1=&lon2=

Re: [R-sig-Geo] anisotropy modeling [SEC=UNCLASSIFIED]

Hi Jon, Thanks for your kind comments. The shape was indeed not expected to change much as the same variogram model (i.e., exp) was used. Cheers, Jin -Original Message- From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-boun...@r-project.org] On Behalf Of Jon Olav Skoien Sent: Tuesd

Re: [R-sig-Geo] merging polygons by indices

On Tue, 16 Aug 2011, majesty wrote: still I am not quiet sure how to achieve this by manipulation of ID. Could you give an example, say I have 5 polygons in my shapefile shape.shp but would like to have 4 polygons, where shape[2] is merged with shape[3]. I am somewhat clueless because ID has to

Re: [R-sig-Geo] merging polygons by indices

still I am not quiet sure how to achieve this by manipulation of ID. Could you give an example, say I have 5 polygons in my shapefile shape.shp but would like to have 4 polygons, where shape[2] is merged with shape[3]. I am somewhat clueless because ID has to be of equal in length to polygons slot

Re: [R-sig-Geo] merging polygons by indices

On Tue, 16 Aug 2011, majesty wrote: Dear R-sig-geo subscribers, I have the following problem, would be grateful for suggestions. I am doing a research on UK output areas, My .shp file contains output areas for UK but there are some duplicates. when I type which(duplicated(shape$output_area)==TR

[R-sig-Geo] merging polygons by indices

Dear R-sig-geo subscribers, I have the following problem, would be grateful for suggestions. I am doing a research on UK output areas, My .shp file contains output areas for UK but there are some duplicates. when I type which(duplicated(shape$output_area)==TRUE) I get indices of 67 duplicates: [

Re: [R-sig-Geo] geo meetup at UseR! 2011

Yes, it is a late-breaking poster (plotKML: a framework for visualization of space-time data). Blue with some Google Earth screenshots on it - that hopefully would appeal to the R-spatial enthusiasts :) See you there, Pierre 2011/8/16 Barry Rowlingson : > On Mon, Aug 15, 2011 at 10:26 PM, Pierre

Re: [R-sig-Geo] geo meetup at UseR! 2011

On Mon, Aug 15, 2011 at 10:26 PM, Pierre Roudier wrote: > Hi Barry, > > I'll be happy to join :) > > I'm sure there's some more spatial useRs out there in Warwick at the moment, The R-sig-geo meetup will begin at Pierre Roudier's poster at the poster session tonight, at 9pm. I dont see it in the

Re: [R-sig-Geo] Error when projecting polygons into Azimuthal Equidistant Projection

On Tue, 16 Aug 2011, bart wrote: Thank you very much for your answer Roger, but the problem is that i want to do this operation for a regular grid around the world. Meaning i will encounter the error in some places. As far as i understand the projection the distance and direction to other loca

Re: [R-sig-Geo] gstat: Variogram with non-Euclidean Distances

Ben and Luca, Thank you very much for your advice and the references. When using non-Euclidean distance I do precisely run into issue of obtaining non-positive definite correlation matrices, so I might have to return to Euclidean distance... Regards, Julia __

Re: [R-sig-Geo] Error when projecting polygons into Azimuthal Equidistant Projection

Thank you very much for your answer Roger, but the problem is that i want to do this operation for a regular grid around the world. Meaning i will encounter the error in some places. As far as i understand the projection the distance and direction to other locations is only valid if lat_0 and l

[R-sig-Geo] model comparisons ruf.fit

Hi, I'm modelling resource use using ruf.fit, but I'm a little confused about the Materns logLik / LS LogLik given by summary(model). It seems that Materns logLik changes if I change the values for theta, but stays fairly stable if I change the predictor variables. LS LogLik stays the same with va

Re: [R-sig-Geo] anisotropy modeling [SEC=UNCLASSIFIED]

Hi Jin, some comments below: On 16-Aug-11 3:38, jin...@ga.gov.au wrote: Hi All, Following Edzer's examples, I have tried to test the effects of the specification of alpha in variogram on the results of fit.variogram as below: v = variogram(zinc~1,meuse) v0 = variogram(zinc~1,meuse, alp

Re: [R-sig-Geo] Segmentation of SpatialLines at Intersection

Hi, You can do this in GRASS with v.clean module : v.clean input=lines1 output=lines2 err=points tool=break type=line You can use spgrass to integrate GRASS into R I hope this helps, Mathieu 2011/8/15 Cam Plouffe > Hi everybody, > > I have a SpatialLines object that I created by using the