提交: 0215a98b26e7b59f4203ce9f60b7e9e4edb9f643 - myprojects-hg-reborn (hg)

Commit MetaInfo

• 修订版: 0215a98b26e7b59f4203ce9f60b7e9e4edb9f643 (tree)
• 时间: 2023-02-07 08:17:57
• 作者: Lorenzo Isella <lorenzo.isella@gmai...>
• Commiter: Lorenzo Isella

Log Message

A function to create artwork.

更改概述

• add: R-codes/artwork_function.R

差异

diff -r 626baeef9bc4 -r 0215a98b26e7 R-codes/artwork_function.R
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/R-codes/artwork_function.R	Tue Feb 07 00:17:57 2023 +0100
@@ -0,0 +1,91 @@
+rm(list=ls())
+
+library(data.table)
+library(scattermore)
+library(viridis)
+
+makeStreams <- function(startpos,m=100,FUN,d=0.1){
+  pos=matrix(nrow=length(startpos), ncol=m)
+  pos2=matrix(nrow=length(startpos), ncol=m)
+  pos[,1] <- startpos
+  pos2[,m] <- startpos
+  for(i in 2:m){
+    v <- FUN(pos[,i-1])
+    pos[,i] <- pos[,i-1] + d * v/Mod(v)
+  }
+  for(i in (m-1):1){
+    v <- FUN(pos2[,i+1])
+    pos2[,i] <- pos2[,i+1] - d * v/Mod(v)
+  }
+  pos <-cbind(pos,pos2) |> as.vector()
+  cbind(Re(pos),Im(pos))
+}
+
+flowFunction <- function(z) {
+  fz = z-zeros[1]
+  for(i in zeros[-1]) fz = fz * (z-i)
+  for(i in poles) fz = fz / (z-i)
+  fz/Mod(fz)
+  }
+
+
+
+set.seed(1211)
+
+## polesandzeros <- seq(-1,1,l=40)-0.5i*seq(-1,1,l=40) +
+##     .7i*seq(-1,1,l=40)^2 - .3i
+
+
+polesandzeros <- seq(-1,1,l=40)-0.5i*seq(-1,1,l=40) ## +
+    ## .7i*seq(-1,1,l=40)^2 - .3i
+
+
+pole = rbinom(39, 1,0.2 ## 0.5
+              )
+poles = polesandzeros[pole==1]
+zeros = polesandzeros[pole==0]
+
+
+
+
+
+
+N=60000
+startZ <- c(
+  0.5*exp(2i * pi * seq(0,1,l=N/6)),
+  0.3*exp(2i * pi * seq(0,1,l=N/6)),
+  0.4*exp(2i * pi * seq(0,1,l=N/6)),
+  0.6*exp(2i * pi * seq(0,1,l=N/6)),
+  0.8*exp(2i * pi * seq(0,1,l=N/6)),
+  1.0*exp(2i * pi * seq(0,1,l=N/6))
+  )
+  
+pos <- makeStreams(startZ,m= 500
+                  ,FUN=flowFunction,d=0.0025 ## 0.005
+                   )
+
+## col=rep(c(hsv(0,  1,0.5,0.4),
+##           hsv(0.0,1,0.5,0.4),
+##           hsv(0.1,1,0.5,0.4),
+##           hsv(0.9,1,0.5,0.4),
+##           hsv(0.8,1,0.5,0.4),
+##           hsv(1,  1,0.5,0.1)
+##           ), each=N/6)
+
+
+col=rep(mako(10)[1:6], each=N/6)
+
+
+col <- rep(col,1000)
+
+png("my_artwork.png", width=5000, height=5000)
+par(mar=0.1*c(1,1,1,1), bg="#dddddd")
+scattermoreplot(pos,size=c(3000,3000),
+                xlim=c(-1,1), 
+                ylim=c(-1,1), 
+                col=col,
+                axes=F)
+dev.off()
+
+
+print("So far so good")