On Tue, 26 Jan 2010, Agustin Lobo wrote:
Hi!
Package StatDA
provides background geoinformation as a set of coordinates:
str(kola.background)
List of 4
$ boundary:'data.frame': 50 obs. of 2 variables:
..$ V1: num [1:50] 388650 388160 386587 384035 383029 ...
..$ V2: num [1:50] 7892400 7881248 7847303 7790797 7769214 ...
$ coast :'data.frame': 6259 obs. of 2 variables:
..$ V1: num [1:6259] 438431 439102 439102 439643 439643 ...
..$ V2: num [1:6259] 7895619 7896495 7896495 7895800 7895542 ...
$ borders :'data.frame': 504 obs. of 2 variables:
..$ V1: num [1:504] 417575 417704 418890 420308 422731 ...
..$ V2: num [1:504] 7612984 7612984 7613293 7614530 7615972 ...
$ lakes :'data.frame': 6003 obs. of 2 variables:
..$ V1: num [1:6003] 547972 546915 NA 547972 547172 ...
..$ V2: num [1:6003] 7815109 7815599 NA 7815109 7813873 ...
is there any spatial function aready availale to convert these coordinates
into Spatial Lines and Spatial Polygons?
This seems to work at least for SpatialLines - for polygons, the rings may
need to be closed:
library(StatDA)
data(kola.background)
xy <- kola.background$boundary
names(xy) <- c("x", "y")
library(maptools)
bdy <- map2SpatialLines(xy)
plot(bdy)
The data are in the legacy S format (like the maps package).
Roger
Thanks
Agus
[[alternative HTML version deleted]]
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
--
Roger Bivand
Economic Geography Section, Department of Economics, Norwegian School of
Economics and Business Administration, Helleveien 30, N-5045 Bergen,
Norway. voice: +47 55 95 93 55; fax +47 55 95 95 43
e-mail: roger.biv...@nhh.no
_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@stat.math.ethz.ch
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
