[R-sig-Geo] how to get the values for corresponding neighbors

2010-05-04 Thread rusers.sh
Hi there, Take the data of columbus(spdep) as an example. Say “i” is one of the polygons, and “j” represents its neighbours. How to get the vlues of “sum(Wij*Xj)”? “Wij” means the weights, which can be get by “nb2listw$weights”. Say “Xj” to be “columbus$PLUMB”. The problem is how to get the

Re: [R-sig-Geo] generate simulation data for a theoretical spatial model

2010-02-03 Thread rusers.sh
://www.leg.ufpr.br/~paulojus On Tue, 2 Feb 2010, rusers.sh wrote: It works. The problem is that it only generates the simulated data based on our observed dataset,e.g. meuse here. I wonder if we can generate the simulated dataset from the user-specified model with covariates included

Re: [R-sig-Geo] generate simulation data for a theoretical spatial model

2010-02-02 Thread rusers.sh
rusers.sh wrote: Hi Tomislav, Thanks for your info on unconditional simulation. For conditional simulations, i still cannot find any useful information. I searched the R site and didnot find the possible method to do conditional simulations. 1. CondSimu(RandomField): trend: Not programmed yet

Re: [R-sig-Geo] generate simulation data for a theoretical spatial model

2010-02-01 Thread rusers.sh
...@spatial-analyst.net Dear rusers.sh, Here are few simple examples of how to simulate (not-normal) distributions and point processes using geoR and spatstat: http://spatial-analyst.net/book/node/388 See also: http://leg.ufpr.br/geoR/geoRdoc/vignette/geoRintro/geoRintrose8.html#x9-120008

[R-sig-Geo] generate simulation data for a theoretical spatial model

2010-01-30 Thread rusers.sh
Hi all, In classical statistics, we always need to generate a theoretical model such as y=a+b1*x1+b2*x2+e to study some new estimation content. I am wondering how to generate the similar spatial dataset for a theoretical model. Say y is response variable, x1 and x2 are explanatory variables. 1.

Re: [R-sig-Geo] complement part of a polygon in another polygon

2010-01-26 Thread rusers.sh
intersect.owin() works for intersection. Say polygon B located inside A. Is there any function in R for us to get the complement part of B in polygon A? 2010/1/26 Alexandre Villers alexandre.viller...@laposte.net Good morning, You can have a look at union.owin() and intersection.owin() in

[R-sig-Geo] Intersection between polygons

2010-01-25 Thread rusers.sh
Hi, Say i have two polygon shape files, A and B. Each one has many polygons. If we overlay the two layers, we should get three parts generally, PART1 only from A, PART2 only from B, and PART3 from both A AND B. In R, how can we get these three parts in shape format? Just point me several key

[R-sig-Geo] error for converting between KML and Shape file

2009-12-10 Thread rusers.sh
Hi all, I try to use the readOGR and writeOGR in rgdal package to do the conversion KML and Shape files, but failed. Following is my codes for a simulated dataset. #library(rgdal);library(maptools);ogrDrivers() #generate a simulated point SHAPE file,its name is pointexamle0 and save in

Re: [R-sig-Geo] error for converting between KML and Shape file

2009-12-10 Thread rusers.sh
Thanks Barry. The problem is i didnot understand the parameters in the functions correctly. Thanks for your detailed explanation. 2009/12/10 Barry Rowlingson b.rowling...@lancaster.ac.uk On Thu, Dec 10, 2009 at 8:17 PM, rusers.sh rusers...@gmail.com wrote: Hi all, I try to use the readOGR

Re: [R-sig-Geo] convert KML file into shape file or CSV file

2009-12-09 Thread rusers.sh
Barry Rowlingson b.rowling...@lancaster.ac.uk On Wed, Dec 9, 2009 at 1:22 AM, rusers.sh rusers...@gmail.com wrote: Hi all, Does anybody know how to import a KML file in R and them save it as a shape or CSV file. Take the following KML file as an example, http://www.nature.com/nature

[R-sig-Geo] convert KML file into shape file or CSV file

2009-12-08 Thread rusers.sh
Hi all, Does anybody know how to import a KML file in R and them save it as a shape or CSV file. Take the following KML file as an example, http://www.nature.com/nature/googleearth/avianflu1.kml I have tried several extentions of ARCGIS,e.g.KML2Shape,KMLCSVConverter,KMLCSVConverter and

Re: [R-sig-Geo] convert points into SHP polygons

2009-11-28 Thread rusers.sh
Perfect. Thanks very much. Also thanks for kind reminding on getting the centroids of polygons. 2009/11/28 Roger Bivand roger.biv...@nhh.no On Fri, 27 Nov 2009, rusers.sh wrote: Dear Robert, My method to generate the Voronoi tessellation maybe not better, You can show me you method to get

[R-sig-Geo] convert points into SHP polygons

2009-11-27 Thread rusers.sh
Hi, I always need to convert between the points and polygons. For converting the polygons to points, we can use the get.Pcent(maptools) to get the centroids of the polygons. But for converting points to polygons (something like GeoDa's tools-shape-point to polygon), i cannot find the

Re: [R-sig-Geo] convert points into SHP polygons

2009-11-27 Thread rusers.sh
#plot(pp) 2009/11/27 rusers.sh rusers...@gmail.com Hi, I always need to convert between the points and polygons. For converting the polygons to points, we can use the get.Pcent(maptools) to get the centroids of the polygons. But for converting points to polygons (something like GeoDa's

Re: [R-sig-Geo] convert points into SHP polygons

2009-11-27 Thread rusers.sh
(xlim=c(0,1),ylim=c(0,1)) 2009/11/27 Robert J. Hijmans r.hijm...@gmail.com I think you want a Voronoi tessalation (diagram): http://en.wikipedia.org/wiki/Voronoi_diagram But what are SHP polygons? Is 'SHP' a class in an R package? Robert On Fri, Nov 27, 2009 at 2:03 PM, rusers.sh rusers

Re: [R-sig-Geo] convert points into SHP polygons

2009-11-27 Thread rusers.sh
Dear Robert, My method to generate the Voronoi tessellation maybe not better, You can show me you method to get it. What i want to get is the the class of SpatialPolygons(sp) based on the generated Voronoi tessellation. Any ideas or suggestions are greatly appreciated. 2009/11/27 rusers.sh