Bresenham-Algorithmus (in Arbeit): Unterschied zwischen den Versionen

Aus Geometrie-Wiki
Wechseln zu: Navigation, Suche
(Algorithmus für Gerade mit der Steigung m <= 1)
(Algorithmus für Gerade mit der Steigung m <= 1)
Zeile 98: Zeile 98:
 
Den ersten Punkt den wir zeichnen ist natürlich der Punkt <math>O (0,0)</math>, denn dies haben wir so festgelegt. Nun müssen wir für den nächsten Punkt entscheiden, ob wir nur ein Pixel (allgemein: einen Plotterschritt) nach rechts zum Punkt <math> T </math> gehen oder ein Pixel nach rechts und nach oben zum Punkt <math> S </math>. Um dies zu entscheiden müssen wir die Abstände, die die jeweiligen Punkte zu der idealen Geraden haben, vergleichen. Der Punkt mit dem kleineren Abstand wird dann angesteuert.<br>
 
Den ersten Punkt den wir zeichnen ist natürlich der Punkt <math>O (0,0)</math>, denn dies haben wir so festgelegt. Nun müssen wir für den nächsten Punkt entscheiden, ob wir nur ein Pixel (allgemein: einen Plotterschritt) nach rechts zum Punkt <math> T </math> gehen oder ein Pixel nach rechts und nach oben zum Punkt <math> S </math>. Um dies zu entscheiden müssen wir die Abstände, die die jeweiligen Punkte zu der idealen Geraden haben, vergleichen. Der Punkt mit dem kleineren Abstand wird dann angesteuert.<br>
 
Hierfür betrachten wir die Differenz der Abstände <math> |d_t|-|d_s|</math>.<br>
 
Hierfür betrachten wir die Differenz der Abstände <math> |d_t|-|d_s|</math>.<br>
 +
Da wir nur die y-Achsenabweichung der gezeichneten Geraden mit dem Auge wahrnehmen, aber der folgende Zusammenhang besteht, können wir mit den reinen Abständen, die wir aus der Hesseform gewinnen rechnen:
 +
 +
zu zeigen ist: <math> d_s >= d_t  \Leftrightarrow s >= t</math><br>
 +
 +
es gilt:<br><br>
 +
<math> d_s >= d_t \ \  d_s = s \cdot sin (\alpha) \ \  d_t = t \cdot sin (alpha) </math><br><br>
 +
<math> \Leftrightarrow d_s = s \cdot sin (\alpha) >= t \cdot sin (\alpha) = d_t </math><br><br>
 +
<math> \Leftrightarrow  s  >= t </math><br><br>
 +
 +
 +
 +
 +
