2D-Systemtheorie

Aus II-Wiki
Zur Navigation springen Zur Suche springen


Anmerkung: Er meinte ja das dies nur eine Auswahl alter Fragen ist, und wir mehr behandelt haben, deswegen stehen ab und zu mal noch ein paar ausgedachte Fragen. ^^

Grundlagen Fouriertransformation, Abtastung

Allgemeines zur Fouriertransformation

  • periodisches Signal: diskrete Verteilung der Signalenergie über der Frequenz (Fourierreihe)
  • aperiodisches Signal: kontinuirliches Frequenzspektrum (Fourierintegral)

1-Dimensional

Fourierreihe

Fourierintegral

Rücktransformation:

2-Dimensional

Fourierreihe

Fourierintegral

Eigenschaften der Fouriertransformation

jegliche Herleitungen zu diesem Abschnitt sind in der Vorlesung oder im Skript zu finden

Linearität:

  • es gilt das Superpositionsprinzip
  • daraus ergibt sich das nur lineare Zusammenhänge beschreibbar sind

Faltung:

Verschiebung:

Drehung:

  • Drehe ich das Signal im Zeitbereich, dreht sich in gleicher Weise das Spektrum mit

Ähnlichkeit:

Parselval'sches Theorem:

Separierbarkeit:

  • wenn sich eine Funktion in ihren Abhängigkeiten separieren lässt, kann die Transformierte ebenfalls separat ermittelt werden

Wichtige Funktionen

Dirac-Stoß:

  • Ausblendeigenschaft:
  • Foriertransformierte:

Liniensingularitäten:

  • als Liniensingularität wird eine Linie in x- oder y-Richtung definiert
  • Transformierte dazu sind:

Projektions-Schnitt-Theorem:

  • Ein Schnitt durch ein 2D-Signal kann durch Multiplikation mit einer Liniensingularität erfolgen:
  • Schnitt entlang der y-Achse:
  • für jedes fest erfolgt eine Integration in -Richtung, d.h. eine Projektion auf die die -Achse

1. Aufgabe

Ermitteln Sie die Ergebnisse der folgenden Faltung !
2DSys-Aufgabe 1.png
2D-ST-Dreieck und Kreis.png 2DSys-Aufgabe 1-Faltung.png

2. Aufgabe

Beschreiben Sie das resultierende Spektrum einer Gitterstruktur (20 µm helle, 20 µm dunkle Streifen im Bild, die um 45 ° zu den Koordinatenachsen gedreht sind) bei Abtastung mit einer CCD-Matrix mit quadratischem Pixelraster von 11 µm * 11 µm Pixelabstand und einer Pixelappertur von 8 µm * 8 µm.
2DSys-Aufgabe 2.png
  • ist schwer zu transformieren, deswegen nutzt man ein Abtastgitter
ohne Thumbnail-Rahmen

3. Aufgabe

Wie kann man ein Schachbrettmuster mathematisch beschreiben? Wie sieht das zugehörige 2-D-Ortsfrequenzspektrum aus?
  • Bei einem Schachbrettmuster ist die Separierbarkeit anwendbar
  • mit als Breite eines Feldes ergibt sich:

//Bild

4. Aufgabe

Beschreiben Sie die Pixelapertur und ihre Wirkung auf das abgetastete Bild bei Interline-Transfer-Matrizen, die im Frame- bzw. Field-Mode betrieben werden!
Pixelraster 11 µm * 11 µm; Pixelapertur = 1 in einem Rechteck H = 6 µm * B = 10 µm außerhalb = 0.
2D-ST-Aufgabe 4.png
  • im Frame Integration Mode werden abwechseln gerade und ungerade Zeilen ausgelesen
  • im Filed Integration Mode werden immer zwei übereinander liegende Pixel in einen Registerplatz geschoben
  • beachte auch interlaced und non-interlaced Variante
  • ändert aber nichts am Ergebnis
  • im Frame Integration Mode bleibt die Pixelappertur (PA) wie gewohnt
  • im Field Integration Mode ensteht durch das ausschieben zweier Pixel in ein Register eine unterbrochene PA
  • diese entspricht einem Rechteck, gefaltet mit zwei Dirac-Stößen
  • im Frequenzbereich entspricht das einer si-Funktion multipliziert mit einer cos-Funktion
  • die si-Funktion bildet hierbei die Einhüllende für den Cosinus

