📑 Circus of Programming Tricks 01. Coordinate Interpolation

tricks_01.html

xxxxxxxxxx

 

1

import numpy as np; from skimage import io,transform

2

file_path='https://raw.githubusercontent.com/'+\

3

          'OlgaBelitskaya/data/main/patterns/'

4

file1,file2='01_01_001.png','01_01_002.png'

5

def preprocess_img(file1,file2,file_path=file_path):

6

    img1,img2=io.imread(file_path+file1),io.imread(file_path+file2)

7

    imgbw1,imgbw2=np.ones(img1.shape[:2]),np.ones(img2.shape[:2])

8

    imgbw1[img1[:,:,1]<int(128)],imgbw2[img2[:,:,1]<int(127)]=0,0

9

    coord1=np.array(np.where(imgbw1<1)).T; coord=coord1

10

    coord2=np.array(np.where(imgbw2<1)).T

11

    if coord1.shape[0]>coord2.shape[0]: coord1,coord2=coord2,coord

12

    while not coord2.shape==coord1.shape:

13

        randi=np.random.randint(

14

            0,coord2.shape[0]-1,coord2.shape[0]-coord1.shape[0])

15

        coord2=np.delete(coord2,list(randi),axis=0)

16

    return coord1,coord2

17

coord1,coord2=preprocess_img(file1,file2)

18

def create_coords(coord1,coord2):

19

    coords=[]

20

    for el in [coord1,coord2]:

21

        fx=el[:,0]; fy=el[:,1]

22

        fx=.025+.95*(fx-fx.min())/(fx.max()-fx.min())

23

        fy=.025+.95*(fy-fy.min())/(fy.max()-fy.min())

24

        el=np.array([[fx[i],fy[i]] for i in range(len(fx))])

25

        coords.append(el)

26

    return np.array(coords)

27

coords=create_coords(coord1,coord2)

28

steps=20; img_size_out=256

29

def interpolate_hypersphere(v1,v2,steps):

30

    v1norm=np.linalg.norm(v1); v2norm=np.linalg.norm(v2)

31

    vectors=[]; v2normalized=v2*(v1norm/v2norm)

32

    for step in range(steps):

33

        interpolated=v1+(v2normalized-v1)*step/(steps-int(1))

34

        interpolated_norm=np.linalg.norm(interpolated)

35

        interpolated_normalized=\

36

        interpolated*(v1norm/interpolated_norm)

37

        vectors.append(interpolated_normalized)

38

    return np.array(vectors)

39

coords_int=np.vstack(

40

    [interpolate_hypersphere(coords[1],coords[0],steps),

41

     interpolate_hypersphere(coords[0],coords[1],steps)])

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages

xxxxxxxxxx

 

1

def create_images(coords_int,img_size_out,img_size=2048):

2

    imgs=[]

3

    for i in range(coords_int.shape[0]):

4

        fx=coords_int[i,:,0]; fy=coords_int[i,:,1]

5

        fx=np.array(np.clip(fx*img_size,0,img_size-1),

6

                    dtype='int32')

7

        fy=np.array(np.clip(fy*img_size,0,img_size-1),

8

                    dtype='int32')

9

        f=np.array([[fx[i],fy[i]] for i in range(len(fx))])

10

        img=np.zeros((img_size,img_size))

11

        for [x,y] in f: img[x,y]=1

12

        imgs.append(transform.resize(

13

            img,(img_size_out,img_size_out)))

14

    return np.array(imgs)

15

imgs=create_images(coords_int,img_size_out)

16

def mplot(imgs,i,fig_size=3):

17

    return matrix_plot(

18

        imgs[i],figsize=(fig_size,fig_size),frame=False) 

19

for i in [0,steps]: mplot(imgs,i).show()

20

animate([mplot(imgs,i) for i in range(2*steps)])

Language: Sage Gap GP HTML Macaulay2 Maxima Octave Python R Singular

---

Messages