<ggb_applet width="1260" height="880"  version="4.2" ggbBase64="UEsDBBQACAgIALCFaEIAAAAAAAAAAAAAAAAWAAAAZ2VvZ2VicmFfamF2YXNjcmlwdC5qc0srzUsuyczPU0hPT/LP88zLLNHQVKiu5QIAUEsHCEXM3l0aAAAAGAAAAFBLAwQUAAgICACwhWhCAAAAAAAAAAAAAAAADAAAAGdlb2dlYnJhLnhtbOVb627bRhb+nT7FQCj2V0TP/dKVW7hZBE3jtEWdXRSLBQJKHEmsJVIlKV+CvMC+x75A9hX2f/eV9swMKYm6OJblJF4XtjQiOdfvzPnO+Uip983VdIIubFGmeXbcIRHuIJsN8iTNRsedeTXs6s43X3/RG9l8ZPtFjIZ5MY2r4w6PaGfZDo4iwlzjNIFeBn07pIR0ZdwXXR6TuNsfYtwVzKohvLjipIPQVZl+leU/xFNbzuKBPRuM7TQ+zQdx5fscV9Xsq6Ojy8vLqBk9yovR0WjUj67KpINg5ll53Kk/fAXdtRpdMl+dYkyOfnl1GrrvpllZxdnAdpBb1Tz9+osnvcs0S/JLdJkm1RhmT5XpoLFNR2NYp9GwqCNXawaLndlBlV7YEtquHPpFV9NZx1eLM3f9SfiEJov1dFCSXqSJLY47OKJMSS2kIYZLzRjXHZQXqc2qujKpBz1quutdpPYy9Os++SF5B1V5PunHrkskDHr3DlFMMXrqChIKCoWU4RIO5zALBQ0FD4UIdXhozkNVHurwUIczsHlapv2JPe4M40kJMKbZsAATLo7L6npi/ZTqE0sEyFNYVpm+hcoMA6wBdziP8VP3kvDiuMZ7ZZ16ZZ3ELeIdIm72vmDIzZv4+buC14cyHCpfEBwKUl/U7s3jJQ9cEbvTijg2ap+Bq2K+E0kjluNSLZ6C4z1VMK4Qm+OSlUFDn3uM2YxIqFxZKuP4KTW7l0oPWecKunw5pCDCD0nF9iHJgUZdLJSvQAtj+X//2hiSHbTM5YgC33ZEyW8YMQxwvwMq3KKbhmtCSeryJhjubVK9o4YMe/WEUDl2destXdlp6abIjCdFRJAA0pAKOEwgYqBQjjwoIgJxAYdEI+lKhZjjC44Y0sjVIwx56hMa3rjnEokE9OVOqkAqiHEkGCKeMDkCFJAnXcCEMqghBBLQyI1O3LBMIi7hgGnEYYKObpWjNAbt4BgGp4gRxFxbohCVSFKkHGUT7phcajd36JQiiZF0TYGzga8DV0MLjZhbDXjBLC/TBbhjO5ktrOJxTLPZvKqxq88PpkmDY5WvVU/ywfm3a2DbuKyaz1AJAtYyLoYA1gqbT3qTuG8nkF2cuX2A0EU8cW7u+x/mWYUWJBPOjYp4Nk4H5ZmtKmhVol/ji/g0ruzVc6hdNhP0Q/tw3rPzwSRN0jj7G2wS14XrEC2juyOvJrprUQ8zyPMiObsuYeugq7/bIj/udCERUBEXVEuiFZOSAGVfh2sEQnfENROMaU2IpnCpHMRu1zONIwiUXEEd11ACh11vvWawCYPbi8Xq4iu7WBMaFc7vVg5elN/mk+WpWZ5m1bN4Vs0Ln67BSIVb10k2mliPr6dfSHwG5/386iwAy0Jfr69ncITDDPqjZ/kkLxB4JRVAeqO67IfS13FTW9TCvg72NXBjqTRZXCeG+hq+7IfS1wLTh6nVSyXNMgluhklLzzfQeWtn+o3j0qh5llanzUGVDs7rpZLQ4If5tA97rt6U7T7JffXZO1rbZr1zW2R2EvZSBsac5/My7O7FDn3Sm5f2p7gan2TJz3YEfvlT7Kixgq5D1eWUEztIp9AwnK/Bi51h/wpTDWcTOypss8SJz5ADtP4qXt3ZG6d9V8+LfPoiu3gNu2Ztqr2jZj29clCkM7c7UR+4+twu91+SljEwfbLaDhZfwioGjnUAyMqB2EHxvBrnhc+BwXUhlqDv//OvLLMFsCVsSOe1V7PClk5NNEZ5ba8q2Dxw4bjzp9/mefXnMxRKP5Cd2Cnky6jy+xj4oOq0W3rXB9uhvP8r8M8i3oQ6K0aB6xu7v9nZKJ7MxrHL2GvcJvE1zHoVSd/hqzxp4+u5rLRFOlzGf/C/VzCEk0kLkgvxcNEMxElR/eT8GvSR01WGcymEoxFqGKYwu2t3GhONmSbaEE4VbNK3QbuFnRmQ2Qkqa4H6eg9Q2d1AxS1C+eygiggzBZILaNgxsdNe1w4+o4GaqXZqR4OOEDtRbYHkSXiBwLNDEMILwr0zQnEGrOF9D0LZLPDXzNpAfWHC8GEG3fmIsTIdzxNljU+zy1YhAEHrluqiSCtqh7NrBHNbvP7yCPAiNV70E+D1/BHgJT4hXi828WqH/mXEfbB46chIbKTWmgIvaSo9eiJikEISSbhhRippPjqW3z8CLGXEjJYSC0aMNpgyjyWLAEYOnE8oZe5f3guYg3w6jbMEZV7+nkLm21mqrhg79kMxccAG0OZVcyEOXdUdbNjFJdEL1OPDKeGAkLyEtku2oBg26jrkgG3XRCCbDMVSEMWFJnxdH1Qg+84zSF+8iKlqueI/fJcmifWiNuin37LQpA776XQ2SQdptbc9ngV7PN+wR38Pe/Q/5CafyiC05llc82yXLDv7ZCi/yJy6ARTWoI4D1P0NqM8GYxA4cJidVzej3manVrs7eYTU3gKu6IficBsIFm3jGx4JDOmm4BR2vr+D6+gGR8ooRgmWzCiNubxH9jmzI3d+zQot0LxBnm0YJHlT3myHsu66gdc1uAfxJci9sBKOnFIShMAblpxR3nIKJ6UMDCzhDaCXVNzgIeJmDwHQJi72LPY8xKrNewDn1s7czZcfs9dFnJXugVpb/B9uxE0CS958wJk2jXhHLyI43Pmh7o79fVEZiYiB4KK1T3cEpDttZhPgZVhxzuBNKs3Uw7Lidh0uWjocvOb2Qlw8lrsb4JHAgIZTQohskjEBJyEzMGBLoEgt97u9IddgrW4Pq3w09zcwY4ISASUUtdpShBIGp7BQHFIvvgvWNuf8WFTjfJRn8WRLwvQ88E28wTeDPRKmwUNJmLpbUtXrHYntW/dMMqIUorW7lcSxEfqh5FaDXVZ5uU9C9fLhZFHA/BFjIHMl04YaUTOFgbNgLa7APIoC9wezkEgJ0BYczkhFwJr3mEbd6A3PduGe7OENyf+rN7CIQ1KLIQZrrTjmn1zQ7XKHZJdZTvdxh9MH5A5SRAp0sxQG/rUQ4akAMFEklJaUGyMVaAv58f1he0Z6uktL2P2SUPuZ05yWjFAMsNXgDAqyEslDVO26GIClwEQQLLCgWvBaVwgFeQ14B8NgEHaAN3x0WeGfF2/X53bDhP/958029A8WFxaC2q49TGReQ0sAGKyMBGSU4drIOmze2w0scksbkztYpHkOWwxWHKe5PTmZ5Jc/2+HEXnlADxN1L3cpueF+LjT8zClt6zGJAbsLygTGWkhSSzhGNQNVLsGrQLeREMS7MiKQ/RtCYK8IyGdvVHQP0YGGu6LO7//ex4Gg9gccSP3h/Ge77lMt3beHmFaPQkwLcBgqtNCKSk2b+MQiCSdhp0CSrBkRdD8trVuY7qGk9aNQ0hyCvgBXwy7bhWSrUdLggBxyAQOZF2jpnV8U2OtOeH2btUUVJ/skqCd3ezz3MRJUChoB8lIDPCUoF8KEe66ERVphJ+EkSDasVJMqgbRzyRIR2DG+/uj5aR1cTzYAH+0XXEcPRqtBW6OJZhI8X0N8VUJuDa+1JoD6QmiuCUgCTiE/1Q85vG434kkw4umGEcf7GXH8kIzo9ByVoNM40VI1jyraEgOkX2NEsKjTI0YwYcgiFXggRtweUkgrpJzcPqSQRxJSpH/epBSDeEJJE6bBE8F7laLuK2gSs/3iNG+Benp7UPkjyX0o54YD7Un/DZkFqAbCNzZEYQzxZucN7+2Ymhamr26PqfmsX8i4r22qFfAJRHCi/YOEGlFKjYA/QhhsWKFul/lsk0f1LaK1Z5d0202j39/vJZjebwomqjGlMF8IjQprZT6TYrrLly+2KSZ+4B2HbfZ4ucsem3cgfn//huxnkTfkQzYxfzyb7AiPuEU7X6J/+CWhL/cIlPizElBansav7S/rP2hY5aX6l2CBlsiSlvCHn2xCJkQU01orgZUWDS1pLEFEuMcUDK99fe8WmLdzkoD4HngflJk8ZLzdzw80FQo0nOIYwkATV53IUJB3ghmo2Pnzg6PVH4P4H2jVv+P++n9QSwcI84pVEq0LAAB3PgAAUEsBAhQAFAAICAgAsIVoQkXM3l0aAAAAGAAAABYAAAAAAAAAAAAAAAAAAAAAAGdlb2dlYnJhX2phdmFzY3JpcHQuanNQSwECFAAUAAgICACwhWhC84pVEq0LAAB3PgAADAAAAAAAAAAAAAAAAABeAAAAZ2VvZ2VicmEueG1sUEsFBgAAAAACAAIAfgAAAEUMAAAAAA==" showResetIcon = "false" enableRightClick = "false" errorDialogsActive = "true" enableLabelDrags = "false" showMenuBar = "false" showToolBar = "false" showToolBarHelp = "false" showAlgebraInput = "false" useBrowserForJS = "true" allowRescaling = "true" />
 +
 +
 +
  
  