5. Aufgabe

Beschreiben Sie die resultierende Pixelapertur bei einem Zeilenscanner (Bewegter Zeilensensor v = 1 m/s)! Pixelabstand 18 µm, Pixelapertur 12 µm * 12 µm, Abtastraster ta=20 µs; Integrationszeit=10 µs.
2DSys-Aufgabe 5.png
  • die Gesammtappertur ergibt sich aus der Pixelappratur gefaltet mit einem Rechteck in y-Richtung der Länge
  • in x-Richtung habe ich also weiterhin ein Rechteck als Pixelappertur (PA)
  • in y-Richtung ergibt sich ein Trapez (Faltung von PA und Rechteck )
  • in 3D quasi ein Prisma

6. Aufgabe

Wie kann man ein Sechseckabtastraster systemtheoretisch beschreiben?
ohne Thumbnail-Rahmen
  • Abtastgitter wird im folgenden mit AG abgekürzt
  • mit und ergibt sich:
ohne Thumbnail-Rahmen

7. Aufgabe

Wie kann man das Abtastverhalten von Einchip-Farbsensoren erfassen? Überlegen Sie eine Beschreibung für das Y-Signal bei Streifensensoren ( R G B R G B ... ), wenn ein unbuntes Bild angeboten wird? Wie wird eine Objektszene, die nur Rot-Anteile enthält, abgebildet?
  • Abtastverhalten des Streifensensors:
  • das Y-Signal ist das Helligkeitssignal = Summe der Farbwerte

Digitale Filter

8. Aufgabe

Beschreiben Sie im Orts- und im Ortsfrequenzbereich die Wirkung eines 3 * 3 Boxoperators, der auf dem im Speicher abgelegten Bild arbeitet. Vergleichen Sie mit einem theoretischen Boxoperator im Originalbild (z.B. unscharfe Abbildung bei quadratischer Blende).
Als Beispiel zur Erklärung soll dieser Tiefpassoperator herhalten:
[1 1 1]
[1 1 1]
[1 1 1]
Dieser wirkt als Weichzeichner auf ein Bild im Speicher. Sprich von einem 3x3 Pixelraster werden alle Werte addiert und bilden den Mittelwert (mit 1/9 Faktor) für das Pixel in der Mitte.
Der Operator ist in x und in y in [1 1 1] separierbar und kann durch 3 Dirac-Stöße beschrieben werden. ( ___|_|_|___, der Mittlere entsprich Position x=0)
Diese Stöße lassen sich im Frequenzbereich als Kosinus darstellen, der um 1/3 nach oben verschoben ist (Verschiebung durch Dirac-Stoß on der Mitte = Mittelwer)
Dieser nach oben verschobene Kosinus ähnelt einer Sinc-Funktion, die bei der Fouriertransformation eines idealen Rechtecks ([1 1 1]x[1 1 1] Boxoperator) entstehen würden
Na eben nicht. Genau das ist ja der Unterschied, bzw. das Problem von dem Boxoperator. Ein mit Offset versehener Kosinus ist deutlich anders als ein Sinc. Der Kosinus lässt theoretisch beliebig hohe Frequenzanteile ohne Dämpfung durch, während diese durch die Sinc-Funktion schon längst auf Null gedämpft worden wären. Praktisch ist die maximale Frequenz in einem Bild natürlich durch die Rasterung begrenzt. Als Übung kannste ja zum Beispiel mal für diese maximale Frequenz den Unterschied bestimmen. Pale 22:55, 21. Sep. 2010 (UTC)
Im Ortsbereich werden die Dirakstöße mit den Spekturm des Bildes gefaltet. Im Frequenzbereich wird die Kosinus-Funktion (~Sinc-Funktion) mit dem Bildspektrum multiplitiert.
(so etwa hat der Prof. mir es erklärt)

9. Aufgabe

Beschreiben Sie die Wirkung eines mehrfach angewandten Boxoperators (2*, 3*, 4*, 5*,...) im Orts- und im Ortsfrequenzbereich !
Die Anwendung eines Boxoperators kann durch eine Faltung im Ortsbereich beschrieben werden. Da die Faltung die Eigenschaft der Assoziativität aufweist, können mehrfache Faltungen auch ausgerechnet werden bevor sie Auf das Bild angewendert werden.
Im Ortsbereich: Rechteck, Dreieck, ..., x-mal so breites gaußänliches Filter
Im Frequenzbereich: Si, Si^2, ..., Si^x sehr schmaler Gauß

