Julia で円をモチーフにした格子模様を描く（5）
plotter.jl を include
https://blog.goo.ne.jp/r-de-r/e/bd71a52a09801335d56f7c47d879bfe3
include("plotter.jl")
function en5(nx=6, ny=5; r=1, factor=0.985, width=600, height=400)
function unit(x, y, r)
deg = 180:0.25:270
xs = zeros(length(deg), 2)
ys = similar(xs)
for (i, θ) = enumerate(deg)
xs[i, 1] = cosd(θ)factor + r
xs[i, 2] = cosd(θ-180)factor
ys[i, 1] = sind(θ)factor + r
ys[i, 2] = sind(θ-180)factor
end
xs = reshape(xs, :, 1)
ys = reshape(ys, :, 1)
xs2 = copy(xs)
ys2 = copy(ys)
for θ = 0:90:270 # degree
for i = 1:length(xs)
xs2[i], ys2[i] = xs[i]cosd(θ) + ys[i]sind(θ), xs[i]sind(θ) - ys[i]cosd(θ)
end
plotpolygon(x .+ xs2, y .+ ys2, lwd=0, col=:palegreen4, fcol=:palegreen4)
end
end
function cherryblossom(x, y, r)
R = 0.6r
deg = -asind(r/2R):0.5:asind(r/2R)
n = floor(Int, 0.9length(deg))
xs = zeros(n);
ys = similar(xs);
for i = 1:n
xs[i] = cosd(deg[i])R - cosd(asind(r/2R))R
ys[i] = sind(deg[i])R + 0.5r
end
xs = vcat(xs, -reverse(xs))
ys = vcat(ys, reverse(ys))
xs[n+1] = 0.5(xs[n] + xs[n+1])
ys[n+1] = 0.85r
xs2 = copy(xs)
ys2 = copy(ys)
randdeg = rand(0:35)
for θ = 36+randdeg:72:324+randdeg # degree
sine, cosine = sincosd(θ)
for i = 1:length(xs)
xs2[i], ys2[i] = xs[i]cosine + ys[i]sine, xs[i]sine - ys[i]cosine
end
plotpolygon(x .+ xs2, y .+ ys2, col=:pink, fcol=:pink)
end
plotcircle(x, y, 0.08r, col=:white, fcol=:white)
end
x1, y1, x2, y2 = 0, r, 2nx, 2ny+r
println("(width, height) = ($(x2 - x1), $(y2 - y1))")
plotbegin(w=width, h=height)
plotlimit(x1, y1, x2, y2)
for y = 1:2ny+1
for x = 0:2nx
if y % 2 != x %2
unit(x, y, r)
cherryblossom(x, y, 0.4r)
end
end
end
plotend()
end
@time en5(5, 4, width=600, height=480)
savefig("en5.png")
