📑 File Exchange Workout 02. Python2R CSV

file_exchange_02.html

Python2R CSV

xxxxxxxxxx

 

1

import numpy as np,pandas as pd; from skimage import io,transform

2

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

3

          'OlgaBelitskaya/data/main/patterns/'

4

file='01_01_001.png'; img_size_out=256; fig_size=4

5

def create_coord(file,file_path=file_path,fig_size=fig_size):

6

    img=io.imread(file_path+file)

7

    matrix_plot(img/255).show(figsize=(fig_size,fig_size))

8

    imgbw=np.ones(img.shape[:2])

9

    imgbw[img[:,:,1]<int(128)]=0

10

    coord=np.array(np.where(imgbw<1)).T

11

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

12

    fx=.05+.9*(fx-fx.min())/(fx.max()-fx.min())

13

    fy=.05+.9*(fy-fy.min())/(fy.max()-fy.min())

14

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

15

    list_plot(coord,color='black',size=1)\

16

    .show(figsize=(fig_size,fig_size),flip_y=True,frame=True)

17

    return coord

18

coord=create_coord(file)

19

def create_image(coord,img_size_out,img_size=4096):

20

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

21

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

22

                dtype='int32')

23

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

24

                dtype='int32')

25

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

26

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

27

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

28

    img=transform.resize(img,(img_size_out,img_size_out))

29

    matrix_plot(img).show(figsize=(fig_size,fig_size))

30

    return img

31

img=create_image(coord,img_size_out)

32

np.savetxt('coord.csv',coord,delimiter=',',

33

           header='x,y',comments='')

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

---

Messages

xxxxxxxxxx

 

1

%%r

2

svg(filename='Rplots.svg',width=5,height=5.5,

3

    onefile=T,family='times',

4

    bg='black',pointsize=10,

5

    antialias=c('default','none','gray','subpixel'))

6

data<-read.csv('coord.csv')

7

matplot(data[,1],1-data[,2],type='p',cex=.3,pch='o',

8

        col='white',axes=FALSE,xlab='',ylab='')

9

dev.off()

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

---

Messages