Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
kurse:efcomputergrafik:knn [2020/03/24 12:12]
Simon Knaus
kurse:efcomputergrafik:knn [2020/03/24 14:55] (current)
Simon Knaus [$k$ Nearest Neighbors]
Line 1: Line 1:
-==== k Nearest Neighbors ====+==== $kNearest Neighbors ==== 
 +Bei $k$ nearest neighbours (kNN) geht es darum, einem dazukommenden Punkt diese Klasse zuzuweisen, welche die nächsten $k$-Punkte mehrheitlich haben.  
 +Der Trainingsdatensatz sind damit alle Punkte, von welchem man die Klasse und Koordinaten kennt. Auf Grund dieser Klassen und Koordinaten wird einem neuen Datenpunkt einzig auf Grund der Koordinaten eine Klasse zugeordnet. 
 + 
 +{{ :kurse:efcomputergrafik:knn:structure.png |}} 
 +=== Bemerkungen ==== 
 +  * Die untenstehenden Codes sind für TigerJython gedacht 
 +  * Die Codes sind nicht unbedingt schnell sondern illustrativ 
 +  * Es werden folgende Daten benötigt: 
 +    * {{kurse:efcomputergrafik:knn:knn_r2_trainingsdaten.zip|Trainingsdaten}} für zweidimensionales kNN 
 +    * {{kurse:efcomputergrafik:knn:knn_r2_testdaten.zip|Testdaten}} für zweidimensionales kNN 
 +    * ZIP-Code Daten 
 +      * Training: {{kurse:efcomputergrafik:knn:knn_rk_trainingsdaten.zip|Bilddateien als ZIP}} 
 +      * Test: {{kurse:efcomputergrafik:knn:knn_rk_testdaten.zip|Bilddateien als ZIP}} 
 +      * (Optional) {{kurse:efcomputergrafik:knn:knn_rk_training_csv.zip|Training als CSV}} 
 +      * (Optional) {{kurse:efcomputergrafik:knn:knn_rk_test_csv.zip|Test als CSV}}
  
 ==== Einführung ==== ==== Einführung ====
Line 14: Line 29:
  
  
-Wer $k$-nearest-neighbour implementiert hat, kann sich überlegen, wie die untenstehende Daten klassifiziert werden sollen: +Die Daten finden sich in einer {{kurse:efcomputergrafik:knn:knn_r2_trainingsdaten.zip|ZIP-Datei}}
-{{ :efinf:blcks2017:bigdata:l2:structure.png?400 |}} +
-Die Daten finden sich in einer neuen {{efinf:blcks2017:bigdata:l2:data2.zip|ZIP-Datein}}+
  
-<hidden>+<hidden Lösung>
 <file python knn.py> <file python knn.py>
 from gpanel import * from gpanel import *
Line 121: Line 134:
 === Ziele === === Ziele ===
  
-  * Klassifizierungsfehler auf {{efinf:blcks2017:bigdata:testdata.zip|Testdaten}} für verschiedene $k$ berechnen und Tabelle erstellen.+  * Klassifizierungsfehler auf {{kurse:efcomputergrafik:knn:knn_r2_testdaten.zip|Testdaten}} für verschiedene $k$ berechnen und Tabelle erstellen.
   * ZIP-Code Problematik verstehen: ZIP-Code -> Ziffer -> 16x16 Bild -> Liste mit 256 Graustufen-Werten -> kNN in $\mathbb{R}^{256}$.   * ZIP-Code Problematik verstehen: ZIP-Code -> Ziffer -> 16x16 Bild -> Liste mit 256 Graustufen-Werten -> kNN in $\mathbb{R}^{256}$.
-  * Einzelne Ziffern als {{efinf:blcks2017:bigdata:digits.zip|Grafikdatei}} einlesen und als 256 Zahlwerte pro Bild als Liste speichern:+  * Einzelne Ziffern als {{kurse:efcomputergrafik:knn:knn_rk_trainingsdaten.zip|Bilddateien als ZIP}} einlesen und als 256 Zahlwerte pro Bild als Liste speichern:
     * Eine Funktion schreiben, die als Argument einen Dateinamen hat und als Rückgabewert eine Liste mit 256 Elementen.     * Eine Funktion schreiben, die als Argument einen Dateinamen hat und als Rückgabewert eine Liste mit 256 Elementen.
     * Diese Funktion auf alle Dateien anwenden (siehe unten) und die Ziffer aus dem Dateinamen in eine Liste von Liste mit 256+1 Elementen speichern     * Diese Funktion auf alle Dateien anwenden (siehe unten) und die Ziffer aus dem Dateinamen in eine Liste von Liste mit 256+1 Elementen speichern
