📑 File Exchange Workout 01. Numpy Interpolation

file_exchange_01.html

Numpy Interpolation

xxxxxxxxxx

 

1

import warnings,numpy as np

2

from skimage import io,transform

3

warnings.filterwarnings('ignore')

4

def load_img(file_path,file_name,img_size):

5

    img=io.imread(file_path+file_name)

6

    img=transform.resize(img,[img_size,img_size])

7

    return img[:,:,:3]/255

8

def interpolate_hypersphere(v1,v2,steps):

9

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

10

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

11

    for step in range(steps):

12

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

13

        interpolated_norm=np.linalg.norm(interpolated)

14

        interpolated_normalized=\

15

        interpolated*(v1norm/interpolated_norm)

16

        vectors.append(interpolated_normalized)

17

    return np.array(vectors)*255

18

@interact

19

def _vector(file1=[1,2,..,5],file2=[11,12,..,15],

20

            steps=list(range(20,100,10)),img_size=[300,350,400]):

21

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

22

              'OlgaBelitskaya/data/main/decors/'

23

    file_name1='00_00_%03d'%(file1)+'.png'

24

    file_name2='00_00_%03d'%(file2)+'.png'

25

    img1=load_img(file_path,file_name1,img_size)

26

    img2=load_img(file_path,file_name2,img_size)

27

    imgs=np.vstack([interpolate_hypersphere(img1,img2,steps),

28

                    interpolate_hypersphere(img2,img1,steps)])

29

    animate([matrix_plot(

30

        imgs[i],frame=False,transparent=True,figsize=4) 

31

             for i in range(2*steps)]).show()

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

---

Messages