Nothing
spodtSpatialLines <- function(object,data){
x<- coordinates(data)[,1]
y<- coordinates(data)[,2]
if (class(object@racine)=="f.spodt"){#root is a leaf
d2<-matrix(ncol=4,nrow=4)
if (!is.null(x) | !is.null(y)){
minX <- min(x) - (max(x) - min(x))/50
maxX <- max(x) + (max(x) - min(x))/50
minY <- min(y) - (max(y) - min(y))/50
maxY <- max(y) + (max(y) - min(y))/50
d2 <- rbind(c(minX,minY,minX,maxY), c(minX,maxY,maxX,maxY), c(maxX,maxY,maxX,minY), c(maxX,minY,minX,minY))
}
colnames(d2)[c(3,4)] <- colnames(d2)[c(1,2)]
rownames(d2)<-1:dim(d2)[1]
}
else{
if((is.na(object@cl.grf)[1])){ #no graft
data.coord <- recursif.carte.sp(object@racine)
ind.na <- which(is.na(data.coord[,1]))
d2<- data.coord[-ind.na,]
if (!is.null(x) | !is.null(y)){
minX <- min(x) - (max(x) - min(x))/50
maxX <- max(x) + (max(x) - min(x))/50
minY <- min(y) - (max(y) - min(y))/50
maxY <- max(y) + (max(y) - min(y))/50
d2 <- rbind(d2, c(minX,minY,minX,maxY), c(minX,maxY,maxX,maxY), c(maxX,maxY,maxX,minY), c(maxX,minY,minX,minY))
}
colnames(d2)[c(3,4)] <- colnames(d2)[c(1,2)]
points.d2 <- rbind(d2[,1:2],d2[,3:4])
points.d2 <- unique(points.d2)
nbLig <- dim(points.d2)[1]
nbLig.d2 <- dim(d2)[1]
d2<-recurs.d2(points.d2,d2)
rownames(d2)<-1:dim(d2)[1]
}
else{ #graft
d2 <- object@brd[,1:4]
points.d2 <- rbind(d2[,1:2],d2[,3:4])
points.d2 <- unique(points.d2)
nbLig <- dim(points.d2)[1]
nbLig.d2 <- dim(d2)[1]
nd2<-recurs.d2(points.d2,d2)
rownames(nd2)<-1:dim(nd2)[1]
d1 <- object@sgmts.grf
d2 <- rm.dupl(d1,nd2)
d1_inv <- cbind(d1[,3:4], d1[,1:2])
d2 <- rm.dupl(d1_inv,d2)
}
}
spodtL <- list()
spodtL <- sapply(1:(dim(d2)[1]), function(i){
x<-rbind(d2[i,1:2],d2[i,3:4])
spodtL[[i]] <- Lines(list(Line(x)),ID=rownames(d2)[i])})
spodtSL<-SpatialLines(spodtL)
proj4string(spodtSL)<-proj4string(data)
return(spodtSL)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.