+
+spirale <- function(g, weigth, center, miny=0.1) {
+ ncoord <- matrix(0, nrow=length(weigth)+1, ncol=2)
+ v.names <- V(g)$name
+ center.name <- v.names[center]
+ first <- which.max(weigth)[1]
+ if (head_of(g, first)$name == center.name) {
+ n.name <- tail_of(g, first)
+ } else {
+ n.name <- head_of(g, first)
+ }
+ n.name <- n.name$name
+ nb <- length(weigth)
+ ncoord[which(v.names==n.name),] <- c(0,1)
+ weigth[first] <- 0
+ rs <- norm.vec(weigth,1, miny)
+ nbt <- nb %/% 50
+ if (nbt == 0) nbt <- 1
+ angler <- ((360 * nbt) / (nb- 1)) * (pi/180)
+ ang <- 90 * (pi/180)
+ rr <- (1-miny) / (nb-1)
+ r <- 1
+ while (max(weigth != 0)) {
+ first <- which.max(weigth)[1]
+ if (head_of(g, first)$name == center.name) {
+ n.name <- tail_of(g, first)
+ } else {
+ n.name <- head_of(g, first)
+ }
+ n.name <- n.name$name
+ #r <- rs[first]
+ r <- r - rr
+ ang <- ang + angler
+ x <- r * cos(ang)
+ y <- r * sin(ang)
+ weigth[first] <- 0
+ ncoord[which(v.names==n.name),] <- c(x,y)
+ }
+ ncoord
+}
+
+spirale3D <- function(g, weigth, center, miny=0.1) {
+ ncoord <- matrix(0, nrow=length(weigth)+1, ncol=3)
+ v.names <- V(g)$name
+ center.name <- v.names[center]
+ first <- which.max(weigth)[1]
+ if (head_of(g, first)$name == center.name) {
+ n.name <- tail_of(g, first)
+ } else {
+ n.name <- head_of(g, first)
+ }
+ n.name <- n.name$name
+ nb <- length(weigth)
+ ncoord[which(v.names==n.name),] <- c(0,0,1)
+ weigth[first] <- 0
+ rs <- norm.vec(weigth,1, miny)
+ nbt <- nb %/% 50
+ if (nbt == 0) nbt <- 1
+ angler <- ((360 * nbt) / (nb- 1)) * (pi/180)
+ theta <- 0
+ phi <- 90 * (pi/180)
+ rr <- (1-miny) / (nb-1)
+ r <- 1
+ while (max(weigth != 0)) {
+ first <- which.max(weigth)[1]
+ if (head_of(g, first)$name == center.name) {
+ n.name <- tail_of(g, first)
+ } else {
+ n.name <- head_of(g, first)
+ }
+ n.name <- n.name$name
+ #r <- rs[first]
+ r <- r - rr
+ theta <- theta + angler
+ phi <- phi + angler/2
+ x <- r * sin(theta) * cos(phi)
+ y <- r * sin(theta) * sin(phi)
+ z <- r * cos(theta)
+ weigth[first] <- 0
+ ncoord[which(v.names==n.name),] <- c(x,y,z)
+ }
+ ncoord
+}
+