imgsize = 500 # размер получаемого изображения
import math
def rhytmics_gen(cycles, X0, Y0, K, G = 0):
"""
"the Seraph's Rhythmics"
cycles - список окружностей состоящий из кортежей (Rn, Fn, Nn)
X0 - координата X главной опорной центральной точки первой окружности
Y0 - координата Y главной опорной центральной точки первой окружности
X1-X6 - координата X окружностей 1 - 6 с центрами в точках X0 - X5
Xn - координата X окружностей n с центрами в точках Xn-1
Y1-Y6 - координата Y окружностей 1 - 6 с центрами в точках Y0 - Y5
Yn - координата Y окружностей n с центрами в точках Yn-1
R1-R6-Rn - величина радиусов окружностей в пикселах
F1-F6-Fn - начальное угловое отклонение радиусов
N1-N6-Nn - количество оборотов радиусов за цикл
K - коэффициент частоты расчета следующей координаты
G - угловой поворот всего рисунка
"""
X = X0
Y = Y0
for i in range(2*K): # это организация циклического построения сдвигов отрезков графика
# корректный код закомментирован
#X = X0
#Y = Y0
for cycle in cycles:
Rn, Fn, Nn = cycle
X += Rn * math.cos(((i + Fn) * math.pi / K) * Nn + G)
Y += Rn * math.sin(((i + Fn) * math.pi / K) * Nn + G)
yield X, Y # отдать Xn, Yn
# создание изображения с черным фоном
from PIL import Image
image = Image.new("RGB", (imgsize, imgsize), (0,0,0))
pix = image.load()
# пример списка окружностей
cycles = [
(0.8, math.pi/10, 5),
(0.5, math.pi/15, -2),
(0.9, math.pi/10, 10),
(0.9, math.pi/2, -30),
]
# получение списка координат
points = list(rhytmics_gen(cycles, X0=0, Y0=0, K=1024))
# вычисление границ изображения
minx = 1.e10 # невозможные границы
miny = 1.e10
maxx = -1.e10
maxy = -1.e10
for p in points:
x, y = p
minx = min(minx, x)
miny = min(miny, y)
maxx = max(maxx, x)
maxy = max(maxy, y)
print "x:", minx, maxx
print "y:", miny, maxy
# коэффициент масштабирования координат
step = imgsize/(max(abs(minx), abs(maxx), abs(miny), abs(maxy))*2)
# отрисовка изображения
for p in points:
x, y = p # координата точки
# вычисление положение точки на изображении с учетом границ
pixx = int((x-minx)*step)
pixy = int((y-miny)*step)
# проверка границ
assert pixx>=0
assert pixy>=0
assert pixx<imgsize
assert pixy<imgsize
pix[pixx, pixy] = (255, 255, 255) # установка белой точки
# сохранение изображения
image.save("rhytmics001.png", "PNG")
