if (require(sp) {
#########################################
##
## Conversion kasc -> SpatialPixelsDataFrame
##
data(puechabon)
toto <- kasc2spixdf(puechabon$kasc)
image(toto)
summary(toto)
#### and conversely
toto <- spixdf2kasc(toto)
image(toto)
hist(toto)
data(meuse.grid)
m <- SpatialPixelsDataFrame(points = meuse.grid[c("x", "y")],
data = meuse.grid)
i <- spixdf2kasc(m)
image(i)
### conversion asc -> SpatialPixelsDataFrame
cuicui <- asc2spixdf(getkasc(toto,1))
image(cuicui)
#########################################
##
## Conversion area -> SpatialRings
##
data(elec88)
ar <- as.area(elec88$area)
plot(ar)
toto <- area2sr(ar)
plot(toto)
#########################################
##
## Conversion SpatialRings -> area
##
## First create an object of class "SpatialRingsDataFrame"
data(ncshp)
nc1 <- as.SpatialRings.Shapes(nc.shp$Shapes, as.character(nc.shp$att.data$FIPS))
plotSpatialRings(nc1)
df <- nc.shp$att.data
rownames(df) <- as.character(nc.shp$att.data$FIPS)
ncSRDF <- SpatialRingsDataFrame(nc1, df)
## and then conversion:
coincoin <- sr2area(ncSRDF)
## please note the warnings
plot(coincoin)
## gets the attributes
haha <- attsr2area(ncSRDF)
area.plot(coincoin, values = df$SID74/df$BIR74)
#########################################
##
## Conversion traj -> SpatialPointsDataFrame
##
data(puechcirc)
plot(puechcirc)
toto <- traj2spdf(puechcirc)
plot(toto)
#########################################
##
## Conversion traj -> SpatialLinesDataFrame
##
data(puechcirc)
plot(puechcirc)
toto <- traj2sldf(puechcirc)
plot(toto)
}
Run the code above in your browser using DataLab