A LGORITMOS DE GRÁFICAS EN

RASTER II
C o m p uta c i ó n G r á fica
 Midpo int Line Algorithm

Bresenham (1965)

Utiliza solamente aritmética de enteros evitando funciones de redondeo.

Permite el cálculo de (xi+1, yi+1) de manera incremental.

Se extiende a aritmética de punto flotante.

Se extiende a círculos.

No se extiende fácilmente a cónicas arbitrarias.

Modificación: midpoint line algorithm (Pitteway,1967), Van Aken (1984).

Para líneas y círculos el algoritmo de Van Aken se reduce a la formulación de Bresenham.

Minimiza el error (distancia) a la primitiva original.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 Midpo int Line Algorithm

Suponemos que la pendiente está entre 0 y 1 (las demás se manejan por

reflexión sobre los ejes principales)

Llamamos al punto extremo inferior izquierdo (x0,y0) y el superior derecho

(x1,y1).
NE
Suponemos haber elegido al pixel P en (xp,yp) y debemos elegir entre el pixel E
Q (east pixel) y el NE (northeast pixel).
M
Sea Q el punto de intersección entre la línea deseada y la línea de la rejilla
P = (xp,yp) E x=xp+1.

Previous Choices Choices Si M, el punto medio entre NE y E, está abajo de la línea, elegimos E, si está
pixel for current for next
pixel pixel arriba elegimos NE.

El error será siempre menor o igual a 1/2.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015
 Midpo int Line Algorithm

Determinar si el punto medio está arriba o abajo de la línea.

Representar la línea con una función implícita con coeficientes a, b y c: F(x,y) = ax + by + c = 0.

Si dy = y1-y0 y dx = x1-x0 podemos escribir y = (dy/dx)x + B.

F(x,y) = dy⦁x - dx⦁y + B⦁dx = 0

a=dy, b=-dx, c= B⦁dx en la forma implícita.

Se puede verificar que F(x,y) es cero sobre la línea, positiva bajo la línea y negativa para puntos
sobre la línea.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 4. iluminar el pixel en (xi , Round(yi )), donde Round(yi ) = ⇥(0.5 + yi )⇤

F (x, y) = ax + by + c = 0
Midpo int LinedyAlgorithm
y= x+B
dx
F (x, y) = dy · x dx · y + B · dx = 0
⇥
1
Necesitamos solamente calcular F (M ) = F xp + 1, yp +
2

y verificar el signo.

Definimos una variable de decisión d=F(xp+1, yp+1/2)

Si d>0, elegimos el pixel NE.

Si d<0, elegimos el pixel E.

Si d=0, podemos elegir cualquiera, por convención elegimos E.