Line 151: Line 164:
 outcsv.close() outcsv.close()
  </file>  </file>
 +=== Lösungen ===
 +<hidden Liste von Bilddateien>
 +<code python pixelist_from_directory.py>
 +
 +import gpanel  #um bilder einzulesen
 +import os  #um Verzeichnisse zu listen
 +import csv #um CSV-Dateien zu lesen.
 +
 +# Pfad zu den Bilddateien
 +digitsdirectory = 'C:/temp/digits/train/'
 +
 +
 +def getPixeListFromFilePath(filepath):
 +    img = gpanel.getImage(filepath)
 +    w = img.getWidth()
 +    h = img.getHeight()
 +    pixellist = []
 +    for y in range(h):
 +        for x in range(w):
 +
 +            # color is ein Objekt mit verschiedenen Attributen, u.a. red, green, blue.
 +            # bei grau sind rot=gruen=blau, d.h., eine Farbe auslesen reicht.
 +            # siehe auch https://docs.oracle.com/javase/7/docs/api/java/awt/Color.html
 +            color = img.getPixelColor(x, y)
 +
 +            # umlegen auf das Intervall [-1,1] zwecks Normalisierung
 +            value = color.red / 255 * 2 - 1
 +            # an liste anhaengen
 +            pixellist.append(value)
 +            
 +    return pixellist
 +
 +# Lese Ziffer aus Dateiname aus.
 +def getDigitFromFileName(filename):
 +    return int(filename.split('_', 3)[2])
 +
 +
 +# leere Liste fuer alle Trainingsdaten der Form [-0.93,0.331,....,0.99,3]
 +
 +trainingset = []
 +
 +# durch alle files im Ziffernverzeichnis loopen
 +for filename in [filename for filename in os.listdir(digitsdirectory) if filename.endswith("gif")]:
 +
 +    # Ziffer auslesen
 +    currdigit = getDigitFromFileName(filename)
 +
 +    # Pixelliste von Datei auslesen
 +    currpixellist = getPixeListFromFilePath(digitsdirectory + filename)
 +
 +    # Der Pixelliste die Ziffer anhaengen
 +    currpixellist.append(currdigit)
 +
 +    # Gesamte Liste dem trainingsset anhaengen.
 +    trainingset.append(currpixellist)
 +
 +# Das Trainingsset kann jetzt verwendet werden
 +# print(trainingsset)
 +
 +</code>
 +</hidden>
 +
 +==== $k$ Nearest Neighbours auf Pixellisten ====
 +
 +=== Ziele ===
 +  * $k$ nearest neighbours in $k$-Dimensionen durchführen
 +  * Einem Bild aus den Testdaten einen Zahlwert auf Grund der Trainingsdaten zuordnen
 +Der Code unten soll als Grundlage für den eigentlichen kNN Klassifizierer gelten. Es muss einizg noch die Funktion ''assignClass'' geschrieben werden. Das Einlesen der Trainings- und Testdaten ist bereits programmiert, ebenso die Umwanldung eines Bildes in eine Pixelliste.
 +
 +<code python stub_knn_digits.py>
 +
 +import gpanel  #um bilder einzulesen
 +import csv #um CSV-Dateien zu lesen.
 +
 +# Pfad zu den Bilddateien
 +digitsdirectory = 'C:/temp/digits/'
 +
 +
 +# Trainingsdaten
 +trainingset = list(csv.reader(open(digitsdirectory + 'digitsdata.csv', 'r'), delimiter=',', quoting=csv.QUOTE_NONNUMERIC))
 +# Testdate
 +testset = list(csv.reader(open(digitsdirectory + 'digitsdatatest.csv', 'r'), delimiter=',', quoting=csv.QUOTE_NONNUMERIC))
 +
 +
 +
 +def getPixeListFromFilePath(filepath):
 +    img = gpanel.getImage(filepath)
 +    w = img.getWidth()
 +    h = img.getHeight()
 +    pixellist = []
 +    for y in range(h):
 +        for x in range(w):
 +
 +            # color is ein Objekt mit verschiedenen Attributen, u.a. red, green, blue.
 +            # bei grau sind rot=gruen=blau, d.h., eine Farbe auslesen reicht.
 +            # siehe auch https://docs.oracle.com/javase/7/docs/api/java/awt/Color.html
 +            color = img.getPixelColor(x, y)
 +
 +            # umlegen auf das Intervall [-1,1] zwecks Normalisierung
 +            value = color.red / 255 * 2 - 1
 +            # an liste anhaengen
 +            pixellist.append(value)
 +            
 +    return pixellist
 +
 +def assignClass(point, k):
 +    # Funktion die den Abstand von point zu den k naechsten
 +    # Nachbarn im trainingset berechnet mit Pythagoras in 
 +    # 16x16=256 Dimensionen
 +    
 +    # gibt die Mehrheitsklasse wieder. Tip: Liste mit Häufigkeiten zurückgeben.
 +    return(assignedclass)
 +
 +## Funktion Testen
 +
 +# Auf einem Testbild
 +testpoint = getPixeListFromFilePath(digitsdirectory + 'test/testimg_993918225911_2_.gif')
 +print(assignClass(testpoint, 100))
 +
 +# Auf Testdaten 
 +for i in range(len(testset)):
 +    print(assignClass(testset[i][0:255], 30))
 +    print(testset[i][256])
 +</code>
 +=== Lösungen === 
 +<hidden Loesung>
 +<code python convertimagetopixellist.py>
 +
 +import gpanel  #um bilder einzulesen
 +import os  #um Verzeichnisse zu listen
 +import csv #um CSV-Dateien zu lesen.
 +
 +# Pfad zu den Bilddateien
 +digitsdirectory = 'C:/temp/digits/'
 +
 +
 +def getPixeListFromFilePath(filepath):
 +    img = gpanel.getImage(filepath)
 +    w = img.getWidth()
 +    h = img.getHeight()
 +    pixellist = []
 +    for y in range(h):
 +        for x in range(w):
 +
 +            # color is ein Objekt mit verschiedenen Attributen, u.a. red, green, blue.
 +            # bei grau sind rot=gruen=blau, d.h., eine Farbe auslesen reicht.
 +            # siehe auch https://docs.oracle.com/javase/7/docs/api/java/awt/Color.html
 +            color = img.getPixelColor(x, y)
 +
 +            # umlegen auf das Intervall [-1,1] zwecks Normalisierung
 +            value = color.red / 255 * 2 - 1
 +            # an liste anhaengen
 +            pixellist.append(value)
 +            
 +    return pixellist
 +
 +# Lese Ziffer aus Dateiname aus.
 +def getDigitFromFileName(filename):
 +    return int(filename.split('_', 3)[2])
 +
 +
 +# leere Liste fuer alle Trainingsdaten der Form [-0.93,0.331,....,0.99,3]
 +
 +trainingset = []
 +
 +# durch alle files im Ziffernverzeichnis loopen
 +for filename in [filename for filename in os.listdir(digitsdirectory) if filename.endswith("gif")]:
 +
 +
 +    # Ziffer auslesen
 +    currdigit = getDigitFromFileName(filename)
 +
 +    # Pixelliste von Datei auslesen
 +    currpixellist = getPixeListFromFilePath(digitsdirectory + filename)
 +
 +    # Der Pixelliste die Ziffer anhaengen
 +    currpixellist.append(currdigit)
 +
 +    # Gesamte Liste dem trainingsset anhaengen.
 +    trainingset.append(currpixellist)
 +
 +# Das Trainingsset kann jetzt verwendet werden
 +# print(trainingsset)
 +
 +# Alternative: Einlesen von Datei
 +# trainingset = list(csv.reader(open('C:/temp/digitsdata.csv', 'r'), delimiter=',', quoting=csv.QUOTE_NONNUMERIC))
 +
 +
 +# Funktion, die einem Punkt eine Klasse auf Grund der k naechsten Nachbarn zuweist.
 +def assignClass(point, k):
 +
 +    # Liste um die Distanzen zu speichern.
 +    distlist = []
 +
 +    # Fuer jeden den Punkt point die Distanzen zu allen Punkten in datalist berechnen.
 +    for i in range(len(trainingset)):
 +        dist = 0
 +
 +        # schlaufe ueber die 256 pixel
 +        for j in range(255):
 +            dist += (point[j] - trainingset[i][j]) ** 2
 +            
 +        distlist.append([trainingset[i][256], sqrt(dist)])
 +
 +    # das waere ein sehr Pythonesquer Weg mit Lambda-Funktionen
 +    # nearest = sorted(distlist,key=lambda result:result[1])
 +
 +    # definiere eine Funktion, welche das zweite Element zurueckgibt.
 +    def sortFunction(item):
 +        return item[1]
 +
 +    # Sortiere die liste! Achtung: Man koennte auch ohne key Arbeiten, wenn Distanz an 1. Stelle waere
 +    nearest = sorted(distlist, key=sortFunction)
 +
 +    # Zaehle k naechsten Klassennummern und entscheide ueber Klasse.
 +    # Liste die die Anzahl Vorkommen von 0,1,2,...,9 enthaelt. An Stelle 0 die Null, an Stelle 1 die Eins etc.
 +    classcounts = [0] * 10
 +
 +    # Loope durch die k naechsten Nachbarn und erhoehe den classcount jeweils um 1 
 +    # wenn eine Ziffer (Klasse) der jeweiligen Ziffer gefunden worden ist.
 +    for j in range(k):
 +        classcounts[int(nearest[j][0])] += 1
 +
 +    # Rueckggabe der Anzahl/Klasse.
 +    # Wenn anstelle der Anzahl/Klasse die eigentlich Ziffer zurueckgegeben werden soll, 
 +    # kann dies mit    classcounts.index(max(classcounts)) passieren.
 +
 +    return classcounts
 +
 +
 +# Teste das Programm
 +
 +testpoint = getPixeListFromFilePath(digitsdirectory + 'testimg_993918225911_2_.gif')
 +print assignClass(testpoint, 100)
 +testpoint = getPixeListFromFilePath(digitsdirectory + 'testimg_899286188295_1_.gif')
 +print assignClass(testpoint, 100)
 +testpoint = getPixeListFromFilePath(digitsdirectory + 'testimg_809996871457_8_.gif')
 +print assignClass(testpoint, 100)
 +testpoint = getPixeListFromFilePath('C:/temp/man1.gif')
 +print assignClass(testpoint, 5)
 +
 +</code>
 +</hidden>
  • kurse/efcomputergrafik/knn.1585048334.txt.gz
  • Last modified: 2020/03/24 12:12
  • by Simon Knaus