Version vom 8. März 2013, 16:55 Uhr

Inhaltsverzeichnis

Allgemeines zum Bresenham Algorithmus

Motivation des Algorithmus

Der Bresenham Algorihtmus, ist ein Verfahren um Geraden möglichst "gut" auf Anzeigegeräten zu zeichnen. Hier heisst "gut", dass die Abweichung zwischen dem gezeichneten und dem gedachten Objekt möglichst gering ist.

Das Problem beim Darstellen einer Strecke auf einem Anzeigegerät ist, dass das erzeugte Bild nur durch endlich viele Punkte aufgebaut ist, dadurch entstehen "Lücken" beim zeichnen. Da unser Auge nur eine endliche Auflösung verarbeiten kann, scheint uns ein Bild auf dem Monitor nicht durch einzelne Punkte aufgebaut zu sein, sondern es entstehen für uns geometrische Formen, die wir mit unserer Vorstellung von mathematischen Objekten in Einklang bringen können.



Eine Gerade wie wir sie auf dem Bildschirm wahrnehmen.
Gerade in normaler Größe



Hier sieht man dieselbe Gerade um das Achtfache vergrößert.


Gerade bis auf die Pixeleben vergrößert

Effektivität des Algorithmus

Dieser Algorithmus ist auch in heutiger Zeit in seiner Geschwindigkeit und Einfachheit ungeschlagen, da die Rechenoperationen sich auf Addition und die Multiplikation mit der Zwei beschränken lässt. Dadurch kann dieser Rechenalgorithmus direkt in die Hardware der Grafikkarten implementiert werden. Potenzen: Obwohl in den Rechentermen Potenzen auftreten, kann man diese umgehen. Hierfür betrachten wir den Ausdruck:

 n^2