La siguiente posición de M y el valor de d dependerá del pixel elegido en la última iteración.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 ular yi = mxi + B para cada F xi(x, y) = ax + by + c = 0
minar el pixel en (xi , Round(yi )), donde Round(yi ) =dy ⇥(0.5 + yi )⇤
y= x+B
yMidpo xint Line Algorithm
dy dx
= +B
F (x, y) = ax + by dx+ c = 0
F (x, y) = dy · x dx · y + B · dx = 0
F (x, y) = dy dy· x dx · y + B · dx = 0 ⇥
y= x+B 1
dx F (M ) = F x ⇥
p + 1, yp +
Si elegimos E, incrementamos M un paso en la dirección
1 x: 2
F (x, y) =Fdy(M· )x = F
dx · xyp++B1,· ydxp⇥+
= 20 ⇥
1 1
dnew = F xp ⇥ + 2, yp + ⇥ = a (xp + 2) + ⇥b yp + + c.
1 12 1 2
dnew = F Fx(M p+ ) 2,
= yFp +xp + 1,= yap(x+p + 2) + b yp + + c.
pero 2 2 2
⇥ ⇥ ⇥
1 1 1
dnew = F xp + d2,old yp= + a (xp=+a1)(x+ p+ b 2)y+p+ b yp ++ c. + c.
2 2 2
⇥
y restando dold de dnew escribimos: 1
dold = a (xp + 1) + b yp + + c.
2
dnew = dold + a = dold + dy

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 F (x, y) = ax + by + c = 0
F (x, y) = ax + by + c = 0
Midpo
y=
dy int Line Algorithm
x+B
dx dy
y= x+B
dx
F (x, y) = dy · x dx · y + B · dx = 0
F (x, y) = dy · x dx · y + B · dx = 0
Si elegimos NE, incrementamos M un paso 1
⇥ en la dirección
⇥ x y un paso en la
F (M ) = FF (M xp) + 1, y + 1
dirección de y: = F pxp + 1, yp +
2 2
⇥ ⇥ ⇥ ⇥
3 3 3 3
dnew = F xdnewp += 2, Fyp +
xp + 2,=ypa+(xp +=2) +pb+ 2)
a (x yp ++b yp + + c. + c.
2 2 2 2
⇥ ⇥
1 1
y restandodold
dold=dea (x
dnew escribimos:
= (x +
+ 1) + b yp +
dp old a p 1) + b y +
+ c. 2
p + c.
2
dnew = dold + a = dold + dy
dnew = dold + a + b = dold + dy dx

Nótese que esta versión del algoritmo funciona para pendientes entre 0 y 1.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 Midpo int Line Algorithm

void MidpointLine( int x0, int y0, int x1, int y1, int value )
{
int dx = x1-x0;
int dy = y1-y0;
int d = 2*dy-dx; /* valor inicial de d */
int incrE = 2*dy; /* incremento para E */
int incrNE = 2*(dy-dx); /* incremento para NE */
int x = x0;
int y = y0;

WritePixel(x,y,value); /* the start pixel */

while( x < x1 ){
if( d <= 0 ){ /* choose E */
d += incrE;
x++;
} else { /* choose NE */
d += incrNE;
x++;
y++;
}
WritePixel(x,y,value); /* the selected pixel */
} /* while */
} /*MidpointLine */

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 Conversión de sca n de un
círculo

Cualquier algoritmo para generar un círculo puede aprovechar de

la simetría de la figura.

Por cada punto calculado, se pueden graficar 8 reflejando en

ejes a 45 grados.

P1 = ( x, y) P5 = (-x,-y)
P2 = ( y, x) P6 = (-y,-x)
P3 = (-y, x) P7 = ( y,-x)
P4 = (-x, y) P8 = ( x,-y)

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 Conversión de sca n de un
círculo

Definir matemáticamente un círculo centrado en el origen:

2 2 2
Polinomio de segundo orden: y = r - x .

Cada coordenada x en el sector de 90° a 45° se

encuentra avanzando r/√2.
2 2
Cada coordenada y se encuentra evaluando √(r -x )
para cada paso de x.

Ineficiente porque par cada punto hay cuadrados y

raíces.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 Conversión de sca n de un
círculo

Definir matemáticamente un círculo centrado en

el origen:

Usando funciones trigonométricas.

x = rcosθ, y=rsinθ.

Avanzamos θ de θ a π/4 y calculamos cada

valor de x y y.

Aún más costoso que el método anterior.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015
 Algoritmo de Bresenha m pa ra
círculos

Aprovechar la simetría del círculo.

Generar puntos solamente en un sector

de 45 grados (de 90 a 45, esto es,
en las direcciones +x y -y).

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 Algoritmo de Bresenha m pa ra
círculos

Dos posibles acciones:

moverse una unidad en dirección x.

moverse una unidad en dirección x y una unidad en dirección -y.

(xi,yi) - coordenada del último pixel elegido al entrar al paso i.

D(T) = distancia del origen al pixel T al cuadrado menos la

distancia del origen al círculo real al cuadrado.

D(S) = distancia del origen al pixel S al cuadrado, menos la

distancia del origen al círculo real al cuadrado.

COMPUTACIÓN GRÁFICA | ALGORITMOS DE RASTER | ENERO-JUNIO 2015 | 28/01/2015

 A LGORITMOS DE
RELLENA DO DE PO LÍGO NOS
C o m p uta c i ó n G r á fica
 Polígono en un plano

• Un polígono en un plano es una serie de segmentos =(s1,s2,...,sn) con

n∈N* y para i=1,...,n el segmento si es de la forma si=[Pi-1,Pi] donde Pi es un
punto del plano para i=0,...,n, con Pn=P0.!
• Los puntos Pi, con i=0,..,n se llaman vértices del polígono .!
• Los segmentos si se llaman arístas del polígono .!
• Para i=2,...,n, el origen Pi-1 del segmento si es igual al extremo del segmento
si-1. !
• El origen P0 del segmento s1 es igual al extremo Pn del segmento sn.
 Polígono en un plano

P2
P4
P1 P3
s3
s2 s4 s5 s2
s4
autointersección
P1 P5 s1 s3
P3

s1 s6
P0 = P4 P2

P0 = P6

• Decimos que un polígono se auto-intersecta si dos segmentos que componen

al polígono y que son no-consecutivos se intersectan.
 Puntos interiores y exteriores
Sea Π=(s1,s2,...,sn) un polígono y sea M un punto en el plano que no forma parte de
ninguna de las aristas del polígono Π.
Consideramos una semirecta Δ que sale de M y no contiene a ninguno de los
vértices del polígono Π.
Decimos que el punto M es interior al polígono Π si la semirecta Δ intersecta un
número impar de segmentos del polígono Π.
Decimos que el punto M es exterior al polígono Π si la semirecta Δ intersecta un
número par de segmentos del polígono Π.

I M

17
 Polígonos convexos

Un polígono se llama convexo si para todos los puntos M y M’ interiores al

polígono Π, el segmento [M,M’] está enteramente compuestos de puntos
interiores al polígono Π.

M M’ M M’
 Principio de rellenado de polígonos

Consideremos un polígono Π del plano, donde los vértices tienen coordenadas

enteras.
El problema de rellenado consiste en encontrar un algoritmo que permita recorrer
todos los puntos de coordenadas enteras del polígono que sean interiores.
D
11
Calcular los segmentos horizontales F
10
incluídos en el polígono. 9

Determinar por cada horizontal cuáles 8

b d
son los pixels al interior del polígono 7 a c
para desplegarlos con su color 6
E
apropiado.
5

4
C
Haciendo eso de manera iterativa para 3

todas las horizontales se rellenará el 2

A
polígono. 1
B
1 2 3 4 5 6 7 8 9 10 11 12 13
 Rellenado de polígonos

El rellenado, para cada línea horizontal se descompone en 3 etapas:

Encontrar las extremidades de los segmentos horizontales incluídos en el polígono:
intersecciones de la línea horizontal con las aristas del polígono.
Ordenar los extremos obtenidos en orden creciente de las coordenadas x.
Desplegar todos los pixels al interior del polígono entre pares de extremidades.
Utilizar regla de paridad para determinar si un pixel se encuentra al interior o al
exterior de un polígono.
 Recorrido de aristas

¿Cómo encontrar los extremos sobre cada horizontal?

rellenar cada arista con el algoritmo del punto medio.

11 Algunos pixels están al exterior del polígono.

10 Problema si dibujamos varios polígonos diferentes con
9 diferentes colores.
8 Podría mostrar polígonos disjuntos como juntos.
7
No vamos a utilizar el algoritmo del punto medio.
6

5 Algoritmo incremental: elegir los pixels a la derecha de la

4
arista ideal si se trata de una arista izquierda y viceversa.
3 Memorizamos los pixels extremo generados por este
2
algoritmo en una tabla.
1

1 2 3 4 5 6 7 8 9 10 11 12 13