Artificial Data & H5 Storing

xxxxxxxxxx

 

1

import warnings; warnings.filterwarnings('ignore')

2

import imageio,numpy as np,pandas as pd

3

import os,h5py,seaborn as sn,pylab as pl

4

from skimage.transform import resize

5

from IPython.display import HTML

6

def randintcoord(img_size_out,img_size=1024):

7

    a=(.5+.1**6*np.random.randint(1,999999))*\

8

      np.random.choice([-1,1],1)[0]

9

    b=np.random.randint(3,10)

10

    c=.1**3*np.random.randint(1,99)*\

11

      np.random.choice([-1,1],1)[0]

12

    t=np.arange(0,12*np.pi,1/(2880*b))

13

    fx=np.sin(t/6)+a*np.sin(b*t)*np.cos(t)-\

14

       c*np.sin(16*b*t)

15

    fy=np.cos(t/6)+a*np.sin(b*t)*np.sin(t)-\

16

       c*np.cos(16*b*t)

17

    fx=.951*(fx-1.051*fx.min())/(fx.max()-fx.min())

18

    fy=.951*(fy-1.051*fy.min())/(fy.max()-fy.min())

19

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

20

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

21

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

22

    img=np.ones((img_size,img_size,3))

23

    randcol=.9-.8*np.random.random(3)

24

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

25

    img=resize(img,(img_size_out,img_size_out))

26

    return img,np.around(a,6),b,np.around(c,3),randcol

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

---

Messages

xxxxxxxxxx

 

1

img_size=256; num_images=10

2

images=np.zeros((num_images,img_size,img_size,3),dtype=np.float32)

3

labels=np.zeros((num_images,),dtype=np.int32)

4

targets=np.zeros((num_images,2),dtype=np.float32)

5

colors=np.zeros((num_images,3),dtype=np.float32)

6

for i in range(num_images):

7

    print('=>',end='',flush=True)

8

    img,a,b,c,col=randintcoord(img_size)

9

    images[i,:,:,:]=img

10

    labels[i],targets[i,0],targets[i,1]=b-3,a,c

11

    colors[i,:]=col

12

h5f='ArtificialImages'+str(img_size)+'.h5'

13

with h5py.File(h5f,'w') as f:

14

    f.create_dataset('images',data=images,compression='gzip')

15

    f.create_dataset('labels',data=labels,compression='gzip')

16

    f.create_dataset('targets',data=targets,compression='gzip')

17

    f.create_dataset('colors',data=colors,compression='gzip')

18

    f.close()

19

print('\nfile size: %s'%list(os.stat(h5f))[6])

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

---

Messages

xxxxxxxxxx

 

1

with h5py.File(h5f,'r') as f:

2

    keys=list(f.keys()); print(keys)

3

    images=np.array(f[keys[1]]); labels=np.array(f[keys[2]])

4

    targets=np.array(f[keys[3]]); colors=np.array(f[keys[0]])

5

    f.close()

6

fig=pl.figure(figsize=(6,2))

7

for i in range(2):

8

    randi=np.random.randint(1,num_images,2)

9

    ax=fig.add_subplot(1,2,i+1)

10

    pl.imshow(images[randi[i]]); pl.axis('off')

11

    ti=str([labels[randi[i]],list(targets[randi[i]])])

12

    pl.title(ti,color=colors[randi[i]])

13

pl.tight_layout(); pl.show()

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

---

Messages

xxxxxxxxxx

 

1

def interpolate_hypersphere(v1,v2,steps):

2

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

3

    v2normalized=v2*(v1norm/v2norm)

4

    vectors=[]

5

    for step in range(1,steps+1):

6

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

7

        interpolated_norm=np.linalg.norm(interpolated)

8

        interpolated_normalized=interpolated*(v1norm/interpolated_norm)

9

        vectors.append(interpolated_normalized)

10

    return np.array(vectors)

11

steps=60; file_name='pic.gif'

12

imgs=np.vstack(

13

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

14

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

15

imageio.mimsave('pic.gif',imgs)

16

file_path='https://sagecell.sagemath.org/kernel/'

17

file_path+=os.getcwd()[14:]+'/files/'

18

s1='<div id="imgs"><img src="'

19

s2='" height="400" width="400"></img></div>'

20

HTML(s1+file_path+file_name+s2)

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

---

Messages