# NOT RUN {
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 kver -> SpatialPolygons
##
hr <- getverticeshr(kernelUD(puechabon$locs[,c("X","Y")], puechabon$locs$Name,
grid=100))
plot(hr)
class(hr)
spo <- kver2spol(hr)
plot(spo)
class(spo)
#########################################
##
## Conversion area -> SpatialPolygons
##
data(elec88)
ar <- as.area(elec88$area)
plot(ar)
toto <- area2spol(ar)
plot(toto)
#########################################
##
## Conversion SpatialPolygons -> area
##
## First create an object of class "SpatialRingsDataFrame"
tutu <- SpatialPolygonsDataFrame(toto, elec88$tab)
## and then conversion:
coincoin <- spol2area(tutu)
plot(coincoin)
## gets the attributes
haha <- attpol2area(tutu)
area.plot(coincoin, values = haha$Waechter)
#########################################
##
## Conversion ltraj -> SpatialPointsDataFrame
##
data(puechcirc)
plot(puechcirc)
toto <- ltraj2spdf(puechcirc)
plot(toto)
#########################################
##
## Conversion ltraj -> SpatialLinesDataFrame
##
toto <- ltraj2sldf(puechcirc)
plot(toto)
}
# }
Run the code above in your browser using DataLab