10. Aufgabe

Beschreiben Sie im Orts- und im Ortsfrequenzbereich die Wirkung eines Gradientenoperators
1 0 -1
1 0 -1
1 0 -1
Da der Operator separierbar ist lassen sich die Wirkeungen in die verschiedenen Richtungen einzeln betrachten
Im Orstbereich:
  • Gradient in x-Richtung [1 0 -1]
  • Summe in y-Richtung [1 1 1]'
Im Frequenzbereich:
  • Hochpass in Kx-Richtung
  • Tiefpass in Ky-Richtung

11. Aufgabe

Beschreiben Sie die Wirkung des Gradientenoperators (1 0 -1) auf einem Bild mit hellem Quadrat (=1) auf dunklem Untergrund (=0) !
 Bild =
    0     0     0     0     0     0     0     0     0
    0     0     0     0     0     0     0     0     0
    0     0     0     0     0     0     0     0     0
    0     0     0     1     1     1     0     0     0
    0     0     0     1     1     1     0     0     0
    0     0     0     1     1     1     0     0     0
    0     0     0     0     0     0     0     0     0
    0     0     0     0     0     0     0     0     0
    0     0     0     0     0     0     0     0     0
 
 OP_Grad =
    1     0    -1
 
 >> convn(Bild,OP_Grad,'valid')
 
 ans =
    x     0     0     0     0     0     0     0     x
    x     0     0     0     0     0     0     0     x
    x     0     0     0     0     0     0     0     x
    x     0    -1    -1     0     1     1     0     x
    x     0    -1    -1     0     1     1     0     x
    x     0    -1    -1     0     1     1     0     x
    x     0     0     0     0     0     0     0     x
    x     0     0     0     0     0     0     0     x
    x     0     0     0     0     0     0     0     x

Rauschen

12. Aufgabe

Erläutern Sie die wichtigsten Rauschquellen bei CCD-Bildaufnahmesystemen!
Rauschart Ort Verteilung Zusammenhang
Photonenrauschen Silizium Binomial-
Dunkelsignalrauschen Silizium Poisson-
Transferrauschen Silizium Poisson-
thermisches Rauschen der FET Ausgangsverstärker der CCD-Matrix Gauß-
thermisches Rauschen der Elektronik Analog-Elektronik Gauß-
Taktdurchgriffe, Phasenjitter, Klemmfehler, Netzbrummen Analog-Elektronik Gauß-
Digitalisierungsrauschen ADU Gleich-
  • Festmusterrauschen
  • Lichtdurchlässigkeit der Optik -> Staub
  • Inhomogeneitäten im Halbleiter
  • optische elemente auf dem Pixel
  • Zeilenklemmrauschen

13. Aufgabe

Wie kann man Rauschen vermindern? Für welche Rauschanteile greift Ihr Vorschlag?
  • thermisches Rauschen: niedrige Temperatur
  • Photonenrauschen, Dunkelsignalrauschen: lässt sich wegmitteln
  • Festmusterrauschen: wenn Muster bekannt ist, kann es abgezogen werden
  • Shading: für jedes Pixel einen Faktor bestimmen, mit dem multipliziert wird. Korrekturbild kann dann zB im Speicher abgelegt werden.
  • Transferineffektivität: Faktor bestimmen der übrig bleibt
  • TP-Verhalten: ist meist bei Transferineffektivität mit drin, wird sonst in ähnlicher Weise behandelt
  • Smear: durch Filter beseitigen
  • bei CCD-Matrizen im Konsumerbereich ist Photonenrauschen bestimmende Größe
  • wenn bessere SNR erreicht werden sollen, muss ich Photonenrauschen vermindern
    • wird durch Zusammenhang zwischen Sättigungselektronenzahl und Rauschen ermöglicht
    • benötigt werden also größere Pixelflächen und insgesamt mehr Si-Fläche
    • bisher nur in wisschenschaftlichen CCD's

14. Aufgabe

