作者josephw (With or Without u)
看板NCTU-STAT95G
標題[分享] 測辨色能力
時間Fri May 23 00:13:06 2008
## 先把這個背景弄進去, 跟之前的一樣不過參照學姐意見改的比較簡單
rgb_plot <- function(r = 0.5, g=0.5, b = 0.5, times = 10)
{
x <- c(-20,20)
y <- c(-20,20)
ax <- min(x)
bx <- max(x)
ay <- min(y)
by <- max(y)
xcorint = 5
ycorint = 5
plot(x, y,type='n',axes=F, xlab="", ylab="")
i <- 1
while(r+i*0.005 < 1 && g+i*0.005 < 1 && b+i*0.005 < 1){
if(r == 0){
if(b == 0)
polygon(c(ax+(i-1)/times-xcorint,ax+(i-1)/times-xcorint,
bx-(i-1)/times+xcorint,bx-(i-1)/times+xcorint),
c(ay+(i-1)/times-ycorint,by-(i-1)/times+ycorint,
by-(i-1)/times+ycorint, ay+(i-1)/times-ycorint),
col=rgb(0,g+i*0.005,0),lty=0)
else if(g == 0)
polygon(c(ax+(i-1)/times-xcorint,ax+(i-1)/times-xcorint,
bx-(i-1)/times+xcorint,bx-(i-1)/times+xcorint),
c(ay+(i-1)/times-ycorint,by-(i-1)/times+ycorint,
by-(i-1)/times+ycorint, ay+(i-1)/times-ycorint),
col=rgb(0,0,b+i*0.005), lty=0)
else
polygon(c(ax+(i-1)/times-xcorint,ax+(i-1)/times-xcorint,
bx-(i-1)/times+xcorint,bx-(i-1)/times+xcorint),
c(ay+(i-1)/times-ycorint,by-(i-1)/times+ycorint,
by-(i-1)/times+ycorint, ay+(i-1)/times-ycorint),
col=rgb(0,g+i*0.005,b+ i*0.005), lty=0)
}
else if(b == 0){
if(g == 0)
polygon(c(ax+(i-1)/times-xcorint,ax+(i-1)/times-xcorint,
bx-(i-1)/times+xcorint,bx-(i-1)/times+xcorint),
c(ay+(i-1)/times-ycorint,by-(i-1)/times+ycorint,
by-(i-1)/times+ycorint5, ay+(i-1)/times-ycorint),
col=rgb(r+i*0.005,0,0), lty=0)
else
polygon(c(ax+(i-1)/times-xcorint,ax+(i-1)/times-xcorint,
bx-(i-1)/times+xcorint,bx-(i-1)/times+xcorint),
c(ay+(i-1)/times-ycorint,by-(i-1)/times+ycorint,
by-(i-1)/times+ycorint, ay+(i-1)/times-ycorint),
col= rgb(r+i*0.005,g+i*0.005,0), lty=0)
}
else if(g == 0)
polygon(c(ax+(i-1)/times-xcorint,ax+(i-1)/times-xcorint,
bx-(i-1)/times+xcorint,bx-(i-1)/times+xcorint),
c(ay+(i-1)/times-ycorint,by-(i-1)/times+ycorint,
by-(i-1)/times+ycorint, ay+(i-1)/times-ycorint),
col= rgb(r+i*0.005,0,b+ i*0.005), lty=0)
else
polygon(c(ax+(i-1)/times-xcorint,ax+(i-1)/times-xcorint,
bx-(i-1)/times+xcorint,bx-(i-1)/times+xcorint),
c(ay+(i-1)/times-ycorint,by-(i-1)/times+ycorint,
by-(i-1)/times+ycorint, ay+(i-1)/times-ycorint),
col= rgb(r+i*0.005,g+i*0.005,b+ i*0.005), lty=0)
i <- i+1
}
}
####Example:
#######################辨色能力(猜一個位數)###################################
rgb_plot(r=0.1,b=0,g=0.1)
par(new=T)
plot(runif(n=10000,min=-2,max=2),runif(n=10000,min=-2,max=2),xlim=c(-5,5),
ylim=c(-5,5),xlab="",ylab="",col=rainbow(100))
for(i in 1:21){
points(0,1.5-(i-1)/10, col=rgb(0.1+i/100,0.7+i/100,0.7+i/100))
points(0.05,1.5-(i-1)/10, col=rgb(0.001+2*i/50,0.7+i/100,0.7+i/300))
points(-0.1,1.5-(i-1)/10, col=rgb(0.1,0.5,1))
}
for(i in 1:11){
points(-0.5+(i-1)/10,-0.49,col=rgb(0.1+i/80,0.7+i/100,0.7+i/80))
points(-0.5+(i-1)/10,-0.5, col=rgb(0.1+2*i/100,0.7+i/200,0.7+i/100))
points(-0.5+(i-1)/20,-0.51, col=rgb(0.1+2*i/100,0.7+i/200,0.7+i/100))
}
for(i in 1:11){
points(-0.1-(i-1)/30,1.5-(i-1)/20,col=rgb(0.1+i/100,0.7+i/100,0.7+i/100))
points(-0.115-(i-1)/30,1.5-(i-1)/20,col=rgb(0.1+i/100,0.7+i/200,0.7+i/150))
}
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.113.191.14
1F：→ josephw:答案是 1喔XD 05/23 00:27
2F：推 Y0SHIKI:聽說有人說是4 :D 05/23 00:35
3F：推 ym7226:這看久眼睛會花吧 有點像7的1吧@@ 05/23 00:57
4F：→ josephw:黃翁賢一次猜中XD 05/23 09:59
5F：推 mrliang:我我我 我說四 05/24 00:03
6F：→ mrliang:結果發現我才是色盲 嗚嗚 嗡嗡對不起 05/24 00:03
7F：→ josephw:XD.....你最近太累了喔!! 05/24 00:24