In [2]:

    

import scipy.io
import numpy
import os


    

In [2]:

    

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

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

    filtro = numpy.where(coords[:,0]<lmax)
    coords = coords[filtro]
    filtro = numpy.where(coords[:,0]>lmin)
    coords = coords[filtro]

    filtro = numpy.where(coords[:,1]>bmin)
    coords = coords[filtro]
    filtro = numpy.where(coords[:,1]<bmax)
    coords = coords[filtro]
    return coords

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


    

    




(78, 5)

In [2]:

    

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

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

    filtro = numpy.where(coords[:,0]<lmax)
    coords = coords[filtro]
    filtro = numpy.where(coords[:,0]>lmin)
    coords = coords[filtro]

    filtro = numpy.where(coords[:,1]>bmin)
    coords = coords[filtro]
    filtro = numpy.where(coords[:,1]<bmax)
    coords = coords[filtro]
    return coords

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


    

    




(254, 5)

In [ ]: