1

私は、通常の順序とは逆の (ラップ) ダイヤモンド スクエア フラクタルを生成する Python スクリプトに取り組んでいます。最終値が依存するポイントのみを生成します。各ポイントの処理方法を決定するために作成した関数は次のとおりです。

def stepinfo(x, y):
    y1 = y-(y&(y-1))
    x1 = x-(x&(x-1))
    sum = x1 + y1
    return [((sum&(sum-1)) == 0), min(y1, x1)]

返される最初の値は、点をひし形または正方形のステップのどちらで生成するかを指定するブール値です。2 番目の値は、ステップのグリッド サイズを指定します。x または y が 0 の場合を除いて、これはうまく機能します。返されるステップ サイズは 0 で、スクリプトはループします。しかし、フラクタルを 2x2 グリッドに並べて配置し、タイルの右下のコピー (つまり、ウィンドウのゼロ軸から離れた場所) をクリックしてポイントを指定すると、フラクタル全体が完全に生成されます。しかし、タイルサイズの倍数をマウス座標に追加したり、それらの線に沿って何かを追加したりして、この結果を「偽造」することはできないようです。

助言がありますか？

完全なスクリプトは次のとおりです (x または y = 0 の場合に再帰を停止することで無限ループを回避しましたが、不足しているポイントでアーティファクトが生成されます)。

MAPEXPONENT = 8

import pygame
import random

MAPSIZE = 2**int(MAPEXPONENT)
MAP = [[[None] for col in range(MAPSIZE)] for row in range(MAPSIZE)]

def point(x, y):
    return MAP[x%MAPSIZE][y%MAPSIZE]

def displace(value, scale):
    displaced = (value + (random.random()-.5) * scale * 12)
    if displaced > 1 or displaced < 0: displaced = 1-displaced%1
    return displaced

def average(values):
    realvalues = []
    for value in values:
        if value <> None: realvalues.append(value)
    if realvalues: return float(sum(realvalues, 0)) / len(realvalues)
    else: return 0

def color_point(x, y, array):
    pixel = point(x, y)[0]
    green = pixel*255
    red = pixel**.5*255
    blue = 255-green
    array[x%MAPSIZE][y%MAPSIZE] = [min(255, max(0, int(red))),
                                   min(255, max(0, int(green))),
                                   min(255, max(0, int(blue)))]

def stepinfo(x, y):
    y1 = y-(y&(y-1))
    x1 = x-(x&(x-1))
    sum = x1 + y1
    return [((sum&(sum-1)) == 0), min(y1, x1)]

def makepoint(x, y, map_array):
    midpoint = point(x, y)
    if midpoint[0] == None:
        step = stepinfo(x, y)
        size = step[1]
        scale = float(size)/MAPSIZE
        heights = []
        if step[0]: points = [[x-size, y-size], [x+size, y-size],
                              [x-size, y+size], [x+size, y+size]]
        else:       points = [[x, y-size],      [x, y+size],
                              [x-size, y],      [x+size, y]]
        for p in points:
            if point(p[0], p[1])[0] == None and p[0] and p[1]:
                makepoint(p[0], p[1], map_array)
            heights.append(point(p[0], p[1])[0])
        midpoint[0] = displace(average(heights), scale)
        color_point(x, y, map_array)
    return midpoint

pygame.init()
screen = pygame.display.set_mode((2*MAPSIZE, 2*MAPSIZE))
map_surface = pygame.surface.Surface((MAPSIZE, MAPSIZE))
brush = 4+MAPSIZE/32

map_surface_array = pygame.surfarray.pixels3d(map_surface)
MAP[0][0][0] = random.random()
color_point(0, 0, map_surface_array)
del map_surface_array

running = True
while running:
    event = pygame.event.wait()
    if event.type == pygame.MOUSEBUTTONDOWN:
        while event.type <> pygame.MOUSEBUTTONUP:
            event = pygame.event.wait()
            x0, y0 = pygame.mouse.get_pos()
            map_surface_array = pygame.surfarray.pixels3d(map_surface)
            for x in range(x0-brush, x0+brush):
                for y in range(y0-brush, y0+brush):
                    if (x-x0)**2+(y-y0)**2 < brush**2:
                        makepoint(x, y, map_surface_array)
            del map_surface_array
            for x, y in [(0,0), (MAPSIZE,0), (0,MAPSIZE), (MAPSIZE,MAPSIZE)]:
                screen.blit(map_surface, (x, y))
            pygame.display.update()
    elif event.type == pygame.QUIT:
        running = False
pygame.quit()
4

1 に答える 1

1

http://gameprogrammer.com/fractal.html#diamondの説明に基づいて、問題は、x==0 または y==0 の場合に正方形とダイヤモンドで間違った決定を下すことだと思います。

正しい決定を次のようにコーディングできると思います。

def stepinfo(x, y):
    y1 = y-(y&(y-1))
    x1 = x-(x&(x-1))
    return [x1==y1, min(y1, x1)]

または次のようにさらに簡単です：

def stepinfo(x, y):
    y1 = y&-y
    x1 = x&-x
    return [x1==y1, min(y1, x1)]
于 2012-01-11T20:12:37.050 に答える