In [1]:

    

import scipy.io
import numpy
import os


def faiquadrato(coords):
    lmax = 10
    lmin = 350
    bmax = 10
    bmin = -10

    filtro = numpy.where(coords[:,1]>bmin)
    coords = coords[filtro]
    filtro = numpy.where(coords[:,1]<bmax)
    coords = coords[filtro]
    
    filtro = numpy.where(coords[:,0]<lmax)
    coordsUp = coords[filtro]
    filtro = numpy.where(coords[:,0]>lmin)
    coordsDown = coords[filtro]
    coords = numpy.concatenate((coordsDown,coordsUp))
    return coords


    

In [8]:

    

percorsoIn = "/home/protoss/Documenti/TESI/thesis/codici/Archiviati/matlabbo/griglia8192.mat"
percorsoOut = "/home/protoss/Documenti/TESI/thesis/codici/Archiviati/matlabbo/quadrato8192Gal.mat"
griglia = scipy.io.loadmat(percorsoIn)['griglia']

image8192 = faiquadrato(griglia)
print(numpy.shape(image8192))
scipy.io.savemat(percorsoOut, {'quadratoGal':image8192})


    

    




[ 96.38397315  96.38373421  96.84570724 ... 276.96105791 276.84897801
 276.38397315]

In [19]:

    

from matplotlib import pyplot

pyplot.scatter(image8192[:,0],image8192[:,1], s = 0.1)


    

    
Out[19]:





<matplotlib.collections.PathCollection at 0x7f397b3c76d8>






    

In [22]:

    

percorsoIn = "/home/protoss/Documenti/TESI/thesis/codici/Archiviati/matlabbo/griglia4096.mat"
percorsoOut = "/home/protoss/Documenti/TESI/thesis/codici/Archiviati/matlabbo/quadrato4096Gal.mat"
griglia = scipy.io.loadmat(percorsoIn)['griglia']
#print(griglia)

image4096 = faiquadrato(griglia)
print(numpy.shape(image4096))
scipy.io.savemat(percorsoOut, {'quadratoGal':image4096})


    

    




(804, 5)

In [23]:

    

from matplotlib import pyplot

pyplot.scatter(image4096[:,0],image4096[:,1], s = 0.1)


    

    
Out[23]:





<matplotlib.collections.PathCollection at 0x7f397b340668>