Description Usage Arguments Examples
This function weights points with fraction of points closest to them among all border points
1 | my.dcirc.measure(x, r1, r2 = 1)
|
x |
is a 0/1 matrix. |
r1 |
is a positive real value of width to measure circumference. |
r2 |
is a positive real value to construct border points |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 | library(igraph)
n <- 2^5
x <- y <- 1:n
xy <- as.matrix(expand.grid(x,y))
xy.val <- rep(0,length(xy[,1]))
n.circle <- 30
ctr <- matrix(sample(20:40,replace=TRUE,n.circle*2),ncol=2)
r <- sample(1:20,replace=TRUE,n.circle)
for(i in 1:n.circle){
tmp.x <- xy[,1] - ctr[i,1]
tmp.y <- xy[,2] - ctr[i,2]
s <- which(tmp.x^2+tmp.y^2 < r[i]^2)
xy.val[s] <- 1
}
xy.mat <- matrix(xy.val,ncol=n)
yyy <- my.dcirc.measure(xy.mat,10,1)
bdr <- yyy$border
bdr <- rbind(bdr,bdr[2,])
diff.bdr <- apply(bdr,2,diff)
#diff.bdr <- diff.bdr[-length(diff.bdr[,1]),]
diff.bdr.1 <- diff.bdr[-length(diff.bdr[,1]),]
diff.bdr.2 <- diff.bdr[-1,]
diff.bdr.1.cplx <- diff.bdr.1[,1] + diff.bdr.1[,2]*1i
diff.bdr.2.cplx <- diff.bdr.2[,1] + diff.bdr.2[,2]*1i
direc <- diff.bdr.2.cplx/diff.bdr.1.cplx
arg.direc <- Arg(direc)
cumsum.arg <- cumsum(arg.direc)
plot(cumsum.arg,type="l")
plot(yyy[[1]])
plot(yyy[[3]])
col <- yyy[[1]]
col <- (max(col)-col)/(max(col)-min(col))
col <- 0.2 + 0.6*col
plot(yyy[[3]],pch=20,col=gray(col))
loc <- cumsum(yyy[[1]])
plot(loc)
plot(loc,cumsum.arg)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.