Was geschieht mit den örtlichen und zeitlichen Rauschanteilen bei der Subtraktion von zwei Bildern der gleichen Szene?
Mittelwert verschwindet, Streuung steig

15. Aufgabe

Erläutern Sie anhand des Zusammenhanges zwischen Bildspektrum und 1-D-Spektrum auf dem Übertragungskanal den Einfluß eines asynchronen Störers!

16. Aufgabe

Wie wirkt ein Tiefpaß im 1-D-Signalweg auf das Spektrum des abgetasteten, übertragenen und im Speicher abgelegten Bildes?
2D-Bildinformationen werden Zeilenweise über den 1D-Signalweg übertragen. Dabei entsteht ein kontinuierlier Fluss von Intensitätswerten (Unter Vernachlässigung des Dunkelsignales und der Zeilenpause).
Ein Tiefpass würde also hohe Frequenzen unterdrücken und somit alle großen Intensitätsunterschiede auf dem 1D-Signalweg Dämpfen.

Reale CCD-Bauelemente - CCD-Kamera

17. Aufgabe

Nennen Sie Abweichungen vom idealen Verhalten bei realen CCD-Bauelelementen bzw. CCD-Kameras. Wie kann man diese Effekte erfassen und wie kann man sie korrigieren?
Sample & Hold Schaltung:
  • Abtast- und Haltestufen führen zu einer charakteristischen Signalbeeinflussung
  • CCD-Signal ist durch Taktsteuerung korreliert
  • es entsteht Signalpegel der Pixelinformation enthält
  • Rauschen welches bei der Ladungsgenerierung wirkt, ist auch korreliert zum Takt
  • thermisches Rauschen verursacht zufällige Signalanteile unkorreliert zum Takt (auch 1/f-Rauschen genannt)
  • unkorrelierte Anteile werden durch CDS (Correlated Duopple Sampling) wesentlich beeinflusst
2D-ST-S&H.png


Idealer ADU:
  • Fehler im Intervall gleichverteilt von
  • Signalleistung ergibt sich aus Moment zweiter Ordnung:
  • führt zu:
  • Rauschleistung:
Realer ADU:
  • zusätzliche Effekte die SNR verschlechtern:
  • Nichtlinearität der globalen Kennlinie
  • differentielle Nichtlinearität (DNL)
  • ist Maß für unterschiedliche Breite der einzelnen Code-Stufen
  • mit
  • DNL führt zu einer integralen NL (Wanderkennlinie):
  • Sprünge in der Kennlinie
  • zur Beschreibung des Verhaltens, tatsächliches SNR bestimmen und effektiv auflösbare Bitbreite berechnen
  • Tests durchführen um Kennlinie zu erhalten: sollten statisch und dynamisch sein
  • bei statischen Testsignalen separate Spannungsquelle verwenden (Batterie)
  • bei dynamischen Test auf statistische Auswertung stützen
  • reale Auflösung:

18. Aufgabe

Welcher Zusammenhang besteht zwischen den in Aufg. 17 gefundenen relevanten Parametern und einer sinnvollen Auflösung eines ADU in einer Kamera ?
  • reale Auflösung:

19. Aufgabe

Wie kann ein in Aufg. 17 entwickelter Korrekturalgorithmus umgesetzt werden (Software, Hardware) ?

Optische Abbildung

20. Aufgabe

Bestimmen Sie die beugungsbedingte Punktverwaschungsfunktion eines Objektives bei Abbildung aus dem Unendlichen ! (f=50 mm, k=4, \lambda=550nm)
Hinweis:

21. Aufgabe

Aus welchen Anteilen setzt sich die Punktverwaschungsfunktion bei einer realen optischen Abbildungsanordnung zusammen ? Wie kann man die Anteile beschreiben ?

22. Aufgabe

Bestimmen Sie die notwendige Defokussierung, wenn damit eine Tiefpaßfilterung (Antialiasing) mit einer 3dB-Grenzfrequenz von 20 LP/mm erreicht werden soll !

23. Aufgabe

Wie kann man das Übertragungsverhalten eines optischen Abbildungssystems messen?

24. Aufgabe

Skizzieren Sie das Grundprinzip einer optischen Anordnung zur Generierung des Betrages des Ortsfrequenzspektrums einer transparenten Struktur!

25. Aufgabe

