２０２２年３月の趣味コーディング

今月趣味で書いたプログラムの、まとめ！

今月はpythonのグラフ描画ライブラリ（matplotlib）を使って、１ツイートでジェネラティブアートを作る遊びに挑戦していました。

きっかけは、今月の始めころに見かけたこちらのツイートでした。

#rstats #maths #creativecoding
library(tidyverse)
seq(-3,3,by=.01) %>%
expand.grid(x=., y=.) %>%
ggplot(aes(x=(1-x-sin(y^2)), y=(1+y-cos(x^2)))) +
geom_point(alpha=.05, shape=20, size=0)+
theme_void()+
coord_polar() pic.twitter.com/JyaMXHgPI9

— Antonio Sánchez Chinchón (@aschinchon) March 2, 2022

何かの言語を使って、１ツイートでアートが作られている様子。いったい何の言語だろうと調べてみたら、統計処理等でよく使われるR言語のようでした。

統計用のグラフ描画ライブラリでアート、という発想がとても新鮮で面白かったので、普段使っているpythonでも同じようなことができないかな？と挑戦してみました。

pythonのグラフ描画ライブラリのmatplotlibを使用して…

import matplotlib.pyplot as P
import numpy as N
F=P.figure(figsize=(9,9))
X,Y=N.meshgrid(*[N.linspace(-5,5,98)]*2)
Q=abs(X)+abs(Y)
P.streamplot(X,Y,N.cos(Q),N.sin(Q),density=4,color=-Q,arrowsize=0)
P.axis(False)https://t.co/qfBCaCqV8z()
F.savefig("img.png")#つぶやきmatplotlib pic.twitter.com/9g66w94HW4

— ryotako (@ryotakob) March 4, 2022

できました！

これは適当なベクトル場に対して流線図を表示したもの。アートっぽくするために、流線の矢印マークは非表示にしています。matplotlibを使って、１ツイートでもジェネラティブアートっぽい画像を作ることができました。

import matplotlib.pyplot as P
F=P.figure(0,(4,4))
s="A"
for _ in s*6:s=s.translate({65:"B-A-B",66:"A+B+A"})
X,Y,q,z=[],[],1,0
for c in s:
if c=="+":q*=.5+.866j
elif c=="-":q*=.5-.866j
else:z+=q;X+=[z.real];Y+=[z.imag]
P.plot(X,Y)https://t.co/qfBCaCqV8z()#つぶやきmatplotlib pic.twitter.com/9RW3EyimwG

— ryotako (@ryotakob) March 6, 2022

さらにその後、L-systemを使ってフラクタル図形（シェルピンスキーのギャスケット）っぽいプロットを作ってみたり、

import matplotlib.pyplot as P
import numpy as N
F=P.figure(0,(4,4))
X,Y=N.meshgrid(*[N.linspace(-4,4,999)]*2)
S,T=abs(X%2-1),abs(Y%2-1)
S,T=N.maximum(S,T),N.minimum(S,T)
Z=N.sin(N.log(S)-N.sin(9*T+7*S))%1
P.contourf(X,Y,Z)
P.axis(1==0)https://t.co/qfBCaCqV8z()#つぶやきmatplotlib pic.twitter.com/Tn0mpsBuwv

— ryotako (@ryotakob) March 7, 2022

非線形関数のコンター図（等高線図）を描くことで模様を作ってみたり、

import matplotlib.pyplot as P
F=P.figure(0,(4,4))
X,Y,Z,S=[],[],[],[]
for i in range(999):
x,y=(i**.4)%4,(i*.7)%4;z=1
for a,b,c in zip(X,Y,Z):z=min(z,(((x-a)**2+(y-b)**2)**.5-c))
X+=[x];Y+=[y];Z+=[z];S+=[z*z*9900]
P.scatter(X,Y,S,Y)https://t.co/qfBCaCqV8z()#つぶやきmatplotlib pic.twitter.com/dEYxzxLpP0

— ryotako (@ryotakob) March 8, 2022

散布図の点の大きさを調整することでサークルパッキングっぽいグラフを作ったりしました！

いくつか作ったところでネタ切れになってきたので現在は休止中。そのうちネタが沸いたらまた挑戦するかもしれない。

以上、今月の趣味コーディングでした！