die zweite Binomische Formel liefert:

 (n - 1)^2 = n^2 - 2 \cdot n + 1 \ \Leftrightarrow n ^ 2 = ( n - 1)^2 + 2 \cdot n - 1

durch rekursives Anwenden dieser Formel lassen sich Potenzen solange vereinfachen bis die Basis Eins erreicht ist und damit die Potenz verschwindet.

Multiplikation mit Zwei:
Eine Multiplikation mit Zwei ist für einen Computer eine elementare Rechenoperation, da die sogenannte Shift Operation (Verrückungsoperation) durchgeführt werden kann.

Wir erinnern uns an die Darstellung von Zahlen im Binärsystem.


Wertigkeit der Stelle 2^7 = 128 2^6 = 64 2^5 = 32 2^4 = 16 2^3 = 8 2^2 = 4 2^1 = 2 2^0 = 1
Binärzahl 0 1 0 0 1 1 0 1

Hier wird die Zahl (im Dezimalsystem) dargestellt

  0 \cdot 2^7 + 1 \cdot 2^6 + 0 \cdot 2^5 + 0 \cdot 2^4 + 1 \cdot 2^3 +  1 \cdot 2^2 + 0\cdot 2^1 + 1 \cdot 2^0 = 77

Multiplizerien wir nun mit zwei, so folgt:

2 \cdot ( 0 \cdot 2^7 + 1 \cdot 2^6 + 0 \cdot 2^5 + 0 \cdot 2^4 + 1 \cdot 2^3 +  1 \cdot 2^2 + 0\cdot 2^1 + 1 \cdot 2^0 )  = 154

 \Leftrightarrow

 0 \cdot 2^8 + 1 \cdot 2^7 + 0 \cdot 2^6 + 0 \cdot 2^5 + 1 \cdot 2^4 +  1 \cdot 2^3 + 0\cdot 2^2 + 1 \cdot 2^1 = 154

Also die Zahl 154 (Dezimal) in Binärdarstellung:

Wertigkeit der Stelle 2^7 = 128 2^6 = 64 2^5 = 32 2^4 = 16 2^3 = 8 2^2 = 4 2^1 = 2 2^0 = 1
Binärzahl 1 0 0 1 1 0 1 0

Vergleichen wir nun die Darstellung der Zahl 77 mit der Darstellung der Zahl 154, so fällt auf, dass durch die Multiplikation mit Zwei jedes Bit um eine Stelle nach links verschoben wurde.

Bresenham-Algortihmus für Geraden

Algorithmus für Gerade mit der Steigung m <= 1