Mit inkohärenter monochromatischer Beleuchtung sollen Gitterstrukturen (Gewebe) vermessen werden. Die Gitterstrukturen befinden sich in der Blendenebene eines Punktabbildungssystems. Die Brennweite beträgt f = 100 mm. Das Abbildungssystem besteht aus einer Tandemoptik.
2DSys-Aufgabe 25.png
Bild: Abbildungssystem zur Vermessung von Gitterstrukturen
Die CCD-Matrix mit 752 * 582 Pixel hat einen Pixelabstand von 11 µm * 11 µm. Zur Erfassung der Ergebnisstruktur sollen die Peakabstände d > 20 Pixel und d < 200 Pixel sein.
Die Beleuchtung erfolgt mit
Welche Gitterkonstanten (Bereich von - bis) können vermessen werden ?
  • Es gilt:
  • Die Transmission des Gitters wird der einfachheit halber als einfache Cos-Funktion angenommen. Ausserdem ist es ausreichend wenn man das Gitter nur in eine Dimension betrachtet, da die andere Dimension analog zu berechnen ist und das Gitter in X und in Y Richtung seperierbar ist. Damit gilt für das Gitter als Blendenfunktion:
    • D ... Gitterkonstante
    • ... eine Dimension des Gitters
  • Eingesetzt in die obige Formel gibt das:
  • Durch Substitution von kann das Integral in das Fourierintegral überführt werden
  • Damit ergibt sich:
  • Dieses Integral liefert nur für Werte von , (Gleichanteil ist in unserer Formel nicht erfasst, muss aber mit betrachtet werden, da dieser in der Realität auftritt) und Werte
  • Damit die Ergebnisstruktur erfasst werden kann muss zwischen diesen Werten aussreichend Platz sein. Daher folgt:
  • Nach Rückgängigmachen der Substitution und umstellen erhält man:
    • ... Wellenlänge
    • ... Brennweite der Linse
    • ... Abstand in m (Pixelanzahl * Pixelbreite)
  • Damit ergibt sich für das Größte Gitter:
  • und für das kleinste Gitter:
  • Damit ergibt sich als Ergebnis: 0.025mm < D < 0.25 mm

Andere Abbildungssysteme

26. Aufgabe

Geben Sie die Abbildungseigenschaften bei Röntgenabbildung (Abbildungsmaßstab, Schwächung) an !
Definition Abbildungsmaßstab
Schwächung beim Materialdurchgang
  • ist material- und wellenlängenabhängig
  • zurückgelegter Weg d und damit die Schwächung sind winkelabhängig
  • Ausbreitung der Röntgenstrahlung in Materie und Vakuum geradlinig
  • kann deshalb durch Zentralprojektion beschrieben werden
  • Abbildungsmaßstab:
  • gilt damit immer

27. Aufgabe

Erläutern Sie das Grundprinzip von Tomografieverfahren !
Aufbau Tomographieverfahren
Prinzip Tomographieverfahren
  • Röntgenstrahlen werden beim Durchgang durch Materie gedämpft (Gleichung siehe Aufgabe 26)
  • logarithmiertes Verhältnis, stellt Projketion der Materialeigenschaften längs des Weges dar:
  • Tomographieverfahren bietet Möglichkeit Information über Materialeigenschaften in bereit zu stellen
  • Vereinfachungen:
    • Strahlungsquellen mit paralleler Ausstrahlung
    • Strahlungsquellen und Empfänger parallel zu einander
    • Strahlungsquellen und Empfänger mehrfach vorhanden oder gemeinsam um Objekt gedreht
    • ausgesendete Strahlung soll auf konstant sein
  • Logarithmus der empfangenen Strahlung entspricht Projektion der Materialeigenschaften entlang des Weges
  • heißt soviel wie, das eine Projektion im Ortsbereich einem Schnitt durch den Ortsfrequenzbereich entspricht
  • Verteilung kann somit aus vielen Projektionen/Schnitten gewonnen werden:
  • Projektionen werden jeweils um den Winkel gedreht und im Ortsfrequenzbereich eingetragen
  • solange durchführen bis Ortsfrequenzbereich hinreichend dicht belegt ist
  • 2D-Rücktransformation liefert dann gewünschte Verteilung in Ebene
  • 3D-Abbildung durch Hinzunahme der Koordinate