Julia set: Difference between revisions

Line 1,451:

[https://imgur.com/a/o3j9Zzn Example output.]

<lang python>"""

Solution from:

This solution is an improved version of an efficient Julia set solver

https://codereview.stackexchange.com/questions/210271/generating-julia-set

from:

'Bauckhage C. NumPy/SciPy Recipes for Image Processing:

Creating Fractal Images. researchgate. net, Feb. 2015.'

"""

from functools import partial

Line 1,493:

Line 1,491:

:return: 2D pixels array of intensities

"""

⚫

im, re = np.ogrid[min_coordinate.imag: max_coordinate.imag: height * 1j,

imaginary_axis, real_axis = np.ogrid[

min_coordinate.~~imag~~: max_coordinate.~~imag~~: ~~height~~ * 1j,

min_coordinate.real: max_coordinate.real: width * 1j]

z = (re + 1j * im).flatten()

⚫

min_coordinate.~~real~~: max_coordinate.~~real~~: ~~width~~ * 1j]

complex_plane = real_axis + 1j * imaginary_axis

~~result~~ = np.~~ones~~(~~complex_plane~~.shape)

live, = np.indices(z.shape) # indexes of pixels that have not escaped

iterations = np.empty_like(z, dtype=int)

for _ in range(iterations_count):

for i in range(iterations_count):

~~mask~~ = ~~np.abs(complex_plane)~~ <= ~~threshold~~

z_live = z[live] = mapping(z[live])

if ~~not~~ ~~mask.any~~():

escaped = abs(z_live) > threshold

iterations[live[escaped]] = i

live = live[~escaped]

if live.size == 0:

break

else:

complex_plane[mask] = mapping(complex_plane[mask])

~~result~~[~~~mask~~] += 1

iterations[live] = iterations_count

return ~~result~~

return iterations.reshape((height, width))

Line 1,521:

Line 1,522:

plt.axis('off')

plt.imshow(image,

cmap='~~nipy_spectral~~',

cmap='nipy_spectral_r',

origin='lower')

plt.show()