Wir wollen uns nun mit der folgenden Problemstellung beschäftigen, wir wollen eine Gerade, welche in Hesseform gegeben ist, mit einem Plotter (z.B. Monitor) zeichnen lassen. Wir wollen für den ersten Teil nur Geraden mit einer Steigung  m <= 1 betrachten.

Sei die Gerade  g : \vec{n_n} \cdot ( \vec{r} - \vec{a} ) = 0 gegeben. Hierbei ist  n_n = \begin{pmatrix} n_x \\ n_y \end{pmatrix} der normierte Normalenvektor,  \vec{a} ein gegebener Ortsvektor eines Punktes der Geraden.
Wir geben den Startpunkt der Geraden auf unserem Plotter die Koordinaten (0,0), d.h. jede Gerade, die wir zeichnen erhält damit ein eigenes Koordinatensystem. Den ersten Punkt den wir zeichnen ist natürlich der Punkt O (0,0), denn dies haben wir so festgelegt. Nun müssen wir für den nächsten Punkt entscheiden, ob wir nur ein Pixel (allgemein: einen Plotterschritt) nach rechts zum Punkt  T gehen oder ein Pixel nach rechts und nach oben zum Punkt  S . Um dies zu entscheiden müssen wir die Abstände, die die jeweiligen Punkte zu der idealen Geraden haben, vergleichen. Der Punkt mit dem kleineren Abstand wird dann angesteuert.
Hierfür betrachten wir die Differenz der Abstände  |d_t|-|d_s|.
Da wir nur die y-Achsenabweichung der gezeichneten Geraden mit dem Auge wahrnehmen, aber der folgende Zusammenhang besteht, können wir mit den reinen Abständen, die wir aus der Hesseform gewinnen rechnen:

zu zeigen ist:  d_s >= d_t  \Leftrightarrow s >= t

es gilt:

 d_s >= d_t \ \  d_s = s \cdot sin (\alpha) \ \  d_t = t \cdot sin (alpha)

 \Leftrightarrow d_s = s \cdot sin (\alpha) >= t \cdot sin (\alpha) = d_t

 \Leftrightarrow  s  >= t









Nun gehen wir nicht mehr von dem ersten Punkt (Ursprung)  O= \begin{pmatrix} 0 \\ 0 \end{pmatrix} aus, sondern von einem allgemeinen Punkt  P = \begin{pmatrix} p_x \\ p_y \end{pmatrix} , dementsprechend lautet dann der Punkt  T = \begin{pmatrix} p_x + 1 \\ p_y \end{pmatrix} und der Punkt  S = \begin{pmatrix} p_x + 1 \\ p_y + 1 \end{pmatrix}




Achtung:
Wie wir unter dem Punkt Hesseform gesehen haben, kommt es bei dem Vorzeichen der Abständen der Punkte auf die Lage der Punkte im Bezug zur Geraden und dem Ursprung an. Da der Punkt  S nicht zwischen Gerade und Ursprung liegt hat, ist der Abstand  d_s positiv, analog ist der Abstand  d_t negativ.

Also folgt mit:

|t| - |s| = -t - s = - (t+s)

= -(\vec{n_n} \cdot (\vec {t} - \vec {a}) - \vec n_n \cdot ( \vec{s} - \vec {a}))

= - ( \vec{n} \cdot ( \vec{t} + \vec{s}))

= \begin{pmatrix} n_x \\ n_y \end{pmatrix}  \cdot \begin{pmatrix} p_x + 1 + p_x + 1 \\ p_y + p_y + 1 \end{pmatrix}

= \begin{pmatrix} n_x \\ n_y \end{pmatrix}  \cdot \begin{pmatrix} 2 \cdot p_x + 2 \\ 2 \cdot p_y + 1 \end{pmatrix}

= - ( n_x \cdot 2 \cdot p_x \ + 2 \cdot n_x + n_y \cdot 2 \cdot p_y + n_y ) =: d

Nun entscheidet das Vorzeichen von  d welcher Punkt angesteuert wird:

Für  d >= 0  der Punkt  S

Für  d < 0  der Punkt  T

Algorithmus für Gerade mit der Steigung m > 1

Wir wollen nun den Fall für Geraden betrachten deren Steigung größer ist als Eins. Ist dieser Fall nun auch betrachtet, steht uns der Algorithmus für alle Geraden zur Verfügung.

Die Vorausetzungen seien wie in dem Fall für die Steigung kleiner Eins.