INQUINAMENTO ACUSTICO

def main():
    # -------------SALVO IL FILE IN MEMORIA NEL PROGRAMMA----------

    INPUT = "input_file.csv"
    OUTPUT = "region_bad.txt"
    PRIMOANNO = 2008
    ULTIMOANNO = 2012
    valori = []

    try:
        with open(INPUT, "r") as file:
            for line in file:
                lista = line.strip().split(";")
                if lista[3].isnumeric():
                    valori.append(lista)

    except IOError:
        print("il programma non è riuscito ad aprire il file contenente i dati")
        exit()

    # ---------------CALCOLO LA MEDIA NAZIONALE ANNUA-------------

    medieNazionali = mediaAnnua(valori, PRIMOANNO, ULTIMOANNO)

    # ------------CALCOLO LA MEDIA REGIONALE ANNUA------------
    medieRegionali = []
    # ---------creo un insieme ausiliario per avere il nome delle regioni----

    regioni = set()
    for line in valori:
        regioni.add(line[2])

    for regione in sorted(regioni):
        valoriRegionali = []
        for dato in valori:
            if dato[2] == regione:
                valoriRegionali.append(dato)

        # creo una tabella regione-media
        mediaRegionale = mediaAnnua(valoriRegionali, PRIMOANNO, ULTIMOANNO)
        listaRegione = [regione, mediaRegionale, len(valoriRegionali)]
        medieRegionali.append(listaRegione)

    # ------------ SALVO IL NUMERO DI ANNI IN CUI LE SINGOLE REGIONI SONO STATE SOPRA E LO NORMALIZZO----------------
    maggNormalizzate = {}
    for valore in medieRegionali:
        maggNormalizzate[valore[0]] = sopraMedia(valore[1], medieNazionali) / valore[2]

    with open(OUTPUT, "w")as file:
        for regione in maggNormalizzate.keys():
            file.write(f"{regione} >>> {maggNormalizzate[regione]}\n")

    # ------------SALVO IL NUMERO DI ANNI IN CUI LE SINGOLE REGIONI SONO STATE SOPRA LA MEDIA E CREO LA LISTA DI REGIONI PIU VIRTUOSE-------------
    minNormalizzate = {}
    for valore in medieRegionali:
        minNormalizzate[valore[0]] = sottoMedia(valore[1], medieNazionali) / valore[2]

    #-----------trovo la regione con la media normalizzata più virtuosa----------------------------
    minorInquinamento = 0
    for regione in minNormalizzate.keys():

        if minNormalizzate[regione] > minorInquinamento:
            minorInquinamento = minNormalizzate[regione]

    #------------------aggiungo alla lista tutte le regioni più virtuose------------------
    migliori = []
    for regione in minNormalizzate.keys():
        if minNormalizzate[regione] == minorInquinamento:
            migliori.append(regione)


    #------------------aggiungo al file di testo le regioni virtuose-------------
    with open(OUTPUT, "a") as file:
        file.write(f"regione/i più virtuosa/e: {migliori}")

    with open(OUTPUT, "r") as file:
        for line in file:
            print(line)


# -------------------FUNZIONI--------------------------
def mediaAnnua(tabella, primoAnno, ultimoAnno):  # definisco la funzione per calcolare le medie annue
    medieAnnue = {}

    for anno in range(primoAnno, ultimoAnno + 1):
        media = 0
        counter = 0
        for linea in tabella:
            if int(linea[3]) == anno:
                media += float(linea[4])
                counter += 1
        media = media / counter
        medieAnnue[anno] = media
    return medieAnnue


def sopraMedia(medieRegione, medieNazione):
    counter = 0
    for anno in sorted(medieRegione.keys()):
        if medieRegione[anno] >= medieNazione[anno]:
            counter += 1
    return counter


def sottoMedia(medieRegione, medieNazione):
    counter = 0
    for anno in sorted(medieRegione.keys()):
        if medieRegione[anno] < medieNazione[anno]:
            counter += 1
    return counter


main()