Re: [R-sig-Geo] Converting kml to shape (readOGR)

2009-01-07 Thread Henk Sierdsema
Hi, Thanks for the all the information on kml's and rebuilding gdal which sets me on a new paths using stuff like VC++, source codes, etc. I started this 'excercise' because we receive more and more kml-files (from Google Earth) of observers with the outlines of plots and/or observations. But

Re: [R-sig-Geo] combining two image files into one database

2009-01-07 Thread Greg Lee
Hello Alok, We really need more details about what you are trying to acheive in order to help you. c(X, Y) concatenates the objects X and Y. The result will be a vector. If X and Y happen to be m x n matrices, the result will be a vector of length 2 * m * n: probably not what you want. You can

[R-sig-Geo] no boundaries between countries by using unionSpatialPolygons

2009-01-07 Thread Hans-Jörg Bibiko
Hi, I have a tiny problem. I'm using Roger Bivand's ( http://spatial.nhh.no/R/etc/TM_WORLD_BORDERS_SIMPL-0.2.RData ) to create a plot showing two countries without boundaries. I did this: library(foreign) library(sp) library(maptools) library(gpclib)

[R-sig-Geo] Wald statistics

2009-01-07 Thread Annachiara Saguatti
Hi, I'd need some information about the Wald statistics which is displayed in the summaries of the lagsarlm and errorsarlm regressions. Can anyone explain me which the null hypothesis is? Thank you very much. Best regards, Annachiara Saguatti [[alternative HTML version deleted]]

Re: [R-sig-Geo] no boundaries between countries by using unionSpatialPolygons

2009-01-07 Thread Roger Bivand
On Wed, 7 Jan 2009, Hans-Jörg Bibiko wrote: Hi, I have a tiny problem. I'm using Roger Bivand's ( http://spatial.nhh.no/R/etc/TM_WORLD_BORDERS_SIMPL-0.2.RData ) to create a plot showing two countries without boundaries. I did this: library(foreign) library(sp) library(maptools)

Re: [R-sig-Geo] Wald statistics

2009-01-07 Thread Roger Bivand
On Wed, 7 Jan 2009, Annachiara Saguatti wrote: Hi, I'd need some information about the Wald statistics which is displayed in the summaries of the lagsarlm and errorsarlm regressions. Can anyone explain me which the null hypothesis is? ?Wald1.sarlm explains - the test reports the values for

Re: [R-sig-Geo] no boundaries between countries by using unionSpatialPolygons

2009-01-07 Thread Hans-Jörg Bibiko
On 07.01.2009, at 15:47, Roger Bivand wrote: Your ID has two unique values: ID [1] (-67.1,-64.6] [-74.4,-71.9] Levels: [-74.4,-71.9] (-71.9,-69.5] (-69.5,-67.1] (-67.1,-64.6] with two country polygons: length(slot(pe_bo_data, polygons)) [1] 2 so you get two output polygons. If you do:

Re: [R-sig-Geo] combining two image files into one database

2009-01-07 Thread Alok K. Bohara, PhD
Thanks Greg for your response: I want to do what the examples in Spatstat do with regard to the point and pixel/image data bei and bei.extra.. (bei.extra contains a list of two image files --elevation and gradient. They plot these two items (and also take log, I believe) In my case:

Re: [R-sig-Geo] combining two image files into one database

2009-01-07 Thread Tomislav Hengl
Dear Alok, Maybe you should consider converting the im data to sp class. Then, you have as many grids in a single SpatialGridDataFrame e.g.: library(maptools) library(rgdal) grids - as(bei.extra[[1]], SpatialGridDataFrame) names(grids)[1] - elev grids$grad - as(bei.extra[[2]],