In [2]:

    

import scipy.io
import numpy
import os


    

In [21]:

    

percorsoFile = "/home/protoss/Documenti/TESI/thesis/codici/matlabbo/griglia4096.mat"

griglia = scipy.io.loadmat(percorsoFile)['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))


    

    




(71, 2)

In [17]:

    

from matplotlib import pyplot
%matplotlib inline
#pyplot.figure(figsize=(100, 30))
a = pyplot.scatter(image4096[:,0],image4096[:,1], s = 1)


    

In [36]:

    

percorsoFile = "/home/protoss/Documenti/TESI/thesis/codici/matlabbo/griglia8192.mat"

griglia = scipy.io.loadmat(percorsoFile)['griglia']

lat = griglia[:,0]
lon = griglia[:,1]

def faiquadrato(l,b):
    lmax = 10
    lmin = -10
    bmax = 10
    bmin = -10

    coords = numpy.stack((lat,lon),1)
    
    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[:,0], coords[:,1]

image8192 = faiquadrato(lat,lon)
print(numpy.shape(image8192))


    

    




(2, 231)

In [23]:

    

from matplotlib import pyplot
%matplotlib inline
#pyplot.figure(figsize=(100, 30))
a = pyplot.scatter(image8192[:,0],image8192[:,1], s = 1)


    

In [ ]: