動かす

• ロゴ(こちら)…(つい、自分もやってしまう)
• ImageJを使う(こちら)

size<-100
z<-matrix(0,size,size)

z[,1:size*0.1]<-1
for(i in (size*0.1+1):(size*0.3)){
 z[(i-size*0.1):(i+size*0.15),i]<-1
}
for(i in (size*0.3+1):(size*0.5)){
 z[(size*0.5-i):(size*0.75-i),i]<-1
}
z[1:size*0.5,(size*0.5+1):(size*0.6)]<-1
z[1:(size*0.5),(size*0.72):(size*0.83)]<-1
image(z)

Niter<-100
for(i in 1:Niter){
 #tmpz<-z*matrix(runif(size^2),size,size)
 tmpz<-z*cos(i/50*2*pi)
persp(tmpz,col="green",theta=i*5,phi=80,shade=0.3,zlim=c(0,1))
 #plot3d(row(z),col(z),z,col="green")

 file.name=paste("test",100+i,".png",sep="")
	png(file=file.name)
	persp(tmpz,col="green",theta=i*5,phi=80,shade=0.3,zlim=c(0,1))

	dev.off()
}